/[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.214 by jmc,
Sun Mar 30 21:44:09 2008 UTC
+revision 1.240 by mlosch,
Sat Sep 25 23:09:54 2010 UTC
 Line 24 
 C     Useful physical values
        Real*8 deg2rad
        PARAMETER ( deg2rad = 2.D0*PI/360.D0           )
- C     Symbolic values
- C     precXXXX :: Used to indicate what precision to use for
- C                dumping model state.
-       INTEGER precFloat32
-       PARAMETER ( precFloat32 = 32 )
-       INTEGER precFloat64
-       PARAMETER ( precFloat64 = 64 )
- C     UNSET_xxx :: Used to indicate variables that have not been given a value
-       Real*8 UNSET_FLOAT8
-       PARAMETER ( UNSET_FLOAT8 = 1.234567D5 )
-       Real*4 UNSET_FLOAT4
-       PARAMETER ( UNSET_FLOAT4 = 1.234567E5 )
-       _RL    UNSET_RL
-       PARAMETER ( UNSET_RL     = 1.234567D5 )
-       _RS    UNSET_RS
-       PARAMETER ( UNSET_RS     = 1.234567E5 )
-       INTEGER UNSET_I
-       PARAMETER ( UNSET_I      = 123456789  )
  C--   COMMON /PARM_C/ Character valued parameters used by the model.
  C     buoyancyRelation :: Flag used to indicate which relation to use to
  C                         get buoyancy.
-Line 57 
 C     sRefFile      :: File containing r
+Line 38 
 C     sRefFile      :: File containing r
  C     rhoRefFile    :: File containing reference density profile rhoRef (1.D)
  C     delRFile      :: File containing vertical grid spacing delR  (1.D array)
  C     delRcFile     :: File containing vertical grid spacing delRc (1.D array)
+ C     hybSigmFile   :: File containing hybrid-sigma vertical coord. coeff. (2x 1.D)
  C     delXFile      :: File containing X-spacing grid definition (1.D array)
  C     delYFile      :: File containing Y-spacing grid definition (1.D array)
  C     horizGridFile :: File containing horizontal-grid definition
-Line 65 
 C     bathyFile       :: File containing
+Line 47 
 C     bathyFile       :: File containing
  C                        is taken from inline function.
  C     topoFile        :: File containing the topography of the surface (unit=m)
  C                        (mainly used for the atmosphere = ground height).
- C     shelfIceFile    :: File containing the topography of the shelfice draught
- C                        (unit=m)
  C     hydrogThetaFile :: File containing initial hydrographic data (3-D)
  C                        for potential temperature.
  C     hydrogSaltFile  :: File containing initial hydrographic data (3-D)
  C                        for salinity.
  C     diffKrFile      :: File containing 3D specification of vertical diffusivity
+ C     viscAhDfile     :: File containing 3D specification of horizontal viscosity
+ C     viscAhZfile     :: File containing 3D specification of horizontal viscosity
+ C     viscA4Dfile     :: File containing 3D specification of horizontal viscosity
+ C     viscA4Zfile     :: File containing 3D specification of horizontal viscosity
  C     zonalWindFile   :: File containing zonal wind data
  C     meridWindFile   :: File containing meridional wind data
  C     thetaClimFile   :: File containing surface theta climataology used
-Line 89 
 C                 m/s when using externa
+Line 73 
 C                 m/s when using externa
  C                 to kg/m2/s by multiplying by rhoConstFresh.
  C     saltFluxFile    :: File containing surface salt flux
  C     pLoadFile       :: File containing pressure loading
- C     eddyTauxFile    :: File containing zonal Eddy stress data
+ C     eddyPsiXFile    :: File containing zonal Eddy streamfunction data
- C     eddyTauyFile    :: File containing meridional Eddy stress data
+ C     eddyPsiYFile    :: File containing meridional Eddy streamfunction data
  C     the_run_name    :: string identifying the name of the model "run"
        COMMON /PARM_C/
       &                buoyancyRelation, eosType,
       &                pickupSuff, mdsioLocalDir,
       &                tRefFile, sRefFile, rhoRefFile,
-      &                delRFile, delRcFile,
+      &                delRFile, delRcFile, hybSigmFile,
       &                delXFile, delYFile, horizGridFile,
-      &                bathyFile, topoFile, shelfIceFile,
+      &                bathyFile, topoFile,
+      &                viscAhDfile, viscAhZfile,
+      &                viscA4Dfile, viscA4Zfile,
       &                hydrogThetaFile, hydrogSaltFile, diffKrFile,
       &                zonalWindFile, meridWindFile, thetaClimFile,
       &                saltClimFile,
-Line 107 
 C     the_run_name    :: string identify
+Line 93 
 C     the_run_name    :: string identify
       &                lambdaThetaFile, lambdaSaltFile,
       &                uVelInitFile, vVelInitFile, pSurfInitFile,
       &                dQdTfile, ploadFile,
-      &                eddyTauxFile, eddyTauyFile,
+      &                eddyPsiXFile, eddyPsiYFile,
       &                the_run_name
        CHARACTER*(MAX_LEN_FNAM) buoyancyRelation
        CHARACTER*(6)  eosType
-Line 118 
 C     the_run_name    :: string identify
+Line 104 
 C     the_run_name    :: string identify
        CHARACTER*(MAX_LEN_FNAM) rhoRefFile
        CHARACTER*(MAX_LEN_FNAM) delRFile
        CHARACTER*(MAX_LEN_FNAM) delRcFile
+       CHARACTER*(MAX_LEN_FNAM) hybSigmFile
        CHARACTER*(MAX_LEN_FNAM) delXFile
        CHARACTER*(MAX_LEN_FNAM) delYFile
        CHARACTER*(MAX_LEN_FNAM) horizGridFile
-       CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile, shelfIceFile
+       CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile
        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile, hydrogSaltFile
        CHARACTER*(MAX_LEN_FNAM) diffKrFile
+       CHARACTER*(MAX_LEN_FNAM) viscAhDfile
+       CHARACTER*(MAX_LEN_FNAM) viscAhZfile
+       CHARACTER*(MAX_LEN_FNAM) viscA4Dfile
+       CHARACTER*(MAX_LEN_FNAM) viscA4Zfile
        CHARACTER*(MAX_LEN_FNAM) zonalWindFile
        CHARACTER*(MAX_LEN_FNAM) meridWindFile
        CHARACTER*(MAX_LEN_FNAM) thetaClimFile
-Line 138 
 C     the_run_name    :: string identify
+Line 129 
 C     the_run_name    :: string identify
        CHARACTER*(MAX_LEN_FNAM) pSurfInitFile
        CHARACTER*(MAX_LEN_FNAM) dQdTfile
        CHARACTER*(MAX_LEN_FNAM) ploadFile
-       CHARACTER*(MAX_LEN_FNAM) eddyTauxFile
+       CHARACTER*(MAX_LEN_FNAM) eddyPsiXFile
-       CHARACTER*(MAX_LEN_FNAM) eddyTauyFile
+       CHARACTER*(MAX_LEN_FNAM) eddyPsiYFile
        CHARACTER*(MAX_LEN_FNAM) lambdaThetaFile
        CHARACTER*(MAX_LEN_FNAM) lambdaSaltFile
        CHARACTER*(MAX_LEN_PREC/2) the_run_name
-Line 164 
 C     nonlinFreeSurf      :: option rela
+Line 155 
 C     nonlinFreeSurf      :: option rela
  C                           =0 Linear free surface ; >0 Non-linear
  C     select_rStar        :: option related to r* vertical coordinate
  C                           =0 (default) use r coord. ; > 0 use r*
+ C     selectNHfreeSurf    :: option for Non-Hydrostatic (free-)Surface formulation:
+ C                           =0 (default) hydrostatic surf. ; > 0 add NH effects.
+ C     selectSigmaCoord    :: option related to sigma vertical coordinate
+ C     selectAddFluid      :: option to add mass source/sink of fluid in the interior
+ C                            (3-D generalisation of oceanic real-fresh water flux)
+ C                           =0 off ; =1 add fluid ; =-1 virtual flux (no mass added)
  C     momForcingOutAB     :: =1: take momentum forcing contribution
  C                           out of (=0: in) Adams-Bashforth time stepping.
  C     tracForcingOutAB    :: =1: take tracer (Temp,Salt,pTracers) forcing contribution
-Line 187 
 C     debugLevel          :: debug level
+Line 184 
 C     debugLevel          :: debug level
       &        writeStatePrec,
       &        writeBinaryPrec, readBinaryPrec,
       &        nonlinFreeSurf, select_rStar,
+      &        selectNHfreeSurf,
+      &        selectSigmaCoord,
+      &        selectAddFluid,
       &        momForcingOutAB, tracForcingOutAB,
       &        tempAdvScheme, tempVertAdvScheme,
       &        saltAdvScheme, saltVertAdvScheme,
 Line 205 
 C     debugLevel          :: debug level
        INTEGER readBinaryPrec
        INTEGER nonlinFreeSurf
        INTEGER select_rStar
+       INTEGER selectNHfreeSurf
+       INTEGER selectSigmaCoord
+       INTEGER selectAddFluid
        INTEGER momForcingOutAB, tracForcingOutAB
        INTEGER tempAdvScheme, tempVertAdvScheme
        INTEGER saltAdvScheme, saltVertAdvScheme
-Line 213 
 C     debugLevel          :: debug level
+Line 216 
 C     debugLevel          :: debug level
        INTEGER monitorSelect
        INTEGER debugLevel
- C
-       INTEGER debLevZero
-       PARAMETER(debLevZero=0)
-       INTEGER debLevA
-       PARAMETER(debLevA=1)
-       INTEGER debLevB
-       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     maskIniTemp    :: apply mask to initial Pot.Temp.
+ C     maskIniSalt    :: apply mask to initial salinity
+ C     checkIniTemp   :: check for points with identically zero initial Pot.Temp.
+ C     checkIniSalt   :: check for points with identically zero initial salinity
  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- Pressure solver related parameters (PARM02)
- C                       on irregular grids)
+ C     useSRCGSolver  :: Set to true to use conjugate gradient
- C     useAreaViscLength :: Set to true to use old scaling for viscous
+ C                       solver with single reduction (only one call of
- C              lengths, e.g., L2=Raz.  May be preferable for cube sphere.
+ C                       s/r mpi_allreduce), default is false
- C     useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms
+ C- Time-stepping params:
  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     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 291 
 C     momDissip_In_AB   :: if False, put
