/[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.92 by heimbach, Tue Jul 8 15:00:26 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         &                mdsioLocalDir
87        CHARACTER*(5) checkPtSuff(maxNoChkptLev)        CHARACTER*(5) checkPtSuff(maxNoChkptLev)
88        CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile        CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile
89        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile
# Line 100  C                        LINEAR, POLY3, Line 101  C                        LINEAR, POLY3,
101        CHARACTER*(MAX_LEN_FNAM) pSurfInitFile        CHARACTER*(MAX_LEN_FNAM) pSurfInitFile
102        CHARACTER*(MAX_LEN_FNAM) dQdTfile        CHARACTER*(MAX_LEN_FNAM) dQdTfile
103        CHARACTER*(MAX_LEN_FNAM) ploadFile        CHARACTER*(MAX_LEN_FNAM) ploadFile
104          CHARACTER*(MAX_LEN_FNAM) mdsioLocalDir
105        CHARACTER*(6) eosType        CHARACTER*(6) eosType
106          CHARACTER*(10) pickupSuff
107    
108  C--   COMMON /PARM_I/ Integer valued parameters used by the model.  C--   COMMON /PARM_I/ Integer valued parameters used by the model.
109  C     cg2dMaxIters        :: Maximum number of iterations in the  C     cg2dMaxIters        :: Maximum number of iterations in the
# Line 121  C     readBinaryPrec      :: Precision u Line 124  C     readBinaryPrec      :: Precision u
124  C     nCheckLev           :: Holds current checkpoint level  C     nCheckLev           :: Holds current checkpoint level
125  C     nonlinFreeSurf      :: option related to non-linear free surface  C     nonlinFreeSurf      :: option related to non-linear free surface
126  C                           =0 Linear free surface ; >0 Non-linear  C                           =0 Linear free surface ; >0 Non-linear
127    C     select_rStar        :: option related to r* vertical coordinate
128    C                           =0 (default) use r coord. ; > 0 use r*
129    
130        COMMON /PARM_I/        COMMON /PARM_I/
131       &        cg2dMaxIters,       &        cg2dMaxIters,
# Line 131  C                           =0 Linear fr Line 136  C                           =0 Linear fr
136       &        numStepsPerPickup,       &        numStepsPerPickup,
137       &        writeStatePrec, nCheckLev,       &        writeStatePrec, nCheckLev,
138       &        writeBinaryPrec, readBinaryPrec,       &        writeBinaryPrec, readBinaryPrec,
139       &        nonlinFreeSurf,       &        nonlinFreeSurf, select_rStar,
140       &        tempAdvScheme, saltAdvScheme, tracerAdvScheme       &        tempAdvScheme, saltAdvScheme, tracerAdvScheme,
141         &        debugLevel
142        INTEGER cg2dMaxIters        INTEGER cg2dMaxIters
143        INTEGER cg2dChkResFreq        INTEGER cg2dChkResFreq
144        INTEGER cg3dMaxIters        INTEGER cg3dMaxIters
# Line 146  C                           =0 Linear fr Line 152  C                           =0 Linear fr
152        INTEGER readBinaryPrec        INTEGER readBinaryPrec
153        INTEGER nCheckLev        INTEGER nCheckLev
154        INTEGER nonlinFreeSurf        INTEGER nonlinFreeSurf
155          INTEGER select_rStar
156        INTEGER tempAdvScheme        INTEGER tempAdvScheme
157        INTEGER saltAdvScheme        INTEGER saltAdvScheme
158        INTEGER tracerAdvScheme        INTEGER tracerAdvScheme
159          INTEGER debugLevel
160    
161    C
162          INTEGER debLevZero
163          PARAMETER(debLevZero=0)
164          INTEGER debLevA
165          PARAMETER(debLevA=1)
166          INTEGER debLevB
167          PARAMETER(debLevB=2)
168    
169  C--   COMMON /PARM_L/ Logical valued parameters used by the model.  C--   COMMON /PARM_L/ Logical valued parameters used by the model.
170  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian
# Line 195  C     tr1Stepping   :: Turns passive tra Line 211  C     tr1Stepping   :: Turns passive tra
211  C     useConstantF  :: Coriolis parameter set to f0  C     useConstantF  :: Coriolis parameter set to f0
212  C     useBetaPlaneF :: Coriolis parameter set to f0 + beta.y  C     useBetaPlaneF :: Coriolis parameter set to f0 + beta.y
213  C     useSphereF    :: Coriolis parameter set to 2.omega.sin(phi)  C     useSphereF    :: Coriolis parameter set to 2.omega.sin(phi)
214    C     useCDscheme   :: use CD-scheme to calculate Coriolis terms.
215  C     useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986)  C     useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986)
216  C     implicitDiffusion :: Turns implicit vertical diffusion on  C     implicitDiffusion :: Turns implicit vertical diffusion on
217  C     implicitViscosity :: Turns implicit vertical viscosity on  C     implicitViscosity :: Turns implicit vertical viscosity on
# Line 207  C     doSaltClimRelax  :: Set true if re Line 224  C     doSaltClimRelax  :: Set true if re
224  C                        climatology is required.  C                        climatology is required.
225  C     periodicExternalForcing :: Set true if forcing is time-dependant  C     periodicExternalForcing :: Set true if forcing is time-dependant
226  C     usingPCoords     :: Set to indicate that we are working in pressure  C     usingPCoords     :: Set to indicate that we are working in pressure
227  C                        coords.  C                        coords. (jmc: is it still used ?)
228  C     usingZCoords     :: Set to indicate that we are working in height  C     usingZCoords     :: Set to indicate that we are working in height
229  C                        coords.  C                        coords. (jmc: is it still used ?)
230    C     useDynP_inEos_Zc :: use the dynamical pressure in EOS (with Z-coord.)
231    C                         this requires specific code for restart & exchange
232  C     setCenterDr    :: set cell Center depth and put Interface at the middle  C     setCenterDr    :: set cell Center depth and put Interface at the middle
233  C     nonHydrostatic :: Using non-hydrostatic terms  C     nonHydrostatic :: Using non-hydrostatic terms
234  C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm  C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm
235  C     globalFiles    :: Selects between "global" and "tiled" files  C     globalFiles    :: Selects between "global" and "tiled" files
236    C     useSingleCpuIO :: On SGI platforms, option globalFiles is either
237    C                       slow (f77) or does not work (f90).  When
238    C                       useSingleCpuIO is set, mdsio_writefield.F
239    C                       outputs from master mpi process only.
240  C     allowFreezing  :: Allows water to freeze and form ice  C     allowFreezing  :: Allows water to freeze and form ice
241  C     groundAtK1  :: put the surface(k=1) at the Lower Boundary (=ground)  C     groundAtK1  :: put the surface(k=1) at the Lower Boundary (=ground)
242        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
# Line 229  C     groundAtK1  :: put the surface(k=1 Line 252  C     groundAtK1  :: put the surface(k=1
252       & momStepping, tempStepping, saltStepping, tr1Stepping,       & momStepping, tempStepping, saltStepping, tr1Stepping,
253       & metricTerms, usingSphericalPolarMTerms, useNHMTerms,       & metricTerms, usingSphericalPolarMTerms, useNHMTerms,
254       & useConstantF, useBetaPlaneF, useSphereF,       & useConstantF, useBetaPlaneF, useSphereF,
255         & useCDscheme,
256       & useEnergyConservingCoriolis, useJamartWetPoints,       & useEnergyConservingCoriolis, useJamartWetPoints,
257       & implicitDiffusion, implicitViscosity,       & implicitDiffusion, implicitViscosity,
258       & multiDimAdvection, forcing_In_AB,       & multiDimAdvection, forcing_In_AB,
259       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,
260       & periodicExternalForcing,       & periodicExternalForcing,
261       & usingPCoords, usingZCoords, setCenterDr,       & usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,
262       & nonHydrostatic, quasiHydrostatic, globalFiles,       & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,
263       & allowFreezing, groundAtK1,       & allowFreezing, groundAtK1,
264       & usePickupBeforeC35, debugMode,       & usePickupBeforeC35, debugMode,
265       & readPickupWithTracer, writePickupWithTracer       & readPickupWithTracer, writePickupWithTracer
# Line 272  C     groundAtK1  :: put the surface(k=1 Line 296  C     groundAtK1  :: put the surface(k=1
296        LOGICAL useConstantF        LOGICAL useConstantF
297        LOGICAL useBetaPlaneF        LOGICAL useBetaPlaneF
298        LOGICAL useSphereF        LOGICAL useSphereF
299          LOGICAL useCDscheme
300        LOGICAL useEnergyConservingCoriolis        LOGICAL useEnergyConservingCoriolis
301        LOGICAL useJamartWetPoints        LOGICAL useJamartWetPoints
302        LOGICAL implicitDiffusion        LOGICAL implicitDiffusion
# Line 284  C     groundAtK1  :: put the surface(k=1 Line 309  C     groundAtK1  :: put the surface(k=1
309        LOGICAL periodicExternalForcing        LOGICAL periodicExternalForcing
310        LOGICAL usingPCoords        LOGICAL usingPCoords
311        LOGICAL usingZCoords        LOGICAL usingZCoords
312          LOGICAL useDynP_inEos_Zc
313        LOGICAL setCenterDr        LOGICAL setCenterDr
314        LOGICAL nonHydrostatic        LOGICAL nonHydrostatic
315        LOGICAL quasiHydrostatic        LOGICAL quasiHydrostatic
316        LOGICAL globalFiles        LOGICAL globalFiles
317          LOGICAL useSingleCpuIO
318        LOGICAL allowFreezing        LOGICAL allowFreezing
319        LOGICAL groundAtK1        LOGICAL groundAtK1
320        LOGICAL usePickupBeforeC35        LOGICAL usePickupBeforeC35
# Line 339  C     f0        :: Reference coriolis pa Line 366  C     f0        :: Reference coriolis pa
366  C                 ( Southern edge f for beta plane )  C                 ( Southern edge f for beta plane )
367  C     beta      :: df/dy ( s^-1.m^-1 )  C     beta      :: df/dy ( s^-1.m^-1 )
368  C     omega     :: Angular velocity ( rad/s )  C     omega     :: Angular velocity ( rad/s )
369    C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)
370  C     viscAh    :: Eddy viscosity coeff. for mixing of  C     viscAh    :: Eddy viscosity coeff. for mixing of
371  C                 momentum laterally ( m^2/s )  C                 momentum laterally ( m^2/s )
372  C     viscAz    :: Eddy viscosity coeff. for mixing of  C     viscAz    :: Eddy viscosity coeff. for mixing of
# Line 457  C      --"-"--  Quadratic  ( linear: 1/s Line 485  C      --"-"--  Quadratic  ( linear: 1/s
485       & endTime, chkPtFreq, pchkPtFreq, dumpFreq,       & endTime, chkPtFreq, pchkPtFreq, dumpFreq,
486       & diagFreq, taveFreq, tave_lastIter, monitorFreq,       & diagFreq, taveFreq, tave_lastIter, monitorFreq,
487       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
488       & cosPower, cAdjFreq, omega,       & cosPower, cAdjFreq, omega, rotationPeriod,
489       & tauThetaClimRelax, lambdaThetaClimRelax,       & tauThetaClimRelax, lambdaThetaClimRelax,
490       & tauSaltClimRelax, lambdaSaltClimRelax,       & tauSaltClimRelax, lambdaSaltClimRelax,
491       & tauTr1ClimRelax, lambdaTr1ClimRelax,       & tauTr1ClimRelax, lambdaTr1ClimRelax,
# Line 551  C      --"-"--  Quadratic  ( linear: 1/s Line 579  C      --"-"--  Quadratic  ( linear: 1/s
579        _RL cosPower        _RL cosPower
580        _RL cAdjFreq        _RL cAdjFreq
581        _RL omega        _RL omega
582          _RL rotationPeriod
583        _RL tauThetaClimRelax        _RL tauThetaClimRelax
584        _RL lambdaThetaClimRelax        _RL lambdaThetaClimRelax
585        _RL tauSaltClimRelax        _RL tauSaltClimRelax
# Line 570  C      --"-"--  Quadratic  ( linear: 1/s Line 599  C      --"-"--  Quadratic  ( linear: 1/s
599        _RL bottomDragLinear        _RL bottomDragLinear
600        _RL bottomDragQuadratic        _RL bottomDragQuadratic
601    
602        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp,  C--   COMMON /PARM_A/ Thermodynamics constants ?
603       &                Lamba_theta        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp
604        _RL HeatCapacity_Cp        _RL HeatCapacity_Cp
       _RL Lamba_theta  
