/[MITgcm]/MITgcm/model/inc/PARAMS.h

Diff of /MITgcm/model/inc/PARAMS.h

Parent Directory | Revision Log | View Revision Graph Revision Graph | View Patch Patch

-revision 1.212 by mlosch,
Tue Feb  5 13:32:18 2008 UTC
+revision 1.216 by jmc,
Sat Apr  5 18:09:30 2008 UTC
 Line 173 
 C     tempVertAdvScheme   :: Temp. Vert.
  C     saltAdvScheme       :: Salt. Horiz.advection scheme selector
  C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector
  C     selectKEscheme      :: Kinetic Energy scheme selector (Vector Inv.)
+ C     selectVortScheme    :: Scheme selector for Vorticity term (Vector Inv.)
  C     monitorSelect       :: select group of variables to monitor
  C                            =1 : dynvars ; =2 : + vort ; =3 : + surface
  C     debugLevel          :: debug level selector: higher -> more writing
-Line 189 
 C     debugLevel          :: debug level
+Line 190 
 C     debugLevel          :: debug level
       &        momForcingOutAB, tracForcingOutAB,
       &        tempAdvScheme, tempVertAdvScheme,
       &        saltAdvScheme, saltVertAdvScheme,
-      &        selectKEscheme,
+      &        selectKEscheme, selectVortScheme,
       &        monitorSelect, debugLevel
        INTEGER cg2dMaxIters
        INTEGER cg2dChkResFreq
-Line 208 
 C     debugLevel          :: debug level
+Line 209 
 C     debugLevel          :: debug level
        INTEGER tempAdvScheme, tempVertAdvScheme
        INTEGER saltAdvScheme, saltVertAdvScheme
        INTEGER selectKEscheme
+       INTEGER selectVortScheme
        INTEGER monitorSelect
        INTEGER debugLevel
-Line 220 
 C
+Line 222 
 C
        PARAMETER(debLevB=2)
  C--   COMMON /PARM_L/ Logical valued parameters used by the model.
+ C- Coordinate + Grid params:
+ C     fluidIsAir       :: Set to indicate that the fluid major constituent
+ C                        is air
+ C     fluidIsWater     :: Set to indicate that the fluid major constituent
+ C                        is water
+ C     usingPCoords     :: Set to indicate that we are working in a pressure
+ C                        type coordinate (p or p*).
+ C     usingZCoords     :: Set to indicate that we are working in a height
+ C                        type coordinate (z or z*)
+ C     useDynP_inEos_Zc :: use the dynamical pressure in EOS (with Z-coord.)
+ C                         this requires specific code for restart & exchange
  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian
  C                          coordinate frame.
  C     usingSphericalPolarGrid :: If TRUE grid generation will be in a
  C                               spherical polar frame.
- C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical
+ C     rotateGrid      :: rotate grid coordinates to geographical coordinates
+ C                        according to Euler angles phiEuler, thetaEuler, psiEuler
  C     usingCurvilinearGrid :: If TRUE, use a curvilinear grid (to be provided)
+ C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical
  C     deepAtmosphere :: deep model (drop the shallow-atmosphere approximation)
- C     no_slip_sides :: Impose "no-slip" at lateral boundaries.
+ C     setInterFDr    :: set Interface depth (put cell-Center at the middle)
+ C     setCenterDr    :: set cell-Center depth (put Interface at the middle)
+ C- Momentum params:
+ C     no_slip_sides  :: Impose "no-slip" at lateral boundaries.
  C     no_slip_bottom :: Impose "no-slip" at bottom boundary.
+ C     useFullLeith   :: Set to true to use full Leith viscosity(may be unstable
+ C                       on irregular grids)
+ C     useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms
+ C     useAreaViscLength :: Set to true to use old scaling for viscous lengths,
+ C                          e.g., L2=Raz.  May be preferable for cube sphere.
  C     momViscosity  :: Flag which turns momentum friction terms on and off.
  C     momAdvection  :: Flag which turns advection of momentum on and off.
  C     momForcing    :: Flag which turns external forcing of momentum on
- C                     and off.
+ C                      and off.
  C     momPressureForcing :: Flag which turns pressure term in momentum equation
  C                          on and off.
  C     metricTerms   :: Flag which turns metric terms on or off.
  C     useNHMTerms   :: If TRUE use non-hydrostatic metric terms.
  C     useCoriolis   :: Flag which turns the coriolis terms on and off.
