make testwhich runs the
testXY
and testXZ
cases in turn. These cases can
be individually run with the commands
make testXYand
make testXZCentral differencing for advection of momentum (
suntans.dat: nonlinear = 2
) is employed and nonhydrostatic pressure
used for the testXZ
case (suntans.dat: nonhydrostatic = 1
). Small
time steps of 0.001
and 0.005
are employed to resolve
eddies over a total time of 40
, corresponding to a rough approximation
for steady state. Viscosity controls the Reynolds number and is set to
yield (suntans.dat: nu = 0.0003125
and
suntans.dat: nu_H = 0.0003125
) for comparisons with
Zang et al. [6]. Results from the simulation can be compared with
these literature data values via the Matlab scripts CompareResultXY.m
and
CompareResultXZ.m
. A more resolved case than the examples (13,500 cells)
is shown in Figure 13, where is the square domain's length,
is lid velocity, and and are the zonal and meridional velocities,
respectively.
0.8
