[MITgcm-support] Rotating tank in Cartesian.

[aslavin at mun.ca]

Hello Martin,

Thank you very much for your interest to my problem. 
Here is my bed profile:
http://img135.imageshack.us/img135/8871/bed.png
Here is initial temperatures:
http://img256.imageshack.us/img256/5267/10826309.png
In cylindrical grid I got 4 vortexes (as it must be), but in Cartesian I
got this after 10 sec (I also added 1 more in center to check what happens in
center):  http://img256.imageshack.us/img256/4277/10sec.png
There are also some unnatural disturbances on the left side of the tank.
And after 100 sec (without center vortex):
http://img135.imageshack.us/img135/1855/100sec.png

Alexander



Quoting Martin Losch <Martin.Losch at awi.de>:

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