+Line 316 
 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     periodicExternalForcing :: Set true if forcing is time-dependant
- 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     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     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)
+ C     periodicExternalForcing :: Set true if forcing is time-dependant
+ C- I/O parameters -
+ C     globalFiles    :: Selects between "global" and "tiled" files.
+ C                       On some platforms with MPI, option globalFiles is either
+ C                       slow or does not work. Use useSingleCpuIO instead.
+ C     useSingleCpuIO :: moved to EEPARAMS.h
+ C     pickupStrictlyMatch :: check and stop if pickup-file do not stricly match
+ 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)
-Line 330 
 C     snapshot_mdsio     :: use mdsio fo
+Line 345 
 C     snapshot_mdsio     :: use mdsio fo
  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,
+      & maskIniTemp, maskIniSalt, checkIniTemp, checkIniSalt,
       & 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,
-      & 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,
       & 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,
-      & balanceEmPmR, balanceQnet, balancePrintMean,
+      & useSRCGSolver
-      & 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 393 
 C                        according to Eu
+Line 410 
 C                        according to Eu
        LOGICAL useCoriolis
        LOGICAL vectorInvariantMomentum
        LOGICAL tempAdvection
+       LOGICAL tempIsActiveTr
        LOGICAL tempForcing
        LOGICAL saltAdvection
