/[MITgcm]/MITgcm/model/inc/PARAMS.h
ViewVC logotype

Diff of /MITgcm/model/inc/PARAMS.h

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.79 by adcroft, Thu Nov 7 21:51:15 2002 UTC revision 1.90 by adcroft, Thu Jun 5 16:03:05 2003 UTC
# Line 73  C     pLoadFile       :: File containing Line 73  C     pLoadFile       :: File containing
73  C     buoyancyRelation :: Flag used to indicate which relation to use to  C     buoyancyRelation :: Flag used to indicate which relation to use to
74  C                        get buoyancy.  C                        get buoyancy.
75  C     eosType         :: choose the equation of state:  C     eosType         :: choose the equation of state:
76  C                        LINEAR, POLY3, UNESCO, JMD95Z, JMD95P  C                        LINEAR, POLY3, UNESCO, JMD95Z, JMD95P, MDJWF, IDEALGAS
77        COMMON /PARM_C/ checkPtSuff,        COMMON /PARM_C/ checkPtSuff,
78       &                bathyFile, topoFile,       &                bathyFile, topoFile,
79       &                hydrogThetaFile, hydrogSaltFile,       &                hydrogThetaFile, hydrogSaltFile,
# Line 82  C                        LINEAR, POLY3, Line 82  C                        LINEAR, POLY3,
82       &                EmPmRfile, surfQfile, surfQswfile,       &                EmPmRfile, surfQfile, surfQswfile,
83       &                uVelInitFile, vVelInitFile, pSurfInitFile,       &                uVelInitFile, vVelInitFile, pSurfInitFile,
84       &                dQdTfile, ploadFile,       &                dQdTfile, ploadFile,
85       &                eosType       &                eosType, pickupSuff
86        CHARACTER*(5) checkPtSuff(maxNoChkptLev)        CHARACTER*(5) checkPtSuff(maxNoChkptLev)
87        CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile        CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile
88        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile
# Line 101  C                        LINEAR, POLY3, Line 101  C                        LINEAR, POLY3,
101        CHARACTER*(MAX_LEN_FNAM) dQdTfile        CHARACTER*(MAX_LEN_FNAM) dQdTfile
102        CHARACTER*(MAX_LEN_FNAM) ploadFile        CHARACTER*(MAX_LEN_FNAM) ploadFile
103        CHARACTER*(6) eosType        CHARACTER*(6) eosType
104          CHARACTER*(10) pickupSuff
105    
106  C--   COMMON /PARM_I/ Integer valued parameters used by the model.  C--   COMMON /PARM_I/ Integer valued parameters used by the model.
107  C     cg2dMaxIters        :: Maximum number of iterations in the  C     cg2dMaxIters        :: Maximum number of iterations in the
# Line 121  C     readBinaryPrec      :: Precision u Line 122  C     readBinaryPrec      :: Precision u
122  C     nCheckLev           :: Holds current checkpoint level  C     nCheckLev           :: Holds current checkpoint level
123  C     nonlinFreeSurf      :: option related to non-linear free surface  C     nonlinFreeSurf      :: option related to non-linear free surface
124  C                           =0 Linear free surface ; >0 Non-linear  C                           =0 Linear free surface ; >0 Non-linear
125    C     select_rStar        :: option related to r* vertical coordinate
126    C                           =0 (default) use r coord. ; > 0 use r*
127    
128        COMMON /PARM_I/        COMMON /PARM_I/
129       &        cg2dMaxIters,       &        cg2dMaxIters,
# Line 131  C                           =0 Linear fr Line 134  C                           =0 Linear fr
134       &        numStepsPerPickup,       &        numStepsPerPickup,
135       &        writeStatePrec, nCheckLev,       &        writeStatePrec, nCheckLev,
136       &        writeBinaryPrec, readBinaryPrec,       &        writeBinaryPrec, readBinaryPrec,
137       &        nonlinFreeSurf,       &        nonlinFreeSurf, select_rStar,
138       &        tempAdvScheme, saltAdvScheme, tracerAdvScheme       &        tempAdvScheme, saltAdvScheme, tracerAdvScheme
139        INTEGER cg2dMaxIters        INTEGER cg2dMaxIters
140        INTEGER cg2dChkResFreq        INTEGER cg2dChkResFreq
# Line 146  C                           =0 Linear fr Line 149  C                           =0 Linear fr
149        INTEGER readBinaryPrec        INTEGER readBinaryPrec
150        INTEGER nCheckLev        INTEGER nCheckLev
151        INTEGER nonlinFreeSurf        INTEGER nonlinFreeSurf
152          INTEGER select_rStar
153        INTEGER tempAdvScheme        INTEGER tempAdvScheme
154        INTEGER saltAdvScheme        INTEGER saltAdvScheme
155        INTEGER tracerAdvScheme        INTEGER tracerAdvScheme
# Line 195  C     tr1Stepping   :: Turns passive tra Line 199  C     tr1Stepping   :: Turns passive tra
199  C     useConstantF  :: Coriolis parameter set to f0  C     useConstantF  :: Coriolis parameter set to f0
200  C     useBetaPlaneF :: Coriolis parameter set to f0 + beta.y  C     useBetaPlaneF :: Coriolis parameter set to f0 + beta.y
201  C     useSphereF    :: Coriolis parameter set to 2.omega.sin(phi)  C     useSphereF    :: Coriolis parameter set to 2.omega.sin(phi)
202    C     useCDscheme   :: use CD-scheme to calculate Coriolis terms.
203  C     useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986)  C     useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986)
204  C     implicitDiffusion :: Turns implicit vertical diffusion on  C     implicitDiffusion :: Turns implicit vertical diffusion on
205  C     implicitViscosity :: Turns implicit vertical viscosity on  C     implicitViscosity :: Turns implicit vertical viscosity on
# Line 207  C     doSaltClimRelax  :: Set true if re Line 212  C     doSaltClimRelax  :: Set true if re
212  C                        climatology is required.  C                        climatology is required.
213  C     periodicExternalForcing :: Set true if forcing is time-dependant  C     periodicExternalForcing :: Set true if forcing is time-dependant
214  C     usingPCoords     :: Set to indicate that we are working in pressure  C     usingPCoords     :: Set to indicate that we are working in pressure
215  C                        coords.  C                        coords. (jmc: is it still used ?)
216  C     usingZCoords     :: Set to indicate that we are working in height  C     usingZCoords     :: Set to indicate that we are working in height
217  C                        coords.  C                        coords. (jmc: is it still used ?)
218    C     useDynP_inEos_Zc :: use the dynamical pressure in EOS (with Z-coord.)
219    C                         this requires specific code for restart & exchange
220  C     setCenterDr    :: set cell Center depth and put Interface at the middle  C     setCenterDr    :: set cell Center depth and put Interface at the middle
221  C     nonHydrostatic :: Using non-hydrostatic terms  C     nonHydrostatic :: Using non-hydrostatic terms
222  C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm  C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm
223  C     globalFiles    :: Selects between "global" and "tiled" files  C     globalFiles    :: Selects between "global" and "tiled" files
224    C     useSingleCpuIO :: On SGI platforms, option globalFiles is either
225    C                       slow (f77) or does not work (f90).  When
226    C                       useSingleCpuIO is set, mdsio_writefield.F
227    C                       outputs from master mpi process only.
228  C     allowFreezing  :: Allows water to freeze and form ice  C     allowFreezing  :: Allows water to freeze and form ice
229  C     groundAtK1  :: put the surface(k=1) at the Lower Boundary (=ground)  C     groundAtK1  :: put the surface(k=1) at the Lower Boundary (=ground)
230        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
# Line 229  C     groundAtK1  :: put the surface(k=1 Line 240  C     groundAtK1  :: put the surface(k=1
240       & momStepping, tempStepping, saltStepping, tr1Stepping,       & momStepping, tempStepping, saltStepping, tr1Stepping,
241       & metricTerms, usingSphericalPolarMTerms, useNHMTerms,       & metricTerms, usingSphericalPolarMTerms, useNHMTerms,
242       & useConstantF, useBetaPlaneF, useSphereF,       & useConstantF, useBetaPlaneF, useSphereF,
243         & useCDscheme,
244       & useEnergyConservingCoriolis, useJamartWetPoints,       & useEnergyConservingCoriolis, useJamartWetPoints,
245       & implicitDiffusion, implicitViscosity,       & implicitDiffusion, implicitViscosity,
246       & multiDimAdvection, forcing_In_AB,       & multiDimAdvection, forcing_In_AB,
247       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,
248       & periodicExternalForcing,       & periodicExternalForcing,
249       & usingPCoords, usingZCoords, setCenterDr,       & usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,
250       & nonHydrostatic, quasiHydrostatic, globalFiles,       & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,
251       & allowFreezing, groundAtK1,       & allowFreezing, groundAtK1,
252       & usePickupBeforeC35, debugMode,       & usePickupBeforeC35, debugMode,
253       & readPickupWithTracer, writePickupWithTracer       & readPickupWithTracer, writePickupWithTracer
# Line 272  C     groundAtK1  :: put the surface(k=1 Line 284  C     groundAtK1  :: put the surface(k=1
284        LOGICAL useConstantF        LOGICAL useConstantF
285        LOGICAL useBetaPlaneF        LOGICAL useBetaPlaneF
286        LOGICAL useSphereF        LOGICAL useSphereF
287          LOGICAL useCDscheme
288        LOGICAL useEnergyConservingCoriolis        LOGICAL useEnergyConservingCoriolis
289        LOGICAL useJamartWetPoints        LOGICAL useJamartWetPoints
290        LOGICAL implicitDiffusion        LOGICAL implicitDiffusion
# Line 284  C     groundAtK1  :: put the surface(k=1 Line 297  C     groundAtK1  :: put the surface(k=1
297        LOGICAL periodicExternalForcing        LOGICAL periodicExternalForcing
298        LOGICAL usingPCoords        LOGICAL usingPCoords
299        LOGICAL usingZCoords        LOGICAL usingZCoords
300          LOGICAL useDynP_inEos_Zc
301        LOGICAL setCenterDr        LOGICAL setCenterDr
302        LOGICAL nonHydrostatic        LOGICAL nonHydrostatic
303        LOGICAL quasiHydrostatic        LOGICAL quasiHydrostatic
304        LOGICAL globalFiles        LOGICAL globalFiles
305          LOGICAL useSingleCpuIO
306        LOGICAL allowFreezing        LOGICAL allowFreezing
307        LOGICAL groundAtK1        LOGICAL groundAtK1
308        LOGICAL usePickupBeforeC35        LOGICAL usePickupBeforeC35
# Line 339  C     f0        :: Reference coriolis pa Line 354  C     f0        :: Reference coriolis pa
354  C                 ( Southern edge f for beta plane )  C                 ( Southern edge f for beta plane )
355  C     beta      :: df/dy ( s^-1.m^-1 )  C     beta      :: df/dy ( s^-1.m^-1 )
356  C     omega     :: Angular velocity ( rad/s )  C     omega     :: Angular velocity ( rad/s )
357    C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)
358  C     viscAh    :: Eddy viscosity coeff. for mixing of  C     viscAh    :: Eddy viscosity coeff. for mixing of
359  C                 momentum laterally ( m^2/s )  C                 momentum laterally ( m^2/s )
360  C     viscAz    :: Eddy viscosity coeff. for mixing of  C     viscAz    :: Eddy viscosity coeff. for mixing of
# Line 457  C      --"-"--  Quadratic  ( linear: 1/s Line 473  C      --"-"--  Quadratic  ( linear: 1/s
473       & endTime, chkPtFreq, pchkPtFreq, dumpFreq,       & endTime, chkPtFreq, pchkPtFreq, dumpFreq,
474       & diagFreq, taveFreq, tave_lastIter, monitorFreq,       & diagFreq, taveFreq, tave_lastIter, monitorFreq,
475       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
476       & cosPower, cAdjFreq, omega,       & cosPower, cAdjFreq, omega, rotationPeriod,
477       & tauThetaClimRelax, lambdaThetaClimRelax,       & tauThetaClimRelax, lambdaThetaClimRelax,
478       & tauSaltClimRelax, lambdaSaltClimRelax,       & tauSaltClimRelax, lambdaSaltClimRelax,
479       & tauTr1ClimRelax, lambdaTr1ClimRelax,       & tauTr1ClimRelax, lambdaTr1ClimRelax,
# Line 551  C      --"-"--  Quadratic  ( linear: 1/s Line 567  C      --"-"--  Quadratic  ( linear: 1/s
567        _RL cosPower        _RL cosPower
568        _RL cAdjFreq        _RL cAdjFreq
569        _RL omega        _RL omega
570          _RL rotationPeriod
571        _RL tauThetaClimRelax        _RL tauThetaClimRelax
572        _RL lambdaThetaClimRelax        _RL lambdaThetaClimRelax
573        _RL tauSaltClimRelax        _RL tauSaltClimRelax
# Line 576  C      --"-"--  Quadratic  ( linear: 1/s Line 593  C      --"-"--  Quadratic  ( linear: 1/s
593        _RL Lamba_theta        _RL Lamba_theta
594        _RL recip_Cp        _RL recip_Cp
595    
 CmlC Equation of State (polynomial coeffients)  
 Cml      COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS  
 Cml      _RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1)  
 CmlC Linear equation of state  
 CmlC     tAlpha    :: Linear EOS thermal expansion coefficient ( 1/degree ).  
 CmlC     sBeta     :: Linear EOS haline contraction coefficient.  
 Cml      COMMON /PARM_EOS_LIN/ tAlpha,sBeta,eosType  
 Cml      _RL tAlpha  
 Cml      _RL sBeta  
   
