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

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

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

revision 1.6 by adcroft, Wed May 20 21:29:31 1998 UTC revision 1.162 by jmc, Sun Jul 31 22:01:15 2005 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3  C  C
4  C     /==========================================================\  
5  C     | PARAMS.h                                                 |  CBOP
6  C     | o Header file defining model "parameters".               |  C     !ROUTINE: PARAMS.h
7  C     |==========================================================|  C     !INTERFACE:
8  C     | The values from the model's standard input file are      |  C     #include PARAMS.h
9  C     | stored into the variables held here. Notes describing    |  
10  C     | the parameters can also be found here.                   |  C     !DESCRIPTION:
11  C     \==========================================================/  C     Header file defining model "parameters".  The values from the
12    C     model standard input file are stored into the variables held
13    C     here. Notes describing the parameters can also be found here.
14    
15    CEOP
16    
17    C     Macros for special grid options
18    #include "PARAMS_MACROS.h"
19    
20  C--   Contants  C--   Contants
 C     nOBands - No. of offline data time bands  
       INTEGER nOBands  
       PARAMETER ( nOBands = 12 )  
21  C     Useful physical values  C     Useful physical values
22        Real*8 PI        Real*8 PI
23        PARAMETER ( PI    = 3.14159265358979323844D0   )        PARAMETER ( PI    = 3.14159265358979323844D0   )
24        Real*8 deg2rad        Real*8 deg2rad
25        PARAMETER ( deg2rad = 2.D0*PI/360.D0           )        PARAMETER ( deg2rad = 2.D0*PI/360.D0           )
26    
27    C     Symbolic values
28    C     precXXXX :: Used to indicate what precision to use for
29    C                dumping model state.
30          INTEGER precFloat32
31          PARAMETER ( precFloat32 = 32 )
32          INTEGER precFloat64
33          PARAMETER ( precFloat64 = 64 )
34    C     UNSET_xxx :: Used to indicate variables that have not been given a value
35          Real*8 UNSET_FLOAT8
36          PARAMETER ( UNSET_FLOAT8 = 1.234567D5 )
37          Real*4 UNSET_FLOAT4
38          PARAMETER ( UNSET_FLOAT4 = 1.234567E5 )
39          _RL    UNSET_RL    
40          PARAMETER ( UNSET_RL     = 1.234567D5 )
41          _RS    UNSET_RS    
42          PARAMETER ( UNSET_RS     = 1.234567E5 )
43          INTEGER UNSET_I
44          PARAMETER ( UNSET_I      = 123456789  )
45    
46    C     Checkpoint data
47          INTEGER maxNoChkptLev
48          PARAMETER ( maxNoChkptLev = 2 )
49    
50  C--   COMMON /PARM_C/ Character valued parameters used by the model.  C--   COMMON /PARM_C/ Character valued parameters used by the model.
51  C     oBandId  - Offline dataset identifiers for different periods.  C     checkPtSuff :: List of checkpoint file suffices
52        COMMON /PARM_C/ oBandId  C     delXFile      :: File containing X-spacing grid definition (1.D array)
53        CHARACTER*3 oBandId(nOBands)  C     delYFile      :: File containing Y-spacing grid definition (1.D array)
54    C     horizGridFile :: File containing horizontal-grid definition
55    C                        (only when using curvilinear_grid)
56    C     bathyFile   :: File containing bathymetry. If not defined bathymetry
57    C                   is taken from inline function.
58    C     topoFile    :: File containing the topography of the surface (unit=m)
59    C                   (mainly used for the atmosphere = ground height).
60    C     hydrogThetaFile :: File containing initial hydrographic data for potential
61    C                       temperature.
62    C     hydrogSaltFile  :: File containing initial hydrographic data for salinity.
63    C     zonalWindFile   :: File containing zonal wind data
64    C     meridWindFile   :: File containing meridional wind data
65    C     thetaClimFile   :: File containing theta climataology used
66    C                       in relaxation term -lambda(theta-theta*)
67    C     saltClimFile    :: File containing salt climataology used
68    C                       in relaxation term -lambda(salt-salt*)
69    C     surfQfile       :: File containing surface heat flux, excluding SW
70    C                        (old version, kept for backward compatibility)
71    C     surfQnetFile    :: File containing surface net heat flux
72    C     surfQswFile     :: File containing surface shortwave radiation
73    C     dQdTfile        :: File containing thermal relaxation coefficient
74    C     EmPmRfile       :: File containing surface fresh water flux
75    C     saltFluxFile    :: File containing surface salt flux
76    C     pLoadFile       :: File containing pressure loading
77    C     eddyTauxFile    :: File containing zonal Eddy stress data
78    C     eddyTauyFile    :: File containing meridional Eddy stress data
79    C     buoyancyRelation :: Flag used to indicate which relation to use to
80    C                        get buoyancy.
81    C     eosType         :: choose the equation of state:
82    C                        LINEAR, POLY3, UNESCO, JMD95Z, JMD95P, MDJWF, IDEALGAS
83    C     the_run_name    :: string identifying the name of the model "run"
84          COMMON /PARM_C/ checkPtSuff,
85         &                delXFile, delYFile, horizGridFile,
86         &                bathyFile, topoFile,
87         &                hydrogThetaFile, hydrogSaltFile,
88         &                zonalWindFile, meridWindFile, thetaClimFile,
89         &                saltClimFile, buoyancyRelation,
90         &                EmPmRfile, saltFluxFile,
91         &                surfQfile, surfQnetFile, surfQswFile,
92         &                lambdaThetaFile, lambdaSaltFile,
93         &                uVelInitFile, vVelInitFile, pSurfInitFile,
94         &                dQdTfile, ploadFile,
95         &                eddyTauxFile, eddyTauyFile,
96         &                eosType, pickupSuff,
97         &                mdsioLocalDir,
98         &                the_run_name
99          CHARACTER*(5) checkPtSuff(maxNoChkptLev)
100          CHARACTER*(MAX_LEN_FNAM) delXFile
101          CHARACTER*(MAX_LEN_FNAM) delYFile
102          CHARACTER*(MAX_LEN_FNAM) horizGridFile
103          CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile
104          CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile
105          CHARACTER*(MAX_LEN_FNAM) hydrogSaltFile
106          CHARACTER*(MAX_LEN_FNAM) zonalWindFile
107          CHARACTER*(MAX_LEN_FNAM) meridWindFile
108          CHARACTER*(MAX_LEN_FNAM) thetaClimFile
109          CHARACTER*(MAX_LEN_FNAM) saltClimFile
110          CHARACTER*(MAX_LEN_FNAM) surfQfile
111          CHARACTER*(MAX_LEN_FNAM) surfQnetFile
112          CHARACTER*(MAX_LEN_FNAM) surfQswFile
113          CHARACTER*(MAX_LEN_FNAM) EmPmRfile
114          CHARACTER*(MAX_LEN_FNAM) saltFluxFile
115          CHARACTER*(MAX_LEN_FNAM) buoyancyRelation
116          CHARACTER*(MAX_LEN_FNAM) uVelInitFile
117          CHARACTER*(MAX_LEN_FNAM) vVelInitFile
118          CHARACTER*(MAX_LEN_FNAM) pSurfInitFile
119          CHARACTER*(MAX_LEN_FNAM) dQdTfile
120          CHARACTER*(MAX_LEN_FNAM) ploadFile
121          CHARACTER*(MAX_LEN_FNAM) eddyTauxFile
122          CHARACTER*(MAX_LEN_FNAM) eddyTauyFile
123          CHARACTER*(MAX_LEN_FNAM) lambdaThetaFile
124          CHARACTER*(MAX_LEN_FNAM) lambdaSaltFile
125          CHARACTER*(MAX_LEN_FNAM) mdsioLocalDir
126          CHARACTER*(MAX_LEN_FNAM) the_run_name
127          CHARACTER*(6) eosType
128          CHARACTER*(10) pickupSuff
129    
130  C--   COMMON /PARM_I/ Integer valued parameters used by the model.  C--   COMMON /PARM_I/ Integer valued parameters used by the model.
131  C     cg2dMaxIters        - Maximum number of iterations in the  C     cg2dMaxIters        :: Maximum number of iterations in the
132  C                           two-dimensional con. grad solver.  C                           two-dimensional con. grad solver.
133  C     cg2dChkResFreq      - Frequency with which to check residual  C     cg2dChkResFreq      :: Frequency with which to check residual
134    C                           in con. grad solver.
135    C     cg2dPreCondFreq     :: Frequency for updating cg2d preconditioner
136    C                            (non-linear free-surf.)
137    C     cg3dMaxIters        :: Maximum number of iterations in the
138    C                           three-dimensional con. grad solver.
139    C     cg3dChkResFreq      :: Frequency with which to check residual
140  C                           in con. grad solver.  C                           in con. grad solver.
141  C     nIter0              - Start time-step number of for this run  C     nIter0              :: Start time-step number of for this run
142  C     nTimeSteps          - Number of timesteps to execute  C     nTimeSteps          :: Number of timesteps to execute
143  C     numStepsPerPickup   - For offline setup. Frequency of pickup  C     numStepsPerPickup   :: For offline setup. Frequency of pickup
144  C                           of flow fields.  C                           of flow fields.
145    C     writeStatePrec      :: Precision used for writing model state.
146    C     writeBinaryPrec     :: Precision used for writing binary files
147    C     readBinaryPrec      :: Precision used for reading binary files
148    C     nCheckLev           :: Holds current checkpoint level
149    C     nonlinFreeSurf      :: option related to non-linear free surface
150    C                           =0 Linear free surface ; >0 Non-linear
151    C     select_rStar        :: option related to r* vertical coordinate
152    C                           =0 (default) use r coord. ; > 0 use r*
153    C     tempAdvScheme       :: Temp. Horiz.Advection scheme selector
154    C     tempVertAdvScheme   :: Temp. Vert. Advection scheme selector
155    C     saltAdvScheme       :: Salt. Horiz.advection scheme selector
156    C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector
157    C     debugLevel          :: debug level selector: higher -> more writing
158    
159        COMMON /PARM_I/        COMMON /PARM_I/
160       &        cg2dMaxIters,       &        cg2dMaxIters,
161       &        cg2dChkResFreq,       &        cg2dChkResFreq, cg2dPreCondFreq,
162       &        nIter0, nTimeSteps,       &        cg3dMaxIters,
163       &        numStepsPerPickup       &        cg3dChkResFreq,
164         &        nIter0, nTimeSteps, nEndIter,
165         &        numStepsPerPickup,
166         &        writeStatePrec, nCheckLev,
167         &        writeBinaryPrec, readBinaryPrec,
168         &        nonlinFreeSurf, select_rStar,
169         &        tempAdvScheme, tempVertAdvScheme,
170         &        saltAdvScheme, saltVertAdvScheme,
171         &        debugLevel
172        INTEGER cg2dMaxIters        INTEGER cg2dMaxIters
173        INTEGER cg2dChkResFreq        INTEGER cg2dChkResFreq
174          INTEGER cg2dPreCondFreq
175          INTEGER cg3dMaxIters
176          INTEGER cg3dChkResFreq
177        INTEGER nIter0        INTEGER nIter0
178        INTEGER nTimeSteps        INTEGER nTimeSteps
179          INTEGER nEndIter
180        INTEGER numStepsPerPickup        INTEGER numStepsPerPickup
181          INTEGER writeStatePrec
182          INTEGER writeBinaryPrec
183          INTEGER readBinaryPrec
184          INTEGER nCheckLev
185          INTEGER nonlinFreeSurf
186          INTEGER select_rStar
187          INTEGER tempAdvScheme, tempVertAdvScheme
188          INTEGER saltAdvScheme, saltVertAdvScheme
189          INTEGER debugLevel
190    
191    C
192          INTEGER debLevZero
193          PARAMETER(debLevZero=0)
194          INTEGER debLevA
195          PARAMETER(debLevA=1)
196          INTEGER debLevB
197          PARAMETER(debLevB=2)
198    
199  C--   COMMON /PARM_L/ Logical valued parameters used by the model.  C--   COMMON /PARM_L/ Logical valued parameters used by the model.
200  C     usingCartesianGrid - If TRUE grid generation will be in a cartesian  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian
201  C                          coordinate frame.  C                          coordinate frame.
202  C     usingSphericalPolarGrid - If TRUE grid generation will be in a  C     usingSphericalPolarGrid :: If TRUE grid generation will be in a
203  C                               spherical polar frame.  C                               spherical polar frame.
204  C     momViscosity  - Flag which turns momentum friction terms on and off.  C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical
205  C     momAdvection  - Flag which turns advection of momentum on and off.  C     no_slip_sides :: Impose "no-slip" at lateral boundaries.
206  C     momForcing    - Flag which turns external forcing of momentum on  C     no_slip_bottom :: Impose "no-slip" at bottom boundary.
207  C                     and off.  C     staggerTimeStep :: enable a Stagger time stepping T,S Rho then U,V
208  C     useCoriolis   - Flag which turns the coriolis terms on and off.  C     momViscosity  :: Flag which turns momentum friction terms on and off.
209  C     tempDiffusion - Flag which turns diffusion of temperature on  C     momAdvection  :: Flag which turns advection of momentum on and off.
210  C                     and off.  C     momForcing    :: Flag which turns external forcing of momentum on
 C     tempAdvection - Flag which turns advection of temperature on  
