[MITgcm-support] Thermal wind in ocean cause too large velocity?
Chun-Yan Zhou
c.zhou at dundee.ac.uk
Mon May 21 07:38:38 EDT 2012
Based on
The temperature difference in horizontal direction can cause large velocity because of thermal wind.
And I did a test where the initial temperature field changed form 20-13 degree along the X direction, no other forces are applied.
Large velocity is expected, but the maximum velocity is 4m/s, and the minimum is -4m/s, most of the velocity is around 1.5m/s. Is the velocity should be so large?
Any suggestion is appreciated.
chunyan
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
&PARM01
tRef=20*20.,
sRef=20*35.,
viscAh=4.E-2,
viscAz=4.E-2,
no_slip_sides=.FALSE.,
no_slip_bottom=.FALSE.,
diffKhT=4.E-2,
diffKzT=4.E-2,
f0=0.E-4,
beta=0.E-11,
tAlpha=2.0E-4,
sBeta =0.,
tempAdvScheme=77,
gravity=10.,
rhoConst=1000.,
rhoNil=1000.,
heatCapacity_Cp=4000.,
#rigidLid=.TRUE.,
implicitFreeSurface=.TRUE.,
#exactConserv=.TRUE.,
eosType='LINEAR',
nonHydrostatic=.TRUE.,
saltStepping=.FALSE.,
staggerTimeStep=.TRUE.,
&
# Elliptic solver parameters
&PARM02
cg2dMaxIters=1000,
cg2dTargetResidual=1.E-9,
cg3dMaxIters=1000,
cg3dTargetResidual=1.E-9,
&
# Time stepping parameters
&PARM03
nIter0=0,
#endTime=43200.,
nTimeSteps=21600,
deltaT=10.,
abEps=0.1,
pChkptFreq=43200.,
chkptFreq=7200.,
dumpFreq=1800.,
monitorFreq=600.,
monitorSelect=1,
monitorFreq=1.,
&
# Gridding parameters
&PARM04
usingCartesianGrid=.TRUE.,
dXspacing=100.,
dYspacing=100.,
delZ=50*20.,
&
# Input datasets
&PARM05
hydrogThetaFile='T.bin',
&
%%%%%%%%%%%gendata.m
var=zeros(100,100,50);
for i=1:100
var(i,:,:)=20-0.07*(i-1);
end
fid=fopen('T.bin','w','b'); fwrite(fid,var,'real*4'); fclose(fid);
Dear all,
The University of Dundee is a registered Scottish Charity, No: SC015096
More information about the MITgcm-support
mailing list