596  C Atmospheric physical parameters (Ideal Gas EOS, ...)  C Atmospheric physical parameters (Ideal Gas EOS, ...)
597  C     atm_po    :: standard reference pressure  C     atm_Po    :: standard reference pressure
598  C     atm_cp    :: specific heat (Cp) of the (dry) air at constant pressure  C     atm_Cp    :: specific heat (Cp) of the (dry) air at constant pressure
599    C     atm_Rd    :: gas constant for dry air
600  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)
601  C     Integr_GeoPot :: option to select the way we integrate the geopotential  C     integr_GeoPot :: option to select the way we integrate the geopotential
602  C                     (still a subject of discussions ...)  C                     (still a subject of discussions ...)
603        COMMON /PARM_ATM/ atm_cp, atm_kappa, atm_po,  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is
604       &                  Integr_GeoPot  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)
605        _RL atm_cp, atm_kappa, atm_po        COMMON /PARM_ATM/ atm_Cp, atm_Rd, atm_kappa, atm_Po,
606        INTEGER Integr_GeoPot       &                  integr_GeoPot, selectFindRoSurf
607          _RL atm_Po, atm_Cp, atm_Rd, atm_kappa
608          INTEGER integr_GeoPot, selectFindRoSurf
609    
610  C Logical flags for selecting packages  C Logical flags for selecting packages
611        LOGICAL useKPP        LOGICAL useKPP
# Line 608  C Logical flags for selecting packages Line 618  C Logical flags for selecting packages
618        LOGICAL useZONAL_FILT        LOGICAL useZONAL_FILT
619        LOGICAL useFLT        LOGICAL useFLT
620        LOGICAL usePTRACERS        LOGICAL usePTRACERS
621          LOGICAL useSBO
622          LOGICAL useSEAICE
623    cswdblk -- add ---
624          LOGICAL useBulkforce
625          LOGICAL useThermSEAICE
626    cswdblk --- end add ---
627        COMMON /PARM_PACKAGES/        COMMON /PARM_PACKAGES/
628       &        useKPP, useGMRedi, useOBCS, useAIM, useECCO,       &        useKPP, useGMRedi, useOBCS, useAIM, useECCO,
629       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
630       &        usePTRACERS       &        usePTRACERS,  useSBO, useSEAICE,
631    cswdblk -- add ---
632         &        useThermSEAICE, useBulkforce
633    cswdblk --- end add ---
Legend:
Removed from v.1.79  
changed lines
  Added in v.1.90
  ViewVC Help
Powered by ViewVC 1.1.22