/[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.85 by dimitri, Tue Feb 18 05:33:53 2003 UTC revision 1.98 by jmc, Sun Nov 23 01:25:55 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 134  C                           =0 (default) Line 137  C                           =0 (default)
137       &        writeStatePrec, nCheckLev,       &        writeStatePrec, nCheckLev,
138       &        writeBinaryPrec, readBinaryPrec,       &        writeBinaryPrec, readBinaryPrec,
139       &        nonlinFreeSurf, select_rStar,       &        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 152  C                           =0 (default) Line 156  C                           =0 (default)
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 198  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 210  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
# Line 221  C     useSingleCpuIO :: On SGI platforms Line 237  C     useSingleCpuIO :: On SGI platforms
237  C                       slow (f77) or does not work (f90).  When  C                       slow (f77) or does not work (f90).  When
238  C                       useSingleCpuIO is set, mdsio_writefield.F  C                       useSingleCpuIO is set, mdsio_writefield.F
239  C                       outputs from master mpi process only.  C                       outputs from master mpi process only.
240  C     allowFreezing  :: Allows water to freeze and form ice  C     allowFreezing  :: Allows surface water to freeze and form ice
241    C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)
242  C     groundAtK1  :: put the surface(k=1) at the Lower Boundary (=ground)  C     groundAtK1  :: put the surface(k=1) at the Lower Boundary (=ground)
243        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
244       & usingCurvilinearGrid,       & usingCurvilinearGrid,
# Line 236  C     groundAtK1  :: put the surface(k=1 Line 253  C     groundAtK1  :: put the surface(k=1
253       & momStepping, tempStepping, saltStepping, tr1Stepping,       & momStepping, tempStepping, saltStepping, tr1Stepping,
254       & metricTerms, usingSphericalPolarMTerms, useNHMTerms,       & metricTerms, usingSphericalPolarMTerms, useNHMTerms,
255       & useConstantF, useBetaPlaneF, useSphereF,       & useConstantF, useBetaPlaneF, useSphereF,
256         & useCDscheme,
257       & useEnergyConservingCoriolis, useJamartWetPoints,       & useEnergyConservingCoriolis, useJamartWetPoints,
258       & implicitDiffusion, implicitViscosity,       & implicitDiffusion, implicitViscosity,
259       & multiDimAdvection, forcing_In_AB,       & multiDimAdvection, forcing_In_AB,
260       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,
261       & periodicExternalForcing,       & periodicExternalForcing,
262       & usingPCoords, usingZCoords, setCenterDr,       & usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,
263       & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,       & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,
264       & allowFreezing, groundAtK1,       & allowFreezing, useOldFreezing, groundAtK1,
265       & usePickupBeforeC35, debugMode,       & usePickupBeforeC35, debugMode,
266       & readPickupWithTracer, writePickupWithTracer       & readPickupWithTracer, writePickupWithTracer
267        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
# Line 279  C     groundAtK1  :: put the surface(k=1 Line 297  C     groundAtK1  :: put the surface(k=1
297        LOGICAL useConstantF        LOGICAL useConstantF
298        LOGICAL useBetaPlaneF        LOGICAL useBetaPlaneF
299        LOGICAL useSphereF        LOGICAL useSphereF
300          LOGICAL useCDscheme
301        LOGICAL useEnergyConservingCoriolis        LOGICAL useEnergyConservingCoriolis
302        LOGICAL useJamartWetPoints        LOGICAL useJamartWetPoints
303        LOGICAL implicitDiffusion        LOGICAL implicitDiffusion
# Line 291  C     groundAtK1  :: put the surface(k=1 Line 310  C     groundAtK1  :: put the surface(k=1
310        LOGICAL periodicExternalForcing        LOGICAL periodicExternalForcing
311        LOGICAL usingPCoords        LOGICAL usingPCoords
312        LOGICAL usingZCoords        LOGICAL usingZCoords
313          LOGICAL useDynP_inEos_Zc
314        LOGICAL setCenterDr        LOGICAL setCenterDr
315        LOGICAL nonHydrostatic        LOGICAL nonHydrostatic
316        LOGICAL quasiHydrostatic        LOGICAL quasiHydrostatic
317        LOGICAL globalFiles        LOGICAL globalFiles
318        LOGICAL useSingleCpuIO        LOGICAL useSingleCpuIO
319        LOGICAL allowFreezing        LOGICAL allowFreezing
320          LOGICAL useOldFreezing
321        LOGICAL groundAtK1        LOGICAL groundAtK1
322        LOGICAL usePickupBeforeC35        LOGICAL usePickupBeforeC35
323        LOGICAL debugMode        LOGICAL debugMode
# Line 324  C           number of iterations for con Line 345  C           number of iterations for con
345  C           drop form 192 -> 134! Need to investigate this further!  C           drop form 192 -> 134! Need to investigate this further!
346  C           For now I have introduced a parameter cg2dpcOffDFac which  C           For now I have introduced a parameter cg2dpcOffDFac which
347  C           defaults to 0.51 but can be set at runtime.  C           defaults to 0.51 but can be set at runtime.
 C     delP      :: Vertical grid spacing ( Pa ).  
 C     delZ      :: Vertical grid spacing ( m  ).  
