/[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.2 by cnh,
Fri Apr 24 02:11:37 1998 UTC
+revision 1.41 by adcroft,
Thu Aug 26 17:47:37 1999 UTC
 Line 4 
 C     /=================================
  C     | PARAMS.h                                                 |
  C     | o Header file defining model "parameters".               |
  C     |==========================================================|
- C     | The values from the model's standard input file are      |
+ C     | The values from the model standard input file are        |
  C     | stored into the variables held here. Notes describing    |
  C     | the parameters can also be found here.                   |
  C     \==========================================================/
+ C     Macros for special grid options
+ #include "PARAMS_MACROS.h"
  C--   Contants
- C     nOBands - No. of offline data time bands
-       INTEGER nOBands
-       PARAMETER ( nOBands = 12 )
  C     Useful physical values
        Real*8 PI
        PARAMETER ( PI    = 3.14159265358979323844D0   )
        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     Checkpoint data
+       INTEGER maxNoChkptLev
+       PARAMETER ( maxNoChkptLev = 2 )
  C--   COMMON /PARM_C/ Character valued parameters used by the model.
- C     oBandId  - Offline dataset identifiers for different periods.
+ C     checkPtSuff - List of checkpoint file suffices
-       COMMON /PARM_C/ oBandId
+ C     bathyFile   - File containing bathymetry. If not defined bathymetry
-       CHARACTER*3 oBandId(nOBands)
+ C                   is taken from inline function.
+ C     hydrogThetaFile - File containing initial hydrographic data for potential
+ C                       temperature.
+ C     hydrogSaltFile  - File containing initial hydrographic data for salinity.
+ C     zonalWindFile   - File containing zonal wind data
+ C     meridWindFile   - File containing meridional wind data
+ C     thetaClimFile   - File containing theta climataology used
+ C                       in relaxation term -lambda(theta-theta*)
+ C     saltClimFile    - File containing salt climataology used
+ C                       in relaxation term -lambda(salt-salt*)
+ C     surfQfile       - File containing surface heat flux
+ C     EmPmRfile       - File containing surface fresh water flux
+ C     buoyancyRelation - Flag used to indicate which relation to use to
+ C                        get buoyancy.
+       COMMON /PARM_C/ checkPtSuff,
+      &                bathyFile, hydrogThetaFile, hydrogSaltFile,
+      &                zonalWindFile, meridWindFile, thetaClimFile,
+      &                saltClimFile, buoyancyRelation,
+      &                EmPmRfile,surfQfile
+       CHARACTER*(5) checkPtSuff(maxNoChkptLev)
+       CHARACTER*(MAX_LEN_FNAM) bathyFile
+       CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile
+       CHARACTER*(MAX_LEN_FNAM) hydrogSaltFile
+       CHARACTER*(MAX_LEN_FNAM) zonalWindFile
+       CHARACTER*(MAX_LEN_FNAM) meridWindFile
+       CHARACTER*(MAX_LEN_FNAM) thetaClimFile
+       CHARACTER*(MAX_LEN_FNAM) saltClimFile
+       CHARACTER*(MAX_LEN_FNAM) surfQfile
+       CHARACTER*(MAX_LEN_FNAM) EmPmRfile
+       CHARACTER*(MAX_LEN_FNAM) buoyancyRelation
  C--   COMMON /PARM_I/ Integer valued parameters used by the model.
  C     cg2dMaxIters        - Maximum number of iterations in the
  C                           two-dimensional con. grad solver.
  C     cg2dChkResFreq      - Frequency with which to check residual
  C                           in con. grad solver.
+ C     cg3dMaxIters        - Maximum number of iterations in the
+ C                           three-dimensional con. grad solver.
+ C     cg3dChkResFreq      - Frequency with which to check residual
+ C                           in con. grad solver.
  C     nIter0              - Start time-step number of for this run
  C     nTimeSteps          - Number of timesteps to execute
  C     numStepsPerPickup   - For offline setup. Frequency of pickup
  C                           of flow fields.
+ C     writeStatePrec      - Precision used for writing model state.
+ C     writeBinaryPrec     - Precision used for writing binary files
+ C     readBinaryPrec      - Precision used for reading binary files
+ C     nCheckLev           - Holds current checkpoint level
+ C     nShap               - "Order" of shapiro filter to apply to
+ C                           model prognositic fields.
+ C                           nShap == 1 => del2
+ C                           nShap == 2 => del4
+ C                           nShap == 3 => del6
+ C                           etc...
        COMMON /PARM_I/
       &        cg2dMaxIters,
       &        cg2dChkResFreq,
-      &        nIter0, nTimeSteps,
+      &        cg3dMaxIters,
-      &        numStepsPerPickup
+      &        cg3dChkResFreq,
+      &        nIter0, nTimeSteps, nEndIter,
+      &        numStepsPerPickup,
+      &        writeStatePrec, nCheckLev,
+      &        writeBinaryPrec, readBinaryPrec,
+      &        nShap
        INTEGER cg2dMaxIters
        INTEGER cg2dChkResFreq
+       INTEGER cg3dMaxIters
+       INTEGER cg3dChkResFreq
        INTEGER nIter0
        INTEGER nTimeSteps
+       INTEGER nEndIter
        INTEGER numStepsPerPickup
+       INTEGER writeStatePrec
+       INTEGER writeBinaryPrec
+       INTEGER readBinaryPrec
+       INTEGER nCheckLev
+       INTEGER nShap
  C--   COMMON /PARM_L/ Logical valued parameters used by the model.
  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     no_slip_sides - Impose "no-slip" at lateral boundaries.
+ C     no_slip_bottom- Impose "no-slip" at bottom boundary.
  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     momPressureForcing - Flag which turns pressure term in momentum equation
+ C                          on and off.
+ C     metricTerms   - Flag which turns metric terms on or off.
+ C     usingSphericalPolarMTerms - If TRUE use spherical polar metric terms.
  C     useCoriolis   - Flag which turns the coriolis terms on and off.
  C     tempDiffusion - Flag which turns diffusion of temperature on
  C                     and off.
-Line 65 
 C     saltAdvection - Flag which turns a
+Line 152 
 C     saltAdvection - Flag which turns a
  C                     and off.
  C     saltForcing   - Flag which turns external forcing of salinit on
  C                     and off.
+ C     implicitFreeSurface - Set to true to use implcit free surface
+ C     rigidLid            - Set to true to use rigid lid
+ 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     implicitDiffusion - Turns implicit vertical diffusion on
+ C     implicitViscosity - Turns implicit vertical viscosity on
+ 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 pressure
+ C                        coords.
+ C     usingZCoords     - Set to indicate that we are working in height
+ C                        coords.
+ C     openBoundaries - Using open-boundaries
+ C     nonHydrostatic - Using non-hydrostatic terms
+ C     globalFiles    - Selects between "global" and "tiled" files
+ C     allowFreezing  - Allows water to freeze and form ice
        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
-      & momViscosity, momAdvection, momForcing, useCoriolis,
+      & no_slip_sides,no_slip_bottom,
-      & tempDiffusion, tempAdvection, tempForcing,
+      & momViscosity, momAdvection, momForcing, useCoriolis,
-      & saltDiffusion, saltAdvection, saltForcing
+      & momPressureForcing,tempDiffusion, tempAdvection, tempForcing,
+      & saltDiffusion, saltAdvection, saltForcing,
+      & implicitFreeSurface, rigidLid,
+      & momStepping, tempStepping, saltStepping,
+      & metricTerms, usingSphericalPolarMTerms,
+      & useConstantF, useBetaPlaneF, useSphereF,
+      & implicitDiffusion, implicitViscosity,
+      & doThetaClimRelax, doSaltClimRelax,
+      & periodicExternalForcing, usingPCoords, usingZCoords,
+      & openBoundaries, nonHydrostatic, globalFiles,
+      & allowFreezing
        LOGICAL usingCartesianGrid
        LOGICAL usingSphericalPolarGrid
+       LOGICAL usingSphericalPolarMTerms
+       LOGICAL no_slip_sides
+       LOGICAL no_slip_bottom
        LOGICAL momViscosity
        LOGICAL momAdvection
        LOGICAL momForcing
+       LOGICAL momPressureForcing
        LOGICAL useCoriolis
        LOGICAL tempDiffusion
        LOGICAL tempAdvection
-Line 81 
 C                     and off.
+Line 205 
 C                     and off.
        LOGICAL saltDiffusion
        LOGICAL saltAdvection
        LOGICAL saltForcing
+       LOGICAL implicitFreeSurface
+       LOGICAL rigidLid
+       LOGICAL momStepping
+       LOGICAL tempStepping
+       LOGICAL saltStepping
+       LOGICAL metricTerms
+       LOGICAL useConstantF
+       LOGICAL useBetaPlaneF
+       LOGICAL useSphereF
+       LOGICAL implicitDiffusion
+       LOGICAL implicitViscosity
+       LOGICAL doThetaClimRelax
+       LOGICAL doSaltClimRelax
+       LOGICAL periodicExternalForcing
+       LOGICAL usingPCoords
+       LOGICAL usingZCoords
+       LOGICAL openBoundaries
+       LOGICAL nonHydrostatic
+       LOGICAL globalFiles
+       LOGICAL allowFreezing
  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
  C     cg2dTargetResidual
  C               - Target residual for cg2d solver.
- C     delZ      - Vertical grid spacing (Pa) - delZ is the distance
+ C     cg3dTargetResidual
- C                 between "w" surfaces.
+ C               - Target residual for cg3d solver.
+ C     cg2dpcOffDFac - Averaging weight for preconditioner off-diagonal.
+ C     Note. 20th May 1998
+ C           I made a weird discovery! In the model paper we argue
+ C           for the form of the preconditioner used here ( see
+ C           A Finite-volume, Incompressible Navier-Stokes Model
+ C           ...., Marshall et. al ). The algebra gives a simple
+ C           0.5 factor for the averaging of ac and aCw to get a
+ C           symmettric pre-conditioner. By using a factor of 0.51
+ C           i.e. scaling the off-diagonal terms in the
+ C           preconditioner down slightly I managed to get the
+ C           number of iterations for convergence in a test case to
+ C           drop form 192 -> 134! Need to investigate this further!
+ C           For now I have introduced a parameter cg2dpcOffDFac which
+ C           defaults to 0.51 but can be set at runtime.
+ C     delP      - Vertical grid spacing ( Pa ).
+ C     delZ      - Vertical grid spacing ( m  ).
+ C     delR      - Vertical grid spacing ( units of r ).
  C     delX      - Separation between cell faces (m) or (deg), depending
  C     delY        on input flags.
  C     gravity   - Accel. due to gravity ( m/s^2 )
+ C     recip_gravity and its inverse
+ C     gBaro     - Accel. due to gravity used in barotropic equation ( m/s^2 )
  C     ronil     - Reference density
- C     startTime - Start time for model
+ C     rhoConst  - Vertically constant reference density
- C     oBandBeg  - Offline bands start and end times.
+ C     startTime - Start time for model ( s )
- C     oBandEnd
- C     oPeriodLength - Lenght of offline period.
  C     phiMin    - Latitude of southern most cell face.
  C     thetaMin  - Longitude of western most cell face (this
  C                 is an "inert" parameter but it is included
  C                 to make geographical references simple.)
- C     rSphere   - Radius of sphere for a spherical polar 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                 ( Southern edge f for beta plane )
  C     beta      - df/dy ( s^-1.m^-1 )
+ C     omega     - Angular velocity ( rad/s )
  C     viscAh    - Eddy viscosity coeff. for mixing of
  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 )