- C     tempAdvection :: Flag which turns advection of temperature on
+ C     use3dCoriolis :: Turns the 3-D coriolis terms (in Omega.cos Phi) on - off
- C                     and off.
+ C     useConstantF  :: Coriolis parameter set to f0
- C     tempForcing   :: Flag which turns external forcing of temperature on
+ C     useBetaPlaneF :: Coriolis parameter set to f0 + beta.y
- C                     and off.
+ C     useSphereF    :: Coriolis parameter set to 2.omega.sin(phi)
- C     saltAdvection :: Flag which turns advection of salinity on
+ C     useCDscheme   :: use CD-scheme to calculate Coriolis terms.
- C                     and off.
+ C     vectorInvariantMomentum :: use Vector-Invariant form (mom_vecinv package)
- C     saltForcing   :: Flag which turns external forcing of salinity on
+ C                                (default = F = use mom_fluxform package)
- C                     and off.
+ C     useJamartWetPoints :: Use wet-point method for Coriolis (Jamart & Ozer 1986)
+ C     useJamartMomAdv :: Use wet-point method for V.I. non-linear term
+ C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term
+ C     highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection)
+ C     useAbsVorticity :: work with f+zeta in Coriolis terms
+ C     upwindShear        :: use 1rst order upwind interp. (V.I., vertical advection)
+ C     momStepping    :: Turns momentum equation time-stepping off
+ C- Temp. & Salt params:
+ C     tempStepping   :: Turns temperature equation time-stepping off
+ C     saltStepping   :: Turns salinity equation time-stepping off
+ C     tempAdvection  :: Flag which turns advection of temperature on and off.
+ C     tempIsActiveTr :: Pot.Temp. is a dynamically active tracer
+ C     tempForcing    :: Flag which turns external forcing of temperature on
+ C                       and off.
+ C     saltAdvection  :: Flag which turns advection of salinity on and off.
+ C     saltIsActiveTr :: Salinity  is a dynamically active tracer
+ C     saltForcing    :: Flag which turns external forcing of salinity on
+ C                       and off.
  C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux
  C                         as a real Fresh Water (=> changes the Sea Level)
  C                         if F, converts P+R-E to salt flux (no SL effect)
- C     useFullLeith   :: Set to true to use full Leith viscosity(may be unstable
+ C- Time-stepping params:
- C                       on irregular grids)
- C     useAreaViscLength :: Set to true to use old scaling for viscous
- C              lengths, e.g., L2=Raz.  May be preferable for cube sphere.
- C     useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms
  C     rigidLid            :: Set to true to use rigid lid
  C     implicitFreeSurface :: Set to true to use implicit free surface
  C     exactConserv        :: Set to true to conserve exactly the total Volume
  C     linFSConserveTr     :: Set to true to correct source/sink of tracer
  C                            at the surface due to Linear Free Surface
  C     uniformLin_PhiSurf  :: Set to true to use a uniform Bo_surf in the
- C                           linear relation Phi_surf = Bo_surf*eta
+ C                            linear relation Phi_surf = Bo_surf*eta
- C     use3Dsolver   :: set to true to use 3-D pressure solver
+ C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm
+ C     nonHydrostatic   :: Using non-hydrostatic algorithm
+ C     use3Dsolver      :: set to true to use 3-D pressure solver
  C     implicitIntGravWave :: treat Internal Gravity Wave implicitly
- C     staggerTimeStep :: enable a Stagger time stepping T,S Rho then U,V
+ C     staggerTimeStep   :: enable a Stagger time stepping U,V (& W) then T,S
- C     momStepping   :: Turns momentum equation time-stepping off
- C     tempStepping  :: Turns temperature equation time-stepping off
- C     saltStepping  :: Turns salinity equation time-stepping off
- C     useConstantF  :: Coriolis parameter set to f0
- C     useBetaPlaneF :: Coriolis parameter set to f0 + beta.y
- C     useSphereF    :: Coriolis parameter set to 2.omega.sin(phi)
- C     use3dCoriolis :: Turns the 3-D coriolis terms (in Omega.cos Phi) on - off
- C     useCDscheme   :: use CD-scheme to calculate Coriolis terms.
- C     useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986)
- C     useJamartMomAdv :: Use wet-point method for V.I. non-linear term
- C     SadournyCoriolis :: use the enstrophy conserving scheme by Sadourny
- C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term
- C     highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection)
- C     upwindShear        :: use 1rst order upwind interp. (V.I., vertical advection)
- C     useAbsVorticity :: work with f+zeta in Coriolis terms
  C     implicitDiffusion :: Turns implicit vertical diffusion on
  C     implicitViscosity :: Turns implicit vertical viscosity on
  C     tempImplVertAdv :: Turns on implicit vertical advection for Temperature
