[MITgcm-support] Does the value of the parameter "viscAh" have a limit? (ideal_2D_oce)
赖燕红
yhlai at pku.edu.cn
Tue Jun 23 20:51:59 EDT 2020
Hi,
When I tried to change the parameter "viscAh" from 12E5(its default value) to 12E6 when usind ideal_2D_oce case, for some reason the simulation proceeds extremely slowly, like it only simulates 5000 years after integrated for 48 hours, but normally it only needs a few hours to complete 10 thousand years simulation if we don't change "viscAh" to 12E6, do you know why or does the value of "viscAh" have a limit?
Here's part of my data file:
# Continuous equation parameters
&PARM01
tRef= 18.8, 16.3, 13.1, 10.4, 8.1, 6.0, 4.5,
3.4, 2.7, 2.2, 1.8, 1.5, 1.2, 1.0, 0.9,
sRef= 15*34.7,
cosPower=1.,
viscAr=1.E-3,
viscAh=12.E6,
no_slip_sides=.FALSE.,
no_slip_bottom=.TRUE.,
bottomDragLinear=1.E-3,
#diffK4T=2.E12,
diffKhT=0.E3,
diffKrNrT=1.E-4, 14*1.E-5,
ivdc_kappa=100.,
implicitDiffusion=.TRUE.,
eosType='LINEAR',
tAlpha=2.E-4,
sBeta =0.E-4,
saltStepping=.FALSE.,
gravity=9.81,
rhonil=1035.,
implicitFreeSurface=.TRUE.,
exactConserv=.TRUE., staggerTimeStep=.TRUE., useCDscheme=.TRUE., readBinaryPrec=64, tempAdvScheme=77, usePickupBeforeC54=.TRUE., &# Elliptic solver parameters &PARM02 cg2dMaxIters=1000,#cg2dTargetResidual=1.E-13, cg2dTargetResWunit=1.E-17, &
...
Thanks in advance!
Yanhong
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20200624/a925aee5/attachment.html>
More information about the MITgcm-support
mailing list