The model domain is decomposed into tiles and within each tile a quasi-regular grid is used. A tile is the basic unit of domain decomposition for parallelization but may be used whether parallelized or not; see section [sec:domain_decomposition] for more details. Although the tiles may be patched together in an unstructured manner (i.e. irregular or non-tessilating pattern), the interior of tiles is a structured grid of quadrilateral cells. The horizontal coordinate system is orthogonal curvilinear meaning we can not necessarily treat the two horizontal directions as separable. Instead, each cell in the horizontal grid is described by the length of it’s sides and it’s area.
The grid information is quite general and describes any of the available coordinates systems, cartesian, spherical-polar or curvilinear. All that is necessary to distinguish between the coordinate systems is to initialize the grid data (descriptors) appropriately.
In the following, we refer to the orientation of quantities on the computational grid using geographic terminology such as points of the compass. This is purely for convenience but should not be confused with the actual geographic orientation of model quantities.
Staggering of horizontal grid descriptors (lengths and areas). The grid lines indicate the tracer cell boundaries and are the reference grid for all panels. a) The area of a tracer cell, A_c, is bordered by the lengths \Delta x_g and \Delta y_g. b) The area of a vorticity cell, A_\zeta, is bordered by the lengths \Delta x_c and \Delta y_c. c) The area of a u cell, A_w, is bordered by the lengths \Delta x_v and \Delta y_f. d) The area of a v cell, A_s, is bordered by the lengths \Delta x_f and \Delta y_u.
:numref:`hgrid-abcd` (a) shows the tracer cell (synonymous with the continuity cell). The length of the southern edge, \Delta x_g, western edge, \Delta y_g and surface area, A_c, presented in the vertical are stored in arrays :varlink:`dxG`, :varlink:`dyG` and :varlink:`rA`. The “g” suffix indicates that the lengths are along the defining grid boundaries. The “c” suffix associates the quantity with the cell centers. The quantities are staggered in space and the indexing is such that dxG(i,j) is positioned to the south of rA(i,j) and dyG(i,j) positioned to the west.
:numref:`hgrid-abcd` (b) shows the vorticity cell. The length of the northern edge, \Delta x_c, eastern edge, \Delta y_c and surface area, A_\zeta, presented in the vertical are stored in arrays :varlink:`dxC`, :varlink:`dyC` and :varlink:`rAz`. The “c” suffix indicates that the lengths are measured between the cell centers and the “\zeta” suffix associates points with the vorticity points. The quantities are staggered in space and the indexing is such that dxC(i,j) is positioned to the north of rAz(i,j) and dyC(i,j) positioned to the east.
:numref:`hgrid-abcd` (c) shows the “u” or western (w) cell. The length of the southern edge, \Delta x_v, eastern edge, \Delta y_f and surface area, A_w, presented in the vertical are stored in arrays :varlink:`dxV`, :varlink:`dyF` and :varlink:`rAw`. The “v” suffix indicates that the length is measured between the v-points, the “f” suffix indicates that the length is measured between the (tracer) cell faces and the “w” suffix associates points with the u-points (w stands for west). The quantities are staggered in space and the indexing is such that dxV(i,j) is positioned to the south of rAw(i,j) and dyF(i,j) positioned to the east.
:numref:`hgrid-abcd` (d) shows the “v” or southern (s) cell. The length of the northern edge, \Delta x_f, western edge, \Delta y_u and surface area, A_s, presented in the vertical are stored in arrays :varlink:`dxF`, :varlink:`dyU` and :varlink:`rAs`. The “u” suffix indicates that the length is measured between the u-points, the “f” suffix indicates that the length is measured between the (tracer) cell faces and the “s” suffix associates points with the v-points (s stands for south). The quantities are staggered in space and the indexing is such that dxF(i,j) is positioned to the north of rAs(i,j) and dyU(i,j) positioned to the west.
S/R :filelink:`INI_CARTESIAN_GRID <model/src/ini_cartesian_grid.F>` , :filelink:`INI_SPHERICAL_POLAR_GRID <model/src/ini_spherical_polar_grid.F>` , :filelink:`INI_CURVILINEAR_GRID <model/src/ini_curvilinear_grid.F>`
A_c , A_\zeta , A_w , A_s : :varlink:`rA`, :varlink:`rAz`, :varlink:`rAw`, :varlink:`rAs` ( :filelink:`GRID.h <model/inc/GRID.h>` )
Lengths and areas appear in the denominator of expressions as much as in the numerator. For efficiency and portability, we pre-calculate the reciprocal of the horizontal grid quantities so that in-line divisions can be avoided.
For each grid descriptor (array) there is a reciprocal named using the prefix
recip_
. This doubles the amount of storage in :filelink:`GRID.h
<model/inc/GRID.h>` but they are all only 2-D descriptors.
S/R :filelink:`INI_MASKS_ETC <model/src/ini_masks_etc.F>`
A_c^{-1} , A_\zeta^{-1} , A_w^{-1} , A_s^{-1} : :varlink:`recip_rA`, :varlink:`recip_rAz`, :varlink:`recip_rAw`, :varlink:`recip_rAs` ( :filelink:`GRID.h <model/inc/GRID.h>` )\Delta x_g^{-1} , \Delta y_g^{-1} : :varlink:`recip_dxG`, :varlink:`recip_dyG` ( :filelink:`GRID.h <model/inc/GRID.h>` )\Delta x_c^{-1} , \Delta y_c^{-1} : :varlink:`recip_dxC`, :varlink:`recip_dyC` ( :filelink:`GRID.h <model/inc/GRID.h>` )\Delta x_f^{-1} , \Delta y_f^{-1} : :varlink:`recip_dxF`, :varlink:`recip_dyF` ( :filelink:`GRID.h <model/inc/GRID.h>` )\Delta x_v^{-1} , \Delta y_u^{-1} : :varlink:`recip_dxV`, :varlink:`recip_dyU` ( :filelink:`GRID.h <model/inc/GRID.h>` )
Cartesian coordinates are selected when the logical flag
:varlink:`usingCartesianGrid` in namelist PARM04
is set to true. The grid
spacing can be set to uniform via scalars :varlink:`dXspacing` and
:varlink:`dYspacing` in namelist PARM04
or to variable resolution by the
vectors :varlink:`DELX` and :varlink:`DELY`. Units are normally meters.
Non-dimensional coordinates can be used by interpreting the
gravitational constant as the Rayleigh number.
Spherical coordinates are selected when the logical flag
:varlink:`usingSphericalPolarGrid` in namelist PARM04
is set to true. The
grid spacing can be set to uniform via scalars :varlink:`dXspacing` and
:varlink:`dYspacing` in namelist PARM04
or to variable resolution by the
vectors :varlink:`DELX` and :varlink:`DELY`. Units of these namelist variables
are alway degrees. The horizontal grid descriptors are calculated from these
namelist variables have units of meters.
Curvilinear coordinates are selected when the logical flag
:varlink:`usingCurvilinearGrid` in namelist PARM04
is set to true. The grid
spacing can not be set via the namelist. Instead, the grid descriptors
are read from data files, one for each descriptor. As for other grids,
the horizontal grid descriptors have units of meters.