/[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.89 by jmc, Thu Apr 17 13:11:31 2003 UTC revision 1.99 by molod, Mon Dec 15 23:00:43 2003 UTC
# 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 224  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 247  C     groundAtK1  :: put the surface(k=1 Line 261  C     groundAtK1  :: put the surface(k=1
261       & periodicExternalForcing,       & periodicExternalForcing,
262       & usingPCoords, usingZCoords, useDynP_inEos_Zc, 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 303  C     groundAtK1  :: put the surface(k=1 Line 317  C     groundAtK1  :: put the surface(k=1
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 330  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 356  C     omega     :: Angular velocity ( ra Line 369  C     omega     :: Angular velocity ( ra
369  C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)  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 435  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 456  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 475  C      --"-"--  Quadratic  ( linear: 1/s Line 478  C      --"-"--  Quadratic  ( linear: 1/s
478       & cosPower, cAdjFreq, omega, rotationPeriod,       & 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 487  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 517  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 573  C      --"-"--  Quadratic  ( linear: 1/s Line 568  C      --"-"--  Quadratic  ( linear: 1/s
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 586  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  C Atmospheric physical parameters (Ideal Gas EOS, ...)  C--   COMMON /PARM_ATM/ Atmospheric physical parameters (Ideal Gas EOS, ...)
591    C     celsius2K :: convert centigrade (Celsius) degree to Kelvin
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 611  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
# Line 619  C Logical flags for selecting packages Line 621  C Logical flags for selecting packages
621        LOGICAL usePTRACERS        LOGICAL usePTRACERS
622        LOGICAL useSBO        LOGICAL useSBO
623        LOGICAL useSEAICE        LOGICAL useSEAICE
624  cswdblk -- add ---        LOGICAL useBulkForce
625        LOGICAL useBulkforce        LOGICAL useThSIce
626        LOGICAL useThermSEAICE        LOGICAL usefizhi
627  cswdblk --- end add ---        LOGICAL usegridalt
628        COMMON /PARM_PACKAGES/        COMMON /PARM_PACKAGES/
629       &        useKPP, useGMRedi, useOBCS, useAIM, useECCO,       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,
630       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
631       &        usePTRACERS,  useSBO, useSEAICE,       &        usePTRACERS,  useSBO, useSEAICE,
632  cswdblk -- add ---       &        useThSIce, useBulkForce, usefizhi,  usegridalt
      &        useThermSEAICE, useBulkforce  
 cswdblk --- end add ---  
Legend:
Removed from v.1.89  
changed lines
  Added in v.1.99
  ViewVC Help
Powered by ViewVC 1.1.22