348  C     delR      :: Vertical grid spacing ( units of r ).  C     delR      :: Vertical grid spacing ( units of r ).
349  C     delRc     :: Vertical grid spacing between cell centers (r unit).  C     delRc     :: Vertical grid spacing between cell centers (r unit).
350  C     delX      :: Separation between cell faces (m) or (deg), depending  C     delX      :: Separation between cell faces (m) or (deg), depending
# Line 347  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 )
 C     viscAz    :: Eddy viscosity coeff. for mixing of  
 C                 momentum vertically ( m^2/s )  
 C     viscAp    :: Eddy viscosity coeff. for mixing of  
 C                 momentum vertically ( Pa^2/s )  
372  C     viscAr    :: Eddy viscosity coeff. for mixing of  C     viscAr    :: Eddy viscosity coeff. for mixing of
373  C                 momentum vertically ( units of r^2/s )  C                 momentum vertically ( units of r^2/s )
374  C     viscA4    :: Biharmonic viscosity coeff. for mixing of  C     viscA4    :: Biharmonic viscosity coeff. for mixing of
375  C                 momentum laterally ( m^4/s )  C                 momentum laterally ( m^4/s )
376  C     diffKhT   :: Laplacian diffusion coeff. for mixing of  C     diffKhT   :: Laplacian diffusion coeff. for mixing of
377  C                 heat laterally ( m^2/s )  C                 heat laterally ( m^2/s )
 C     diffKzT   :: Laplacian diffusion coeff. for mixing of  
 C                 heat vertically ( m^2/s )  
 C     diffKpT   :: Laplacian diffusion coeff. for mixing of  
 C                 heat vertically ( Pa^2/s )  
378  C     diffKrT   :: Laplacian diffusion coeff. for mixing of  C     diffKrT   :: Laplacian diffusion coeff. for mixing of
379  C                 heat vertically ( units of r^2/s )  C                 heat vertically ( units of r^2/s )
380  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of
381  C                 heat laterally ( m^4/s )  C                 heat laterally ( m^4/s )
382  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of
383  C                 salt laterally ( m^2/s )  C                 salt laterally ( m^2/s )
 C     diffKzS   :: Laplacian diffusion coeff. for mixing of  
 C                 salt vertically ( m^2/s )  
 C     diffKpS   :: Laplacian diffusion coeff. for mixing of  
 C                 salt vertically ( Pa^2/s )  
