ECMWF logo:

The ECMWF Model Grids


Introduction

The ECMWF data held at the BADC is archived at a number of grid resolutions. The aim of this page is to explain the difference between the following grid types:
  1. Spectral Resolution
  2. Gaussian Grids
  3. Regular Latitude-Longitude Grids

Information is also provided on the interpolation methods used to convert spectral data to a regular grid.


Spectral Resolution

The ECMWF Integrated Forecasting System (IFS) uses a spectral model with model levels calculated in the vertical axis. Spectral models calculate their parameters using spherical harmonics as opposed to grid points. This is computationally more efficient than a regular grid-point model. Spectral model resolution is typically expressed with the number of levels (e.g. T159L60 - where 'T' is the spectral resolution and 'L' is the number of model levels) , although the number of levels is often omitted. As the Operational model is improved, the resolution increases. The current Operational model runs at T511, equivalent to a grid size of approximately 40 km.

Only upper air data is stored in spectral form so the complementary surface data sets are provided on a Gaussian grid. The ECMWF data archived at the BADC is stored in two spectral resolutions, T159 and T106.

T159 resolution

The ERA-40 Re-analysis and Operational Analyses data are stored at a resolution of T159. The surface data for this data is provided on an N80 Reduced Gaussian grid.

T106 resolution

The ERA-15 Re-analysis project and Operational Analyses are both stored at a resolution of T106. The surface data for the ERA-15 data set is provided on an N80 Reduced Gaussian grid whilst the surface data for the Operational data set is on an N80 Full Gaussian grid.


Gaussian Grids

The ECMWF provides surface fields on an N80 Gaussian grid to be used alongside the spectral upper air fields (on model levels). Gaussian grids are used in the model calculations in preference to regular grids because they provide superior numerical accuracy. The ECMWF has used an N80 Gaussian grid in its Operational spectral model since 1 May 1985.

There are two types of N80 Gaussian grid, the full and reduced Gaussian grid. The BADC holds ECMWF data stored on both types of Gaussian grid. The ECMWF web site contains information on the different types of Gaussian grid.

N80 Full Gaussian Grid

The N80 Full Gaussian grid is equivalent to a 1.125° grid along its longitudinal axis. This means that it has 320 grid points encircling the earth at each latitude. The grid is symmetric about the Equator with 80 latitudinal points in each hemisphere (see Table 1). The ECMWF surface Operational data held at the BADC is stored on an N80 Full Gaussian grid.

N80 Reduced Gaussian Grid

The N80 Reduced Gaussian grid is more complicated than the N80 Full Gaussian grid. The reduced grid has the same number of latitudinal grid points but the number of longitudinal grid points varies. The number of longitudinal grid points decreases from the equator to the poles as shown in Table 1, with 320 grid points around the equatorial region but only 18 grid points at 89.1416°N and 89.1416°S. This reduces storage time and computing resources when running the model. The ECMWF ERA-15 Re-analysis project and the ERA-40 Re-analysis project surface data is available in the N80 Reduced Gaussian grid.

Table 1. Grid points on the "full" and "reduced" N80 Gaussian grids.
"Reduced" N80 Gaussian grid"Full" N80 Gaussian grid
No. from poleLatitude (degrees)No. of points on a Lat. circleNo. of points on a Lat. circle
189.141618320
288.029425320
386.910836320
485.790640320
584.669945320
683.548954320
782.427860320
881.306664320
980.185372320
1079.064072320
1177.942680320
1276.821290320
1375.699896320
1474.5784100320
1573.4570108320
1672.3356120320
1771.2141120320
1870.0927128320
1968.9712135320
2067.8498144320
2166.7283144320
2265.6069150320
2364.4854160320
2463.3639160320
2562.2425180320
2661.1210180320
2759.9995180320
2858.8780192320
2957.7566192320
3056.6351200320
3155.5136200320
3254.3921216320
3353.2707216320
3452.1492216320
3551.0277225320
3649.9062225320
3748.7847240320
3847.6632240320
3946.5418240320
4045.4203256320
4144.2988256320
4243.1773256320
4342.0558256320
4440.9343288320
4539.8129288320
4638.6914288320
4737.5699288320
4836.4484288320
4935.3269288320
5034.2054288320
5133.0839288320
5231.9624288320
5330.8410300320
5429.7195300320
5528.5980300320
5627.4765300320
5726.3550320320
5825.2335320320
5924.1120320320
6022.9905320320
6121.8690320320
6220.7476320320
6319.6261320320
6418.5046320320
6517.3831320320
6616.2616320320
6715.1401320320
6814.0186320320
6912.8971320320
7011.7756320320
7110.6542320320
729.53270320320
738.41120320320
747.28970320320
756.16820320320
765.04670320320
773.92520320320
782.80370320320
791.68220320320
800.56070320320


Regular Latitude-Longitude Grids

The simplest type of grid is the latitude-longitude regular grid. This represents the Earth's surface or upper air layers as a number of grid points at regular intervals (defined in degrees).

1° x 1° grid

Some of the ERA-40 Re-analysis data is archived on a regular 1° x 1° grid.

1.125° x 1.125° grid

The ECMWF Operational data set is archived on a regular 1.125° x 1.125° grid for upper air model levels.

2.5° x 2.5° grid

There are two ECMWF data sets held on a regular 2.5° x 2.5° grid. These are the ERA-15 Re-analysis and the Operational data sets. Both surface and upper air data (on pressure levels) are archived at this resolution.


Interpolation methods

The ECMWF archives its upper air data in spectral form and its surface or single level data on a Gaussian grid. In order to convert the data to a regular latitude-longitude grid the BADC has used interpolation software. The method used to interpolate ECMWF data has varied depending on the dataset. Although the interpolation methods produce similar outputs, information on the methods may be relevant to users performing sensitivity tests using ECMWF data from more than data set.

The ERA-15 Re-analysis and Operational data was interpolated using bi-linear interpolation routines provided by the ECMWF (see the ECMWF web pages for details).

The same software and method of interpolation was used to convert ERA-40 data from the N80 Reduced Gaussian grid to a 1° x 1° regular grid are described here.

The BADC also makes the XCONV/CONVSH software package available from the BADC web pages. XCONV is Graphical User Interface that allows the user to visualise, subset and convert data files. CONVSH is a command line utility that allows access to the same features as XCONV via a shell prompt. Some users will download XCONV/CONVSH to interpolate data themselves.

XCONV/CONVSH also uses a form of bi-linear interpolation to convert from one grid resolution to another. Users of the data should be aware that data values will vary slightly depending on the method of interpolation applied. The bi-linear interpolation used for ERA-40 may produce occasional spurious results such as very slightly negative values in precipitation fields. Users who wish to access the data in its original model form should use the spectral and N80 Gaussian gridded data instead of the interpolated regular gridded data.

To ensure that the two interpolation methods did not produce significantly different output the BADC ran a comparison of data sets interpolated by the ECMWF routines and the XCONV/CONVSH package. Figures 1 and 2 show two data sets interpolated from one original spectral data set of temperature on model level 60 (nearest the Earth's surface) on 1 January 1989. The comparison shows that the interpolation routines produce almost identical results.


Figure 1. Level 60 air temperatures (01/01/1989) interpolated by XCONV/CONVSH package.

Figure 2. Level 60 air temperatures (01/01/1989) interpolated by ECMWF routines.

The difference between the two interpolated data sets was less than 0.0012K (around 0.0005%) throughout the globe. This confirms that the there is minimal difference between the XCONV/CONVSH interpolation routine and the ECMWF's own interpolation routines.