Versions Compared

Old Version 7

changes.mady.by.user Paul Dando

Saved on

compared with

New Version 8

changes.mady.by.user Paul Dando

Saved on

Key

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

...

  • The error occurs whenever the grid spacing Dj=360/Nj at a specific latitude line, j, with Nj longitude points has a remainder smaller than 10 microdegrees (1/100000 of a degree).
  • The error accumulates linearly from 0° to 360° along the line of latitude.
  • The error is minimal at grid points to the east and maximal at those to the west of the 0° meridian.
  • The error is larger the greater the number of longitude points (Nj) along the line of latitude.

The When the error occurs, the incorrect computation of the longitude values leads to two potential issues:

  1. The longitude points in the source grid are computed incorrectly.  This leads to an incorrect computation of interpolation weights and hence to incorrect values at some points of the output grid.
  2. The nearest grid points used for the interpolation may be incorrectly identified due to a numerical 'shift' of the input grid cell.

...

  • Interpolations to rotated latitude-longitude grids are unaffected
  • Transformations from spherical harmonic components to grid point fields are unaffected.
  • Interpolations of wave (WAM, ENS-WAM etc) fields are unaffected.
  • Wave ?

...

For which cases is the problem fixed ?

The issue is fixed in EMOSLIB cycle 000430 and newer for the following cases:

  • interpolations from original reduced  or octahedral reduced Gaussian grids to global, unrotated regular Gaussian grids (gridType=regular_gg);
  • interpolations from original reduced  or octahedral reduced Gaussian grids to global, unrotated latitude-longitude grids (gridType=regular_ll).

In which software versions is the problem fixed ?

The issue is fixed in EMOSLIB cycle 000430 and newer.  This is used by:

  • MARS client to be updated in January February 2016 (until then it is available with "mars -t")
  • Metview version 4.6.1 and newer
  • Product generation Disseminated products for IFS cycle 41r2 (to be implemented in Spring 2016)implementation planned for March 2016)

For the Disseminated products from IFS cycle 41r2, the problem is also fixed for interpolations to sub-areas of latitude-longitude grids.

For which cases does the problem still exist ?

The issue has not been fixed for the following cases

  • interpolations from original reduced  or octahedral reduced Gaussian grids to sub-areas of regular Gaussian grids (gridType=regular_gg);
  • interpolations from original reduced  or octahedral reduced Gaussian grids to sub-areas of unrotated regular latitude-longitude grids (gridType=regular_ll);
  • interpolations from original reduced  or octahedral reduced Gaussian grids to sub-areas of rotated latitude-longitude grids (gridType=rotated_ll).

Examples