[MITgcm-support] runtime error

[Patrick Heimbach]

Jan and Martin,

not entirely sure, but the array size 200 is associated to
MNC_MAX_FID (not MNC_MAX_ID as written in error message)
and is set in MNC_SIZE.h

I think (though no doc on this) that it's a buffer that keeps
track of the mnc file name counters to ensure that it generates
a new filename counter every time it needs one
(see S/R MNC_GET_NEXT_EMPTY_IND in mnc_utils.F)

So you could just up this number to expected MNC_MAX_FID
if you think that so many different names should have been created.
If not, wait for answer by the author of that package :o)

-Patrick



On Mar 2, 2007, at 5:19 AM, Jan Saynisch wrote:

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

---
Dr Patrick Heimbach | heimbach at mit.edu | http://www.mit.edu/~heimbach
MIT | EAPS, 54-1518 | 77 Massachusetts Ave | Cambridge, MA 02139, USA
FON: +1-617-253-5259 | FAX: +1-617-253-4464 | SKYPE: patrick.heimbach