[MITgcm-support] runtime error
Martin Losch
Martin.Losch at awi.de
Fri Mar 2 06:19:09 EST 2007
Hi,
one may add, that this happens, when Jan is trying to write his 50th
pickup files (pickup_cd is written properly, but pickup.
0000008640.t001.nc is not.)
Martin
On 2 Mar 2007, at 11:19, 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
More information about the MITgcm-support
mailing list