[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