Versions Compared

Old Version 12

changes.mady.by.user Paul Dando

Saved on

compared with

New Version 13

changes.mady.by.user Paul Dando

Saved on

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.

...

Section

Gaussian grids

Naming convention

The Gaussian grids are defined by the quadrature points used to facilitate the accurate numerical computation of the integrals involved in the Fourier and Legendre transforms. The grids are labelled by N where N is the number of latitude line lines between the pole and the equator.  For example, for the N640 Gaussian grid, there are 640 lines of latitude between the pole and the equator giving 1280 latitude lines in total.

The grid points in latitude, θk, are given by the zeros of the Legendre polynomial or order N:  PN0(μk=sinθk) = 0.  A consequence of this is that a Gaussian grid has:

  • latitude lines which are not equally spaced;
  • no latitude points at the poles;
  • no line of latitude at the equator.

Regular (or full) Gaussian grid

A regular Gaussian grid has the following characteristics:

  • there are 4N longitude point points along each latitude circle;
  • the longitudinal resolution in latitude-longitude is given by 90o/N;
  • the points get closer together (i.e. more crowded) as the latitude increases towards the poles;
  • the total number of grid points is 8N2.

Reduced (or quasi-regular) Gaussian grid

A reduced Gaussian grid:

  • has the same number of latitude lines (4N) as the corresponding regular Gaussian grid;
  • the has a decreasing number of longitude points decreases towards the poles;
  • has a quasi-regular grid spacing in distance at each latitude;

  • provides a uniform CFL (Courant–Friedrichs–Lewy) condition.

Up to and including IFS cycle 41r1, ECMWF has used a standard reduced Gaussian grid.  This has 4N longitude points at the latitude nearest to the equator, with the number of longitude points reducing in blocks as the latitudes approach the poles.

With the horizontal resolution increase at IFS cycle 42r1, ECMWF introduces a slightly different form of the reduced Gaussian grid which is referred to as the octahedral reduced Gaussian grid or, more simply, the octahedral grid.

Section

Relationship between spectral truncation and grid point resolution - linear, quadratic and cubic grids

The relationship between the spectral resolution, governed by the truncation number T, and the grid resolution depends on the number of points on the grid at the equator, 4N, at which the shortest wavelength field is sampled:

linear:  each wavelength is sampled by 2 grid points → 4N L = 2(TL + 1)

quadratic: each wavelength is sampled by 3 grid points → 4N Q = 3(TQ + 1)

cubic: each wavelength is sampled by 4 grid points → 4N C = 4(TC + 1)

Until the implementation of IFS cycle 18r5 on 1 April 1998, the IFS used a quadratic grid. The introduction of the two-time level semi-Lagrangian numerical scheme at IFS cycle 18r5 made possible the use of a Linear Gaussian Grid reflected by the TL notation.  The linear grid has been used since then, up to and including IFS cycle 41r1. At IFS cycle 42r1, the cubic grid will be used and will be indicated by the TC.

Note
iconfalse

By changing the number of points at which the shortest wavelength is sampled on the grid, the effective grid point resolution can be increased while keeping the spectral truncation number constant.

Grid point computations, which are often non-linear, benefit from the higher resolution of the grid produced with cubic sampling.  Such a grid has no aliasing, less numerical diffusion and more realistic surface fields.

Section

Increasing the horizontal resolution

The relationship between the spectral truncation, T, the sampling (linear, quadratic or cubic) and the grid point resolution, N, allows for three possible approaches for increasing the horizontal resolution:

  1. increase the number of wave numbers in the spectral representation but keep the same Gaussian grid by reducing the sampling frequency (i.e.g, quadratic → linear, increase T with TQ →TL );
  2. increase the number of wave numbers in the spectral representation and the Gaussian grid resolution but keep the sampling frequency constant ;(i.e. increase T and increase N).
  3. keep the number of wave numbers in the spectral representation constant and increase the Gaussian grid resolution by increasing the sampling frequency (i.e., TL →TC, and increase N).
Note
iconfalse

The horizontal resolution upgrade at IFS cycle 42r1 will be is achieved by:

  • increasing the Gaussian grid resolution by keeping the spectral truncation constant but using cubic sampling;
  • introducing a new form of the reduced Gaussian grid - the octahedral grid.

This new grid will be referred to simply as a cubic octahedral grid.  This will be indicated with the notation TCO.

Section
Column

What is the octahedral grid ?

The octahedral grid has been inspired by the Collignon Projection of the sphere onto a an octahedron. It is a form of reduced Gaussian grid with the same number of latitude lines located at the same latitude values as those of a standard Gaussian grid but with the number of longitude points at each latitude circle computed according to the formula:

Mathdisplay
\begin{eqnarray*}
\mbox{N}_{lat}(lat_N) & = & 20 \\
\mbox{N}_{lat}(lat_i) & = & \mbox{N}_{lat}(lat_{i+1}) + 4,  \mbox{ for } i=\mbox{N} - 1,\ldots,1
\end{eqnarray*}

In other words, there are 20 longitude points at the latitude circle closest to the poles with the number of points increasing continuously by 4 at each latitude towards the equator. This is in contrast to the standard reduced Gaussian grid where there are 'jumps' between blocks of latitudes with a constant number of longitude points (a restriction which was imposed by the Fast Fourier Transform routines being used prior to IFS cycle 41r2).

As a consequence, the zonal resolution of the octahedral grid varies more with latitude than for the standard reduced Gaussian grid. This can be seen the figures to the right.

Note in particular that the octahedral grid has 4N 4N + 16 longitude points at the latitude circle closest to the equator whereas the standard reduced Gaussian grid has 4N 4N longitude points at the latitude circle closest to the equator.

There are also fewer total grid points in the octahedral grid compared to the standard reduced Gaussian grid.For example, the N1280 octahedral grid has about 20% fewer grid points than the equivalent N1280 linear grid. Generally, there are 4xNx(N+8) grid points in the octahedral grid or resolution N.

The octahedral grid has been shown to improve the calculation of local derivatives in grid point space.

Column

Figure 1

Comparison of the resolution for the reduced Gaussian grids with latitude.  The full red and yellow curves show the resolution for the standard reduced Gaussian grids at N1024 and N1280.  Note that the resolution remains more or less constant at 10km and 8km, respectively, as the latitude varies.  The corresponding curves for the N1024 and N1280 octahedral grids are shown by the dashed and full blue lines, respectively.  Note that the resolution for the N1280 octahedral grid varies from about 8 km at the equator, increasing to about 10 km at 70oN and 70oS before decreasing again towards the poles. The resolution of the N640 standard reduced Gaussian grid used at IFS cycle 41r1 is at about 16 km. Also shown are the regular Gaussian grids at N1024 (black dashed line) and N1280 (black full line).

Figure 2

Comparison of the zonal variation in resolution for the N1280 standard reduced Gaussian grid (left) with the corresponding octahedral grid (right) on the surface of the sphere