/[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.91 by jmc, Thu Jun 12 18:08:04 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 570  C      --"-"--  Quadratic  ( linear: 1/s Line 587  C      --"-"--  Quadratic  ( linear: 1/s
587        _RL bottomDragLinear        _RL bottomDragLinear
588        _RL bottomDragQuadratic        _RL bottomDragQuadratic
589    
590        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp,  C--   COMMON /PARM_A/ Thermodynamics constants ?
591       &                Lamba_theta        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp
592        _RL HeatCapacity_Cp        _RL HeatCapacity_Cp
       _RL Lamba_theta  
593        _RL recip_Cp        _RL recip_Cp
594    
595  CmlC Equation of State (polynomial coeffients)  C--   COMMON /PARM_ATM/ Atmospheric physical parameters (Ideal Gas EOS, ...)
596  Cml      COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS  C     celsius2K :: convert centigrade (Celsius) degree to Kelvin
597  Cml      _RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1)  C     atm_Po    :: standard reference pressure
598  CmlC Linear equation of state  C     atm_Cp    :: specific heat (Cp) of the (dry) air at constant pressure
599  CmlC     tAlpha    :: Linear EOS thermal expansion coefficient ( 1/degree ).  C     atm_Rd    :: gas constant for dry air
 CmlC     sBeta     :: Linear EOS haline contraction coefficient.  
 Cml      COMMON /PARM_EOS_LIN/ tAlpha,sBeta,eosType  
 Cml      _RL tAlpha  
 Cml      _RL sBeta  
   
 C Atmospheric physical parameters (Ideal Gas EOS, ...)  
 C     atm_po    :: standard reference pressure  
 C     atm_cp    :: specific heat (Cp) of the (dry) air at constant pressure  
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/
606        INTEGER Integr_GeoPot       &            celsius2K,
607         &            atm_Cp, atm_Rd, atm_kappa, atm_Po,
608         &            integr_GeoPot, selectFindRoSurf
609          _RL celsius2K
610          _RL atm_Po, atm_Cp, atm_Rd, atm_kappa
611          INTEGER integr_GeoPot, selectFindRoSurf
612    
613  C Logical flags for selecting packages  C Logical flags for selecting packages
614        LOGICAL useKPP        LOGICAL useKPP
615        LOGICAL useGMRedi        LOGICAL useGMRedi
616        LOGICAL useOBCS        LOGICAL useOBCS
617        LOGICAL useAIM        LOGICAL useAIM
618          LOGICAL useLand
619        LOGICAL useGrdchk        LOGICAL useGrdchk
620        LOGICAL useECCO        LOGICAL useECCO
621        LOGICAL useSHAP_FILT        LOGICAL useSHAP_FILT
622        LOGICAL useZONAL_FILT        LOGICAL useZONAL_FILT
623        LOGICAL useFLT        LOGICAL useFLT
624        LOGICAL usePTRACERS        LOGICAL usePTRACERS
625          LOGICAL useSBO
626          LOGICAL useSEAICE
627    cswdblk -- add ---
628          LOGICAL useBulkforce
629          LOGICAL useThermSEAICE
630    cswdblk --- end add ---
631        COMMON /PARM_PACKAGES/        COMMON /PARM_PACKAGES/
632       &        useKPP, useGMRedi, useOBCS, useAIM, useECCO,       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,
633       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
634       &        usePTRACERS       &        usePTRACERS,  useSBO, useSEAICE,
635    cswdblk -- add ---
636         &        useThermSEAICE, useBulkforce
637    cswdblk --- end add ---
Legend:
Removed from v.1.79  
changed lines
  Added in v.1.91
  ViewVC Help
Powered by ViewVC 1.1.22