+       LOGICAL saltIsActiveTr
        LOGICAL saltForcing
+       LOGICAL maskIniTemp
+       LOGICAL maskIniSalt
+       LOGICAL checkIniTemp
+       LOGICAL checkIniSalt
        LOGICAL useRealFreshWaterFlux
        LOGICAL useFullLeith
        LOGICAL useStrainTensionVisc
-Line 405 
 C                        according to Eu
+Line 428 
 C                        according to Eu
        LOGICAL exactConserv
        LOGICAL linFSConserveTr
        LOGICAL uniformLin_PhiSurf
+       LOGICAL quasiHydrostatic
+       LOGICAL nonHydrostatic
        LOGICAL use3Dsolver
        LOGICAL implicitIntGravWave
        LOGICAL staggerTimeStep
-Line 433 
 C                        according to Eu
+Line 458 
 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 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 periodicExternalForcing
+       LOGICAL globalFiles
        LOGICAL pickupStrictlyMatch
        LOGICAL usePickupBeforeC54
        LOGICAL startFromPickupAB2
-       LOGICAL dumpInitAndLast
-       LOGICAL debugMode
        LOGICAL pickup_read_mdsio, pickup_write_mdsio
        LOGICAL pickup_write_immed, writePickupAtEnd
        LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio
        LOGICAL outputTypesInclusive
