[MITgcm-support] Too short time step ?!
Menemenlis, Dimitris (3248)
Dimitris.Menemenlis at jpl.nasa.gov
Wed Jul 7 11:51:14 EDT 2010
For horizontal viscosity, you can also try the Leith scheme, which attempts to choose a viscosity that is just right for damping grid scale noise while maintaining a realistic energy cascade that extends as close as possible to your grid's Nyquist wavenumber:
viscAh=0.,
viscAhGrid=0.,
viscC2Leith=0.,
viscA4=0.,
useFullLeith=.TRUE.
viscC4Leith=2.0,
viscC4Leithd=2.0,
On Jul 6, 2010, at 9:11 PM, Matthew Mazloff wrote:
> Hi Zhang,
>
> You can turn up viscAh and also add some diffusivity to stabalize the
> model run.
>
> For your domain something like
> viscAz =1.E-4,
> viscAh =1.E1,
> viscA4 =1.E9,
> diffKhT=1.E0,
> diffKzT=1.E-5,
> diffKhS=1.E0,
> diffKzS=1.E-5,
> bottomDragQuadratic=1.E-3,
>
> deltaTmom = 1200.0,
> deltaTtracer= 1200.0,
>
> may be better
>
> -Matt
>
>
>
>
> On Jul 6, 2010, at 7:40 PM, 张正光 wrote:
>
>> 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',
>> &
>> <ATT00001.txt>
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
More information about the MITgcm-support
mailing list