[MITgcm-support] runtime error

Patrick Heimbach heimbach at MIT.EDU
Fri Mar 2 08:23:04 EST 2007


Oops, I just did injustice in claiming "no doc on this".
Many of the MNC parameters are actually documented in the manual,
see section 7.2
http://mitgcm.org/r2_web_testing/latest/online_documents/node259.html
(but not complete, the MNC_MAX_FID is missing).
-p.



On Mar 2, 2007, at 8:08 AM, Patrick Heimbach wrote:

> 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
>
>
> _______________________________________________
> 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





More information about the MITgcm-support mailing list