+ 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     diffKhT   - Laplacian diffusion coeff. for mixing of
  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 )
+ C     diffKrT   - Laplacian diffusion coeff. for mixing of
+ C                 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     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 )
+ C     diffKrS   - Laplacian diffusion coeff. for mixing of
+ C                 salt vertically ( units of r^2/s )
  C     diffK4S   - Biharmonic diffusion coeff. for mixing of
  C                 salt laterally ( m^4/s )
- C     delt      - Timestep ( s )
+ C     deltaT    - Default timestep ( s )
- C     tauCD     - CD scheme coupling timescale ( 1/s )
+ C     deltaTClock  - Timestep used as model "clock". This determines the
- C     rCD       - CD scheme normalised coupling parameter ( 0-1 )
+ C                    IO frequencies and is used in tagging output. It can
- C     startTime - Starting time for this integration ( s ).
+ C                    be totally different to the dynamical time. Typically
- C     endTime   - Ending time for this integration ( s ).
+ C                    it will be the deep-water timestep for accelerated runs.
- C     chkPtFreq - Frequency of check pointing ( s ).
+ C                    Frequency of checkpointing and dumping of the model state
- C     dumpFreq  - Frequency with which model state is written to
+ C                    are referenced to this clock. ( s )
- C                 post-processing files ( s ).
+ C     deltaTMom    - Timestep for momemtum equations ( s )
-       COMMON /PARM_R/ cg2dTargetResidual, delZ, delX, delY,
+ C     deltaTtracer - Timestep for tracer equations ( s )
-      & deltaTmom, deltaTtracer, abeps, startTime, oBandBeg, oBandEnd,
+ C     freesurfFac  - Parameter to turn implicit free surface term on or off
-      & oPeriodLength, phiMin, thetaMin, rSphere, f0, fCori, beta,
+ C                    freesurfac = 1. uses implicit free surface
-      & viscAh, viscAz, viscA4, diffKhT, diffKzT, diffK4T,
+ C                    freesurfac = 0. uses rigid lid
-      & diffKhS, diffKzS, diffK4S, delT, tauCD, rCD,
+ C     hFacMin      - Minimum fraction size of a cell (affects hFacC etc...)
-      & tAlpha, sBeta, gravity, rhonil, tRef, sRef,
+ C     hFacMinDz    - Minimum dimesional size of a cell (affects hFacC etc..., m)
-      & endTime, chkPtFreq, dumpFreq
+ C     hFacMinDp    - Minimum dimesional size of a cell (affects hFacC etc..., Pa)
-       REAL cg2dTargetResidual
+ C     hFacMinDr    - Minimum dimesional size of a cell (affects hFacC etc..., units of r)
-       REAL delZ(Nz)
+ C     tauCD        - CD scheme coupling timescale ( 1/s )
-       REAL delX(Nx)
+ C     rCD          - CD scheme normalised coupling parameter ( 0-1 )
-       REAL delY(Ny)
+ C     GMmaxslope    - max. slope allowed in GM/Redi tensor
-       REAL deltaTmom
+ C     GMlength      - Length to use in Visbeck et al. formula for K (m)
-       REAL deltaTtracer
+ C     GMalpha       - alpha to use in Visbeck et al. formula for K
-       REAL abeps
+ C     GMdepth       - Depth over which to integrate Richardson # (Visbeck et al.)
-       REAL oBandBeg(nOBands),oBandEnd(nOBands)
+ C     GMkbackground - background value of GM/Redi diffusion coefficient ( m^2/s )
-       REAL oPeriodLength
+ C     GMmaxval      - max. value of KapGM allowed in GM/Redi scheme ( m^2/s )
-       REAL phiMin
+ C     startTime     - Starting time for this integration ( s ).
-       REAL thetaMin
+ C     endTime       - Ending time for this integration ( s ).
-       REAL rSphere
+ C     chkPtFreq     - Frequency of rolling check pointing ( s ).
-       REAL f0
+ C     pChkPtFreq    - Frequency of permanent check pointing ( s ).
-       _RL  beta
+ C     dumpFreq      - Frequency with which model state is written to
-       REAL viscAh
+ C                     post-processing files ( s ).
-       REAL viscAz
+ C     afFacMom      - Advection of momentum term tracer parameter
-       REAL viscA4
+ C     vfFacMom      - Momentum viscosity tracer parameter
-       REAL diffKhT
+ C     pfFacMom      - Momentum pressure forcing tracer parameter
-       REAL diffKzT
+ C     cfFacMom      - Coriolis term tracer parameter
-       REAL diffK4T
+ C     foFacMom      - Momentum forcing tracer parameter
-       REAL diffKhS
+ C     mtFacMom      - Metric terms tracer parameter
-       REAL diffKzS
+ C     cosPower      - Power of cosine of latitude to multiply viscosity
-       REAL diffK4S
+ C     cAdjFreq      - Frequency of convective adjustment
-       REAL delt
+ C
-       REAL tauCD
+ C     taveFreq      - Frequency with which time-averaged model state is written to
-       REAL rCD
+ C                     post-processing files ( s ).
-       REAL tAlpha
+ C     tauThetaClimRelax - Relaxation to climatology time scale ( s ).
-       REAL sBeta
+ C     lambdaThetaClimRelax - Inverse time scale for relaxation ( 1/s ).
-       REAL gravity
+ C     tauSaltClimRelax - Relaxation to climatology time scale ( s ).
-       REAL rhonil
+ C     lambdaSaltClimRelax - Inverse time scale for relaxation ( 1/s ).
-       REAL tRef(Nz)
+ C     externForcingPeriod - Is the period of which forcing varies (eg. 1 month)
-       REAL sRef(Nz)
+ C     externForcingCycle - Is the repeat time of the forcing (eg. 1 year)
-       real Fcori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
+ C                          (note: externForcingCycle must be an integer
-       REAL startTime
+ C                           number times externForcingPeriod)
-       REAL endTime
+ C     horiVertRatio      - Ratio on units in vertical to units in horizontal.
-       REAL chkPtFreq
+ C     recip_horiVertRatio  ( 1 if horiz in m and vertical in m ).
-       REAL dumpFreq
+ C                          ( g*rho if horiz in m and vertical in Pa ).
-       COMMON /PARM_A/ HeatCapacity_Cp,
+ C     latFFTFiltLo       - Low latitude for FFT filtering of latitude
+ C                          circles ( see filter*.F )
+       COMMON /PARM_R/ cg2dTargetResidual, cg2dpcOffDFac,
+      & cg3dTargetResidual,
+      & delP, delZ, delR, delX, delY,
+      & deltaT,deltaTmom, deltaTtracer, deltaTClock,abeps, startTime,
+      & phiMin, thetaMin, rSphere, recip_RSphere, f0, fCori, beta,
+      & viscAh,  viscAz,  viscA4,  viscAr,
+      & diffKhT, diffKzT, diffK4T, diffKrT,
+      & diffKhS, diffKzS, diffK4S, diffKrS,
+      & delT, tauCD, rCD, freeSurfFac, hFacMin, hFacMinDz,
+      & GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,GMmaxval,
+      & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
+      & recip_rhoConst, rhoConst, tRef, sRef,
+      & endTime, chkPtFreq, pchkPtFreq, dumpFreq, taveFreq,
+      & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
+      & cosPower,
+      & cAdjFreq, omega, tauThetaClimRelax, lambdaThetaClimRelax,
+      & tauSaltClimRelax, lambdaSaltClimRelax,
+      & externForcingCycle, externForcingPeriod,
+      & viscAp, diffKpT, diffKpS, hFacMinDr, hFacMinDp,
+      & theta_S, specVol_S, horiVertRatio, recip_horiVertRatio,
+      & latFFTFiltLo, ivdc_kappa
+       _RL cg2dTargetResidual
+       _RL cg3dTargetResidual
+       _RL cg2dpcOffDFac
+       _RL delZ(Nr)
+       _RL delP(Nr)
+       _RL delR(Nr)
+       _RL delX(Nx)
+       _RL delY(Ny)
+       _RL deltaT
+       _RL deltaTClock
+       _RL deltaTmom
+       _RL deltaTtracer
+       _RL abeps
+       _RL phiMin
+       _RL thetaMin
+       _RL rSphere
+       _RL recip_RSphere
+       _RL f0
+       _RL freeSurfFac
+       _RL hFacMin
+       _RL hFacMinDz
+       _RL hFacMinDp
+       _RL hFacMinDr
+       _RL beta
+       _RL viscAh
+       _RL viscAz
+       _RL viscAp
+       _RL viscAr
+       _RL viscA4
+       _RL diffKhT
+       _RL diffKrT
+       _RL diffKzT
+       _RL diffKpT
+       _RL diffK4T
+       _RL diffKhS
+       _RL diffKrS
+       _RL diffKzS
+       _RL diffKpS
+       _RL diffK4S
+       _RL delt
+       _RL tauCD
+       _RL rCD
+       _RL GMmaxslope
+       _RL GMlength
+       _RL GMalpha
+       _RL GMdepth
+       _RL GMkbackground
+       _RL GMmaxval
+       _RL gravity
+       _RL recip_gravity
+       _RL gBaro
+       _RL rhonil
+       _RL recip_rhonil
+       _RL rhoConst
+       _RL recip_rhoConst
+       _RL specVol_S(Nr)
+       _RL tRef(Nr)
+       _RL theta_S(Nr)
+       _RL sRef(Nr)
+       _RS Fcori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)
+       _RL startTime
+       _RL endTime
+       _RL chkPtFreq
+       _RL pChkPtFreq
+       _RL dumpFreq
+       _RL taveFreq
+       _RL afFacMom
+       _RL vfFacMom
+       _RL pfFacMom
+       _RL cfFacMom
+       _RL foFacMom
+       _RL mTFacMom
+       _RL cosPower
+       _RL cAdjFreq
+       _RL omega
+       _RL tauThetaClimRelax
+       _RL lambdaThetaClimRelax
+       _RL tauSaltClimRelax
+       _RL lambdaSaltClimRelax
+       _RL externForcingCycle
+       _RL externForcingPeriod
+       _RL horiVertRatio
+       _RL recip_horiVertRatio
+       _RL latFFTFiltLo
+       _RL ivdc_kappa
+       COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp,
       &                Lamba_theta
