/[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.91 by jmc, Thu Jun 12 18:08:04 2003 UTC revision 1.102 by adcroft, Sat Feb 7 16:27:13 2004 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, pickupSuff       &                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        CHARACTER*(10) pickupSuff
107    
# Line 135  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 153  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 203  C     useCDscheme   :: use CD-scheme to Line 215  C     useCDscheme   :: use CD-scheme to
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
218    C     tempImplVertAdv :: Turns on implicit vertical advection for Temperature
219    C     saltImplVertAdv :: Turns on implicit vertical advection for Salinity
220    C     momImplVertAdv  :: Turns on implicit vertical advection for Momentum
221  C     multiDimAdvection :: Flag that enable multi-dimension advection  C     multiDimAdvection :: Flag that enable multi-dimension advection
222  C     forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution  C     forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution
223  C                      out off Adams-Bashforth time stepping.  C                      out off Adams-Bashforth time stepping.
# Line 225  C     useSingleCpuIO :: On SGI platforms Line 240  C     useSingleCpuIO :: On SGI platforms
240  C                       slow (f77) or does not work (f90).  When  C                       slow (f77) or does not work (f90).  When
241  C                       useSingleCpuIO is set, mdsio_writefield.F  C                       useSingleCpuIO is set, mdsio_writefield.F
242  C                       outputs from master mpi process only.  C                       outputs from master mpi process only.
243  C     allowFreezing  :: Allows water to freeze and form ice  C     allowFreezing  :: Allows surface water to freeze and form ice
244    C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)
245  C     groundAtK1  :: put the surface(k=1) at the Lower Boundary (=ground)  C     groundAtK1  :: put the surface(k=1) at the Lower Boundary (=ground)
246        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
247       & usingCurvilinearGrid,       & usingCurvilinearGrid,
# Line 243  C     groundAtK1  :: put the surface(k=1 Line 259  C     groundAtK1  :: put the surface(k=1
259       & useCDscheme,       & useCDscheme,
260       & useEnergyConservingCoriolis, useJamartWetPoints,       & useEnergyConservingCoriolis, useJamartWetPoints,
261       & implicitDiffusion, implicitViscosity,       & implicitDiffusion, implicitViscosity,
262         & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,
263       & multiDimAdvection, forcing_In_AB,       & multiDimAdvection, forcing_In_AB,
264       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,
265       & periodicExternalForcing,       & periodicExternalForcing,
266       & usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,       & usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,
267       & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,       & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,
268       & allowFreezing, groundAtK1,       & allowFreezing, useOldFreezing, groundAtK1,
269       & usePickupBeforeC35, debugMode,       & usePickupBeforeC35, debugMode,
270       & readPickupWithTracer, writePickupWithTracer       & readPickupWithTracer, writePickupWithTracer
271        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
# Line 289  C     groundAtK1  :: put the surface(k=1 Line 306  C     groundAtK1  :: put the surface(k=1
306        LOGICAL useJamartWetPoints        LOGICAL useJamartWetPoints
307        LOGICAL implicitDiffusion        LOGICAL implicitDiffusion
308        LOGICAL implicitViscosity        LOGICAL implicitViscosity
309          LOGICAL tempImplVertAdv
310          LOGICAL saltImplVertAdv
311          LOGICAL momImplVertAdv
312        LOGICAL multiDimAdvection        LOGICAL multiDimAdvection
313        LOGICAL forcing_In_AB        LOGICAL forcing_In_AB
314        LOGICAL doThetaClimRelax        LOGICAL doThetaClimRelax
# Line 304  C     groundAtK1  :: put the surface(k=1 Line 324  C     groundAtK1  :: put the surface(k=1
324        LOGICAL globalFiles        LOGICAL globalFiles
325        LOGICAL useSingleCpuIO        LOGICAL useSingleCpuIO
326        LOGICAL allowFreezing        LOGICAL allowFreezing
327          LOGICAL useOldFreezing
328        LOGICAL groundAtK1        LOGICAL groundAtK1
329        LOGICAL usePickupBeforeC35        LOGICAL usePickupBeforeC35
330        LOGICAL debugMode        LOGICAL debugMode
# Line 331  C           number of iterations for con Line 352  C           number of iterations for con
352  C           drop form 192 -> 134! Need to investigate this further!  C           drop form 192 -> 134! Need to investigate this further!
353  C           For now I have introduced a parameter cg2dpcOffDFac which  C           For now I have introduced a parameter cg2dpcOffDFac which
354  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  ).  
355  C     delR      :: Vertical grid spacing ( units of r ).  C     delR      :: Vertical grid spacing ( units of r ).
356  C     delRc     :: Vertical grid spacing between cell centers (r unit).  C     delRc     :: Vertical grid spacing between cell centers (r unit).
357  C     delX      :: Separation between cell faces (m) or (deg), depending  C     delX      :: Separation between cell faces (m) or (deg), depending
# Line 357  C     omega     :: Angular velocity ( ra Line 376  C     omega     :: Angular velocity ( ra
376  C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)  C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)
377  C     viscAh    :: Eddy viscosity coeff. for mixing of  C     viscAh    :: Eddy viscosity coeff. for mixing of
378  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 )  
379  C     viscAr    :: Eddy viscosity coeff. for mixing of  C     viscAr    :: Eddy viscosity coeff. for mixing of
380  C                 momentum vertically ( units of r^2/s )  C                 momentum vertically ( units of r^2/s )
381  C     viscA4    :: Biharmonic viscosity coeff. for mixing of  C     viscA4    :: Biharmonic viscosity coeff. for mixing of
382  C                 momentum laterally ( m^4/s )  C                 momentum laterally ( m^4/s )
383    C     viscAhMax :: Maximum eddy viscosity coeff. for mixing of
384    C                 momentum laterally ( m^2/s )
385    C     viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of
386    C                 momentum laterally ( m^4/s )
387    C     viscAhGrid:: non-dimensional grid-size dependent viscosity
388    C     viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity
389  C     diffKhT   :: Laplacian diffusion coeff. for mixing of  C     diffKhT   :: Laplacian diffusion coeff. for mixing of
390  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 )  
391  C     diffKrT   :: Laplacian diffusion coeff. for mixing of  C     diffKrT   :: Laplacian diffusion coeff. for mixing of
392  C                 heat vertically ( units of r^2/s )  C                 heat vertically ( units of r^2/s )
393  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of
394  C                 heat laterally ( m^4/s )  C                 heat laterally ( m^4/s )
395  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of
396  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 )  
397  C     diffKrS   :: Laplacian diffusion coeff. for mixing of  C     diffKrS   :: Laplacian diffusion coeff. for mixing of
398  C                 salt vertically ( units of r^2/s )  C                 salt vertically ( units of r^2/s )
399  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of
# Line 436  C     tauThetaClimRelax :: Relaxation to Line 449  C     tauThetaClimRelax :: Relaxation to
449  C     lambdaThetaClimRelax :: Inverse time scale for relaxation ( 1/s ).  C     lambdaThetaClimRelax :: Inverse time scale for relaxation ( 1/s ).
450  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).
451  C     lambdaSaltClimRelax :: Inverse time scale for relaxation ( 1/s ).  C     lambdaSaltClimRelax :: Inverse time scale for relaxation ( 1/s ).
452    C     latBandClimRelax :: latitude band where Relaxation to Clim. is applied,
453    C                         i.e. where |yC| <= latBandClimRelax
454  C     externForcingPeriod :: Is the period of which forcing varies (eg. 1 month)  C     externForcingPeriod :: Is the period of which forcing varies (eg. 1 month)
455  C     externForcingCycle :: Is the repeat time of the forcing (eg. 1 year)  C     externForcingCycle :: Is the repeat time of the forcing (eg. 1 year)
456  C                          (note: externForcingCycle must be an integer  C                          (note: externForcingCycle must be an integer
# Line 457  C     bottomDragLinear   :: Drag coeffic Line 472  C     bottomDragLinear   :: Drag coeffic
472  C      --"-"--  Quadratic  ( linear: 1/s, quadratic: 1/m )  C      --"-"--  Quadratic  ( linear: 1/s, quadratic: 1/m )
473        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
474       & cg2dpcOffDFac, cg3dTargetResidual,       & cg2dpcOffDFac, cg3dTargetResidual,
475       & delP, delZ, delR, delRc, delX, delY,       & delR, delRc, delX, delY,
476       & deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock,       & deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock,
477       & abeps, startTime,       & abeps, startTime,
478       & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,       & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,
479       & fCori, fCoriG, fCoriCos,       & fCori, fCoriG, fCoriCos,
480       & viscAh,  viscAz,  viscA4,  viscAr, viscAstrain, viscAtension,       & viscAh,  viscAhMax, viscAhGrid,
481       & diffKhT, diffKzT, diffK4T, diffKrT,       & viscA4,  viscA4Max, viscA4Grid,
482       & diffKhS, diffKzS, diffK4S, diffKrS,       & viscAr, viscAstrain, viscAtension,
483         & diffKhT, diffK4T, diffKrT,
484         & diffKhS, diffK4S, diffKrS,
485       & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,       & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,
486       & hFacMin, hFacMinDz, hFacInf, hFacSup,       & hFacMin, hFacMinDz, hFacInf, hFacSup,
487       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
# Line 476  C      --"-"--  Quadratic  ( linear: 1/s Line 493  C      --"-"--  Quadratic  ( linear: 1/s
493       & cosPower, cAdjFreq, omega, rotationPeriod,       & cosPower, cAdjFreq, omega, rotationPeriod,
494       & tauThetaClimRelax, lambdaThetaClimRelax,       & tauThetaClimRelax, lambdaThetaClimRelax,
495       & tauSaltClimRelax, lambdaSaltClimRelax,       & tauSaltClimRelax, lambdaSaltClimRelax,
496       & tauTr1ClimRelax, lambdaTr1ClimRelax,       & tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax,
497       & externForcingCycle, externForcingPeriod,       & externForcingCycle, externForcingPeriod,
498       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,
499       & viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp,       & hFacMinDr, hFacMinDp,
500       & horiVertRatio, recip_horiVertRatio,       & horiVertRatio, recip_horiVertRatio,
501       & ivdc_kappa, Ro_SeaLevel,       & ivdc_kappa, Ro_SeaLevel,
502       & bottomDragLinear,bottomDragQuadratic       & bottomDragLinear,bottomDragQuadratic
# Line 488  C      --"-"--  Quadratic  ( linear: 1/s Line 505  C      --"-"--  Quadratic  ( linear: 1/s
505        _RL cg2dTargetResWunit        _RL cg2dTargetResWunit
506        _RL cg3dTargetResidual        _RL cg3dTargetResidual
507        _RL cg2dpcOffDFac        _RL cg2dpcOffDFac
       _RL delZ(Nr)  
       _RL delP(Nr)  
508        _RL delR(Nr)        _RL delR(Nr)
509        _RL delRc(Nr+1)        _RL delRc(Nr+1)
510        _RL delX(Nx)        _RL delX(Nx)
# Line 516  C      --"-"--  Quadratic  ( linear: 1/s Line 531  C      --"-"--  Quadratic  ( linear: 1/s
531        _RL hFacSup        _RL hFacSup
532        _RL beta        _RL beta
533        _RL viscAh        _RL viscAh
534          _RL viscAhMax
535          _RL viscAhGrid
536        _RL viscAstrain        _RL viscAstrain
537        _RL viscAtension        _RL viscAtension
       _RL viscAz  
       _RL viscAp  
538        _RL viscAr        _RL viscAr
539        _RL viscA4        _RL viscA4
540          _RL viscA4Max
541          _RL viscA4Grid
542        _RL diffKhT        _RL diffKhT
543        _RL diffKrT        _RL diffKrT
       _RL diffKzT  
       _RL diffKpT  
544        _RL diffK4T        _RL diffK4T
545        _RL diffKhS        _RL diffKhS
546        _RL diffKrS        _RL diffKrS
       _RL diffKzS  
       _RL diffKpS  
547        _RL diffK4S        _RL diffK4S
548        _RL delt        _RL delt
549        _RL tauCD        _RL tauCD
# Line 574  C      --"-"--  Quadratic  ( linear: 1/s Line 587  C      --"-"--  Quadratic  ( linear: 1/s
587        _RL lambdaSaltClimRelax        _RL lambdaSaltClimRelax
588        _RL tauTr1ClimRelax        _RL tauTr1ClimRelax
589        _RL lambdaTr1ClimRelax        _RL lambdaTr1ClimRelax
590          _RL latBandClimRelax
591        _RL externForcingCycle        _RL externForcingCycle
592        _RL externForcingPeriod        _RL externForcingPeriod
593        _RL convertFW2Salt        _RL convertFW2Salt
# Line 598  C     atm_Po    :: standard reference pr Line 612  C     atm_Po    :: standard reference pr
612  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
613  C     atm_Rd    :: gas constant for dry air  C     atm_Rd    :: gas constant for dry air
614  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)
615    C     atm_Rq    :: water vapour specific volume anomaly relative to dry air
616    C                  (e.g. typical value = (29/18 -1) 10^-3 with q [g/kg])
617  C     integr_GeoPot :: option to select the way we integrate the geopotential  C     integr_GeoPot :: option to select the way we integrate the geopotential
618  C                     (still a subject of discussions ...)  C                     (still a subject of discussions ...)
619  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is
620  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)
621        COMMON /PARM_ATM/        COMMON /PARM_ATM/
622       &            celsius2K,       &            celsius2K,
623       &            atm_Cp, atm_Rd, atm_kappa, atm_Po,       &            atm_Cp, atm_Rd, atm_kappa, atm_Rq, atm_Po,
624       &            integr_GeoPot, selectFindRoSurf       &            integr_GeoPot, selectFindRoSurf
625        _RL celsius2K        _RL celsius2K
626        _RL atm_Po, atm_Cp, atm_Rd, atm_kappa        _RL atm_Po, atm_Cp, atm_Rd, atm_kappa, atm_Rq
627        INTEGER integr_GeoPot, selectFindRoSurf        INTEGER integr_GeoPot, selectFindRoSurf
628    
629  C Logical flags for selecting packages  C Logical flags for selecting packages
# Line 624  C Logical flags for selecting packages Line 640  C Logical flags for selecting packages
640        LOGICAL usePTRACERS        LOGICAL usePTRACERS
641        LOGICAL useSBO        LOGICAL useSBO
642        LOGICAL useSEAICE        LOGICAL useSEAICE
643  cswdblk -- add ---        LOGICAL useBulkForce
644        LOGICAL useBulkforce        LOGICAL useThSIce
645        LOGICAL useThermSEAICE        LOGICAL usefizhi
646  cswdblk --- end add ---        LOGICAL usegridalt
647          LOGICAL usediagnostics
648        COMMON /PARM_PACKAGES/        COMMON /PARM_PACKAGES/
649       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,
650       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
651       &        usePTRACERS,  useSBO, useSEAICE,       &        usePTRACERS,  useSBO, useSEAICE, useThSIce, useBulkForce,
652  cswdblk -- add ---       &        usefizhi,  usegridalt, usediagnostics
      &        useThermSEAICE, useBulkforce  
 cswdblk --- end add ---  
Legend:
Removed from v.1.91  
changed lines
  Added in v.1.102
  ViewVC Help
Powered by ViewVC 1.1.22