+       LOGICAL dumpInitAndLast
+       LOGICAL debugMode
        LOGICAL inAdMode, inAdTrue, inAdFalse, inAdExact
+       LOGICAL useSRCGSolver
-       LOGICAL balanceEmPmR
-       LOGICAL balanceQnet
-       LOGICAL balancePrintMean
-       LOGICAL rotateGrid
  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
  C     cg2dTargetResidual
-Line 498 
 C     rhoConst  :: Vertically constant r
+Line 512 
 C     rhoConst  :: Vertically constant r
  C     rhoFacC   :: normalized (by rhoConst) reference density at cell-Center
  C     rhoFacF   :: normalized (by rhoConst) reference density at cell-interFace
  C     rhoConstFresh :: Constant reference density for fresh water (rain)
+ C     rho1Ref   :: reference vertical profile for density
  C     tRef      :: reference vertical profile for potential temperature
  C     sRef      :: reference vertical profile for salinity/specific humidity
  C     phiRef    :: reference potential (pressure/rho, geopotential) profile
-Line 518 
 C     rUnit2mass :: units conversion fac
+Line 533 
 C     rUnit2mass :: units conversion fac
  C                :: from vertical r-coordinate unit to mass per unit area [kg/m2].
  C                :: z-coord: = rhoConst  ( [m] * rho = [kg/m2] ) ;
  C                :: p-coord: = 1/gravity ( [Pa] /  g = [kg/m2] ) ;
- C     phiMin    :: Latitude of southern most cell face.
+ C     xgOrigin   :: Origin of the X-axis (Cartesian Grid) / Longitude of Western
- C     thetaMin  :: Longitude of western most cell face (this
+ C                :: most cell face (Lat-Lon grid) (Note: this is an "inert"
- C                 is an "inert" parameter but it is included
+ C                :: parameter but it makes geographical references simple.)
- C                 to make geographical references simple.)
+ C     ygOrigin   :: Origin of the Y-axis (Cartesian Grid) / Latitude of Southern
- C     rSphere   :: Radius of sphere for a spherical polar grid ( m ).
+ C                :: most face (Lat-Lon grid).
+ C     rSphere    :: Radius of sphere for a spherical polar grid ( m ).
  C     recip_rSphere  :: Reciprocal radius of sphere ( m ).
- C     f0        :: Reference coriolis parameter ( 1/s )
+ C     f0         :: Reference coriolis parameter ( 1/s )
- C                 ( Southern edge f for beta plane )
+ C                   ( Southern edge f for beta plane )
- C     beta      :: df/dy ( s^-1.m^-1 )
+ C     beta       :: df/dy ( s^-1.m^-1 )
- C     omega     :: Angular velocity ( rad/s )
+ C     omega      :: Angular velocity ( rad/s )
  C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)
- C     viscAr    :: Eddy viscosity coeff. for mixing of
+ C     viscArNr   :: vertical profile of Eddy viscosity coeff.
- C                 momentum vertically ( units of r^2/s )
+ C                   for vertical mixing of momentum ( units of r^2/s )
- C     viscAh    :: Eddy viscosity coeff. for mixing of
+ C     viscAh     :: Eddy viscosity coeff. for mixing of
- C                 momentum laterally ( m^2/s )
+ C                   momentum laterally ( m^2/s )
- C     viscAhW   :: Eddy viscosity coeff. for mixing of vertical
+ C     viscAhW    :: Eddy viscosity coeff. for mixing of vertical
- C                 momentum laterally, no effect for hydrostatic
+ C                   momentum laterally, no effect for hydrostatic
- C                 model, defaults to viscAh if unset ( m^2/s )
+ C                   model, defaults to viscAh if unset ( m^2/s )
- C                 Not used if variable horiz. viscosity is used.
+ C                   Not used if variable horiz. viscosity is used.
- C     viscA4    :: Biharmonic viscosity coeff. for mixing of
+ C     viscA4     :: Biharmonic viscosity coeff. for mixing of
- C                 momentum laterally ( m^4/s )
+ C                   momentum laterally ( m^4/s )
- C     viscA4W   :: Biharmonic viscosity coeff. for mixing of vertical
+ C     viscA4W    :: Biharmonic viscosity coeff. for mixing of vertical
- C                 momentum laterally, no effect for hydrostatic
+ C                   momentum laterally, no effect for hydrostatic
- C                 model, defaults to viscA4 if unset ( m^2/s )
+ C                   model, defaults to viscA4 if unset ( m^2/s )
- C                 Not used if variable horiz. viscosity is used.
+ C                   Not used if variable horiz. viscosity is used.
- C     viscAhD   :: Eddy viscosity coeff. for mixing of momentum laterally
+ C     viscAhD    :: Eddy viscosity coeff. for mixing of momentum laterally
- C                  (act on Divergence part) ( m^2/s )
+ C                   (act on Divergence part) ( m^2/s )
- C     viscAhZ   :: Eddy viscosity coeff. for mixing of momentum laterally
+ C     viscAhZ    :: Eddy viscosity coeff. for mixing of momentum laterally
- C                  (act on Vorticity  part) ( m^2/s )
+ C                   (act on Vorticity  part) ( m^2/s )
- C     viscA4D   :: Biharmonic viscosity coeff. for mixing of momentum laterally
+ C     viscA4D    :: Biharmonic viscosity coeff. for mixing of momentum laterally
- C                  (act on Divergence part) ( m^4/s )
+ C                   (act on Divergence part) ( m^4/s )
- C     viscA4Z   :: Biharmonic viscosity coeff. for mixing of momentum laterally
+ C     viscA4Z    :: Biharmonic viscosity coeff. for mixing of momentum laterally
- C                  (act on Vorticity  part) ( m^4/s )
+ C                   (act on Vorticity  part) ( m^4/s )
  C     viscC2leith  :: Leith non-dimensional viscosity factor (grad(vort))
  C     viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div))
  C     viscC4leith  :: Leith non-dimensional viscosity factor (grad(vort))
-Line 613 
 C     implicSurfPress :: parameter of th
+Line 629 
 C     implicSurfPress :: parameter of th
  C                     Implicit part of Surface Pressure Gradient ( 0-1 )
  C     implicDiv2Dflow :: parameter of the Crank-Nickelson time stepping :
  C                     Implicit part of barotropic flow Divergence ( 0-1 )
