[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