The current version of CANDIE uses the Arakawa A grid formulation. It
was developed by Dan Wright. The main reason for choosing the A grid
CANDIE as our standard model code is based on the finding that the
C-grid model does not reproduce well the bottom trapped baroclinic
topographic Rossby waves in a periodic channel with a uniformly
slopping bottom. The governing equations and time discretization used
on the A grid CANDIE are identical to the C grid (see Sections A.1 and
A.2). The main differences between the A-grid and C-grid formations
are solely associated with the spatial discretization of horizontal
velocity components u and v. On the A grid u and v are
defined at the center of each cell (i.e., the p points, see Figure
1). On the C grid, on the other hand, u and v are defined on the
cell faces normal to the x and y directions, respectively.
The A grid formulation overcomes the difficulty of the C-grid in estimating the Coriolis terms since the state variables u and v are defined at the same location. The conservation equations for heat, salt and momentum, however, each require knowledge of the velocity components on the cell faces, i.e. on the C grid. Thus interpolations of the horizontal velocity components on the A grid onto the C grid are still required. The trial velocity components on the C grid can be calculated from the A-grid velocity components. Using two-point averaging, for example, we have:
The surface pressure correction can then be calculated from the depth integrated horizontal divergence of over each water column. The barotropic velocity corrections can readily be calculated from through (33) and (34). Interpolations of back to the p-points yield the total horizontal velocity components on the A grid at the end of the time step: