|
|
|
Next: 6.3.6 Tapering: Danabasoglu and
Up: 6.3 Gent/McWiliams/Redi SGS Eddy
Previous: 6.3.4 Variable
Contents
Subsections
Experience with the GFDL model showed that the GM scheme has to be
matched to the convective parameterization. This was originally
expressed in connection with the introduction of the KPP boundary
layer scheme (Large et al., 97) but in fact, as subsequent experience
with the MIT model has found, is necessary for any convective
parameterization.
Figure 6.1:
Taper functions used in GKW99 and DM95.
|
Figure 6.2:
Effective slope as a function of ``true'' slope using Cox
slope clipping, GKW91 limiting and DM95 limiting.
|
Deep convection sites and the mixed layer are indicated by
homogenized, unstable or nearly unstable stratification. The slopes in
such regions can be either infinite, very large with a sign reversal
or simply very large. From a numerical point of view, large slopes
lead to large variations in the tensor elements (implying large bolus
flow) and can be numerically unstable. This was first recognized by
Cox, 1987, who implemented ``slope clipping'' in the isopycnal mixing
tensor. Here, the slope magnitude is simply restricted by an upper
limit:
|
|
|
(6.14) |
|
|
|
(6.15) |
|
|
|
(6.16) |
|
|
|
(6.17) |
Notice that this algorithm assumes stable stratification through the
``min'' function. In the case where the fluid is well stratified (
) then the slopes evaluate to:
|
(6.18) |
while in the limited regions (
) the slopes become:
|
(6.19) |
so that the slope magnitude is limited
.
The slope clipping scheme is activated in the model by setting GM_taper_scheme = 'clipping' in data.gmredi.
Even using slope clipping, it is normally the case that the vertical
diffusion term (with coefficient
) is large and must be time-stepped using an
implicit procedure (see section on discretisation and code later).
Fig. 6.3 shows the mixed layer depth resulting from
a) using the GM scheme with clipping and b) no GM scheme (horizontal
diffusion). The classic result of dramatically reduced mixed layers is
evident. Indeed, the deep convection sites to just one or two points
each and are much shallower than we might prefer. This, it turns out,
is due to the over zealous re-stratification due to the bolus transport
parameterization. Limiting the slopes also breaks the adiabatic nature
of the GM/Redi parameterization, re-introducing diabatic fluxes in
regions where the limiting is in effect.
The tapering scheme used in Gerdes et al., 1999, ([42])
addressed two issues with the clipping method: the introduction of
large vertical fluxes in addition to convective adjustment fluxes is
avoided by tapering the GM/Redi slopes back to zero in
low-stratification regions; the adjustment of slopes is replaced by a
tapering of the entire GM/Redi tensor. This means the direction of
fluxes is unaffected as the amplitude is scaled.
The scheme inserts a tapering function, , in front of the
GM/Redi tensor:
|
(6.20) |
where is the maximum slope you want allowed. Where the
slopes,
then
and the tensor is un-tapered
but where
then scales down the tensor so
that the effective vertical diffusivity term
.
The GKW tapering scheme is activated in the model by setting GM_taper_scheme = 'gkw91' in data.gmredi.
Next: 6.3.6 Tapering: Danabasoglu and
Up: 6.3 Gent/McWiliams/Redi SGS Eddy
Previous: 6.3.4 Variable
Contents
mitgcm-support@dev.mitgcm.org
Copyright © 2002
Massachusetts Institute of Technology |
|
|