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

Abhisek Chakraborty abhisek.sac at gmail.com
Tue Jun 12 12:05:49 EDT 2012


Respected Dr. Martin,

Many Thanks for your prompt and guiding reply. In the mean time I
adjusted model time step and ran the global model for 2 years without
any interruption.


regards,
Abhisek


On 6/12/12, Martin Losch <Martin.Losch at awi.de> wrote:
> 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
>
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support
>



More information about the MITgcm-support mailing list