[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