+ C     implicitNHPress :: parameter of the Crank-Nickelson time stepping :
+ C                     Implicit part of Non-Hydrostatic Pressure Gradient ( 0-1 )
  C     hFacMin      :: Minimum fraction size of a cell (affects hFacC etc...)
- C     hFacMinDz    :: Minimum dimesional size of a cell (affects hFacC etc..., m)
+ C     hFacMinDz    :: Minimum dimensional size of a cell (affects hFacC etc..., m)
- C     hFacMinDp    :: Minimum dimesional size of a cell (affects hFacC etc..., Pa)
+ C     hFacMinDp    :: Minimum dimensional size of a cell (affects hFacC etc..., Pa)
- C     hFacMinDr    :: Minimum dimesional size of a cell (affects hFacC etc..., units of r)
+ C     hFacMinDr    :: Minimum dimensional size of a cell (-> hFacC etc..., r units)
  C     hFacInf      :: Threshold (inf and sup) for fraction size of surface cell
- C     hFacSup        that control vanishing and creating levels
+ C     hFacSup          that control vanishing and creating levels
- C     tauCD        :: CD scheme coupling timescale ( 1/s )
+ C     tauCD         :: CD scheme coupling timescale ( s )
- C     rCD          :: CD scheme normalised coupling parameter ( 0-1 )
+ C     rCD           :: CD scheme normalised coupling parameter (= 1 - deltaT/tauCD)
+ C     epsAB_CD      :: Adams-Bashforth-2 stabilizing weight used in CD scheme
  C     baseTime      :: model base time (time origin) = time @ iteration zero
  C     startTime     :: Starting time for this integration ( s ).
  C     endTime       :: Ending time for this integration ( s ).
-Line 656 
 C     convertFW2Salt :: salinity, used t
+Line 675 
 C     convertFW2Salt :: salinity, used t
  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     temp_addMass :: temperature of addMass array
+ C     salt_addMass :: salinity of addMass array
  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     hMixCriteria:: criteria for mixed-layer diagnostic
+ C     dRhoSmall   :: parameter for mixed-layer diagnostic
+ C     hMixSmooth  :: Smoothing parameter for mixed-layer diag (default=0=no smoothing)
  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     rSigmaBnd   :: vertical position (in r-unit) of r/sigma transition (Hybrid-Sigma)
  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)
-Line 681 
 C     psiEuler      :: Euler angle, rota
+Line 705 
 C     psiEuler      :: Euler angle, rota
       & delR, delRc, delX, delY,
       & deltaT, deltaTmom, dTtracerLev, deltaTfreesurf, deltaTClock,
       & abEps, alph_AB, beta_AB,
-      & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,
+      & xgOrigin, ygOrigin, rSphere, recip_RSphere, f0, beta,
       & viscAh, viscAhW, viscAhMax,
       & viscAhGrid, viscAhGridMax, viscAhGridMin,
       & viscC2leith, viscC2leithD,
-Line 690 
 C     psiEuler      :: Euler angle, rota
+Line 714 
 C     psiEuler      :: Euler angle, rota
       & viscA4, viscA4W, viscA4Max,
       & viscA4Grid, viscA4GridMax, viscA4GridMin,
       & viscAhReMax, viscA4ReMax,
-      & viscC4leith, viscC4leithD, viscAr,
+      & viscC4leith, viscC4leithD, viscArNr,
       & diffKhT, diffK4T, diffKrNrT,
       & diffKhS, diffK4S, diffKrNrS,
       & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,
       & BL79LatVary,
       & diffKrBLEQsurf, diffKrBLEQdeep, diffKrBLEQscl, diffKrBLEQHo,
-      & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,
+      & tauCD, rCD, epsAB_CD,
+      & freeSurfFac, implicSurfPress, implicDiv2Dflow, implicitNHPress,
       & hFacMin, hFacMinDz, hFacInf, hFacSup,
       & gravity, recip_gravity, gBaro,
       & rhonil, recip_rhonil, rhoConst, recip_rhoConst,
       & rhoFacC, recip_rhoFacC, rhoFacF, recip_rhoFacF,
