The vorticity is defined at an -point ( corner point) as follows:
The horizontal divergence is defined at a -point. It is given by:
Note that although the vorticity has the same discrete expression in - and -coordinates, its physical meaning is not identical. is a pseudo vorticity along -surfaces (only pseudo because are still defined along geopotential surfaces, but are not necessarily defined at the same depth).
The vorticity and divergence at the before step are used in the computation of the horizontal diffusion of momentum. Note that because they have been calculated prior to the Asselin filtering of the before velocities, the before vorticity and divergence arrays must be included in the restart file to ensure perfect restartability. The vorticity and divergence at the now time step are used for the computation of the nonlinear advection and of the vertical velocity respectively.
The sea surface height is given by :
The vertical velocity is computed by an upward integration of the horizontal divergence starting at the bottom, taking into account the change of the thickness of the levels :
In the case of a non-linear free surface (key_ vvl), the top vertical velocity is , as changes in the divergence of the barotropic transport are absorbed into the change of the level thicknesses, re-orientated downward. In the case of a linear free surface, the time derivative in (6.4) disappears. The upper boundary condition applies at a fixed level . The top vertical velocity is thus equal to the divergence of the barotropic transport ( the first term in the right-hand-side of (6.3)).
Note also that whereas the vertical velocity has the same discrete expression in - and -coordinates, its physical meaning is not the same: in the second case, is the velocity normal to the -surfaces. Note also that the -axis is re-orientated downwards in the FORTRAN code compared to the indexing used in the semi-discrete equations such as (6.4) (see §4.1.3).
Gurvan Madec and the NEMO Team
NEMO European Consortium2017-02-17