[MITgcm-support] centrifugal force

[quentin jamet]

Dear MITgcm support team

I try to simulate the dynamic of an anticyclonic eddy. To correctly
initialize the eddy, I use a three-way balance between Coriolis force,
centrifugal force and horizontal pressure gradient.
The introduction of centrifugal force in the calcul of initial velocities
leads to uncoherent results becaus strong (~O(10-3 - 10-5 m/s) initial
vertical velocities are generated at the begin of the run (i.e.
W.0000000000.001.001.data').

How can I do to correctly initialize my eddy with concidering centrifugal
force?

Thanks You

Quentin Jamet

# Model parameters
# Continuous equation parameters
 &PARM01
#tRef=20.,10.,8.,6.,
#sRef=10.,10.,10.,10.,
 viscAz=1.E-4,
 viscAh=1.0,
 no_slip_sides=.FALSE.,
 no_slip_bottom=.FALSE.,
#diffKhT=4.E2,
#diffKzT=1.E-2,
#rotationPeriod=86400.,
 f0=1.E-4,
 beta=0.E-11,
 tAlpha=2.E-4,
 sBeta =0.,
 gravity=9.81,
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 eosType='LINEAR',
 readBinaryPrec=64,
 &
# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=1000,
 cg2dTargetResidual=1.E-13,
 &
# Time stepping parameters
 &PARM03
 startTime=0.,
#endTime=360000.,
 endTime=180.,
 deltaTmom=180.0,
 deltaTtracer=180.0,
 abEps=0.1,
 pChkptFreq=0.0,
 chkptFreq=8640000.0,
#dumpFreq=1800.,
 dumpFreq=180.,
 monitorFreq=1.,
 &
# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 usingSphericalPolarGrid=.FALSE.,
#ygOrigin=35.,
 delX=150*2.E3,
 delY=150*2.E3,
 delZ=200*4.,
#delZ=500.,500.,500.,500.,
#delZfile='dz.bin',
 &
 &PARM05
 bathyFile='topog.bin',
 hydrogThetaFile='theta2.bin',
#zonalWindFile='tauX_cst.bin',
#meridWindFile='windy.sin_x',
 vVelInitFile='vVelCyclo.bin',
 uVelInitFile='uVelCyclo.bin',
 pSurfInitFile='ini_eta.bin',
 &
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20120229/aa2eb5b6/attachment.htm>