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

[Martin Losch]

Hi Abhisek,

this is a very common problem (and can be found multiple times in the support archives): You cannot use the parameters of coarse model with a different (higher) resolution. Please read up on CFL criteria, both for advection and viscosity/diffusion, and adjust your time step accordingly. The information can be found in the manual, but I'd actually recommend a textbook about basics of numerical modelling.

Martin


On Jun 11, 2012, at 5:41 PM, Abhisek Chakraborty wrote:

> 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
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support