211  C                     and off.  C                     and off.
212  C     tempForcing   - Flag which turns external forcing of temperature on  C     momPressureForcing :: Flag which turns pressure term in momentum equation
213    C                          on and off.
214    C     metricTerms   :: Flag which turns metric terms on or off.
215    C     usingSphericalPolarMTerms :: If TRUE use spherical polar metric terms.
216    C     useNHMTerms :: If TRUE use non-hydrostatic metric terms.
217    C     useCoriolis   :: Flag which turns the coriolis terms on and off.
218    C     tempAdvection :: Flag which turns advection of temperature on
219  C                     and off.  C                     and off.
220  C     saltDiffusion - Flag which turns diffusion of salinit on  C     tempForcing   :: Flag which turns external forcing of temperature on
221  C                     and off.  C                     and off.
222  C     saltAdvection - Flag which turns advection of salinit on  C     saltAdvection :: Flag which turns advection of salinity on
223  C                     and off.  C                     and off.
224  C     saltForcing   - Flag which turns external forcing of salinit on  C     saltForcing   :: Flag which turns external forcing of salinity on
225  C                     and off.  C                     and off.
226    C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux
227    C                         as a real Fresh Water (=> changes the Sea Level)
228    C                         if F, converts P+R-E to salt flux (no SL effect)
229    C     useFullLeith   :: Set to true to use full Leith viscosity (may be unstable
230    C                       on irregular grids)
231    C     useAnisotropicViscAGridMax :: Set to true to use Alistair's latest
232    C              anisotropic length scale.  It is used only for maximum viscosity
233    C              calculations.  Alistair recommends a value of viscA*GridMax=.25
234    C     rigidLid            :: Set to true to use rigid lid
235    C     implicitFreeSurface :: Set to true to use implcit free surface
236    C     exactConserv        :: Set to true to conserve exactly the total Volume
237    C     uniformLin_PhiSurf  :: Set to true to use a uniform Bo_surf in the
238    C                           linear relation Phi_surf = Bo_surf*eta
239    C     momStepping   :: Turns momentum equation time-stepping off
240    C     tempStepping  :: Turns temperature equation time-stepping off
241    C     saltStepping  :: Turns salinity equation time-stepping off
242    C     useConstantF  :: Coriolis parameter set to f0
243    C     useBetaPlaneF :: Coriolis parameter set to f0 + beta.y
244    C     useSphereF    :: Coriolis parameter set to 2.omega.sin(phi)
245    C     useCDscheme   :: use CD-scheme to calculate Coriolis terms.
246    C     useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986)
247    C     useJamartMomAdv :: Use wet-point method for V.I. non-linear term
248    C     SadournyCoriolis :: use the enstrophy conserving scheme by Sadourny
249    C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term
250    C     highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection)
251    C     upwindShear        :: use 1rst order upwind interp. (V.I., vertical advection)
252    C     useAbsVorticity :: work with f+zeta in Coriolis terms
253    C     implicitDiffusion :: Turns implicit vertical diffusion on
254    C     implicitViscosity :: Turns implicit vertical viscosity on
255    C     tempImplVertAdv :: Turns on implicit vertical advection for Temperature
256    C     saltImplVertAdv :: Turns on implicit vertical advection for Salinity
257    C     momImplVertAdv  :: Turns on implicit vertical advection for Momentum
258    C     multiDimAdvection :: Flag that enable multi-dimension advection
259    C     useMultiDimAdvec  :: True if multi-dim advection is used at least once
260    C     forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution
261    C                      out off Adams-Bashforth time stepping.
262    C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup
263    C     doThetaClimRelax :: Set true if relaxation to temperature
264    C                        climatology is required.
265    C     doSaltClimRelax  :: Set true if relaxation to salinity
266    C                        climatology is required.
267    C     periodicExternalForcing :: Set true if forcing is time-dependant
268    C     usingPCoords     :: Set to indicate that we are working in a pressure
269    C                        type coordinate (p or p*).
270    C     usingZCoords     :: Set to indicate that we are working in a height
271    C                        type coordinate (z or z*)
272    C     fluidIsAir       :: Set to indicate that the fluid major constituent
273    C                        is air
274    C     fluidIsWater     :: Set to indicate that the fluid major constituent
275    C                        is water
276    C     useDynP_inEos_Zc :: use the dynamical pressure in EOS (with Z-coord.)
277    C                         this requires specific code for restart & exchange
278    C     setCenterDr    :: set cell Center depth and put Interface at the middle
279    C     nonHydrostatic :: Using non-hydrostatic terms
280    C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm
281    C     globalFiles    :: Selects between "global" and "tiled" files
282    C     useSingleCpuIO :: On SGI platforms, option globalFiles is either
283    C                       slow (f77) or does not work (f90).  When
284    C                       useSingleCpuIO is set, mdsio_writefield.F
285    C                       outputs from master mpi process only.
286    C     allowFreezing  :: Allows surface water to freeze and form ice
287    C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)
288    C     pickup_write_mdsio :: use mdsio to write pickups
289    C     pickup_read_mdsio  :: use mdsio to read  pickups
290    C     pickup_write_immed :: echo the pickup immediately (for conversion)
291    C     timeave_mdsio      :: use mdsio for timeave output
292    C     snapshot_mdsio     :: use mdsio for "snapshot" (dumpfreq/diagfreq) output
293    C     monitor_stdio      :: use stdio for monitor output
294    C     calendarDumps :: When set, approximate months (30-31 days) and years (360-372 days)
295    C                      for parameters chkPtFreq, pChkPtFreq, taveFreq, SEAICE_taveFreq,
296    C                      KPP_taveFreq, and freq in pkg/diagnostics are converted to exact
297    C                      calendar months and years.  Requires pkg/cal.
298        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
299       & momViscosity, momAdvection, momForcing, useCoriolis,       & usingCurvilinearGrid, usingCylindricalGrid,
300       & tempDiffusion, tempAdvection, tempForcing,       & no_slip_sides,no_slip_bottom,
301       & saltDiffusion, saltAdvection, saltForcing       & staggerTimeStep,
302         & momViscosity, momAdvection, momForcing, useCoriolis,
303         & momPressureForcing, vectorInvariantMomentum,
304         & tempAdvection, tempForcing,
305         & saltAdvection, saltForcing,
306         & useRealFreshWaterFlux,
307         & useFullLeith, useAnisotropicViscAGridMax,
308         & rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf,
309         & momStepping, tempStepping, saltStepping,
310         & metricTerms, usingSphericalPolarMTerms, useNHMTerms,
311         & useConstantF, useBetaPlaneF, useSphereF,
312         & useCDscheme,
313         & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,
314         & SadournyCoriolis, upwindVorticity, highOrderVorticity,
315         & useAbsVorticity, upwindShear,
316         & implicitDiffusion, implicitViscosity,
317         & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,
318         & multiDimAdvection, useMultiDimAdvec, forcing_In_AB,
319         & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,
320         & periodicExternalForcing,
321         & fluidIsAir, fluidIsWater,
322         & usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,
323         & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,
324         & allowFreezing, useOldFreezing,
325         & usePickupBeforeC35, usePickupBeforeC54, startFromPickupAB2,
326         & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,
327         & timeave_mdsio, snapshot_mdsio, monitor_stdio,
328         & outputTypesInclusive, debugMode,
329         & inAdMode, inAdTrue, inAdFalse, inAdExact,
330         & calendarDumps
331    
332        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
333        LOGICAL usingSphericalPolarGrid        LOGICAL usingSphericalPolarGrid
334          LOGICAL usingCylindricalGrid
335          LOGICAL usingCurvilinearGrid
336          LOGICAL usingSphericalPolarMTerms
337          LOGICAL useNHMTerms
338          LOGICAL no_slip_sides
339          LOGICAL no_slip_bottom
340          LOGICAL staggerTimeStep
341        LOGICAL momViscosity        LOGICAL momViscosity
342        LOGICAL momAdvection        LOGICAL momAdvection
343        LOGICAL momForcing        LOGICAL momForcing
344          LOGICAL momPressureForcing
345        LOGICAL useCoriolis        LOGICAL useCoriolis
346        LOGICAL tempDiffusion        LOGICAL vectorInvariantMomentum
347        LOGICAL tempAdvection        LOGICAL tempAdvection
348        LOGICAL tempForcing        LOGICAL tempForcing
       LOGICAL saltDiffusion  
