Next: spm.F
Up: The Proudman Oceanographic Laboratory
Previous: setopenbc.F
Subsections
sigmaset.F
subroutine sigmaset
This subroutine calculates the vertical coordinate sigma
defined by
, where is the
equilibrium water depth, is the sea surface elevation
and is the cartesian vertical coordinate. The water column
at each grid point is divided into equallyspaced
levels between the sea bed () and the sea
surface (). On these levels, the vertical velocity,
eddy viscosity and diffusivity and turbulent kinetic energy are
defined ( points in the figure). The horizontal
u and v currents and scalars, such as temperature and salinity,
are defined at the midpoints of these levels (u points
in the figure) and have indices
, where levels
and lie below the sea bed and above the sea
surface, respectively, to facilitate the use of flux boundary
conditions.
Figure 5:
Positions of levels for and u points.

 none
 k index of vertical level.
 dsc, sds variables for use in calculating the
sigma levels.
 sigo, sig, ds, dsu, dsup
 none
 Calculate values at levels where the vertical
velocity, eddy viscosity and diffusivity and turbulent kinetic
energy are defined  array sigo.
 Calculate values at levels where the horizontal
u and v currents and scalars such as temperature and salinity
are defined  array sig  midway between the sigo levels.
 none
 b3dgrid
 lpgrad  true if using pressure gradient
calculations on a horizontal plane.
 zvarybc  true if using depthvarying boundary
condition.
 SCOORD  use horizontally varying vertical
coordinate.
 DEBUG, DEBUG_B3DRUN, DEBUG_MODEL  output variables
for debugging.
 none
Next: spm.F
Up: The Proudman Oceanographic Laboratory
Previous: setopenbc.F
The AMMP Project
20050420