[MITgcm-support] Problem with Global Ocean 1deg run

Abhisek Chakraborty abhisek.sac at gmail.com
Mon Jun 11 11:41:38 EDT 2012


Dear All,

I am trying to run MITgcm for (1deg*1deg) Global oceans with 32
verticals (360*160*32). However, at the very beginning I am getting an
error
"ABNORMAL END: S/R MON_SOLUTION, stops due to EXTREME Pot.Temp" with
an error message:

"(PID.TID 0001.0001) SOLUTION IS HEADING OUT OF BOUNDS: tMin,tMax=
-6.844E+02  1.603E+03
(PID.TID 0001.0001)   exceeds allowed range (monSolutionMaxRange=  1.000E+03)
(PID.TID 0001.0001) MON_SOLUTION: STOPPING CALCULATION at Iter=         4
(PID.TID 0001.0001) *** ERROR *** S/R ALL_PROC_DIE: ending the run"

My "data" file looks like:
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef = 32*20.,
 sRef = 32*35.,
 viscAr=1.E-3,
 viscAh=1.E4,
 diffKhT=1.E2,
 diffKrT=1.E-5,
 diffKhS=1.E2,
 diffKrS=1.E-5,
 rhonil=1035.,
 rhoConstFresh=1000.,
 eosType = 'JMD95Z',
 ivdc_kappa=100.,
 implicitDiffusion=.TRUE.,
 allowFreezing=.TRUE.,
 exactConserv=.TRUE.,
 useRealFreshWaterFlux=.TRUE.,
 useCDscheme=.TRUE.,
# turn on looped cells
 hFacMin=.05,
 hFacMindr=50.,
# set precision of data files
 readBinaryPrec=32,
 &

# Elliptic solver parameters
 &PARM02
 cg2dMaxIters=500,
 cg2dTargetResidual=1.E-13,
 &

# Time stepping parameters
 &PARM03
 nIter0=       0,
 nTimeSteps = 40,
# 100 years of integration will yield a reasonable flow field
# startTime =          0.,
# endTime   = 3110400000.,
 deltaTmom = 3600.,
 tauCD =     321428.,
 deltaTtracer= 86400.,
 deltaTClock = 86400.,
 deltaTfreesurf= 86400.,
 abEps = 0.1,
 pChkptFreq= 1728000.,
 dumpFreq=   864000.,
 taveFreq=   864000.,
 monitorFreq=1.,
# 2 months restoring timescale for temperature
 tauThetaClimRelax=  5184000.,
# 6 months restoring timescale for salinity
 tauSaltClimRelax = 15552000.,
 periodicExternalForcing=.TRUE.,
 externForcingPeriod=2592000.,
 externForcingCycle=31104000.,
 &

# Gridding parameters
 &PARM04
 usingSphericalPolarGrid=.TRUE.,
 delR= 10.,10.,20.,25.,25.,25.,25.,50.,50.,50.,
       50.,100.,100.,100.,100.,100.,100.,100.,100.,100.,100.,
       100.,100.,250.,250.,500.,500.,500.,500.,500.,500.,500.,
 ygOrigin=-79.5,
 dySpacing=1.,
 dxSpacing=1.,
 &

# Input datasets
 &PARM05
 bathyFile=      'bathymetry.bin',
 hydrogThetaFile='lev_t.bin',
 hydrogSaltFile= 'lev_s.bin',
 zonalWindFile=  'taux.bin',
 meridWindFile=  'tauy.bin',
thetaClimFile=  'sst.bin',
 saltClimFile=   'sss.bin',
 surfQFile=      'qnet.bin',
 the_run_name=   'global_oce_latlon',
# fresh water flux is turned on, comment next line to it turn off
# (maybe better with surface salinity restoring)
 EmPmRFile=      'emp.bin',
&"


Can anybody please guide me to overcome this error?

Thanking in advance,

Abhisek



More information about the MITgcm-support mailing list