[MITgcm-support] MITgcm in the netcdf offline mode
Karsten Friis
kfriis at meteo.psu.edu
Fri Oct 1 11:14:54 EDT 2004
Hi,
I try to get the MITgcm running in its netCDF offline mode. It is running
that way on fishtank.mit.edu. The same setup is not running on the Penn
State machine diatom.met.psu.edu. The compiling to the mitgcmuv works fine,
but starting the binary has some problems by reading in the different data*
files [namely:
data data.gchem data.gmredi data.mnc data.off data.pkg data.ptracers].
To read in the data.mnc file I had to exchange the namelist terminator from
& to /, now the program is stuck with the next file [data] which it is
tries to read in:
(PID.TID 0000.0001) MNC_READPARMS: finished reading data.mnc
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) // Model parameter file "data"
(PID.TID 0000.0001) // =======================================================
(PID.TID 0000.0001) ># ====================
(PID.TID 0000.0001) ># | Model parameters |
(PID.TID 0000.0001) ># ====================
(PID.TID 0000.0001) >#
(PID.TID 0000.0001) ># Continuous equation parameters
(PID.TID 0000.0001) > &PARM01
(PID.TID 0000.0001) > tRef=15*20.,
(PID.TID 0000.0001) > sRef=15*35.,
(PID.TID 0000.0001) > viscA4=0.,
(PID.TID 0000.0001) > viscAh=5.E5,
(PID.TID 0000.0001) > diffKhT=0.E3,
(PID.TID 0000.0001) > diffKhS=0.E3,
(PID.TID 0000.0001) > viscAz=1.E-3,
(PID.TID 0000.0001) > diffKzT=5.E-5,
(PID.TID 0000.0001) > diffKzS=5.E-5,
(PID.TID 0000.0001) > beta=1.E-11,
(PID.TID 0000.0001) > f0=1.e-4,
(PID.TID 0000.0001) > tAlpha=2.E-4,
(PID.TID 0000.0001) > sBeta =7.4E-4,
(PID.TID 0000.0001) > gravity=9.81,
(PID.TID 0000.0001) > gBaro=9.81,
(PID.TID 0000.0001) > rhonil=1035.,
(PID.TID 0000.0001) > rigidLid=.FALSE.,
(PID.TID 0000.0001) > implicitFreeSurface=.TRUE.,
(PID.TID 0000.0001) > eosType='POLY3',
(PID.TID 0000.0001) > implicitDiffusion=.TRUE.,
(PID.TID 0000.0001) > implicitViscosity=.TRUE.,
(PID.TID 0000.0001) > ivdc_kappa=100.,
(PID.TID 0000.0001) > tempAdvScheme = 2
(PID.TID 0000.0001) > saltAdvScheme = 2
(PID.TID 0000.0001) > tempStepping = .TRUE.
(PID.TID 0000.0001) > saltStepping = .TRUE.
(PID.TID 0000.0001) > tr1Stepping = .FALSE.
(PID.TID 0000.0001) > multiDimAdvection = .FALSE.
(PID.TID 0000.0001) > allowFreezing=.TRUE.,
(PID.TID 0000.0001) > useCDscheme=.FALSE.,
(PID.TID 0000.0001) >#usePickupBeforeC54=.TRUE.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Elliptic solver parameters
(PID.TID 0000.0001) > &PARM02
(PID.TID 0000.0001) > cg2dMaxIters=1000,
(PID.TID 0000.0001) > cg2dTargetResidual=1.E-13,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Time stepping parameters
(PID.TID 0000.0001) > &PARM03
(PID.TID 0000.0001) >#startTime = 0.,
(PID.TID 0000.0001) > nIter0 = 4248000,
(PID.TID 0000.0001) > nTimeSteps = 720,
(PID.TID 0000.0001) > deltaTmom = 900.0,
(PID.TID 0000.0001) >#tauCD = 321428.,
(PID.TID 0000.0001) > deltaTtracer= 43200.0,
(PID.TID 0000.0001) > deltaTClock = 43200.0,
(PID.TID 0000.0001) > abEps = 0.1,
(PID.TID 0000.0001) >#cAdjFreq = -1,
(PID.TID 0000.0001) > pChkptFreq=3110400000.,
(PID.TID 0000.0001) > chkptFreq= 3110400000.,
(PID.TID 0000.0001) > dumpFreq= 3110400000.,
(PID.TID 0000.0001) > taveFreq= 2592000.,
(PID.TID 0000.0001) > monitorFreq=311040000.,
(PID.TID 0000.0001) > tauThetaClimRelax = 5184000.0,
(PID.TID 0000.0001) > tauSaltClimRelax = 7776000.0,
(PID.TID 0000.0001) > periodicExternalForcing=.TRUE.,
(PID.TID 0000.0001) > externForcingPeriod=2592000.,
(PID.TID 0000.0001) > externForcingCycle=31104000.,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Gridding parameters
(PID.TID 0000.0001) > &PARM04
(PID.TID 0000.0001) > usingCartesianGrid=.FALSE.,
(PID.TID 0000.0001) > usingSphericalPolarGrid=.TRUE.,
(PID.TID 0000.0001) > delZ= 50., 70., 100., 140., 190.,
(PID.TID 0000.0001) > 240., 290., 340., 390., 440.,
(PID.TID 0000.0001) > 490., 540., 590., 640., 690.,
(PID.TID 0000.0001) > phiMin=-90.,
(PID.TID 0000.0001) > delX=128*2.8125,
(PID.TID 0000.0001) > delY=64*2.8125,
(PID.TID 0000.0001) > /
(PID.TID 0000.0001) >
(PID.TID 0000.0001) ># Input datasets
(PID.TID 0000.0001) > &PARM05
(PID.TID 0000.0001) > bathyFile= 'input/depth_g77.bin',
(PID.TID 0000.0001) > hydrogThetaFile='input/lev_clim_temp.bin',
(PID.TID 0000.0001) > hydrogSaltFile= 'input/lev_clim_salt.bin',
(PID.TID 0000.0001) > zonalWindFile= 'input/tren_taux.bin',
(PID.TID 0000.0001) > meridWindFile= 'input/tren_tauy.bin',
(PID.TID 0000.0001) > thetaClimFile= 'input/lev_monthly_temp.bin',
(PID.TID 0000.0001) > saltClimFile= 'input/lev_monthly_salt.bin',
(PID.TID 0000.0001) > surfQFile= 'input/shi_qnet.bin',
(PID.TID 0000.0001) > EmPmRFile= 'input/shi_empmr_year.bin',
(PID.TID 0000.0001) > /
(PID.TID 0000.0001)
(PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM01
Skipping namelist "PARM02": seeking namelist "PARM01".
Skipping namelist "PARM03": seeking namelist "PARM01".
Skipping namelist "PARM04": seeking namelist "PARM01".
Skipping namelist "PARM05": seeking namelist "PARM01".
namelist read start : end of file
apparent state: unit 11 named /tmp/tmp.FbdLY9X
last format: list io
lately reading sequential formatted external IO
Abort (core dumped)
Exchanging the name list terminator & by / doesn't help really but leads to
the following stop:
...
(PID.TID 0000.0001) > EmPmRFile= 'input/shi_empmr_year.bin',
(PID.TID 0000.0001) > /
(PID.TID 0000.0001)
(PID.TID 0000.0001) S/R INI_PARMS ; starts to read PARM01
namelist read start : end of file
apparent state: unit 11 named /tmp/tmp.FmLrVki
last format: list io
lately reading sequential formatted external IO
I am not quite sure if this is a step forward, do you have an idea what to
do next?
-best regards
karsten
email: kfriis at meteo.psu.edu
fax: ++814-865-3663
Mailing address:
Pennsylvania State University
Department of Meteorology
Karsten Friis, PhD
College of Earth and Mineral Sciences
503 Walker Building
University Park, PA 16802-5013
USA
More information about the MITgcm-support
mailing list