Next: 3.15.4 Code Configuration
Up: 3.15 Centennial Time Scale
Previous: 3.15.2 Overview
Contents
Subsections
3.15.3 Discrete Numerical Configuration
The model is configured in hydrostatic form. The domain is discretised with
a uniform grid spacing in latitude and longitude of
, so
that there are ninety grid cells in the and forty in the
direction (Arctic polar regions are not
included in this experiment). Vertically the
model is configured with twenty layers with the following thicknesses
(here the numeric subscript indicates the model level index number, ).
The implicit free surface form of the pressure equation described in Marshall et. al
[39] is employed. A Laplacian operator, , provides viscous
dissipation. Thermal and haline diffusion is also represented by a Laplacian operator.
Wind-stress momentum inputs are added to the momentum equations for both
the zonal flow, and the meridional flow , according to equations
(3.98) and (3.99).
Thermodynamic forcing inputs are added to the equations for
potential temperature, , and salinity, , according to equations
(3.100) and (3.101).
This produces a set of equations solved in this configuration as follows:
|
|
|
(3.100) |
|
|
|
(3.101) |
|
|
0 |
(3.102) |
|
|
|
(3.103) |
|
|
|
(3.104) |
|
|
|
(3.105) |
where and are the and components of the
flow vector . The suffices indicate surface and
interior model levels respectively. As described in
MITgcm Numerical Solution Procedure 2, the time
evolution of potential temperature, , equation is solved prognostically.
The total pressure, , is diagnosed by summing pressure due to surface
elevation and the hydrostatic pressure.
3.15.3.1 Numerical Stability Criteria
The Laplacian dissipation coefficient, , is set to
.
This value is chosen to yield a Munk layer width [1],
|
|
|
(3.106) |
of
km. This is greater than the model
resolution in mid-latitudes , ensuring that the frictional
boundary layer is well resolved.
The model is stepped forward with a
time step
secs. With this time step the stability
parameter to the horizontal Laplacian friction [1]
|
|
|
(3.107) |
evaluates to 0.012, which is well below the 0.3 upper limit
for stability.
The vertical dissipation coefficient, , is set to
. The associated stability limit
|
|
|
(3.108) |
evaluates to
which is again well below
the upper limit.
The values of and are also used for the horizontal ()
and vertical () diffusion coefficients for temperature respectively.
The numerical stability for inertial oscillations
[1]
|
|
|
(3.109) |
evaluates to , which is well below the upper
limit for stability.
The advective CFL [1] for a extreme maximum
horizontal flow
speed of
|
|
|
(3.110) |
evaluates to
. This is well below the stability
limit of 0.5.
The stability parameter for internal gravity waves
[1]
|
|
|
(3.111) |
evaluates to
. This is well below the linear
stability limit of 0.25.
Next: 3.15.4 Code Configuration
Up: 3.15 Centennial Time Scale
Previous: 3.15.2 Overview
Contents
mitgcm-support@dev.mitgcm.org
Copyright © 2002
Massachusetts Institute of Technology |
|
|