[MITgcm-support] Problem with tmin, tmax and monSolutionMaxRange

[Takaya Uchida]

Hi Siwen,

This usually occurs when there’s a numerical instability within your simulation.
It seems that your domain is realistic so I would suggest checking the bathymetry and/or boundary conditions of your domain as a starter.
Saving the model outputs at higher frequency also helps in identifying where the instability is propagating from.

- Takaya






On Oct 11, 2023, at 5:07 PM, Siwen Chen <siwen.chen at mail.mcgill.ca> wrote:

Hi
I'm using the MITgcm for a realistic simulation. I ran the model, and it went well at first, but after 5 steps, it ended run automatically and gives the following error:

(PID.TID 0024.0001) SOLUTION IS HEADING OUT OF BOUNDS: tMin,tMax= -6.713E+08  9.737E+08
(PID.TID 0024.0001)   exceeds allowed range (monSolutionMaxRange=  1.000E+03)
(PID.TID 0024.0001) MON_SOLUTION: STOPPING CALCULATION at Iter=     22711
(PID.TID 0024.0001) *** ERROR *** S/R ALL_PROC_DIE: ending the run

The size.H file is set as follows:

 PARAMETER (
     &           sNx =  60,
     &           sNy =  20,
     &           OLx =   3,
     &           OLy =   3,
     &           nSx =   1,
     &           nSy =   1,
     &           nPx =   4,
     &           nPy =   16,
     &           Nx  = sNx*nSx*nPx,
     &           Ny  = sNy*nSy*nPy,
     &           Nr  =   55)

C     MAX_OLX :: Set to the maximum overlap region size of any array
C     MAX_OLY    that will be exchanged. Controls the sizing of exch
C                routine buffers.
      INTEGER MAX_OLX
      INTEGER MAX_OLY
      PARAMETER ( MAX_OLX = OLx,
     &            MAX_OLY = OLy )

And I did't notice where appears the tMin and tMax. And the data file is set as:

# Model parameters
# Continuous equation parameters
 &PARM01
# ---------------------
# I/O parameters
 readBinaryPrec=32,
 writeBinaryPrec=32,
 useSingleCpuIO=.TRUE.,
# ---------------------
# activate partial cells
 hFacMinDr=5.,
 hFacMin=0.1,
# ---------------------
# free surface parameters
 implicitFreeSurface=.TRUE.,
 useRealFreshWaterFlux=.TRUE.
 exactConserv=.TRUE.,
# ---------------------
# Coriolis parameters at 60deg S for beta plane
 f0=-1.405e-4,
 beta=1.145e-11,
 selectCoriScheme=1,
# ---------------------
# Density and equation of state
 eosType='JMD95Z',
 tRef= 55*0.,
 sRef= 55*34.8,
 rhoConst=1035.,
 rhoConstFresh=1000.,
# ---------------------
# advection scheme
 tempAdvScheme=80,
 saltAdvScheme=80,
 staggerTimeStep=.TRUE.,
 vectorInvariantMomentum=.TRUE.,
# ---------------------
# Viscosity
# viscA4 maximum is ~8.0E13
 viscA4=4.6E10
 viscAr=1.E-5,
 implicitViscosity=.TRUE.,
 bottomDragQuadratic = 5E-2,
# ---------------------
# Diffusivity and convection
 diffKhT=0.,
 diffKrT=1.E-5,
 diffKhS=0.,
 diffKrS=1.E-5,
 implicitDiffusion=.TRUE.,
 &

# Elliptic solver parameters
 &PARM02
 cg2dTargetResidual=1.E-17,
 cg2dMaxIters=1000,
 &

# Time stepping parameters
 &PARM03

# for the longer (50 years) run set pickup frequency to  1 year and
# rolling pickup to 1 month, monitor monthly
 deltaT=600.0,
 nIter0=0,
 nTimeSteps=518400,
# pChkptFreq=31104000.,
# chkptFreq=2592000.,
# monitorFreq=2592000.,
# general settings
 dumpFreq=0.,
 dumpInitAndLast=.FALSE.,
 abEps=0.1,
 forcing_In_AB=.FALSE.,
 momDissip_In_AB=.FALSE.,
 &

# Gridding parameters
 &PARM04
 usingCartesianGrid=.TRUE.,
 delX=240*10.E3,
 delY=320*10.E3,
 delR=  2.12950489,   2.4308905 ,   2.77490323,   3.1675587 ,   3.61571496,
        4.1271872 ,   4.71087665,   5.37691441,   6.13682095,   7.00368101,
        7.99233291,   9.11956927,  10.40434342,  11.86797199,  13.5343178 ,
       15.42992871,  17.58409559,  20.02877577,  22.79830608,  25.92880133,
       29.4571029 ,  33.41911167,  37.8473202 ,  42.7673766 ,  48.19359034,
       54.12347726,  60.5317637 ,  67.36473409,  74.53630909,  81.92757038,
       89.39130539,  96.76222993, 103.87196047, 110.56607025, 116.71949717,
      122.24691006, 127.10620628, 131.29541339, 134.84493639, 137.80768112,
      140.2493299 , 142.24020421, 143.84931689, 145.14057911, 146.17078598,
      146.98889983, 147.63617797, 148.14678645, 148.54864422, 148.86433209,
      149.11196865, 149.30600265, 149.45790119, 149.57673112, 149.66964071,
 interDiffKr_pCell=.TRUE.,
 interViscAr_pCell=.TRUE.,
 &

# I/O datasets
 &PARM05
 bathyFile='TESTSETUP-SOMECONFIG-bathymetry_shelf.bin',
 hydrogThetaFile='TESTSETUP-SOMECONFIG-woa18_temperature_3D.bin',
 hydrogSaltFile='TESTSETUP-SOMECONFIG-woa18_salinity_3D.bin',
 &

What can I do to avoid such error?

Thanks very much!

Best,
Siwen
_______________________________________________
MITgcm-support mailing list
MITgcm-support at mitgcm.org<mailto:MITgcm-support at mitgcm.org>
https://urldefense.com/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support__;!!PhOWcWs!yHCfZujDzmxH8Vi6b8FsgXAwwiltfxUZFDKJe7CDxNI4-Akf6VjXiXulghsHJKJL0MOY-9U0mRHlc3PadA-A6jcB$

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20231011/e2332e07/attachment-0001.html>