-Line 290 
 C     momDissip_In_AB   :: if False, put
+Line 313 
 C     momDissip_In_AB   :: if False, put
  C                          out off Adams-Bashforth time stepping.
  C     doAB_onGtGs       :: if the Adams-Bashforth time stepping is used, always
  C                          apply AB on tracer tendencies (rather than on Tracer)
- C     pickupStrictlyMatch :: check and stop if pickup-file do not stricly match
+ C- Other forcing params -
- C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup
+ C     balanceEmPmR    :: substract global mean of EmPmR at every time step
- C     usePickupBeforeC54 :: start from old-pickup files, generated with code from
+ C     balanceQnet     :: substract global mean of Qnet at every time step
- C                           before checkpoint-54a, Jul 06, 2004.
+ C     balancePrintMean:: print substracted global means to STDOUT
  C     doThetaClimRelax :: Set true if relaxation to temperature
  C                        climatology is required.
  C     doSaltClimRelax  :: Set true if relaxation to salinity
  C                        climatology is required.
+ C     allowFreezing  :: Allows surface water to freeze and form ice
+ C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)
  C     periodicExternalForcing :: Set true if forcing is time-dependant
- C     usingPCoords     :: Set to indicate that we are working in a pressure
+ C- I/O parameters -
- C                        type coordinate (p or p*).
- C     usingZCoords     :: Set to indicate that we are working in a height
- C                        type coordinate (z or z*)
- C     fluidIsAir       :: Set to indicate that the fluid major constituent
- C                        is air
- C     fluidIsWater     :: Set to indicate that the fluid major constituent
- C                        is water
- C     useDynP_inEos_Zc :: use the dynamical pressure in EOS (with Z-coord.)
- C                         this requires specific code for restart & exchange
- C     setInterFDr    :: set Interface depth (put cell-Center at the middle)
- C     setCenterDr    :: set cell-Center depth (put Interface at the middle)
- C     nonHydrostatic :: Using non-hydrostatic terms
- C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm
  C     globalFiles    :: Selects between "global" and "tiled" files
  C     useSingleCpuIO :: On SGI platforms, option globalFiles is either
  C                       slow (f77) or does not work (f90).  When
  C                       useSingleCpuIO is set, mdsio_writefield.F
  C                       outputs from master mpi process only.
- C     allowFreezing  :: Allows surface water to freeze and form ice
+ C     pickupStrictlyMatch :: check and stop if pickup-file do not stricly match
- C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)
+ C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup
+ C     usePickupBeforeC54 :: start from old-pickup files, generated with code from
+ C                           before checkpoint-54a, Jul 06, 2004.
  C     pickup_write_mdsio :: use mdsio to write pickups
  C     pickup_read_mdsio  :: use mdsio to read  pickups
  C     pickup_write_immed :: echo the pickup immediately (for conversion)
+ C     writePickupAtEnd   :: write pickup at the last timestep
  C     timeave_mdsio      :: use mdsio for timeave output
  C     snapshot_mdsio     :: use mdsio for "snapshot" (dumpfreq/diagfreq) output
  C     monitor_stdio      :: use stdio for monitor output
  C     dumpInitAndLast :: dumps model state to files at Initial (nIter0)
  C                        & Last iteration, in addition multiple of dumpFreq iter.
- C     balanceEmPmR    :: substract global mean of EmPmR at every time step
- C     balanceQnet     :: substract global mean of Qnet at every time step
- C     balancePrintMean:: print substracted global means to STDOUT
- C     rotateGrid      :: rotate grid coordinates to geographical coordinates
- C                        according to Euler angles phiEuler, thetaEuler, psiEuler
-       COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
+       COMMON /PARM_L/
+      & fluidIsAir, fluidIsWater,
+      & usingPCoords, usingZCoords, useDynP_inEos_Zc,
+      & usingCartesianGrid, usingSphericalPolarGrid, rotateGrid,
       & usingCurvilinearGrid, usingCylindricalGrid,
       & deepAtmosphere, setInterFDr, setCenterDr,
