[MITgcm-support] Rotating tank in Cartesian.

[Martin Losch]

You data files looks OK to me (at first glance).
What do you mean by "wrong in 4 directions", space-time?
Please be more specific, maybe post a figure of your result somewhere online (please don't send large email attachments to this list).

Martin

On Jul 12, 2010, at 10:03 PM, aslavin at mun.ca wrote:

> 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.,
> &
> 
> 
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support