349        LOGICAL saltAdvection        LOGICAL saltAdvection
350        LOGICAL saltForcing        LOGICAL saltForcing
351          LOGICAL useRealFreshWaterFlux
352          LOGICAL useFullLeith, useAnisotropicViscAGridMax
353          LOGICAL rigidLid
354          LOGICAL implicitFreeSurface
355          LOGICAL exactConserv
356          LOGICAL uniformLin_PhiSurf
357          LOGICAL momStepping
358          LOGICAL tempStepping
359          LOGICAL saltStepping
360          LOGICAL metricTerms
361          LOGICAL useConstantF
362          LOGICAL useBetaPlaneF
363          LOGICAL useSphereF
364          LOGICAL useCDscheme
365          LOGICAL useEnergyConservingCoriolis
366          LOGICAL useJamartWetPoints
367          LOGICAL useJamartMomAdv
368          LOGICAL SadournyCoriolis
369          LOGICAL upwindVorticity
370          LOGICAL highOrderVorticity
371          LOGICAL useAbsVorticity
372          LOGICAL upwindShear
373          LOGICAL implicitDiffusion
374          LOGICAL implicitViscosity
375          LOGICAL tempImplVertAdv
376          LOGICAL saltImplVertAdv
377          LOGICAL momImplVertAdv
378          LOGICAL multiDimAdvection
379          LOGICAL useMultiDimAdvec
380          LOGICAL forcing_In_AB
381          LOGICAL doThetaClimRelax
382          LOGICAL doSaltClimRelax
383          LOGICAL doTr1ClimRelax
384          LOGICAL periodicExternalForcing
385          LOGICAL fluidIsAir
386          LOGICAL fluidIsWater
387          LOGICAL usingPCoords
388          LOGICAL usingZCoords
389          LOGICAL useDynP_inEos_Zc
390          LOGICAL setCenterDr
391          LOGICAL nonHydrostatic
392          LOGICAL quasiHydrostatic
393          LOGICAL globalFiles
394          LOGICAL useSingleCpuIO
395          LOGICAL allowFreezing
396          LOGICAL useOldFreezing
397          LOGICAL usePickupBeforeC35
398          LOGICAL usePickupBeforeC54
399          LOGICAL startFromPickupAB2
400          LOGICAL debugMode
401          LOGICAL pickup_read_mdsio, pickup_write_mdsio
402          LOGICAL pickup_write_immed
403          LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio
404          LOGICAL outputTypesInclusive
405          LOGICAL inAdMode, inAdTrue, inAdFalse, inAdExact
406          LOGICAL calendarDumps
407    
408  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
409  C     cg2dTargetResidual  C     cg2dTargetResidual
410  C               - Target residual for cg2d solver.  C          :: Target residual for cg2d solver; no unit (RHS normalisation)
411  C     delZ      - Vertical grid spacing ( m ) - delZ is the distance  C     cg2dTargetResWunit
412  C                 between "w" surfaces.  C          :: Target residual for cg2d solver; W unit (No RHS normalisation)
413  C     delX      - Separation between cell faces (m) or (deg), depending  C     cg3dTargetResidual
414    C               :: Target residual for cg3d solver.
415    C     cg2dpcOffDFac :: Averaging weight for preconditioner off-diagonal.
416    C     Note. 20th May 1998
417    C           I made a weird discovery! In the model paper we argue
418    C           for the form of the preconditioner used here ( see
419    C           A Finite-volume, Incompressible Navier-Stokes Model
420    C           ...., Marshall et. al ). The algebra gives a simple
421    C           0.5 factor for the averaging of ac and aCw to get a
422    C           symmettric pre-conditioner. By using a factor of 0.51
423    C           i.e. scaling the off-diagonal terms in the
424    C           preconditioner down slightly I managed to get the
425    C           number of iterations for convergence in a test case to
426    C           drop form 192 -> 134! Need to investigate this further!
427    C           For now I have introduced a parameter cg2dpcOffDFac which
428    C           defaults to 0.51 but can be set at runtime.
429    C     delR      :: Vertical grid spacing ( units of r ).
430    C     delRc     :: Vertical grid spacing between cell centers (r unit).
431    C     delX      :: Separation between cell faces (m) or (deg), depending
432  C     delY        on input flags.  C     delY        on input flags.
433  C     gravity   - Accel. due to gravity ( m/s^2 )  C     gravity   :: Accel. due to gravity ( m/s^2 )
434  C     ronil     - Reference density  C     recip_gravity and its inverse
435  C     startTime - Start time for model ( s )  C     gBaro     :: Accel. due to gravity used in barotropic equation ( m/s^2 )
436  C     phiMin    - Latitude of southern most cell face.  C     rhoNil    :: Reference density for the linear equation of state
437  C     thetaMin  - Longitude of western most cell face (this  C     rhoConst  :: Vertically constant reference density
438    C     rhoConstFresh :: Constant reference density for fresh water (rain)
439    C     tRef      :: reference vertical profile for potential temperature
440    C     sRef      :: reference vertical profile for salinity/specific humidity
441    C     phiMin    :: Latitude of southern most cell face.
442    C     thetaMin  :: Longitude of western most cell face (this
443  C                 is an "inert" parameter but it is included  C                 is an "inert" parameter but it is included
444  C                 to make geographical references simple.)  C                 to make geographical references simple.)
445  C     rSphere   - Radius of sphere for a spherical polar grid ( m ).  C     rSphere   :: Radius of sphere for a spherical polar grid ( m ).
446  C     f0        - Reference coriolis parameter ( 1/s )  C     recip_RSphere  :: Reciprocal radius of sphere ( m ).
447    C     f0        :: Reference coriolis parameter ( 1/s )
448  C                 ( Southern edge f for beta plane )  C                 ( Southern edge f for beta plane )
449  C     beta      - df/dy ( s^-1.m^-1 )  C     beta      :: df/dy ( s^-1.m^-1 )
450  C     viscAh    - Eddy viscosity coeff. for mixing of  C     omega     :: Angular velocity ( rad/s )
451    C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)
452    C     viscAh    :: Eddy viscosity coeff. for mixing of
453    C                 momentum laterally ( m^2/s )
454    C     viscAhW   :: Eddy viscosity coeff. for mixing of vertical
455    C                 momentum laterally, no effect for hydrostatic
456    C                 model, defaults to viscAh if unset ( m^2/s )
457    C     viscAr    :: Eddy viscosity coeff. for mixing of
458    C                 momentum vertically ( units of r^2/s )
459    C     viscA4    :: Biharmonic viscosity coeff. for mixing of
460    C                 momentum laterally ( m^4/s )
461    C     viscA4W   :: Biharmonic viscosity coeff. for mixing of vertical
462    C                 momentum laterally, no effect for hydrostatic
463    C                 model, defaults to viscA4 if unset ( m^2/s )
464    C     viscAhD   :: Eddy viscosity coeff. for mixing of momentum laterally
465    C                  (act on Divergence part) ( m^2/s )
466    C     viscAhZ   :: Eddy viscosity coeff. for mixing of momentum laterally
467    C                  (act on Vorticity  part) ( m^2/s )
468    C     viscA4D   :: Biharmonic viscosity coeff. for mixing of momentum laterally
469    C                  (act on Divergence part) ( m^4/s )
470    C     viscA4Z   :: Biharmonic viscosity coeff. for mixing of momentum laterally
471    C                  (act on Vorticity  part) ( m^4/s )
472    C     viscC2leith :: Leith non-dimensional viscosity factor (grad(vort))
473    C     viscC2leithD :: Modified Leith non-dimensional viscosity factor (grad(div))
474    C     viscC2smag  :: Smagorinsky non-dimensional viscosity factor
475    C     viscAhMax :: Maximum eddy viscosity coeff. for mixing of
476  C                 momentum laterally ( m^2/s )  C                 momentum laterally ( m^2/s )
477  C     viscAz    - Eddy viscosity coeff. for mixing of  C     viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ...
478  C                 momentum vertically ( m^2/s )  C     viscAhGridMin::  in terms of non-dimensional grid-size dependent viscosity
479  C     viscA4    - Biharmonic viscosity coeff. for mixing of  C     viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of
480  C                 momentum laterally ( m^4/s )  C                 momentum laterally ( m^4/s )
481  C     diffKhT   - Laplacian diffusion coeff. for mixing of  C     viscAhGrid:: non-dimensional grid-size dependent viscosity
482    C     viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity
483    C     viscA4GridMax:: maximum and minimum biharmonic viscosity coefficients ...
484    C     viscA4GridMin::  in terms of non-dimensional grid-size dependent viscosity
485    C     viscC4leith :: Leith non-dimensional viscosity factor (grad(vort))
486    C     viscC4leithD :: Modified Leith non-dimensional viscosity factor (grad(div))
487    C     diffKhT   :: Laplacian diffusion coeff. for mixing of
488  C                 heat laterally ( m^2/s )  C                 heat laterally ( m^2/s )
489  C     diffKzT   - Laplacian diffusion coeff. for mixing of  C     diffKrNrT :: vertical profile of Laplacian diffusion coeff.
490  C                 heat vertically ( m^2/s )  C                 for mixing of heat vertically ( units of r^2/s )
491  C     diffK4T   - Biharmonic diffusion coeff. for mixing of  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of
492  C                 heat laterally ( m^4/s )  C                 heat laterally ( m^4/s )
493  C     diffKhS  -  Laplacian diffusion coeff. for mixing of  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of
494  C                 salt laterally ( m^2/s )  C                 salt laterally ( m^2/s )
495  C     diffKzS   - Laplacian diffusion coeff. for mixing of  C     diffKrNrS :: vertical profile of Laplacian diffusion coeff.
496  C                 salt vertically ( m^2/s )  C                 for mixing of salt vertically ( units of r^2/s ),
497  C     diffK4S   - Biharmonic diffusion coeff. for mixing of  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of
498  C                 salt laterally ( m^4/s )  C                 salt laterally ( m^4/s )
499  C     deltaT    - Default timestep ( s )  C     diffKrBL79surf :: T/S surface diffusivity (m^2/s) Bryan and Lewis, 1979
500  C     deltaTMom    - Timestep for momemtum equations ( s )  C     diffKrBL79deep :: T/S deep diffusivity (m^2/s) Bryan and Lewis, 1979
501  C     deltaTtracer - Timestep for tracer equations ( s )  C     diffKrBL79scl  :: depth scale for arctan fn (m) Bryan and Lewis, 1979
502  C     tauCD     - CD scheme coupling timescale ( 1/s )  C     diffKrBL79Ho   :: depth offset for arctan fn (m) Bryan and Lewis, 1979
503  C     rCD       - CD scheme normalised coupling parameter ( 0-1 )  C     deltaT    :: Default timestep ( s )
504  C     GMmaxslope  - max. slope allowed in GM/Redi tensor  C     deltaTClock  :: Timestep used as model "clock". This determines the
505  C     GMlength  - Length to use in Visbeck et al. formula for K  C                    IO frequencies and is used in tagging output. It can
506  C     GMalpha   - alpha to use in Visbeck et al. formula for K  C                    be totally different to the dynamical time. Typically
507  C     GMdepth   - Depth over which to integrate Richardson # (Visbeck et al.)  C                    it will be the deep-water timestep for accelerated runs.
508  C     GMbackground - background value of GM/Redi coefficient  C                    Frequency of checkpointing and dumping of the model state
509  C     startTime - Starting time for this integration ( s ).  C                    are referenced to this clock. ( s )
510  C     endTime   - Ending time for this integration ( s ).  C     deltaTMom    :: Timestep for momemtum equations ( s )
511  C     chkPtFreq - Frequency of check pointing ( s ).  C     dTtracerLev  :: Timestep for tracer equations ( s ), function of level k
512  C     dumpFreq  - Frequency with which model state is written to  C     deltaTfreesurf :: Timestep for free-surface equation ( s )
513  C                 post-processing files ( s ).  C     freesurfFac  :: Parameter to turn implicit free surface term on or off
514        COMMON /PARM_R/ cg2dTargetResidual, delZ, delX, delY, deltaT,  C                    freesurfac = 1. uses implicit free surface
515       & deltaTmom, deltaTtracer, abeps, startTime, phiMin, thetaMin,  C                    freesurfac = 0. uses rigid lid
516       & rSphere, f0, fCori, beta, viscAh, viscAz, viscA4, diffKhT, diffKzT,  C     abEps        :: Adams-Bashforth-2 stabilizing weight
517       & diffK4T, diffKhS, diffKzS, diffK4S, delT, tauCD, rCD,  C     alph_AB      :: Adams-Bashforth-3 primary factor
518       & GMmaxslope,GMlength,GMalpha,GMdepth,GMkbackground,  C     beta_AB      :: Adams-Bashforth-3 secondary factor
519       & gravity, rhonil, tRef, sRef,  C     implicSurfPress :: parameter of the Crank-Nickelson time stepping :
520       & endTime, chkPtFreq, dumpFreq  C                     Implicit part of Surface Pressure Gradient ( 0-1 )
521    C     implicDiv2Dflow :: parameter of the Crank-Nickelson time stepping :
522    C                     Implicit part of barotropic flow Divergence ( 0-1 )
523    C     hFacMin      :: Minimum fraction size of a cell (affects hFacC etc...)
524    C     hFacMinDz    :: Minimum dimesional size of a cell (affects hFacC etc..., m)
525    C     hFacMinDp    :: Minimum dimesional size of a cell (affects hFacC etc..., Pa)
526    C     hFacMinDr    :: Minimum dimesional size of a cell (affects hFacC etc..., units of r)
527    C     hFacInf      :: Threshold (inf and sup) for fraction size of surface cell
528    C     hFacSup        that control vanishing and creating levels
529    C     tauCD        :: CD scheme coupling timescale ( 1/s )
530    C     rCD          :: CD scheme normalised coupling parameter ( 0-1 )
531    C     baseTime      :: model base time (time origin) = time @ iteration zero
532    C     startTime     :: Starting time for this integration ( s ).
533    C     endTime       :: Ending time for this integration ( s ).
534    C     chkPtFreq     :: Frequency of rolling check pointing ( s ).
535    C     pChkPtFreq    :: Frequency of permanent check pointing ( s ).
536    C     dumpFreq      :: Frequency with which model state is written to
537    C                     post-processing files ( s ).
538    C     diagFreq      :: Frequency with which model writes diagnostic output
539    C                     of intermediate quantities.
540    C     afFacMom      :: Advection of momentum term tracer parameter
541    C     vfFacMom      :: Momentum viscosity tracer parameter
542    C     pfFacMom      :: Momentum pressure forcing tracer parameter
543    C     cfFacMom      :: Coriolis term tracer parameter
544    C     foFacMom      :: Momentum forcing tracer parameter
545    C     mtFacMom      :: Metric terms tracer parameter
546    C     cosPower      :: Power of cosine of latitude to multiply viscosity
547    C     cAdjFreq      :: Frequency of convective adjustment
548    C
549    C     taveFreq      :: Frequency with which time-averaged model state
550    C                      is written to post-processing files ( s ).
551    C     tave_lastIter :: (for state variable only) fraction of the last time
552    C                      step (of each taveFreq period) put in the time average.
553    C                      (fraction for 1rst iter = 1 - tave_lastIter)
554    C     tauThetaClimRelax :: Relaxation to climatology time scale ( s ).
555    C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).
556    C     latBandClimRelax :: latitude band where Relaxation to Clim. is applied,
557    C                         i.e. where |yC| <= latBandClimRelax
558    C     externForcingPeriod :: Is the period of which forcing varies (eg. 1 month)
559    C     externForcingCycle :: Is the repeat time of the forcing (eg. 1 year)
560    C                          (note: externForcingCycle must be an integer
561    C                           number times externForcingPeriod)
562    C     convertFW2Salt :: salinity, used to convert Fresh-Water Flux to Salt Flux
563    C                       (use model surface (local) value if set to -1)
564    C     temp_EvPrRn :: temperature of Rain & Evap.
565    C     salt_EvPrRn :: salinity of Rain & Evap.
566    C        (notes: a) tracer content of Rain/Evap only used if both
567    C                     NonLin_FrSurf & useRealFreshWater are set.
568    C                b) use model surface (local) value if set to UNSET_RL)
569    C     horiVertRatio      :: Ratio on units in vertical to units in horizontal.
570    C     recip_horiVertRatio  ( 1 if horiz in m and vertical in m ).
571    C                          ( g*rho if horiz in m and vertical in Pa ).
572    C     Ro_SeaLevel        :: standard position of Sea-Level in "R" coordinate, used as
573    C                          starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)
574    C     bottomDragLinear   :: Drag coefficient built in to core dynamics
575    C      --"-"--  Quadratic  ( linear: 1/s, quadratic: 1/m )
576          COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
577         & cg2dpcOffDFac, cg3dTargetResidual,
578         & delR, delRc, delX, delY,
579         & deltaT, deltaTmom, dTtracerLev, deltaTfreesurf, deltaTClock,
580         & abEps, alph_AB, beta_AB,
581         & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,
582         & viscAh, viscAhW, viscAhMax,
583         & viscAhGrid, viscAhGridMax, viscAhGridMin,
584         & viscC2leith, viscC2leithD,
585         & viscC2smag,
586         & viscAhD, viscAhZ, viscA4D, viscA4Z,
587         & viscA4, viscA4W,
588         & viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin,
589         & viscC4leith, viscC4leithD, viscAr, viscAstrain, viscAtension,
590         & diffKhT, diffK4T, diffKrNrT,
591         & diffKhS, diffK4S, diffKrNrS,
592         & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,
593         & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,
594         & hFacMin, hFacMinDz, hFacInf, hFacSup,
595         & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
596         & recip_rhoConst, rhoConst,
597         & rhoConstFresh, convertEmP2rUnit, tRef, sRef,
598         & baseTime, startTime, endTime,
599         & chkPtFreq, pchkPtFreq, dumpFreq, adjDumpFreq,
600         & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq,
601         & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
602         & cosPower, cAdjFreq, omega, rotationPeriod,
603         & tauThetaClimRelax,
604         & tauSaltClimRelax,
605         & tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax,
606         & externForcingCycle, externForcingPeriod,
607         & convertFW2Salt, temp_EvPrRn, salt_EvPrRn,
608         & hFacMinDr, hFacMinDp,
609         & horiVertRatio, recip_horiVertRatio,
610         & ivdc_kappa, Ro_SeaLevel,
611         & bottomDragLinear,bottomDragQuadratic,nh_Am2,
612         & tCylIn, tCylOut
613    
614          _RL nh_Am2
615        _RL cg2dTargetResidual        _RL cg2dTargetResidual
616        _RL delZ(Nz)        _RL cg2dTargetResWunit
617          _RL cg3dTargetResidual
618          _RL cg2dpcOffDFac
619          _RL delR(Nr)
620          _RL delRc(Nr+1)
621        _RL delX(Nx)        _RL delX(Nx)
622        _RL delY(Ny)        _RL delY(Ny)
623        _RL deltaT        _RL deltaT
624          _RL deltaTClock
625        _RL deltaTmom        _RL deltaTmom
626        _RL deltaTtracer        _RL dTtracerLev(Nr)
627        _RL abeps        _RL deltaTfreesurf
628          _RL abEps, alph_AB, beta_AB
629        _RL phiMin        _RL phiMin
630        _RL thetaMin        _RL thetaMin
631        _RL rSphere        _RL rSphere
632          _RL recip_RSphere
633        _RL f0        _RL f0
634          _RL freeSurfFac
635          _RL implicSurfPress
636          _RL implicDiv2Dflow
637          _RL hFacMin
638          _RL hFacMinDz
639          _RL hFacMinDp
640          _RL hFacMinDr
641          _RL hFacInf
642          _RL hFacSup
643        _RL beta        _RL beta
644        _RL viscAh        _RL viscAh
645        _RL viscAz        _RL viscAhW
646          _RL viscAhD
647          _RL viscAhZ
648          _RL viscAhMax
649          _RL viscAhGrid
650          _RL viscAhGridMax
651          _RL viscAhGridMin
652          _RL viscC2leith
653          _RL viscC2leithD
654          _RL viscC2smag
655          _RL viscAstrain
656          _RL viscAtension
657          _RL viscAr
658        _RL viscA4        _RL viscA4
659          _RL viscA4W
660          _RL viscA4D
661          _RL viscA4Z
662          _RL viscA4Max
663          _RL viscA4Grid, viscA4GridMax, viscA4GridMin
664          _RL viscC4leith
665          _RL viscC4leithD
666        _RL diffKhT        _RL diffKhT
667        _RL diffKzT        _RL diffKrNrT(Nr)
668        _RL diffK4T        _RL diffK4T
669        _RL diffKhS        _RL diffKhS
670        _RL diffKzS        _RL diffKrNrS(Nr)
671        _RL diffK4S        _RL diffK4S
672          _RL diffKrBL79surf
673          _RL diffKrBL79deep
674          _RL diffKrBL79scl
675          _RL diffKrBL79Ho
676        _RL delt        _RL delt
677        _RL tauCD        _RL tauCD
678        _RL rCD        _RL rCD
       _RL GMmaxslope  
       _RL GMlength  
       _RL GMalpha  
       _RL GMdepth  
       _RL GMkbackground  