-      & no_slip_sides,no_slip_bottom,
+      & no_slip_sides, no_slip_bottom,
-      & momViscosity, momAdvection, momForcing, useCoriolis,
+      & useFullLeith, useStrainTensionVisc, useAreaViscLength,
-      & momPressureForcing, vectorInvariantMomentum,
+      & momViscosity, momAdvection, momForcing,
-      & tempAdvection, tempForcing,
+      & momPressureForcing, metricTerms, useNHMTerms,
-      & saltAdvection, saltForcing,
+      & useCoriolis, use3dCoriolis,
+      & useConstantF, useBetaPlaneF, useSphereF,
+      & useCDscheme, vectorInvariantMomentum,
+      & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,
+      & upwindVorticity, highOrderVorticity,
+      & useAbsVorticity, upwindShear,
+      & momStepping, tempStepping, saltStepping,
+      & tempAdvection, tempIsActiveTr, tempForcing,
+      & saltAdvection, saltIsActiveTr, saltForcing,
       & useRealFreshWaterFlux,
-      & useFullLeith, useStrainTensionVisc,
-      & useAreaViscLength,
       & rigidLid, implicitFreeSurface, exactConserv, linFSConserveTr,
       & uniformLin_PhiSurf,
+      & quasiHydrostatic, nonHydrostatic,
       & use3Dsolver, implicitIntGravWave, staggerTimeStep,
-      & momStepping, tempStepping, saltStepping,
-      & metricTerms, useNHMTerms,
-      & useConstantF, useBetaPlaneF, useSphereF, use3dCoriolis,
-      & useCDscheme,
-      & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,
-      & SadournyCoriolis, upwindVorticity, highOrderVorticity,
-      & useAbsVorticity, upwindShear,
       & implicitDiffusion, implicitViscosity,
       & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,
       & multiDimAdvection, useMultiDimAdvec,
       & momDissip_In_AB, doAB_onGtGs,
-      & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,
+      & balanceEmPmR, balanceQnet, balancePrintMean,
-      & periodicExternalForcing,
+      & doThetaClimRelax, doSaltClimRelax,
-      & fluidIsAir, fluidIsWater,
-      & usingPCoords, usingZCoords, useDynP_inEos_Zc,
-      & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,
       & allowFreezing, useOldFreezing,
+      & periodicExternalForcing,
+      & globalFiles, useSingleCpuIO,
       & pickupStrictlyMatch, usePickupBeforeC54, startFromPickupAB2,
       & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,
+      & writePickupAtEnd,
       & timeave_mdsio, snapshot_mdsio, monitor_stdio,
       & outputTypesInclusive, dumpInitAndLast, debugMode,
-      & inAdMode, inAdTrue, inAdFalse, inAdExact,
+      & inAdMode, inAdTrue, inAdFalse, inAdExact
-      & balanceEmPmR, balanceQnet, balancePrintMean,
-      & rotateGrid
+       LOGICAL fluidIsAir
+       LOGICAL fluidIsWater
+       LOGICAL usingPCoords
+       LOGICAL usingZCoords
+       LOGICAL useDynP_inEos_Zc
        LOGICAL usingCartesianGrid
-       LOGICAL usingSphericalPolarGrid
+       LOGICAL usingSphericalPolarGrid, rotateGrid
        LOGICAL usingCylindricalGrid
        LOGICAL usingCurvilinearGrid
        LOGICAL deepAtmosphere
-Line 390 
 C                        according to Eu
+Line 406 
 C                        according to Eu
        LOGICAL useCoriolis
        LOGICAL vectorInvariantMomentum
        LOGICAL tempAdvection
+       LOGICAL tempIsActiveTr
        LOGICAL tempForcing
        LOGICAL saltAdvection
+       LOGICAL saltIsActiveTr
        LOGICAL saltForcing
        LOGICAL useRealFreshWaterFlux
        LOGICAL useFullLeith
-Line 402 
 C                        according to Eu
+Line 420 
 C                        according to Eu
        LOGICAL exactConserv
        LOGICAL linFSConserveTr
        LOGICAL uniformLin_PhiSurf
+       LOGICAL quasiHydrostatic
+       LOGICAL nonHydrostatic
        LOGICAL use3Dsolver
        LOGICAL implicitIntGravWave
        LOGICAL staggerTimeStep
-Line 417 
 C                        according to Eu
+Line 437 
 C                        according to Eu
        LOGICAL useEnergyConservingCoriolis
        LOGICAL useJamartWetPoints
        LOGICAL useJamartMomAdv
-       LOGICAL SadournyCoriolis
        LOGICAL upwindVorticity
        LOGICAL highOrderVorticity
        LOGICAL useAbsVorticity
-Line 431 
 C                        according to Eu
+Line 450 
 C                        according to Eu
        LOGICAL useMultiDimAdvec
        LOGICAL momDissip_In_AB
        LOGICAL doAB_onGtGs
+       LOGICAL balanceEmPmR
+       LOGICAL balanceQnet
+       LOGICAL balancePrintMean
        LOGICAL doThetaClimRelax
        LOGICAL doSaltClimRelax
-       LOGICAL doTr1ClimRelax
+       LOGICAL allowFreezing
+       LOGICAL useOldFreezing
        LOGICAL periodicExternalForcing
-       LOGICAL fluidIsAir
-       LOGICAL fluidIsWater
-       LOGICAL usingPCoords
-       LOGICAL usingZCoords
-       LOGICAL useDynP_inEos_Zc
-       LOGICAL nonHydrostatic
-       LOGICAL quasiHydrostatic
        LOGICAL globalFiles
        LOGICAL useSingleCpuIO
-       LOGICAL allowFreezing
-       LOGICAL useOldFreezing
        LOGICAL pickupStrictlyMatch
        LOGICAL usePickupBeforeC54
        LOGICAL startFromPickupAB2
-       LOGICAL dumpInitAndLast
-       LOGICAL debugMode
        LOGICAL pickup_read_mdsio, pickup_write_mdsio
-       LOGICAL pickup_write_immed
+       LOGICAL pickup_write_immed, writePickupAtEnd
        LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio
        LOGICAL outputTypesInclusive
+       LOGICAL dumpInitAndLast
+       LOGICAL debugMode
        LOGICAL inAdMode, inAdTrue, inAdFalse, inAdExact
-       LOGICAL balanceEmPmR
-       LOGICAL balanceQnet
-       LOGICAL balancePrintMean
-       LOGICAL rotateGrid
  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
  C     cg2dTargetResidual
  C          :: Target residual for cg2d solver; no unit (RHS normalisation)
-Line 527 
 C                 ( Southern edge f for
+Line 535 
 C                 ( Southern edge f for
  C     beta      :: df/dy ( s^-1.m^-1 )
  C     omega     :: Angular velocity ( rad/s )
  C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)
+ C     viscAr    :: Eddy viscosity coeff. for mixing of
+ C                 momentum vertically ( units of r^2/s )
  C     viscAh    :: Eddy viscosity coeff. for mixing of
  C                 momentum laterally ( m^2/s )
  C     viscAhW   :: Eddy viscosity coeff. for mixing of vertical
  C                 momentum laterally, no effect for hydrostatic
  C                 model, defaults to viscAh if unset ( m^2/s )
  C                 Not used if variable horiz. viscosity is used.
- C     viscAr    :: Eddy viscosity coeff. for mixing of
- C                 momentum vertically ( units of r^2/s )
  C     viscA4    :: Biharmonic viscosity coeff. for mixing of
  C                 momentum laterally ( m^4/s )
  C     viscA4W   :: Biharmonic viscosity coeff. for mixing of vertical
-Line 549 
 C     viscA4D   :: Biharmonic viscosity
+Line 557 
 C     viscA4D   :: Biharmonic viscosity
  C                  (act on Divergence part) ( m^4/s )
  C     viscA4Z   :: Biharmonic viscosity coeff. for mixing of momentum laterally
  C                  (act on Vorticity  part) ( m^4/s )
- C     viscC2leith :: Leith non-dimensional viscosity factor (grad(vort))
+ C     viscC2leith  :: Leith non-dimensional viscosity factor (grad(vort))
  C     viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div))
- C     viscC2smag  :: Smagorinsky non-dimensional viscosity factor (harmonic)
+ C     viscC4leith  :: Leith non-dimensional viscosity factor (grad(vort))
- C     viscC4smag  :: Smagorinsky non-dimensional viscosity factor (biharmonic)
+ C     viscC4leithD :: Modified Leith non-dimensional viscosity factor (grad(div))
- C     viscAhMax :: Maximum eddy viscosity coeff. for mixing of
+ C     viscC2smag   :: Smagorinsky non-dimensional viscosity factor (harmonic)
- C                   momentum laterally ( m^2/s )
+ C     viscC4smag   :: Smagorinsky non-dimensional viscosity factor (biharmonic)
- C     viscAhReMax :: Maximum gridscale Reynolds number for eddy viscosity
+ C     viscAhMax    :: Maximum eddy viscosity coeff. for mixing of
- C                   coeff. for mixing of momentum laterally (non-dim)
+ C                    momentum laterally ( m^2/s )
+ C     viscAhReMax  :: Maximum gridscale Reynolds number for eddy viscosity
+ C                     coeff. for mixing of momentum laterally (non-dim)
+ C     viscAhGrid   :: non-dimensional grid-size dependent viscosity
  C     viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ...
  C     viscAhGridMin::  in terms of non-dimensional grid-size dependent visc.
- C     viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of
+ C     viscA4Max    :: Maximum biharmonic viscosity coeff. for mixing of
- C                 momentum laterally ( m^4/s )
+ C                     momentum laterally ( m^4/s )
- C     viscA4ReMax :: Maximum Gridscale Reynolds number for
+ C     viscA4ReMax  :: Maximum Gridscale Reynolds number for
- C                  biharmonic viscosity coeff. momentum laterally (non-dim)
+ C                     biharmonic viscosity coeff. momentum laterally (non-dim)
- C     viscAhGrid:: non-dimensional grid-size dependent viscosity
+ C     viscA4Grid   :: non-dimensional grid-size dependent bi-harmonic viscosity
- C     viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity
  C     viscA4GridMax:: maximum and minimum biharmonic viscosity coefficients ...
  C     viscA4GridMin::  in terms of non-dimensional grid-size dependent viscosity
- C     viscC4leith :: Leith non-dimensional viscosity factor (grad(vort))
- C     viscC4leithD :: Modified Leith non-dimensional viscosity factor (grad(div))
  C     diffKhT   :: Laplacian diffusion coeff. for mixing of
  C                 heat laterally ( m^2/s )
  C     diffKrNrT :: vertical profile of Laplacian diffusion coeff.
  C                 for mixing of heat vertically ( units of r^2/s )
  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of
  C                 heat laterally ( m^4/s )
  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of
  C                 salt laterally ( m^2/s )
  C     diffKrNrS :: vertical profile of Laplacian diffusion coeff.
  C                 for mixing of salt vertically ( units of r^2/s ),
  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of
  C                 salt laterally ( m^4/s )
  C     diffKrBL79surf :: T/S surface diffusivity (m^2/s) Bryan and Lewis, 1979
-Line 637 
 C     mtFacMom      :: Metric terms trac
+Line 645 
 C     mtFacMom      :: Metric terms trac
  C     cosPower      :: Power of cosine of latitude to multiply viscosity
  C     cAdjFreq      :: Frequency of convective adjustment
  C
  C     taveFreq      :: Frequency with which time-averaged model state
  C                      is written to post-processing files ( s ).
  C     tave_lastIter :: (for state variable only) fraction of the last time
  C                      step (of each taveFreq period) put in the time average.
  C                      (fraction for 1rst iter = 1 - tave_lastIter)
  C     tauThetaClimRelax :: Relaxation to climatology time scale ( s ).
  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).
