[MITgcm-support] What is the demand for bathyFile ?

Brainf brainf at 163.com
Tue Oct 31 21:03:01 EST 2006


Dear mitgcm,
    I want to simulate my 0.25*0.25 degree (100x80x24 grid) example base on Global Ocean Simulation at 4 Resolution, but something wrong happend.  Does some demand for the bathyFile ? Example the bathyfile must smooth enjoyably or some other demand?
   I generate my bathy file from ETOP5, and the corresponding hydrogThetaFile ,hydrogSaltFile from Levitus. 
    Of course ,I edit  the SIZE.h file and compiled it successed
...............//my size.h//......... 
     PARAMETER (
     &           sNx =  50,
     &           sNy =  40,
     &           OLx =   2,
     &           OLy =   2,
     &           nSx =   2,
     &           nSy =   1,
     &           nPx =   1,
     &           nPy =   1,
     &           Nx  = sNx*nSx*nPx,
     &           Ny  = sNy*nSy*nPy,
     &           Nr  =  24) 
.....///...............................

But I run it when "%MON time_tsnumber = 3", it say errors:
   SOLUTION IS HEADING OUT OF BOUNDS: tMin,tMax=  -1.900E+00   1.558E+04
MON_SOLUTION: STOPPING CALCULATION
 S/R EEDIE: Only     0 threads have completed,     1 are expected for this confi
guration!

 Possibly you have different setenv PARALLEL and nThreads?
STOP MON_SOLUTION: STOPPED DUE TO EXTREME VALUES OF SOLUTION statement executed

///////////// and my data file//////////////
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
 &PARM01
 tRef = 24*20.,
 sRef = 24*35.,
 viscAr=1.E-3,
 viscAh=5.E5,
 diffKhT=0.0, 
 diffKrT=3.E-5,
 diffKhS=0.0,  
 diffKrS=3.E-5,
 rhonil=1035.,
 gravity=9.81,
 eosType = 'JMD95Z',
 ivdc_kappa=100.,
 implicitDiffusion=.TRUE.,
 useOldFreezing=.TRUE.,
 useRealFreshWaterFlux=.TRUE.,
 useCDscheme=.TRUE.,
 useNHMTerms=.TRUE.,
# turn on looped cells
 hFacMin=.05,
 hFacMindr=50.,
# set precision of data files
 readBinaryPrec=32,
 &

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

# Time stepping parameters
 &PARM03
 nIter0 =      0,
 nTimeSteps = 1000,
# 100 years of integration will yield a reasonable flow field
# startTime  =          0.,
# endTime    = 3110400000.,
 deltaTmom = 1200.0,
 tauCD =     321428.,
 deltaTtracer= 1200.0,
 deltaTClock = 1200.0,
# if you are using a version later than checkpoint45d on the main branch
# you can uncomment the following line and increase the time step 
# deltaTtracer and deltaTClock to 172800.0 as well to speed up the 
# asynchronous time stepping
# deltaTfreesurf = 172800.0,
 abEps = 0.1,
 pChkptFreq= 311040000.,
# dumpFreq=   311040000.,
# dumpFreq=   86400.,
 dumpFreq=   43200.,
 taveFreq=   311040000.,
 taveFreq=   864000.,
#monitorFreq=31104000.,
 monitorFreq=1.,
# 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= 10., 10., 10., 20., 25., 
       25., 25., 25., 50., 50., 
       50., 100., 100., 100., 100.,
       100., 100., 100., 100., 100.,
       100., 100., 150., 250., 
 phiMin=5.,
 thetaMin=105.,
 dySpacing=0.25,
 dxSpacing=0.25,
 &

# Input datasets
 &PARM05
 bathyFile=      'fc_dixing100x80x24.bin',
 hydrogThetaFile='fc_t_100x80x24.bin',
 hydrogSaltFile= 'fc_s_100x80x24.bin',
 &
 ////////////////////////////end///////////

I will upload my  bathyFile at next mail if you want test it.

Best regards. 

Brainf
brainf at 163.com
2006-11-01





More information about the MITgcm-support mailing list