Next: 6.3.4 Variable
Up: 6.3 Gent/McWiliams/Redi SGS Eddy
Previous: 6.3.2 GM parameterization
Contents
Griffies notes that the discretisation of bolus velocities involves
multiple layers of differencing and interpolation that potentially
lead to noisy fields and computational modes. He pointed out that the
bolus flux can be re-written in terms of a non-divergent flux and a
skew-flux:
The first vector is non-divergent and thus has no effect on the tracer
field and can be dropped. The remaining flux can be written:
|
(6.10) |
where
|
(6.11) |
is an anti-symmetric tensor.
This formulation of the GM parameterization involves fewer derivatives
than the original and also involves only terms that already appear in
the Redi mixing scheme. Indeed, a somewhat fortunate cancellation
becomes apparent when we use the GM parameterization in conjunction
with the Redi isoneutral mixing scheme:
|
(6.12) |
In the instance that
then
|
(6.13) |
which differs from the variable Laplacian diffusion tensor by only
two non-zero elements in the -row.
Next: 6.3.4 Variable
Up: 6.3 Gent/McWiliams/Redi SGS Eddy
Previous: 6.3.2 GM parameterization
Contents
mitgcm-support@dev.mitgcm.org
Copyright © 2002
Massachusetts Institute of Technology |
|
|