[MITgcm-support] runtime error

[Jan Saynisch]

Hello

I've encountered the following error after approx fifty years of model 
run on a Linux cluster:

(PID.TID 0000.0001) *** ERROR *** MNC_GET_NEXT_EMPTY_IND: array size 
200 exceeded--try increasing MNC_MAX_ID
(PID.TID 0000.0001) *** ERROR *** MNC_GET_NEXT_EMPTY_IND: occurred 
within the 'mnc_f_names' array
ABNORMAL END: S/R MNC_GET_NEXT_EMPTY_IND
3408.440u 533.780s 5:40:02.32 19.3%     0+0k 0+0io 354pf+0w

What goes wrong here?
Thank you - Jan

My configuration is as follows (packages & SIZE.h & data):

  &PACKAGES
  useGMRedi=.TRUE.,
#usePTRACERS=.TRUE.,
  useMNC=.TRUE.,
  useSBO=.TRUE.,
  &
********************************************************************

       PARAMETER (
      &           sNx =  45,
      &           sNy =  40,
      &           OLx =   2,
      &           OLy =   2,
      &           nSx =   2,
      &           nSy =   1,
      &           nPx =   1,
      &           nPy =   1,
      &           Nx  = sNx*nSx*nPx,
      &           Ny  = sNy*nSy*nPy,
      &           Nr  =  15)

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 buufers.
       INTEGER MAX_OLX
       INTEGER MAX_OLY
       PARAMETER ( MAX_OLX = OLx,
      &            MAX_OLY = OLy )
*************************************************************************

# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
  &PARM01
  tRef = 15*20.,
  sRef = 15*35.,
  viscAr=1.E-3,
  viscAh=5.E5,
  diffKhT=0.0,
  diffKrT=3.E-5,
  diffKhS=0.0,
  diffKrS=3.E-5,
  rhonil=1035.,
  rotationPeriod=86400.,
  gravity=9.81,
  eosType ='JMD95P',
  rigidLid=.FALSE.,
  implicitFreeSurface=.TRUE.,
  ivdc_kappa=100.,
  implicitDiffusion=.TRUE.,
  useOldFreezing=.TRUE.,
  exactConserv = .true.,
  useRealFreshWaterFlux=.TRUE.,
# turn on non-linear free surface
  nonlinFreeSurf=4,
  hFacInf=0.2,
  hFacSup=2.0,
# end
  useCDscheme=.TRUE.,
  useNHMTerms=.TRUE.,
# turn on looped cells
  hFacMin=.05,
  hFacMindr=50.,
# set precision of data files
  readBinaryPrec=32,
  writeBinaryPrec=64,
  &

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

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

# Gridding parameters
  &PARM04
  usingCartesianGrid=.FALSE.,
  usingSphericalPolarGrid=.TRUE.,
  delR= 50., 70., 100., 140., 190.,
        240., 290., 340., 390., 440.,
        490., 540., 590., 640., 690.,
  phiMin=-80.,
  dySpacing=4.,
  dxSpacing=4.,
  &

# Input datasets
  &PARM05
  bathyFile=      'bathymetry.bin',
  hydrogThetaFile='lev_t.bin',
  hydrogSaltFile= 'lev_s.bin',
  zonalWindFile=  'trenberth_taux.bin',
  meridWindFile=  'trenberth_tauy.bin',
  thetaClimFile=  'lev_sst.bin',
  saltClimFile=   'lev_sss.bin',
  surfQFile=      'ncep_qnet.bin',
  the_run_name=   'global_ocean.90x40x15',
# fresh water flux is turned off, uncomment next line to turn on
# (not recommened together with surface salinity restoring)
  EmPmRFile=      'ncep_emp.bin',
  &