[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