Yin-Yang grid

When running GEM on a global grid the grid used is a Yin-Yang grid. The Yin-Yang grid is constructed of two slightly overlapping LAM grids. To cover the whole globe, the two LAM grids are organized like the two parts of a tennis or baseball but with to rectangular parts.

Each LAM grid covers 270°x90° and can (or rather should) be enlarged by an overlap respecting the same ratio between X and Y (3 to 1) by using the parameter 'Grd_overlap' (in degrees). The overlap is needed because the two LAM grids will give different answers in the overlap region and the larger the overlap region the easier it is for the solver to find the common solution. The minimum overlap should be 3 to 8 times the grid resolution (delta Y degrees).

If a Grd_overlap (defined by degrees) is given, then each grid is expanded as such:

    (270° + (6*Grd_overlap)) x (90°+ (2*Grd_overlap))

For example: if Grd_overlap=2° then, each LAM grid covers 282°x94°.

When written in an RPN file the two LAM grid are stored as one grid, with one LAM "on top" of the other - see image below.

If need be the Yin-Yang grid can get split into the two LAM grids with the command: yydecode.Abs
And the two LAM grids can get put back together with the command: yyencode.Abs

For more information about the Yin-Yang grid in GEM have a look at the following publication:

    Qaddouri, A., & Lee, V. (2011). The Canadian Global Environmental Multiscale model on the Yin‐Yang grid system. Quarterly Journal of the Royal Meteorological Society. Royal Meteorological Society (Great Britain), 137(660), 1913–1926. https://doi.org/10.1002/qj.873

Example for Yin-Yang grid centered over Africa at :    Grd_xlat1 = 24.5   , Grd_xlon1 = 22.

An overlap option is provided because the two LAM grids will give different answers in the overlap region and the larger the overlap region, the easier it is for the solver to find the common solution. Note that the overlap is 3 times larger on the left and right border than at the top and bottom of each of the two LAM grids! This is done to keep the aspect ratio of each LAM grid at 3:1 (minus 1 point).

When projected on a regular latlon grid the two LAM grids above would look like this: