[MITgcm-support] Rotating tank in Cartesian.

[Martin Losch]

Alexander,

are you sure,  you are plotting the right thing? You plots look very much like you have accidentally exchanged U and V-velocities?
u=rdmds(U,Inf);
v=rdmds(V,Inf);
uc = .5*(u([end 1:end],:,:)+u);
vc = .5*(v(:,[end 1:end],:)+v); 
quiver(uc,vc)

I don't know where the velocities near the tank edge come from, neither the "stripes" in the vector field? If you give me your code and input directories (preferably as tar files to download), I can try to reproduce the result.

Martin



On Jul 16, 2010, at 5:41 PM, aslavin at mun.ca wrote:

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