[MITgcm-support] Too short time step ?!

张正光 forget53 at gmail.com
Tue Jul 6 22:40:08 EDT 2010


Dear mitgcm-support team,
        Hi, I'm Zhang.
        I'm modeling the overflow, now I'm running a idealized case, but I
have to set my time step very short (dt=120s) in order to run my case,
otherwise, it will go crazy or shut down.
        I'v checked the "Numerical Stability Criteria", which told me I can
use a much longer time step.
        I'm confused, and my case need to be integrated for 30-50year, so
use dt=120s will cost too much computer time. What can I do to use a longer
time step?
        My grid size: dx=25km, dy=50km, dz=100m (uniform grid)
        and I think give you my "data" file maybe better:

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef=40*20.,
 sRef=40*35.,
 viscA4=0.0E4,
 viscAh=1.E-2,
 viscAz=1.E-3,
 no_slip_sides=.TRUE.,
 no_slip_bottom=.TRUE.,
 diffK4T=0.E4,
 diffKhT=0.E-2,
 diffKzT=0.E-3,
 diffK4S=0.E4,
 diffKhS=0.E-2,
 diffKzS=0.E-5,
 f0=0.e-5,
 beta=0.E-11,
 rhonil=1035.,
 rhoConstFresh=1000.,
 eosType = 'JMD95Z',
 rigidLid=.FALSE.,
 implicitFreeSurface=.TRUE.,
 hFacMin=0.05,
 nonHydrostatic=.TRUE.,
 readBinaryPrec=64,
#- not safe to use globalFiles in multi-processors runs
#globalFiles=.TRUE.,
 bottomDragLinear=0.E-4,
 bottomDragQuadratic=0.003,
 tempAdvScheme=33,
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=300,
 cg2dTargetResidual=1.E-13,
 cg3dMaxIters=400,
 cg3dTargetResidual=1.E-8,
 &

# Time stepping parameters
 &PARM03
 niter0=0,
 nTimeSteps=10000,
#nTimeSteps=20000,
 deltaT=600.0,
 abEps=0.01,
 pChkptFreq=120000000.0,
 chkptFreq=600000.0,
 dumpFreq=600000.0,
#dumpFreq=120000.0,
 taveFreq=600000.0,
 tauThetaClimRelax=100000.0,
 monitorFreq=6000.,
# staggerTimeStep=.TRUE.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 usingSphericalPolarGrid=.FALSE.,
 delXfile='dx.bin',
 delYfile='dy.bin',
# dYspacing=200.,
 delZ=40*100.0,
 &

# Input datasets
 &PARM05
 bathyFile='topog.slope',
 hydrogThetaFile='T.init',
 surfQfile='Qnet.forcing',
 thetaClimFile='SST.bin',
 diffKrFile='diffK.bin'
#hydrogThetaFile='T.pickup',
#hydrogSaltFile='S.pikup',
#uVelInitFile='U.pickup',
#vVelInitFile='V.pickup',
#pSurfInitFile='Eta.pickup',
 &
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20100707/20ca82e6/attachment.htm>


More information about the MITgcm-support mailing list