/[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.93 by heimbach, Fri Jul 18 20:06:45 2003 UTC revision 1.110 by adcroft, Wed May 26 14:50:09 2004 UTC
# Line 215  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 237  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 255  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 301  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 316  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 343  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 369  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     viscC2leith :: Leith non-dimensional viscosity factor
384    C     viscAhMax :: Maximum eddy viscosity coeff. for mixing of
385    C                 momentum laterally ( m^2/s )
386    C     viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of
387    C                 momentum laterally ( m^4/s )
388    C     viscAhGrid:: non-dimensional grid-size dependent viscosity
389    C     viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity
390    C     viscC4leith :: Leith non-dimensional viscosity factor
391  C     diffKhT   :: Laplacian diffusion coeff. for mixing of  C     diffKhT   :: Laplacian diffusion coeff. for mixing of
392  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 )  
393  C     diffKrT   :: Laplacian diffusion coeff. for mixing of  C     diffKrT   :: Laplacian diffusion coeff. for mixing of
394  C                 heat vertically ( units of r^2/s )  C                 heat vertically ( units of r^2/s )
395  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of
396  C                 heat laterally ( m^4/s )  C                 heat laterally ( m^4/s )
397  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of
398  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 )  
399  C     diffKrS   :: Laplacian diffusion coeff. for mixing of  C     diffKrS   :: Laplacian diffusion coeff. for mixing of
400  C                 salt vertically ( units of r^2/s )  C                 salt vertically ( units of r^2/s )
401  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of
# Line 448  C     tauThetaClimRelax :: Relaxation to Line 451  C     tauThetaClimRelax :: Relaxation to
451  C     lambdaThetaClimRelax :: Inverse time scale for relaxation ( 1/s ).  C     lambdaThetaClimRelax :: Inverse time scale for relaxation ( 1/s ).
452  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).
453  C     lambdaSaltClimRelax :: Inverse time scale for relaxation ( 1/s ).  C     lambdaSaltClimRelax :: Inverse time scale for relaxation ( 1/s ).
454    C     latBandClimRelax :: latitude band where Relaxation to Clim. is applied,
455    C                         i.e. where |yC| <= latBandClimRelax
456  C     externForcingPeriod :: Is the period of which forcing varies (eg. 1 month)  C     externForcingPeriod :: Is the period of which forcing varies (eg. 1 month)
457  C     externForcingCycle :: Is the repeat time of the forcing (eg. 1 year)  C     externForcingCycle :: Is the repeat time of the forcing (eg. 1 year)
458  C                          (note: externForcingCycle must be an integer  C                          (note: externForcingCycle must be an integer
# Line 469  C     bottomDragLinear   :: Drag coeffic Line 474  C     bottomDragLinear   :: Drag coeffic
474  C      --"-"--  Quadratic  ( linear: 1/s, quadratic: 1/m )  C      --"-"--  Quadratic  ( linear: 1/s, quadratic: 1/m )
475        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
476       & cg2dpcOffDFac, cg3dTargetResidual,       & cg2dpcOffDFac, cg3dTargetResidual,
477       & delP, delZ, delR, delRc, delX, delY,       & delR, delRc, delX, delY,
478       & deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock,       & deltaT, deltaTmom, deltaTtracer, deltaTfreesurf, deltaTClock,
479       & abeps, startTime,       & abeps, startTime,
480       & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,       & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,
481       & fCori, fCoriG, fCoriCos,       & fCori, fCoriG, fCoriCos,
482       & viscAh,  viscAz,  viscA4,  viscAr, viscAstrain, viscAtension,       & viscAh,  viscAhMax, viscAhGrid, viscC2leith,
483       & diffKhT, diffKzT, diffK4T, diffKrT,       & viscA4,  viscA4Max, viscA4Grid, viscC4leith,
484       & diffKhS, diffKzS, diffK4S, diffKrS,       & viscAr, viscAstrain, viscAtension,
485         & diffKhT, diffK4T, diffKrT,
486         & diffKhS, diffK4S, diffKrS,
487       & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,       & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,
488       & hFacMin, hFacMinDz, hFacInf, hFacSup,       & hFacMin, hFacMinDz, hFacInf, hFacSup,
489       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
490       & recip_rhoConst, rhoConst,       & recip_rhoConst, rhoConst,
491       & rhoConstFresh, convertEmP2rUnit, tRef, sRef,       & rhoConstFresh, convertEmP2rUnit, tRef, sRef,
492       & endTime, chkPtFreq, pchkPtFreq, dumpFreq,       & endTime, chkPtFreq, pchkPtFreq, dumpFreq, adjDumpFreq,
493       & diagFreq, taveFreq, tave_lastIter, monitorFreq,       & diagFreq, taveFreq, tave_lastIter, monitorFreq,
494       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
495       & cosPower, cAdjFreq, omega, rotationPeriod,       & cosPower, cAdjFreq, omega, rotationPeriod,
496       & tauThetaClimRelax, lambdaThetaClimRelax,       & tauThetaClimRelax, lambdaThetaClimRelax,
497       & tauSaltClimRelax, lambdaSaltClimRelax,       & tauSaltClimRelax, lambdaSaltClimRelax,
498       & tauTr1ClimRelax, lambdaTr1ClimRelax,       & tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax,
499       & externForcingCycle, externForcingPeriod,       & externForcingCycle, externForcingPeriod,
500       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn, trac_EvPrRn,
501       & viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp,       & hFacMinDr, hFacMinDp,
502       & horiVertRatio, recip_horiVertRatio,       & horiVertRatio, recip_horiVertRatio,
503       & ivdc_kappa, Ro_SeaLevel,       & ivdc_kappa, Ro_SeaLevel,
504       & bottomDragLinear,bottomDragQuadratic       & bottomDragLinear,bottomDragQuadratic
# Line 500  C      --"-"--  Quadratic  ( linear: 1/s Line 507  C      --"-"--  Quadratic  ( linear: 1/s
507        _RL cg2dTargetResWunit        _RL cg2dTargetResWunit
508        _RL cg3dTargetResidual        _RL cg3dTargetResidual
509        _RL cg2dpcOffDFac        _RL cg2dpcOffDFac
       _RL delZ(Nr)  
       _RL delP(Nr)  
510        _RL delR(Nr)        _RL delR(Nr)
511        _RL delRc(Nr+1)        _RL delRc(Nr+1)
512        _RL delX(Nx)        _RL delX(Nx)
# Line 528  C      --"-"--  Quadratic  ( linear: 1/s Line 533  C      --"-"--  Quadratic  ( linear: 1/s
533        _RL hFacSup        _RL hFacSup
534        _RL beta        _RL beta
535        _RL viscAh        _RL viscAh
536          _RL viscAhMax
537          _RL viscAhGrid
538          _RL viscC2leith
539        _RL viscAstrain        _RL viscAstrain
540        _RL viscAtension        _RL viscAtension
       _RL viscAz  
       _RL viscAp  
541        _RL viscAr        _RL viscAr
542        _RL viscA4        _RL viscA4
543          _RL viscA4Max
544          _RL viscA4Grid
545          _RL viscC4leith
546        _RL diffKhT        _RL diffKhT
547        _RL diffKrT        _RL diffKrT
       _RL diffKzT  
       _RL diffKpT  
548        _RL diffK4T        _RL diffK4T
549        _RL diffKhS        _RL diffKhS
550        _RL diffKrS        _RL diffKrS
       _RL diffKzS  
       _RL diffKpS  
551        _RL diffK4S        _RL diffK4S
552        _RL delt        _RL delt
553        _RL tauCD        _RL tauCD
# Line 566  C      --"-"--  Quadratic  ( linear: 1/s Line 571  C      --"-"--  Quadratic  ( linear: 1/s
571        _RL chkPtFreq        _RL chkPtFreq
572        _RL pChkPtFreq        _RL pChkPtFreq
573        _RL dumpFreq        _RL dumpFreq
574          _RL adjDumpFreq
575        _RL diagFreq        _RL diagFreq
576        _RL taveFreq        _RL taveFreq
577        _RL tave_lastIter        _RL tave_lastIter
# Line 586  C      --"-"--  Quadratic  ( linear: 1/s Line 592  C      --"-"--  Quadratic  ( linear: 1/s
592        _RL lambdaSaltClimRelax        _RL lambdaSaltClimRelax
593        _RL tauTr1ClimRelax        _RL tauTr1ClimRelax
594        _RL lambdaTr1ClimRelax        _RL lambdaTr1ClimRelax
595          _RL latBandClimRelax
596        _RL externForcingCycle        _RL externForcingCycle
597        _RL externForcingPeriod        _RL externForcingPeriod
598        _RL convertFW2Salt        _RL convertFW2Salt
# Line 610  C     atm_Po    :: standard reference pr Line 617  C     atm_Po    :: standard reference pr
617  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
618  C     atm_Rd    :: gas constant for dry air  C     atm_Rd    :: gas constant for dry air
619  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)  C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)
620    C     atm_Rq    :: water vapour specific volume anomaly relative to dry air
621    C                  (e.g. typical value = (29/18 -1) 10^-3 with q [g/kg])
622  C     integr_GeoPot :: option to select the way we integrate the geopotential  C     integr_GeoPot :: option to select the way we integrate the geopotential
623  C                     (still a subject of discussions ...)  C                     (still a subject of discussions ...)
624  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is
625  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)
626        COMMON /PARM_ATM/        COMMON /PARM_ATM/
627       &            celsius2K,       &            celsius2K,
628       &            atm_Cp, atm_Rd, atm_kappa, atm_Po,       &            atm_Cp, atm_Rd, atm_kappa, atm_Rq, atm_Po,
629       &            integr_GeoPot, selectFindRoSurf       &            integr_GeoPot, selectFindRoSurf
630        _RL celsius2K        _RL celsius2K
631        _RL atm_Po, atm_Cp, atm_Rd, atm_kappa        _RL atm_Po, atm_Cp, atm_Rd, atm_kappa, atm_Rq
632        INTEGER integr_GeoPot, selectFindRoSurf        INTEGER integr_GeoPot, selectFindRoSurf
633    
634  C Logical flags for selecting packages  C Logical flags for selecting packages
# Line 636  C Logical flags for selecting packages Line 645  C Logical flags for selecting packages
645        LOGICAL usePTRACERS        LOGICAL usePTRACERS
646        LOGICAL useSBO        LOGICAL useSBO
647        LOGICAL useSEAICE        LOGICAL useSEAICE
648        LOGICAL useBulkforce        LOGICAL useBulkForce
649        LOGICAL useThermSEAICE        LOGICAL useThSIce
650          LOGICAL usefizhi
651          LOGICAL usegridalt
652          LOGICAL usediagnostics
653          LOGICAL useEBM
654        COMMON /PARM_PACKAGES/        COMMON /PARM_PACKAGES/
655       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,
656       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
657       &        usePTRACERS,  useSBO, useSEAICE,       &        usePTRACERS,  useSBO, useSEAICE, useThSIce, useBulkForce,
658       &        useThermSEAICE, useBulkforce       &        usefizhi,  usegridalt, usediagnostics, useEBM
659    
660    C     Run-time flags for early-initialization of MNC
661          LOGICAL
662         &     useMNC,
663         &     mnc_use_indir, mnc_use_outdir, mnc_outdir_date,
664         &     mnc_echo_gvtypes, mnc_pickup_write, mnc_pickup_read,
665         &     mnc_use_for_mon
666          CHARACTER*(MAX_LEN_FNAM) mnc_outdir_str
667          CHARACTER*(MAX_LEN_FNAM) mnc_indir_str
668          COMMON /PARM_MNC/
669         &     useMNC,
670         &     mnc_use_indir, mnc_use_outdir, mnc_outdir_date,
671         &     mnc_echo_gvtypes, mnc_pickup_write, mnc_pickup_read,
672         &     mnc_outdir_str, mnc_indir_str,
673         &     mnc_use_for_mon
674    
675    
676    CEH3 ;;; Local Variables: ***
677    CEH3 ;;; mode:fortran ***
678    CEH3 ;;; End: ***
Legend:
Removed from v.1.93  
changed lines
  Added in v.1.110
  ViewVC Help
Powered by ViewVC 1.1.22