[MITgcm-support] Can not run a regional model for long time
Sabri MUTLU (MAM-ÇTÜE)
sabri.mutlu at tubitak.gov.tr
Sun Sep 18 09:08:19 EDT 2022
Hi all,
As a newbie in this field, I need help for the following problem, which I have been trying to deal with for a long time.
I built a regional model for my thesis study. When I try to run it for 12 or more days, it gives "NaN" values after 8 or 9 days later in the model.
Could you please examine the following model configuration and give me an idea for possible solution?
(I had to choose Orlanski method since don't have any boundary condition to feed the model.)
Thank you in advance.
Sabri
&PARM01
viscAh=2.E-2,
viscAz=1.5E-5,
no_slip_sides=.TRUE.,
no_slip_bottom=.TRUE.,
diffKhT=2.E-2,
diffKhS=2.E-2,
diffKzT=2.E-5,
diffKzS=2.E-5,
f0=0.948E-4,
implicitFreeSurface=.TRUE.,
exactConserv=.TRUE.,
implicitDiffusion=.TRUE.,
implicitViscosity=.TRUE.,
eosType = 'JMD95Z',
tRef=22.2803,22.3383,22.3646,22.4220,22.5469,22.5081,21.0406,19.6075,18.3619,16.7902,15.8310,15.4016,15.3625,15.3935,15.4334,15.4465,15.4508,15.4710,15.5226,
15.5428,15.5485,15.5589,15.5829,15.5860,15.5859,15.5767,15.5572,15.5477,15.5357,15.5068,15.4782,15.4640,15.4337,15.3816,15.3480,15.2998,15.2423,15.2127,
15.1271,15.0639,15.0086,14.9722,14.9229,14.8937,14.8470,14.8073,14.7838,14.7677,14.7392,14.7210,14.7037,14.7036,14.7035,14.7034,14.7033,
sRef=24.3775,24.4042,24.4730,24.5782,24.7042,24.8407,25.8279,26.4026,28.0163,30.1542,31.7531,33.9951,37.0132,37.5629,37.8303,37.9592,38.0011,38.1433,38.2744,
38.3184,38.3548,38.4269,38.5535,38.6043,38.6738,38.7411,38.7479,38.7496,38.7510,38.7512,38.7469,38.7529,38.7569,38.7644,38.7666,38.7683,38.7698,38.7696,
38.7635,38.7563,38.7513,38.7476,38.7427,38.7399,38.7354,38.7326,38.7305,38.7280,38.7254,38.7236,38.7220,38.7220,38.7220,38.7220,38.7220,
staggerTimeStep = .TRUE.,
tempAdvScheme = 33,
saltAdvScheme = 33,
useSingleCpuIO=.TRUE.,
bottomDragQuadratic = 0.02,
select_rStar = 0,
&
&PARM02
cg2dMaxIters=300,
cg2dTargetResidual=1.E-13,
cg3dMaxIters=20,
cg3dTargetResidual=1.E-8,
&
---------
&PARM03
startTime = 0.,
endTime = 1036800.,
deltaT=5,
dumpFreq=10800.,
monitorFreq=3600.,
monitorSelect=2,
periodicExternalForcing=.TRUE.,
externForcingPeriod=10800.,
externForcingCycle=1036800.,
&
&PARM04
usingSphericalPolarGrid=.TRUE.,
delz = 32*2,22*8,10,
delX = 528*0.0013,
delY = 164*0.001,
ygOrigin=40.656,
xgOrigin=29.2612,
&
&PARM05
bathyFile='bathy_v6_528-164.bin',
zonalWindFile='windU_2309-0510.bin',
meridWindFile='windV_2309-0510.bin',
&
# Open-boundaries
&OBCS_PARM01
OB_Iwest=164*1,
useOrlanskiWest=.TRUE.,
useOBCSbalance=.TRUE.,
useOBCSsponge=.FALSE.,
OBCS_balanceFacW = -1.,
&
# Packages
&PACKAGES
useOBCS=.TRUE.,
useKPP =.TRUE.,
&
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20220918/53dda429/attachment.html>
More information about the MITgcm-support
mailing list