Set up a grid

GEM offers the possibilities of running global uniform (up to version 4), global YinYang (starting version 4) and limited area (LAM) grids. Have a look at the grid types.

The grid is specified in the file 'gem_settings.nml' in the namelist 'grid'. This GEM grid specification parameters link explains how to define the parameters. Alternatively, see below.

You have to set variables for the

• grid type
• grid rotation (optional)
• number of grid points and grid size
• location of limited area (just for LAM grids)
• FFT criteria

Grid type

First thing to do is to set 'Grd_typ_S' to the grid type you want:

Grid rotation

Next specify whether your grid is to be rotated or not (if not just skip this paragraph).
Four variables are used to define the rotation of the grid. They define two points:

• The first point will be the new center of your grid (at the grid(!) equator and the grid(!) 180° longitude).
• The second point is any other point on the grid(!) equator east of the center and will therefore set the rotation of the grid equator wrt the real equator in the first point.
  

Both points are specified in latitudes and longitudes of the real globe. For example:

To create a grid that is shifted but not rotated in the center point Grd_xlat2 and Grd_xlon2 should get set the following way:

   Grd_xlat2 = 0.,  
   Grd_xlon2 = Grd_xlon1 + 90.       =>  You need to set it to the value 'Grd_xlon1 + 90.', do not actually put 'Grd_xlon1' here!!!

Number of grid points and grid size

Global uniform grid
It is enough to specify the total number of points. The model will then calculate size of the grid boxes automatically.

Global YinYang grid

It is enough to specify the core number of points in y-direction as well as the "overlap" in degrees.

LAM grid
Specify the total number of points and the size of the grid boxes in degrees.
If at all possible, insure that the total number of grid points in x-direction matches the FFT criteria.

Location of limited area (just for LAM grids)

The only thing left now when setting up a LAM grid is to specify where the limited area is located on the rotated grid.

Click here to find out how to set blending and pilot area (halo).

FFT criteria

To fulfill the FFT criteria the number of grid points in x-direction, excluding the pilot zone, must be a multiple of 2, 3, and 5. Use the findfft GEMCLIM script to find the possible numbers of grid points fulfilling the FFT criteria.
Again you need to insure that the model environment is set.

Then you can use findfft:

This command will then print all "good" values you can use for 'Grd_ni' and 'Grd_nj' (can be different) between the range you specified with 'min' and 'max'.

grille

The command 'grille' reads the namelist 'grid' from an input file and :

    1. creates 3 files with the corresponding tictacs of the:
         a) model grid (tape1_model) - this file is needed as input to GenPhysX!!!
         b) core grid (tape1_core)
         c) free_grid (tape1_free)

    2. opens 'xrec' with the domain specified in 'grid' - only if requested - see table below.

These are the keys that can get used with the command:

After having specified a grid in the file 'gem_settings.nml' you can have a look at the grid you specified using the grille script.

In the directory in which you have your file 'gem_settings.nml' simply type :

  grille -xrec

This command will open 'xrec' with three fields of your grid: