2.11.4. Horizontal grid

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.

hgrid-abcd

Figure 2.8 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\).

Figure 2.8 (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 dxG, dyG and 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.

Figure 2.8 (b) shows the vorticity cell. The length of the southern edge, \(\Delta x_c\), western edge, \(\Delta y_c\) and surface area, \(A_\zeta\), presented in the vertical are stored in arrays dxC, dyC and rAz. The “z” 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.

Figure 2.8 (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 dxV, dyF and 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.

Figure 2.8 (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 dxF, dyU and 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 INI_CARTESIAN_GRID , INI_SPHERICAL_POLAR_GRID , INI_CURVILINEAR_GRID

\(A_c , A_\zeta , A_w , A_s\) : rA, rAz, rAw, rAs ( GRID.h )
\(\Delta x_g , \Delta y_g\) : dxG, dyG ( GRID.h )
\(\Delta x_c , \Delta y_c\) : dxC, dyC ( GRID.h )
\(\Delta x_f , \Delta y_f\) : dxF, dyF ( GRID.h )
\(\Delta x_v , \Delta y_u\) : dxV, dyU ( GRID.h )

2.11.4.1. Reciprocals of horizontal grid descriptors

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 GRID.h but they are all only 2-D descriptors.

S/R INI_MASKS_ETC

\(A_c^{-1} , A_\zeta^{-1} , A_w^{-1} , A_s^{-1}\) : recip_rA, recip_rAz, recip_rAw, recip_rAs ( GRID.h )
\(\Delta x_g^{-1} , \Delta y_g^{-1}\) : recip_dxG, recip_dyG ( GRID.h )
\(\Delta x_c^{-1} , \Delta y_c^{-1}\) : recip_dxC, recip_dyC ( GRID.h )
\(\Delta x_f^{-1} , \Delta y_f^{-1}\) : recip_dxF, recip_dyF ( GRID.h )
\(\Delta x_v^{-1} , \Delta y_u^{-1}\) : recip_dxV, recip_dyU ( GRID.h )

2.11.4.2. Cartesian coordinates

Cartesian coordinates are selected when the logical flag usingCartesianGrid in namelist PARM04 is set to true. The grid spacing can be set to uniform via scalars dXspacing and dYspacing in namelist PARM04 or to variable resolution by the vectors DELX and DELY. Units are normally meters. Non-dimensional coordinates can be used by interpreting the gravitational constant as the Rayleigh number.

2.11.4.3. Spherical-polar coordinates

Spherical coordinates are selected when the logical flag usingSphericalPolarGrid in namelist PARM04 is set to true. The grid spacing can be set to uniform via scalars dXspacing and dYspacing in namelist PARM04 or to variable resolution by the vectors DELX and DELY. Units of these namelist variables are alway degrees. The horizontal grid descriptors are calculated from these namelist variables have units of meters.

2.11.4.4. Curvilinear coordinates

Curvilinear coordinates are selected when the logical flag 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.