-      & rhoConstFresh, tRef, sRef, phiRef, dBdrRef,
+      & rhoConstFresh, rho1Ref, tRef, sRef, phiRef, dBdrRef,
       & rVel2wUnit, wUnit2rVel, mass2rUnit, rUnit2mass,
       & baseTime, startTime, endTime,
       & chkPtFreq, pChkPtFreq, dumpFreq, adjDumpFreq,
       & 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,
+      & temp_addMass, salt_addMass, hFacMinDr, hFacMinDp,
-      & ivdc_kappa, hMixCriteria, Ro_SeaLevel,
+      & ivdc_kappa, hMixCriteria, dRhoSmall, hMixSmooth,
+      & Ro_SeaLevel, rSigmaBnd,
       & sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2,
       & smoothAbsFuncRange,
       & tCylIn, tCylOut,
-Line 733 
 C     psiEuler      :: Euler angle, rota
+Line 758 
 C     psiEuler      :: Euler angle, rota
        _RL dTtracerLev(Nr)
        _RL deltaTfreesurf
        _RL abEps, alph_AB, beta_AB
-       _RL phiMin
+       _RL xgOrigin
-       _RL thetaMin
+       _RL ygOrigin
        _RL rSphere
        _RL recip_rSphere
        _RL f0
        _RL freeSurfFac
        _RL implicSurfPress
        _RL implicDiv2Dflow
+       _RL implicitNHPress
        _RL hFacMin
        _RL hFacMinDz
        _RL hFacMinDp
-Line 748 
 C     psiEuler      :: Euler angle, rota
+Line 774 
 C     psiEuler      :: Euler angle, rota
        _RL hFacInf
        _RL hFacSup
        _RL beta
-       _RL viscAr
+       _RL viscArNr(Nr)
        _RL viscAh
        _RL viscAhW
        _RL viscAhD
-Line 784 
 C     psiEuler      :: Euler angle, rota
+Line 810 
 C     psiEuler      :: Euler angle, rota
        _RL diffKrBLEQdeep
        _RL diffKrBLEQscl
        _RL diffKrBLEQHo
-       _RL delt
+       _RL tauCD, rCD, epsAB_CD
-       _RL tauCD
-       _RL rCD
        _RL gravity
        _RL recip_gravity
        _RL gBaro
-Line 795 
 C     psiEuler      :: Euler angle, rota
+Line 819 
 C     psiEuler      :: Euler angle, rota
        _RL rhoFacC(Nr),   recip_rhoFacC(Nr)
        _RL rhoFacF(Nr+1), recip_rhoFacF(Nr+1)
        _RL rhoConstFresh
+       _RL rho1Ref(Nr)
        _RL tRef(Nr)
        _RL sRef(Nr)
        _RL phiRef(2*Nr+1)
-Line 825 
 C     psiEuler      :: Euler angle, rota
+Line 850 
 C     psiEuler      :: Euler angle, rota
        _RL rotationPeriod
        _RL tauThetaClimRelax
        _RL tauSaltClimRelax
-       _RL tauTr1ClimRelax
-       _RL lambdaTr1ClimRelax
        _RL latBandClimRelax
        _RL externForcingCycle
        _RL externForcingPeriod
        _RL convertFW2Salt
        _RL temp_EvPrRn
        _RL salt_EvPrRn
+       _RL temp_addMass
+       _RL salt_addMass
        _RL ivdc_kappa
        _RL hMixCriteria
+       _RL dRhoSmall
+       _RL hMixSmooth
        _RL Ro_SeaLevel
+       _RL rSigmaBnd
        _RL sideDragFactor
        _RL bottomDragLinear
        _RL bottomDragQuadratic
-Line 871 
 C             derived from the orography
+Line 899 
 C             derived from the orography
        INTEGER integr_GeoPot, selectFindRoSurf
  C Logical flags for selecting packages