605        _RL recip_Cp        _RL recip_Cp
606    
607  CmlC Equation of State (polynomial coeffients)  C--   COMMON /PARM_ATM/ Atmospheric physical parameters (Ideal Gas EOS, ...)
608  Cml      COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS  C     celsius2K :: convert centigrade (Celsius) degree to Kelvin
609  Cml      _RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1)  C     atm_Po    :: standard reference pressure
610  CmlC Linear equation of state  C     atm_Cp    :: specific heat (Cp) of the (dry) air at constant pressure
611  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  
612  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)
613  C     Integr_GeoPot :: option to select the way we integrate the geopotential  C     integr_GeoPot :: option to select the way we integrate the geopotential
614  C                     (still a subject of discussions ...)  C                     (still a subject of discussions ...)
615        COMMON /PARM_ATM/ atm_cp, atm_kappa, atm_po,  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is
616       &                  Integr_GeoPot  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)
617        _RL atm_cp, atm_kappa, atm_po        COMMON /PARM_ATM/
618        INTEGER Integr_GeoPot       &            celsius2K,
619         &            atm_Cp, atm_Rd, atm_kappa, atm_Po,
620         &            integr_GeoPot, selectFindRoSurf
621          _RL celsius2K
622          _RL atm_Po, atm_Cp, atm_Rd, atm_kappa
623          INTEGER integr_GeoPot, selectFindRoSurf
624    
625  C Logical flags for selecting packages  C Logical flags for selecting packages
626        LOGICAL useKPP        LOGICAL useKPP
627        LOGICAL useGMRedi        LOGICAL useGMRedi
628        LOGICAL useOBCS        LOGICAL useOBCS
629        LOGICAL useAIM        LOGICAL useAIM
630          LOGICAL useLand
631        LOGICAL useGrdchk        LOGICAL useGrdchk
632        LOGICAL useECCO        LOGICAL useECCO
633        LOGICAL useSHAP_FILT        LOGICAL useSHAP_FILT
634        LOGICAL useZONAL_FILT        LOGICAL useZONAL_FILT
635        LOGICAL useFLT        LOGICAL useFLT
636        LOGICAL usePTRACERS        LOGICAL usePTRACERS
637          LOGICAL useSBO
638          LOGICAL useSEAICE
639    cswdblk -- add ---
640          LOGICAL useBulkforce
641          LOGICAL useThermSEAICE
642    cswdblk --- end add ---
643        COMMON /PARM_PACKAGES/        COMMON /PARM_PACKAGES/
644       &        useKPP, useGMRedi, useOBCS, useAIM, useECCO,       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,
645       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
646       &        usePTRACERS       &        usePTRACERS,  useSBO, useSEAICE,
647    cswdblk -- add ---
648         &        useThermSEAICE, useBulkforce
649    cswdblk --- end add ---
Legend:
Removed from v.1.79  
changed lines
  Added in v.1.92
  ViewVC Help
Powered by ViewVC 1.1.22