/[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.73 by mlosch, Wed Jul 31 16:38:30 2002 UTC revision 1.95 by jmc, Wed Oct 15 22:57:43 2003 UTC
# Line 72  C     EmPmRfile       :: File containing Line 72  C     EmPmRfile       :: File containing
72  C     pLoadFile       :: File containing pressure loading  C     pLoadFile       :: File containing pressure loading
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:
76    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 79  C                        get buoyancy. Line 81  C                        get buoyancy.
81       &                saltClimFile, buoyancyRelation,       &                saltClimFile, buoyancyRelation,
82       &                EmPmRfile, surfQfile, surfQswfile,       &                EmPmRfile, surfQfile, surfQswfile,
83       &                uVelInitFile, vVelInitFile, pSurfInitFile,       &                uVelInitFile, vVelInitFile, pSurfInitFile,
84       &                dQdTfile, ploadFile       &                dQdTfile, ploadFile,
85         &                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 97  C                        get buoyancy. Line 101  C                        get buoyancy.
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
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 117  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 127  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 142  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 162  C     momPressureForcing :: Flag which t Line 182  C     momPressureForcing :: Flag which t
182  C                          on and off.  C                          on and off.
183  C     metricTerms   :: Flag which turns metric terms on or off.  C     metricTerms   :: Flag which turns metric terms on or off.
184  C     usingSphericalPolarMTerms :: If TRUE use spherical polar metric terms.  C     usingSphericalPolarMTerms :: If TRUE use spherical polar metric terms.
185    C     useNHMTerms :: If TRUE use non-hydrostatic metric terms.
186  C     useCoriolis   :: Flag which turns the coriolis terms on and off.  C     useCoriolis   :: Flag which turns the coriolis terms on and off.
187  C     tempDiffusion :: Flag which turns diffusion of temperature on  C     tempDiffusion :: Flag which turns diffusion of temperature on
188  C                     and off.  C                     and off.
# Line 190  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 202  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
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 221  C     groundAtK1  :: put the surface(k=1 Line 250  C     groundAtK1  :: put the surface(k=1
250       & useRealFreshWaterFlux,       & useRealFreshWaterFlux,
251       & rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf,       & rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf,
252       & momStepping, tempStepping, saltStepping, tr1Stepping,       & momStepping, tempStepping, saltStepping, tr1Stepping,
253       & metricTerms, usingSphericalPolarMTerms,       & 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, globalFiles,       & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,
263       & allowFreezing, groundAtK1,       & allowFreezing, groundAtK1,
264       & usePickupBeforeC35, debugMode,       & usePickupBeforeC35, debugMode,
265       & readPickupWithTracer, writePickupWithTracer       & readPickupWithTracer, writePickupWithTracer
# Line 237  C     groundAtK1  :: put the surface(k=1 Line 267  C     groundAtK1  :: put the surface(k=1
267        LOGICAL usingSphericalPolarGrid        LOGICAL usingSphericalPolarGrid
268        LOGICAL usingCurvilinearGrid        LOGICAL usingCurvilinearGrid
269        LOGICAL usingSphericalPolarMTerms        LOGICAL usingSphericalPolarMTerms
270          LOGICAL useNHMTerms
271        LOGICAL no_slip_sides        LOGICAL no_slip_sides
272        LOGICAL no_slip_bottom        LOGICAL no_slip_bottom
273        LOGICAL staggerTimeStep        LOGICAL staggerTimeStep
# Line 265  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 277  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
316        LOGICAL globalFiles        LOGICAL globalFiles
317          LOGICAL useSingleCpuIO
318        LOGICAL allowFreezing        LOGICAL allowFreezing
319        LOGICAL groundAtK1        LOGICAL groundAtK1
320        LOGICAL usePickupBeforeC35        LOGICAL usePickupBeforeC35
# Line 317  C     delY        on input flags. Line 352  C     delY        on input flags.
352  C     gravity   :: Accel. due to gravity ( m/s^2 )  C     gravity   :: Accel. due to gravity ( m/s^2 )
353  C     recip_gravity and its inverse  C     recip_gravity and its inverse
354  C     gBaro     :: Accel. due to gravity used in barotropic equation ( m/s^2 )  C     gBaro     :: Accel. due to gravity used in barotropic equation ( m/s^2 )
355  C     ronil     :: Reference density  C     rhoNil    :: Reference density for the linear equation of state
356  C     rhoConst  :: Vertically constant reference density  C     rhoConst  :: Vertically constant reference density
357    C     rhoConstFresh :: Constant reference density for fresh water (rain)
358  C     startTime :: Start time for model ( s )  C     startTime :: Start time for model ( s )
359  C     phiMin    :: Latitude of southern most cell face.  C     phiMin    :: Latitude of southern most cell face.
360  C     thetaMin  :: Longitude of western most cell face (this  C     thetaMin  :: Longitude of western most cell face (this
# Line 330  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 402  C     mtFacMom      :: Metric terms trac Line 439  C     mtFacMom      :: Metric terms trac
439  C     cosPower      :: Power of cosine of latitude to multiply viscosity  C     cosPower      :: Power of cosine of latitude to multiply viscosity
440  C     cAdjFreq      :: Frequency of convective adjustment  C     cAdjFreq      :: Frequency of convective adjustment
441  C  C
442  C     taveFreq      :: Frequency with which time-averaged model state is written to  C     taveFreq      :: Frequency with which time-averaged model state
443  C                     post-processing files ( s ).  C                      is written to post-processing files ( s ).
444    C     tave_lastIter :: (for state variable only) fraction of the last time
445    C                      step (of each taveFreq period) put in the time average.
446    C                      (fraction for 1rst iter = 1 - tave_lastIter)
447  C     tauThetaClimRelax :: Relaxation to climatology time scale ( s ).  C     tauThetaClimRelax :: Relaxation to climatology time scale ( s ).
448  C     lambdaThetaClimRelax :: Inverse time scale for relaxation ( 1/s ).  C     lambdaThetaClimRelax :: Inverse time scale for relaxation ( 1/s ).
449  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).
450  C     lambdaSaltClimRelax :: Inverse time scale for relaxation ( 1/s ).  C     lambdaSaltClimRelax :: Inverse time scale for relaxation ( 1/s ).
451    C     latBandClimRelax :: latitude band where Relaxation to Clim. is applied,
452    C                         i.e. where |yC| <= latBandClimRelax
453  C     externForcingPeriod :: Is the period of which forcing varies (eg. 1 month)  C     externForcingPeriod :: Is the period of which forcing varies (eg. 1 month)
454  C     externForcingCycle :: Is the repeat time of the forcing (eg. 1 year)  C     externForcingCycle :: Is the repeat time of the forcing (eg. 1 year)
455  C                          (note: externForcingCycle must be an integer  C                          (note: externForcingCycle must be an integer
# Line 433  C      --"-"--  Quadratic  ( linear: 1/s Line 475  C      --"-"--  Quadratic  ( linear: 1/s
475       & deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock,       & deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock,
476       & abeps, startTime,       & abeps, startTime,
477       & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,       & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,
478       & fCori, fCoriG,       & fCori, fCoriG, fCoriCos,
479       & viscAh,  viscAz,  viscA4,  viscAr, viscAstrain, viscAtension,       & viscAh,  viscAz,  viscA4,  viscAr, viscAstrain, viscAtension,
480       & diffKhT, diffKzT, diffK4T, diffKrT,       & diffKhT, diffKzT, diffK4T, diffKrT,
481       & diffKhS, diffKzS, diffK4S, diffKrS,       & diffKhS, diffKzS, diffK4S, diffKrS,
482       & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,       & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,
483       & hFacMin, hFacMinDz, hFacInf, hFacSup,       & hFacMin, hFacMinDz, hFacInf, hFacSup,
484       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
485       & recip_rhoConst, rhoConst, tRef, sRef,       & recip_rhoConst, rhoConst,
486         & rhoConstFresh, convertEmP2rUnit, tRef, sRef,
487       & endTime, chkPtFreq, pchkPtFreq, dumpFreq,       & endTime, chkPtFreq, pchkPtFreq, dumpFreq,
488       & diagFreq, taveFreq, monitorFreq,       & diagFreq, taveFreq, tave_lastIter, monitorFreq,
489       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
490       & cosPower, cAdjFreq, omega,       & cosPower, cAdjFreq, omega, rotationPeriod,
491       & tauThetaClimRelax, lambdaThetaClimRelax,       & tauThetaClimRelax, lambdaThetaClimRelax,
492       & tauSaltClimRelax, lambdaSaltClimRelax,       & tauSaltClimRelax, lambdaSaltClimRelax,
493       & tauTr1ClimRelax, lambdaTr1ClimRelax,       & tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax,
494       & externForcingCycle, externForcingPeriod,       & externForcingCycle, externForcingPeriod,
495       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,
496       & viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp,       & viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp,
# Line 513  C      --"-"--  Quadratic  ( linear: 1/s Line 556  C      --"-"--  Quadratic  ( linear: 1/s
556        _RL recip_rhonil        _RL recip_rhonil
557        _RL rhoConst        _RL rhoConst
558        _RL recip_rhoConst        _RL recip_rhoConst
559          _RL rhoConstFresh
560          _RL convertEmP2rUnit
561        _RL tRef(Nr)        _RL tRef(Nr)
562        _RL sRef(Nr)        _RL sRef(Nr)
563        _RS fCori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS fCori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
564        _RS fCoriG(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RS fCoriG(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
565          _RS fCoriCos(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
566        _RL startTime        _RL startTime
567        _RL endTime        _RL endTime
568        _RL chkPtFreq        _RL chkPtFreq
# Line 524  C      --"-"--  Quadratic  ( linear: 1/s Line 570  C      --"-"--  Quadratic  ( linear: 1/s
570        _RL dumpFreq        _RL dumpFreq
571        _RL diagFreq        _RL diagFreq
572        _RL taveFreq        _RL taveFreq
573          _RL tave_lastIter
574        _RL monitorFreq        _RL monitorFreq
575        _RL afFacMom        _RL afFacMom
576        _RL vfFacMom        _RL vfFacMom
# Line 534  C      --"-"--  Quadratic  ( linear: 1/s Line 581  C      --"-"--  Quadratic  ( linear: 1/s
581        _RL cosPower        _RL cosPower
582        _RL cAdjFreq        _RL cAdjFreq
583        _RL omega        _RL omega
584          _RL rotationPeriod
585        _RL tauThetaClimRelax        _RL tauThetaClimRelax
586        _RL lambdaThetaClimRelax        _RL lambdaThetaClimRelax
587        _RL tauSaltClimRelax        _RL tauSaltClimRelax
588        _RL lambdaSaltClimRelax        _RL lambdaSaltClimRelax
589        _RL tauTr1ClimRelax        _RL tauTr1ClimRelax
590        _RL lambdaTr1ClimRelax        _RL lambdaTr1ClimRelax
591          _RL latBandClimRelax
592        _RL externForcingCycle        _RL externForcingCycle
593        _RL externForcingPeriod        _RL externForcingPeriod
594        _RL convertFW2Salt        _RL convertFW2Salt
# Line 553  C      --"-"--  Quadratic  ( linear: 1/s Line 602  C      --"-"--  Quadratic  ( linear: 1/s
602        _RL bottomDragLinear        _RL bottomDragLinear
603        _RL bottomDragQuadratic        _RL bottomDragQuadratic
604    
605        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp,  C--   COMMON /PARM_A/ Thermodynamics constants ?
606       &                Lamba_theta        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp
607        _RL HeatCapacity_Cp        _RL HeatCapacity_Cp
       _RL Lamba_theta  
608        _RL recip_Cp        _RL recip_Cp
609    
610  C Equation of State (polynomial coeffients)  C--   COMMON /PARM_ATM/ Atmospheric physical parameters (Ideal Gas EOS, ...)
611        COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS  C     celsius2K :: convert centigrade (Celsius) degree to Kelvin
612        _RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1)  C     atm_Po    :: standard reference pressure
613  C Linear equation of state  C     atm_Cp    :: specific heat (Cp) of the (dry) air at constant pressure
614  C     tAlpha    :: Linear EOS thermal expansion coefficient ( 1/degree ).  C     atm_Rd    :: gas constant for dry air
 C     sBeta     :: Linear EOS haline contraction coefficient.  
       COMMON /PARM_EOS_LIN/ tAlpha,sBeta,eosType  
       _RL tAlpha  
       _RL sBeta  
       character*(6) eosType  
   
 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  
615  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)
616  C     Integr_GeoPot :: option to select the way we integrate the geopotential  C     atm_Rq    :: water vapour specific volume anomaly relative to dry air
617    C                  (e.g. typical value = (29/18 -1) 10^-3 with q [g/kg])
618    C     integr_GeoPot :: option to select the way we integrate the geopotential
619  C                     (still a subject of discussions ...)  C                     (still a subject of discussions ...)
620        COMMON /PARM_ATM/ atm_cp, atm_kappa, atm_po,  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is
621       &                  Integr_GeoPot  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)
622        _RL atm_cp, atm_kappa, atm_po        COMMON /PARM_ATM/
623        INTEGER Integr_GeoPot       &            celsius2K,
624         &            atm_Cp, atm_Rd, atm_kappa, atm_Rq, atm_Po,
625         &            integr_GeoPot, selectFindRoSurf
626          _RL celsius2K
627          _RL atm_Po, atm_Cp, atm_Rd, atm_kappa, atm_Rq
628          INTEGER integr_GeoPot, selectFindRoSurf
629    
630  C Logical flags for selecting packages  C Logical flags for selecting packages
631        LOGICAL useKPP        LOGICAL useKPP
632        LOGICAL useGMRedi        LOGICAL useGMRedi
633        LOGICAL useOBCS        LOGICAL useOBCS
634        LOGICAL useAIM        LOGICAL useAIM
635          LOGICAL useLand
636        LOGICAL useGrdchk        LOGICAL useGrdchk
637        LOGICAL useECCO        LOGICAL useECCO
638        LOGICAL useSHAP_FILT        LOGICAL useSHAP_FILT
639        LOGICAL useZONAL_FILT        LOGICAL useZONAL_FILT
640        LOGICAL useFLT        LOGICAL useFLT
641        LOGICAL usePTRACERS        LOGICAL usePTRACERS
642          LOGICAL useSBO
643          LOGICAL useSEAICE
644          LOGICAL useBulkforce
645          LOGICAL useThermSEAICE
646        COMMON /PARM_PACKAGES/        COMMON /PARM_PACKAGES/
647       &        useKPP, useGMRedi, useOBCS, useAIM, useECCO,       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,
648       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
649       &        usePTRACERS       &        usePTRACERS,  useSBO, useSEAICE,
650         &        useThermSEAICE, useBulkforce
Legend:
Removed from v.1.73  
changed lines
  Added in v.1.95
  ViewVC Help
Powered by ViewVC 1.1.22