[MITgcm-support] Global_ocean query
Amitabh Mitra
amitabhmitra1 at rediffmail.com
Thu Mar 17 05:05:20 EST 2005
Dear Sir,
I have run the global_ocean.90x40x15 experiment
with a changed resolution of 1x1. (80S to 80N )
Earlier the model was blowing up when I was putting
the value of viscAh = 5.E5 but when I changed it to
viscAh = 0.3E5 in the data file the model executed
well with a Normal End.
The file Size.h , Data file and the Output file
are attached with the mail.
I want to know if the output generated is valid.
How can I check its validity ?
Also how can I choose the value of viscAh and viscAr
in the Data file for a given experiment ?
Thanking You for the kind support,
With regards,
Amitabh Mitra
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20050317/3addb3f1/attachment.htm>
-------------- next part --------------
C $Header: /u/gcmpack/MITgcm/verification/global_ocean.90x40x15/code/SIZE.h,v 1.4 2003/12/10 16:25:57 adcroft Exp $
C $Name: checkpoint53 $
C
C /==========================================================\
C | SIZE.h Declare size of underlying computational grid. |
C |==========================================================|
C | The design here support a three-dimensional model grid |
C | with indices I,J and K. The three-dimensional domain |
C | is comprised of nPx*nSx blocks of size sNx along one axis|
C | nPy*nSy blocks of size sNy along another axis and one |
C | block of size Nz along the final axis. |
C | Blocks have overlap regions of size OLx and OLy along the|
C | dimensions that are subdivided. |
C \==========================================================/
C Voodoo numbers controlling data layout.
C sNx - No. X points in sub-grid.
C sNy - No. Y points in sub-grid.
C OLx - Overlap extent in X.
C OLy - Overlat extent in Y.
C nSx - No. sub-grids in X.
C nSy - No. sub-grids in Y.
C nPx - No. of processes to use in X.
C nPy - No. of processes to use in Y.
C Nx - No. points in X for the total domain.
C Ny - No. points in Y for the total domain.
C Nr - No. points in Z for full process domain.
INTEGER sNx
INTEGER sNy
INTEGER OLx
INTEGER OLy
INTEGER nSx
INTEGER nSy
INTEGER nPx
INTEGER nPy
INTEGER Nx
INTEGER Ny
INTEGER Nr
PARAMETER (
& sNx = 180,
& sNy = 160,
& OLx = 2,
& OLy = 2,
& nSx = 2,
& nSy = 1,
& nPx = 1,
& nPy = 1,
& Nx = sNx*nSx*nPx,
& Ny = sNy*nSy*nPy,
& Nr = 18)
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 )
-------------- next part --------------
# ====================
# | Model parameters |
# ====================
#
# Continuous equation parameters
&PARM01
tRef = 18*20.,
sRef = 18*35.,
viscAr=1.E-3,
viscAh=0.3E5,
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 = 20,
# 100 years of integration will yield a reasonable flow field
# startTime = 0.,
# endTime = 3110400000.,
deltaTmom = 1200.0,
tauCD = 321428.,
deltaTtracer= 43200.0,
deltaTClock = 43200.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= 864000.,
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.,
phiMin=-80.,
dySpacing=1.,
dxSpacing=1.,
&
# Input datasets
&PARM05
bathyFile= 'bathymetry.bin',
hydrogThetaFile='lev_t1.bin',
hydrogSaltFile= 'lev_s1.bin',
zonalWindFile= 'trenberth_taux1.bin',
meridWindFile= 'trenberth_tauy1.bin',
thetaClimFile= 'lev_sst1.bin',
saltClimFile= 'lev_sss1.bin',
surfQFile= 'ncep_qnet.bin',
# fresh water flux is turned off, uncomment next line to turn on
# (not recommened together with surface salinity restoring)
# EmPmRFile= 'ncep_emp.bin',
&
-------------- next part --------------
SGI_2100 10% mpirun -np 1 ./mitgcmuv
S/R EXTERNAL_FIELDS_LOAD: Reading new data 0.E+0, 0
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 0.0000E+00 -6.7200E-01 3.3977E+01 -1.306
2E-03 2.2001E-02 6.3990E+01 0.0000E+00 12 1 5.0000E-01 5.0000E-01
time,fu0,fu1,fu = 0.0000E+00 -2.6089E-03 -3.6054E-06 -1.3062E-03 5.0000E-01
5.0000E-01
cg2d: Sum(rhs),rhsMax = -9.43689570931383E-16 6.25759651411038E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 4.3200E+04 -6.4810E-01 3.3980E+01 -1.262
8E-03 2.3000E-02 6.4303E+01 0.0000E+00 12 1 5.1667E-01 4.8333E-01
time,fu0,fu1,fu = 4.3200E+04 -2.6089E-03 -3.6054E-06 -1.2628E-03 5.1667E-01
4.8333E-01
cg2d: Sum(rhs),rhsMax = 1.24622534514174E-14 8.20129536863549E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 8.6400E+04 -6.2419E-01 3.3984E+01 -1.219
4E-03 2.4000E-02 6.4616E+01 0.0000E+00 12 1 5.3333E-01 4.6667E-01
time,fu0,fu1,fu = 8.6400E+04 -2.6089E-03 -3.6054E-06 -1.2194E-03 5.3333E-01
4.6667E-01
cg2d: Sum(rhs),rhsMax = -8.38495939348149E-14 7.10786477693076E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 1.2960E+05 -6.0029E-01 3.3987E+01 -1.176
0E-03 2.5000E-02 6.4929E+01 0.0000E+00 12 1 5.5000E-01 4.5000E-01
time,fu0,fu1,fu = 1.2960E+05 -2.6089E-03 -3.6054E-06 -1.1760E-03 5.5000E-01
4.5000E-01
cg2d: Sum(rhs),rhsMax = -3.85819848292002E-13 5.51159196363458E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 1.7280E+05 -5.7639E-01 3.3990E+01 -1.132
6E-03 2.6000E-02 6.5241E+01 0.0000E+00 12 1 5.6667E-01 4.3333E-01
time,fu0,fu1,fu = 1.7280E+05 -2.6089E-03 -3.6054E-06 -1.1326E-03 5.6667E-01
4.3333E-01
cg2d: Sum(rhs),rhsMax = -9.00002294912383E-13 3.88138580152627E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 2.1600E+05 -5.5248E-01 3.3994E+01 -1.089
1E-03 2.7000E-02 6.5554E+01 0.0000E+00 12 1 5.8333E-01 4.1667E-01
time,fu0,fu1,fu = 2.1600E+05 -2.6089E-03 -3.6054E-06 -1.0891E-03 5.8333E-01
4.1667E-01
cg2d: Sum(rhs),rhsMax = -1.32316380074826E-12 3.47090721467047E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 2.5920E+05 -5.2858E-01 3.3997E+01 -1.045
7E-03 2.8000E-02 6.5867E+01 0.0000E+00 12 1 6.0000E-01 4.0000E-01
time,fu0,fu1,fu = 2.5920E+05 -2.6089E-03 -3.6054E-06 -1.0457E-03 6.0000E-01
4.0000E-01
cg2d: Sum(rhs),rhsMax = -1.73971947958762E-12 3.49876028393096E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 3.0240E+05 -5.0468E-01 3.4001E+01 -1.002
3E-03 2.9000E-02 6.6180E+01 0.0000E+00 12 1 6.1667E-01 3.8333E-01
time,fu0,fu1,fu = 3.0240E+05 -2.6089E-03 -3.6054E-06 -1.0023E-03 6.1667E-01
3.8333E-01
cg2d: Sum(rhs),rhsMax = -1.55830903736387E-12 3.48255840815390E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 3.4560E+05 -4.8077E-01 3.4004E+01 -9.588
7E-04 3.0000E-02 6.6493E+01 0.0000E+00 12 1 6.3333E-01 3.6667E-01
time,fu0,fu1,fu = 3.4560E+05 -2.6089E-03 -3.6054E-06 -9.5887E-04 6.3333E-01
3.6667E-01
cg2d: Sum(rhs),rhsMax = -1.29851684960158E-12 3.81930848390533E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 3.8880E+05 -4.5687E-01 3.4007E+01 -9.154
5E-04 3.0999E-02 6.6806E+01 0.0000E+00 12 1 6.5000E-01 3.5000E-01
time,fu0,fu1,fu = 3.8880E+05 -2.6089E-03 -3.6054E-06 -9.1545E-04 6.5000E-01
3.5000E-01
cg2d: Sum(rhs),rhsMax = -1.64401825486493E-12 3.47391159117064E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 4.3200E+05 -4.3297E-01 3.4011E+01 -8.720
3E-04 3.1999E-02 6.7118E+01 0.0000E+00 12 1 6.6667E-01 3.3333E-01
time,fu0,fu1,fu = 4.3200E+05 -2.6089E-03 -3.6054E-06 -8.7203E-04 6.6667E-01
3.3333E-01
cg2d: Sum(rhs),rhsMax = -8.74411654194773E-13 6.48870960518738E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 4.7520E+05 -4.0906E-01 3.4014E+01 -8.286
1E-04 3.2999E-02 6.7431E+01 0.0000E+00 12 1 6.8333E-01 3.1667E-01
time,fu0,fu1,fu = 4.7520E+05 -2.6089E-03 -3.6054E-06 -8.2861E-04 6.8333E-01
3.1667E-01
cg2d: Sum(rhs),rhsMax = -1.24433796599988E-12 5.41338069110558E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 5.1840E+05 -3.8516E-01 3.4017E+01 -7.851
8E-04 3.3999E-02 6.7744E+01 0.0000E+00 12 1 7.0000E-01 3.0000E-01
time,fu0,fu1,fu = 5.1840E+05 -2.6089E-03 -3.6054E-06 -7.8518E-04 7.0000E-01
3.0000E-01
cg2d: Sum(rhs),rhsMax = -7.14983627858601E-13 1.05918014829022E+01
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 5.6160E+05 -3.6126E-01 3.4021E+01 -7.417
6E-04 3.4999E-02 6.8057E+01 0.0000E+00 12 1 7.1667E-01 2.8333E-01
time,fu0,fu1,fu = 5.6160E+05 -2.6089E-03 -3.6054E-06 -7.4176E-04 7.1667E-01
2.8333E-01
cg2d: Sum(rhs),rhsMax = -8.23341395062016E-13 8.99124789778651E+00
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 6.0480E+05 -3.3735E-01 3.4024E+01 -6.983
4E-04 3.5999E-02 6.8370E+01 0.0000E+00 12 1 7.3333E-01 2.6667E-01
time,fu0,fu1,fu = 6.0480E+05 -2.6089E-03 -3.6054E-06 -6.9834E-04 7.3333E-01
2.6667E-01
cg2d: Sum(rhs),rhsMax = -3.57047724719450E-13 1.44485639578020E+01
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 6.4800E+05 -3.1345E-01 3.4027E+01 -6.549
2E-04 3.6999E-02 6.8682E+01 0.0000E+00 12 1 7.5000E-01 2.5000E-01
time,fu0,fu1,fu = 6.4800E+05 -2.6089E-03 -3.6054E-06 -6.5492E-04 7.5000E-01
2.5000E-01
cg2d: Sum(rhs),rhsMax = -2.61124455391837E-13 1.22537992960919E+01
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 6.9120E+05 -2.8955E-01 3.4031E+01 -6.115
0E-04 3.7998E-02 6.8995E+01 0.0000E+00 12 1 7.6667E-01 2.3333E-01
time,fu0,fu1,fu = 6.9120E+05 -2.6089E-03 -3.6054E-06 -6.1150E-04 7.6667E-01
2.3333E-01
cg2d: Sum(rhs),rhsMax = -1.24789067967868E-13 2.70498307318316E+01
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 7.3440E+05 -2.6564E-01 3.4034E+01 -5.680
8E-04 3.8998E-02 6.9308E+01 0.0000E+00 12 1 7.8333E-01 2.1667E-01
time,fu0,fu1,fu = 7.3440E+05 -2.6089E-03 -3.6054E-06 -5.6808E-04 7.8333E-01
2.1667E-01
cg2d: Sum(rhs),rhsMax = -6.32827124036339E-14 3.63142289857903E+01
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 7.7760E+05 -2.4174E-01 3.4038E+01 -5.246
6E-04 3.9998E-02 6.9621E+01 0.0000E+00 12 1 8.0000E-01 2.0000E-01
time,fu0,fu1,fu = 7.7760E+05 -2.6089E-03 -3.6054E-06 -5.2466E-04 8.0000E-01
2.0000E-01
cg2d: Sum(rhs),rhsMax = 2.81996648254790E-13 3.23252016866515E+01
time,SST,SSS,fu,fv,Q,E-P,i0,i1,a,b = 8.2080E+05 -2.1784E-01 3.4041E+01 -4.812
4E-04 4.0998E-02 6.9934E+01 0.0000E+00 12 1 8.1667E-01 1.8333E-01
time,fu0,fu1,fu = 8.2080E+05 -2.6089E-03 -3.6054E-06 -4.8124E-04 8.1667E-01
1.8333E-01
cg2d: Sum(rhs),rhsMax = 7.50510764646606E-14 7.84673516208490E+01
hallo 1
STOP NORMAL END
STOP
More information about the MITgcm-support
mailing list