679        _RL gravity        _RL gravity
680          _RL recip_gravity
681          _RL gBaro
682        _RL rhonil        _RL rhonil
683        _RL tRef(Nz)        _RL recip_rhonil
684        _RL sRef(Nz)        _RL rhoConst
685        _RL Fcori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)        _RL recip_rhoConst
686          _RL rhoConstFresh
687          _RL convertEmP2rUnit
688          _RL tRef(Nr)
689          _RL sRef(Nr)
690          _RL baseTime
691        _RL startTime        _RL startTime
692        _RL endTime        _RL endTime
693        _RL chkPtFreq        _RL chkPtFreq
694          _RL pChkPtFreq
695        _RL dumpFreq        _RL dumpFreq
696          _RL adjDumpFreq
697          _RL diagFreq
698          _RL taveFreq
699          _RL tave_lastIter
700          _RL monitorFreq
701          _RL adjMonitorFreq
702          _RL afFacMom
703          _RL vfFacMom
704          _RL pfFacMom
705          _RL cfFacMom
706          _RL foFacMom
707          _RL mTFacMom
708          _RL cosPower
709          _RL cAdjFreq
710          _RL omega
711          _RL rotationPeriod
712          _RL tauThetaClimRelax
713          _RL tauSaltClimRelax
714          _RL tauTr1ClimRelax
715          _RL lambdaTr1ClimRelax
716          _RL latBandClimRelax
717          _RL externForcingCycle
718          _RL externForcingPeriod
719          _RL convertFW2Salt
720          _RL temp_EvPrRn
721          _RL salt_EvPrRn
722          _RL horiVertRatio
723          _RL recip_horiVertRatio
724          _RL ivdc_kappa
725          _RL Ro_SeaLevel
726          _RL bottomDragLinear
727          _RL bottomDragQuadratic
728          _RL tCylIn
729          _RL tCylOut
730    
731        COMMON /PARM_A/ HeatCapacity_Cp,  C--   COMMON /PARM_A/ Thermodynamics constants ?
732       &                Lamba_theta        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp
733        _RL HeatCapacity_Cp        _RL HeatCapacity_Cp
734        _RL Lamba_theta        _RL recip_Cp
735    
736    C--   COMMON /PARM_ATM/ Atmospheric physical parameters (Ideal Gas EOS, ...)
737    C     celsius2K :: convert centigrade (Celsius) degree to Kelvin
738    C     atm_Po    :: standard reference pressure
739    C     atm_Cp    :: specific heat (Cp) of the (dry) air at constant pressure
740    C     atm_Rd    :: gas constant for dry air
741    C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)
742    C     atm_Rq    :: water vapour specific volume anomaly relative to dry air
743    C                  (e.g. typical value = (29/18 -1) 10^-3 with q [g/kg])
744    C     integr_GeoPot :: option to select the way we integrate the geopotential
745    C                     (still a subject of discussions ...)
746    C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is
747    C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)
748          COMMON /PARM_ATM/
749         &            celsius2K,
750         &            atm_Cp, atm_Rd, atm_kappa, atm_Rq, atm_Po,
751         &            integr_GeoPot, selectFindRoSurf
752          _RL celsius2K
753          _RL atm_Po, atm_Cp, atm_Rd, atm_kappa, atm_Rq
754          INTEGER integr_GeoPot, selectFindRoSurf
755    
756    C Logical flags for selecting packages
757          LOGICAL useOPPS
758          LOGICAL usePP81
759          LOGICAL useMY82
760          LOGICAL useGGL90
761          LOGICAL useKPP
762          LOGICAL useGMRedi
763          LOGICAL useOBCS
764          LOGICAL useAIM
765          LOGICAL useLand
766          LOGICAL useGrdchk
767          LOGICAL useECCO
768          LOGICAL useSHAP_FILT
769          LOGICAL useZONAL_FILT
770          LOGICAL useFLT
771          LOGICAL usePTRACERS
772          LOGICAL useGCHEM
773          LOGICAL useSBO
774          LOGICAL useSEAICE
775          LOGICAL useBulkForce
776          LOGICAL useThSIce
777          LOGICAL usefizhi
778          LOGICAL usegridalt
779          LOGICAL usediagnostics
780          LOGICAL useEBM
781          LOGICAL useMNC
782          LOGICAL useMATRIX
783          LOGICAL useRunClock
784          COMMON /PARM_PACKAGES/
785         &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,
786         &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
787         &        usePTRACERS, useGCHEM,
788         &        useSBO, useSEAICE, useThSIce, useBulkForce,
789         &        usefizhi, usegridalt, usediagnostics, useEBM, useMNC,
790         &        usePP81, useMY82, useOPPS, useGGL90, useMATRIX,
791         &        useRunClock
792    
793  C Equation of State (polynomial coeffients)  CEH3 ;;; Local Variables: ***
794        COMMON /PARM_EOS_NL/ eosC,eosSig0,eosRefT,eosRefS  CEH3 ;;; mode:fortran ***
795        _RL eosC(Nz+1,9),eosSig0(Nz+1),eosRefT(Nz+1),eosRefS(Nz+1)  CEH3 ;;; End: ***
 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  
       _RL tAlpha  
       _RL sBeta  
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.162
  ViewVC Help
Powered by ViewVC 1.1.22