-Line 652 
 C                          (note: extern
+Line 660 
 C                          (note: extern
  C                           number times externForcingPeriod)
  C     convertFW2Salt :: salinity, used to convert Fresh-Water Flux to Salt Flux
  C                       (use model surface (local) value if set to -1)
  C     temp_EvPrRn :: temperature of Rain & Evap.
  C     salt_EvPrRn :: salinity of Rain & Evap.
  C        (notes: a) tracer content of Rain/Evap only used if both
  C                     NonLin_FrSurf & useRealFreshWater are set.
  C                b) use model surface (local) value if set to UNSET_RL)
  C     hMixCrit    :: criteria for mixed-layer diagnostic
  C     ivdc_kappa  :: implicit vertical diffusivity for convection [m^2/s]
  C     Ro_SeaLevel :: standard position of Sea-Level in "R" coordinate, used as
  C                    starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)
  C     sideDragFactor     :: side-drag scaling factor (used only if no_slip_sides)
  C                           (default=2: full drag ; =1: gives half-slip BC)
  C     bottomDragLinear    :: Linear    bottom-drag coefficient (units of [r]/s)
  C     bottomDragQuadratic :: Quadratic bottom-drag coefficient (units of [r]/m)
-Line 685 
 C     psiEuler      :: Euler angle, rota
+Line 693 
 C     psiEuler      :: Euler angle, rota
       & viscC2leith, viscC2leithD,
       & viscC2smag, viscC4smag,
       & viscAhD, viscAhZ, viscA4D, viscA4Z,
-      & viscA4, viscA4W,
+      & viscA4, viscA4W, viscA4Max,
-      & viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin,
+      & viscA4Grid, viscA4GridMax, viscA4GridMin,
-      & viscAhRemax, viscA4Remax,
+      & viscAhReMax, viscA4ReMax,
       & viscC4leith, viscC4leithD, viscAr,
       & diffKhT, diffK4T, diffKrNrT,
       & diffKhS, diffK4S, diffKrNrS,
-Line 706 
 C     psiEuler      :: Euler angle, rota
+Line 714 
 C     psiEuler      :: Euler angle, rota
       & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq,
       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
       & cosPower, cAdjFreq, omega, rotationPeriod,
-      & tauThetaClimRelax, tauSaltClimRelax,
+      & tauThetaClimRelax, tauSaltClimRelax, latBandClimRelax,
-      & tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax,
       & externForcingCycle, externForcingPeriod,
       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn,
       & hFacMinDr, hFacMinDp,
-Line 746 
 C     psiEuler      :: Euler angle, rota
+Line 753 
 C     psiEuler      :: Euler angle, rota
        _RL hFacInf
        _RL hFacSup
        _RL beta
+       _RL viscAr
        _RL viscAh
        _RL viscAhW
        _RL viscAhD
        _RL viscAhZ
        _RL viscAhMax
        _RL viscAhReMax
-       _RL viscAhGrid
+       _RL viscAhGrid, viscAhGridMax, viscAhGridMin
-       _RL viscAhGridMax
-       _RL viscAhGridMin
        _RL viscC2leith
        _RL viscC2leithD
        _RL viscC2smag
-       _RL viscC4smag
-       _RL viscAr
        _RL viscA4
        _RL viscA4W
        _RL viscA4D
-Line 769 
 C     psiEuler      :: Euler angle, rota
+Line 773 
 C     psiEuler      :: Euler angle, rota
        _RL viscA4Grid, viscA4GridMax, viscA4GridMin
        _RL viscC4leith
        _RL viscC4leithD
+       _RL viscC4smag
        _RL diffKhT
        _RL diffKrNrT(Nr)
        _RL diffK4T
-Line 825 
 C     psiEuler      :: Euler angle, rota
+Line 830 
 C     psiEuler      :: Euler angle, rota
        _RL rotationPeriod
        _RL tauThetaClimRelax
        _RL tauSaltClimRelax
-       _RL tauTr1ClimRelax
-       _RL lambdaTr1ClimRelax
        _RL latBandClimRelax
        _RL externForcingCycle
        _RL externForcingPeriod
-Line 859 
 C     atm_kappa :: kappa = R/Cp (R: cons
+Line 862 
 C     atm_kappa :: kappa = R/Cp (R: cons
  C     atm_Rq    :: water vapour specific volume anomaly relative to dry air
  C                  (e.g. typical value = (29/18 -1) 10^-3 with q [g/kg])
  C     integr_GeoPot :: option to select the way we integrate the geopotential
  C                     (still a subject of discussions ...)
  C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is
  C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)
        COMMON /PARM_ATM/

 Legend:



Removed from v.1.212
 


changed lines


 
Added in v.1.216
 Legend:



Removed from v.1.212
 


changed lines


 
Added in v.1.216
-Removed from v.1.212
+Added in v.1.216

	ViewVC Help
Powered by ViewVC 1.1.22