384  C     diffKrS   :: Laplacian diffusion coeff. for mixing of  C     diffKrS   :: Laplacian diffusion coeff. for mixing of
385  C                 salt vertically ( units of r^2/s )  C                 salt vertically ( units of r^2/s )
386  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of
# Line 428  C     tauThetaClimRelax :: Relaxation to Line 436  C     tauThetaClimRelax :: Relaxation to
436  C     lambdaThetaClimRelax :: Inverse time scale for relaxation ( 1/s ).  C     lambdaThetaClimRelax :: Inverse time scale for relaxation ( 1/s ).
437  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).
438  C     lambdaSaltClimRelax :: Inverse time scale for relaxation ( 1/s ).  C     lambdaSaltClimRelax :: Inverse time scale for relaxation ( 1/s ).
439    C     latBandClimRelax :: latitude band where Relaxation to Clim. is applied,
440    C                         i.e. where |yC| <= latBandClimRelax
441  C     externForcingPeriod :: Is the period of which forcing varies (eg. 1 month)  C     externForcingPeriod :: Is the period of which forcing varies (eg. 1 month)
442  C     externForcingCycle :: Is the repeat time of the forcing (eg. 1 year)  C     externForcingCycle :: Is the repeat time of the forcing (eg. 1 year)
443  C                          (note: externForcingCycle must be an integer  C                          (note: externForcingCycle must be an integer
# Line 449  C     bottomDragLinear   :: Drag coeffic Line 459  C     bottomDragLinear   :: Drag coeffic
459  C      --"-"--  Quadratic  ( linear: 1/s, quadratic: 1/m )  C      --"-"--  Quadratic  ( linear: 1/s, quadratic: 1/m )
460        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
461       & cg2dpcOffDFac, cg3dTargetResidual,       & cg2dpcOffDFac, cg3dTargetResidual,
462       & delP, delZ, delR, delRc, delX, delY,       & delR, delRc, delX, delY,
463       & deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock,       & deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock,
464       & abeps, startTime,       & abeps, startTime,
465       & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,       & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,
466       & fCori, fCoriG, fCoriCos,       & fCori, fCoriG, fCoriCos,
467       & viscAh,  viscAz,  viscA4,  viscAr, viscAstrain, viscAtension,       & viscAh,  viscA4,  viscAr, viscAstrain, viscAtension,
468       & diffKhT, diffKzT, diffK4T, diffKrT,       & diffKhT, diffK4T, diffKrT,
469       & diffKhS, diffKzS, diffK4S, diffKrS,       & diffKhS, diffK4S, diffKrS,
470       & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,       & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,
471       & hFacMin, hFacMinDz, hFacInf, hFacSup,       & hFacMin, hFacMinDz, hFacInf, hFacSup,
472       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
# Line 465  C      --"-"--  Quadratic  ( linear: 1/s Line 475  C      --"-"--  Quadratic  ( linear: 1/s
475       & endTime, chkPtFreq, pchkPtFreq, dumpFreq,       & endTime, chkPtFreq, pchkPtFreq, dumpFreq,
476       & diagFreq, taveFreq, tave_lastIter, monitorFreq,       & diagFreq, taveFreq, tave_lastIter, monitorFreq,
477       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
478       & cosPower, cAdjFreq, omega,       & cosPower, cAdjFreq, omega, rotationPeriod,
479       & tauThetaClimRelax, lambdaThetaClimRelax,       & tauThetaClimRelax, lambdaThetaClimRelax,
480       & tauSaltClimRelax, lambdaSaltClimRelax,       & tauSaltClimRelax, lambdaSaltClimRelax,
481       & tauTr1ClimRelax, lambdaTr1ClimRelax,       & tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax,
482       & externForcingCycle, externForcingPeriod,       & externForcingCycle, externForcingPeriod,
483       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,
484       & viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp,       & hFacMinDr, hFacMinDp,
485       & horiVertRatio, recip_horiVertRatio,       & horiVertRatio, recip_horiVertRatio,
486       & ivdc_kappa, Ro_SeaLevel,       & ivdc_kappa, Ro_SeaLevel,
487       & bottomDragLinear,bottomDragQuadratic       & bottomDragLinear,bottomDragQuadratic
# Line 480  C      --"-"--  Quadratic  ( linear: 1/s Line 490  C      --"-"--  Quadratic  ( linear: 1/s
490        _RL cg2dTargetResWunit        _RL cg2dTargetResWunit
491        _RL cg3dTargetResidual        _RL cg3dTargetResidual
492        _RL cg2dpcOffDFac        _RL cg2dpcOffDFac
       _RL delZ(Nr)  
       _RL delP(Nr)  
493        _RL delR(Nr)        _RL delR(Nr)
494        _RL delRc(Nr+1)        _RL delRc(Nr+1)
495        _RL delX(Nx)        _RL delX(Nx)
# Line 510  C      --"-"--  Quadratic  ( linear: 1/s Line 518  C      --"-"--  Quadratic  ( linear: 1/s
518        _RL viscAh        _RL viscAh
519        _RL viscAstrain        _RL viscAstrain
520        _RL viscAtension        _RL viscAtension
       _RL viscAz  
       _RL viscAp  
521        _RL viscAr        _RL viscAr
522        _RL viscA4        _RL viscA4
523        _RL diffKhT        _RL diffKhT
524        _RL diffKrT        _RL diffKrT
       _RL diffKzT  
       _RL diffKpT  
525        _RL diffK4T        _RL diffK4T
526        _RL diffKhS        _RL diffKhS
527        _RL diffKrS        _RL diffKrS
       _RL diffKzS  
       _RL diffKpS  
528        _RL diffK4S        _RL diffK4S
529        _RL delt        _RL delt
530        _RL tauCD        _RL tauCD
# Line 559  C      --"-"--  Quadratic  ( linear: 1/s Line 561  C      --"-"--  Quadratic  ( linear: 1/s
561        _RL cosPower        _RL cosPower
562        _RL cAdjFreq        _RL cAdjFreq
563        _RL omega        _RL omega
564          _RL rotationPeriod
565        _RL tauThetaClimRelax        _RL tauThetaClimRelax
566        _RL lambdaThetaClimRelax        _RL lambdaThetaClimRelax
567        _RL tauSaltClimRelax        _RL tauSaltClimRelax
568        _RL lambdaSaltClimRelax        _RL lambdaSaltClimRelax
569        _RL tauTr1ClimRelax        _RL tauTr1ClimRelax
570        _RL lambdaTr1ClimRelax        _RL lambdaTr1ClimRelax
571          _RL latBandClimRelax
572        _RL externForcingCycle        _RL externForcingCycle
573        _RL externForcingPeriod        _RL externForcingPeriod
574        _RL convertFW2Salt        _RL convertFW2Salt
# Line 578  C      --"-"--  Quadratic  ( linear: 1/s Line 582  C      --"-"--  Quadratic  ( linear: 1/s
582        _RL bottomDragLinear        _RL bottomDragLinear
583        _RL bottomDragQuadratic        _RL bottomDragQuadratic
584    
585        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp,  C--   COMMON /PARM_A/ Thermodynamics constants ?
586       &                Lamba_theta        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp
587        _RL HeatCapacity_Cp        _RL HeatCapacity_Cp
       _RL Lamba_theta  
588        _RL recip_Cp        _RL recip_Cp
589    
590  CmlC Equation of State (polynomial coeffients)  C--   COMMON /PARM_ATM/ Atmospheric physical parameters (Ideal Gas EOS, ...)
591  Cml      COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS  C     celsius2K :: convert centigrade (Celsius) degree to Kelvin
 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  
   
 C Atmospheric physical parameters (Ideal Gas EOS, ...)  
592  C     atm_Po    :: standard reference pressure  C     atm_Po    :: standard reference pressure
593  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
594  C     atm_Rd    :: gas constant for dry air  C     atm_Rd    :: gas constant for dry air
595  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)
596    C     atm_Rq    :: water vapour specific volume anomaly relative to dry air
597    C                  (e.g. typical value = (29/18 -1) 10^-3 with q [g/kg])
598  C     integr_GeoPot :: option to select the way we integrate the geopotential  C     integr_GeoPot :: option to select the way we integrate the geopotential
599  C                     (still a subject of discussions ...)  C                     (still a subject of discussions ...)
600  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is
601  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)
602        COMMON /PARM_ATM/ atm_Cp, atm_Rd, atm_kappa, atm_Po,        COMMON /PARM_ATM/
603       &                  integr_GeoPot, selectFindRoSurf       &            celsius2K,
604        _RL atm_Po, atm_Cp, atm_Rd, atm_kappa       &            atm_Cp, atm_Rd, atm_kappa, atm_Rq, atm_Po,
605         &            integr_GeoPot, selectFindRoSurf
606          _RL celsius2K
607          _RL atm_Po, atm_Cp, atm_Rd, atm_kappa, atm_Rq
608        INTEGER integr_GeoPot, selectFindRoSurf        INTEGER integr_GeoPot, selectFindRoSurf
609    
610  C Logical flags for selecting packages  C Logical flags for selecting packages
# Line 613  C Logical flags for selecting packages Line 612  C Logical flags for selecting packages
612        LOGICAL useGMRedi        LOGICAL useGMRedi
613        LOGICAL useOBCS        LOGICAL useOBCS
614        LOGICAL useAIM        LOGICAL useAIM
615          LOGICAL useLand
616        LOGICAL useGrdchk        LOGICAL useGrdchk
617        LOGICAL useECCO        LOGICAL useECCO
618        LOGICAL useSHAP_FILT        LOGICAL useSHAP_FILT
619        LOGICAL useZONAL_FILT        LOGICAL useZONAL_FILT
620        LOGICAL useFLT        LOGICAL useFLT
621        LOGICAL usePTRACERS        LOGICAL usePTRACERS
622          LOGICAL useSBO
623        LOGICAL useSEAICE        LOGICAL useSEAICE
624  cswdblk -- add ---        LOGICAL useBulkForce
625        LOGICAL useBulkforce        LOGICAL useThSIce
       LOGICAL useThermSEAICE  
 cswdblk --- end add ---  
626        COMMON /PARM_PACKAGES/        COMMON /PARM_PACKAGES/
627       &        useKPP, useGMRedi, useOBCS, useAIM, useECCO,       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,
628       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
629       &        usePTRACERS, useSEAICE,       &        usePTRACERS,  useSBO, useSEAICE,
630  cswdblk -- add ---       &        useThSIce, useBulkForce
      &        useThermSEAICE, useBulkforce  
 cswdblk --- end add ---  
Legend:
Removed from v.1.85  
changed lines
  Added in v.1.98
  ViewVC Help
Powered by ViewVC 1.1.22