[MITgcm-support] Rotating tank in Cartesian.
aslavin at mun.ca
aslavin at mun.ca
Mon Jul 12 16:03:12 EDT 2010
Dear colleagues,
My simulations based on Rotating tank example. It works well in Cylindrical
Grid. I decided to convert into Cartesian Grid (usingCartesianGrid=.TRUE.). But
I got problem. There is no problem with grid area itself, the problem with
results obtained. They are not same as in Cylindrical and they are mistaken.
For example if I must have a vortex (after initial hot temperature circle) I
got circle, but direction of the velocities are wrong in 4 directions.
Maybe someone can help me with my problem? Thank you in advance.
Alexander
Below my data file:
# Continuous equation parameters
&PARM01
tRef=20*20.0,
sRef=20*35.0,
viscAh=5.0E-6,
viscAz=5.0E-6,
no_slip_sides=.FALSE.,
no_slip_bottom=.FALSE.,
diffKhT=2.5E-6,
diffKzT=2.5E-6,
diffKhS=1.0E-6,
diffKzS=1.0E-6,
#beta=0.E-11
f0=2,
eosType='LINEAR',
sBeta =0.,
gravity=9.81,
rhoConst=1000.0,
rhoNil=1000.0,
#heatCapacity_Cp=3900.0,
rigidLid=.TRUE.,
implicitFreeSurface=.FALSE.,
nonHydrostatic=.TRUE.,
readBinaryPrec=32,
&
# Elliptic solver parameters
&PARM02
cg2dMaxIters=1000,
cg2dTargetResidual=1.E-7,
cg3dMaxIters=10,
cg3dTargetResidual=1.E-9,
&
# Time stepping parameters
&PARM03
nIter0=0,
nTimeSteps=100,
#nTimeSteps=36000000,
deltaT=0.1,
abEps=0.1,
pChkptFreq=0.0,
#chkptFreq=2.0,
dumpFreq=0.0,
monitorSelect=2,
monitorFreq=0.1,
&
# Gridding parameters
&PARM04
usingCartesianGrid=.TRUE.,
#usingCylindricalGrid=.FALSE.,
#usingCurvilinearGrid=.FALSE.,
delX=120*0.01,
delY=120*0.01,
delZ=20*0.01,
#ygOrigin=0.07,
&
# Input datasets
&PARM05
hydrogThetaFile='TVorCART.bin',
bathyFile='BEDCart.bin',
#vVelInitFile = 'Vstr.bin',
#tCylIn = 20.,
#tCylOut = 20.,
&
More information about the MITgcm-support
mailing list