-       REAL HeatCapacity_Cp
+       _RL HeatCapacity_Cp
-       REAL Lamba_theta
+       _RL Lamba_theta
+       _RL recip_Cp
+ C Equation of State (polynomial coeffients)
+       COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS
+       _RL eosC(9,Nr+1),eosSig0(Nr+1),eosRefT(Nr+1),eosRefS(Nr+1)
+ C Linear equation of state
+ C     tAlpha    - Linear EOS thermal expansion coefficient ( 1/degree ).
+ C     sBeta     - Linear EOS haline contraction coefficient.
+       COMMON /PARM_EOS_LIN/ tAlpha,sBeta,eosType
+       _RL tAlpha
+       _RL sBeta
+       character*(6) eosType
+ C These are input arrays (of integers) that contain the *absolute*
+ C computational index of an open-boundary (OB) point.
+ C A zero (0) element means there is no corresponding OB in that column/row.
+ C The computational coordinate refers to "tracer" cells.
+ C For a northern/southern OB, the OB V point is to the south/north.
+ C For an eastern/western OB, the OB U point is to the west/east.
+ C eg.
+ C     OB_Jnorth(3)=34  means that:
+ C          T( 3 ,34) is a an OB point
+ C          U(3:4,34) is a an OB point
+ C          V( 4 ,34) is a an OB point
+ C while
+ C     OB_Jsouth(3)=1  means that:
+ C          T( 3 ,1) is a an OB point
+ C          U(3:4,1) is a an OB point
+ C          V( 4 ,2) is a an OB point
+ C
+ C For convenience, negative values for Jnorth/Ieast refer to
+ C points relative to the Northern/Eastern edges of the model
+ C eg. OB_Jnorth(3)=-1  means that the point (3,Ny-1) is a northern O-B.
+ C
+       COMMON /PARM_OB/
+      & OB_Jnorth,OB_Jsouth,OB_Ieast,OB_Iwest
+       INTEGER OB_Jnorth(Nx)
+       INTEGER OB_Jsouth(Nx)
+       INTEGER OB_Ieast(Ny)
+       INTEGER OB_Iwest(Ny)

 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.41
 Legend:



Removed from v.1.2
 


changed lines


 
Added in v.1.41
-Removed from v.1.2
+Added in v.1.41

	ViewVC Help
Powered by ViewVC 1.1.22