+       LOGICAL useGAD
+       LOGICAL useOBCS
+       LOGICAL useSHAP_FILT
+       LOGICAL useZONAL_FILT
        LOGICAL useOPPS
        LOGICAL usePP81
        LOGICAL useMY82
        LOGICAL useGGL90
        LOGICAL useKPP
-       LOGICAL useGAD
        LOGICAL useGMRedi
-       LOGICAL useOBCS
+       LOGICAL useDOWN_SLOPE
-       LOGICAL useAIM
-       LOGICAL useLand
        LOGICAL useCAL
        LOGICAL useEXF
+       LOGICAL useBulkForce
        LOGICAL useEBM
+       LOGICAL useCheapAML
        LOGICAL useGrdchk
+       LOGICAL useSMOOTH
        LOGICAL useECCO
-       LOGICAL useSHAP_FILT
+       LOGICAL useSBO
-       LOGICAL useZONAL_FILT
        LOGICAL useFLT
        LOGICAL usePTRACERS
        LOGICAL useGCHEM
        LOGICAL useRBCS
        LOGICAL useOffLine
        LOGICAL useMATRIX
-       LOGICAL useSBO
        LOGICAL useSEAICE
+       LOGICAL useSALT_PLUME
        LOGICAL useShelfIce
+       LOGICAL useICEFRONT
        LOGICAL useThSIce
        LOGICAL useATM2d
-       LOGICAL useBulkForce
+       LOGICAL useAIM
-       LOGICAL usefizhi
+       LOGICAL useLand
-       LOGICAL usegridalt
+       LOGICAL useFizhi
+       LOGICAL useGridAlt
        LOGICAL useDiagnostics
-       LOGICAL useMNC
        LOGICAL useREGRID
+       LOGICAL useLayers
+       LOGICAL useMNC
        LOGICAL useRunClock
        LOGICAL useEMBED_FILES
+       LOGICAL useOASIS
        LOGICAL useMYPACKAGE
-       LOGICAL useSALT_PLUME
        COMMON /PARM_PACKAGES/
+      &        useGAD, useOBCS, useSHAP_FILT, useZONAL_FILT,
       &        useOPPS, usePP81, useMY82, useGGL90, useKPP,
-      &        useGAD, useGMRedi, useOBCS, useAIM, useLand,
+      &        useGMRedi, useDOWN_SLOPE,
-      &        useCAL, useEXF, useEBM, useGrdchk, useECCO,
+      &        useCAL, useEXF, useBulkForce, useEBM, useCheapAML,
-      &        useSHAP_FILT, useZONAL_FILT, useFLT,
+      &        useGrdchk, useSMOOTH, useECCO, useSBO, useFLT,
       &        usePTRACERS, useGCHEM, useRBCS, useOffLine, useMATRIX,
-      &        useSBO, useSEAICE, useShelfIce,
+      &        useSEAICE, useSALT_PLUME, useShelfIce,
-      &        useThSIce, useATM2D, useBulkForce,
+      &        useICEFRONT, useThSIce,
-      &        usefizhi, usegridalt, useDiagnostics, useMNC, useREGRID,
+      &        useATM2D, useAIM, useLand, useFizhi, useGridAlt,
-      &        useRunClock, useEMBED_FILES, useMYPACKAGE, useSALT_PLUME
+      &        useDiagnostics, useREGRID, useLayers, useMNC,
+      &        useRunClock, useEMBED_FILES,
+      &        useOASIS,
+      &        useMYPACKAGE
  CEH3 ;;; Local Variables: ***
  CEH3 ;;; mode:fortran ***
  CEH3 ;;; End: ***

 Legend:



Removed from v.1.214
 


changed lines


 
Added in v.1.240
 Legend:



Removed from v.1.214
 


changed lines


 
Added in v.1.240
-Removed from v.1.214
+Added in v.1.240

	ViewVC Help
Powered by ViewVC 1.1.22