/[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.4 by cnh, Mon Apr 27 04:24:22 1998 UTC revision 1.170 by jmc, Mon Oct 24 22:38:49 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     selectKEscheme      :: Kinetic Energy scheme selector (Vector Inv.)
158    C     debugLevel          :: debug level selector: higher -> more writing
159    
160        COMMON /PARM_I/        COMMON /PARM_I/
161       &        cg2dMaxIters,       &        cg2dMaxIters,
162       &        cg2dChkResFreq,       &        cg2dChkResFreq, cg2dPreCondFreq,
163       &        nIter0, nTimeSteps,       &        cg3dMaxIters,
164       &        numStepsPerPickup       &        cg3dChkResFreq,
165         &        nIter0, nTimeSteps, nEndIter,
166         &        numStepsPerPickup,
167         &        writeStatePrec, nCheckLev,
168         &        writeBinaryPrec, readBinaryPrec,
169         &        nonlinFreeSurf, select_rStar,
170         &        tempAdvScheme, tempVertAdvScheme,
171         &        saltAdvScheme, saltVertAdvScheme,
172         &        selectKEscheme,
173         &        debugLevel
174        INTEGER cg2dMaxIters        INTEGER cg2dMaxIters
175        INTEGER cg2dChkResFreq        INTEGER cg2dChkResFreq
176          INTEGER cg2dPreCondFreq
177          INTEGER cg3dMaxIters
178          INTEGER cg3dChkResFreq
179        INTEGER nIter0        INTEGER nIter0
180        INTEGER nTimeSteps        INTEGER nTimeSteps
181          INTEGER nEndIter
182        INTEGER numStepsPerPickup        INTEGER numStepsPerPickup
183          INTEGER writeStatePrec
184          INTEGER writeBinaryPrec
185          INTEGER readBinaryPrec
186          INTEGER nCheckLev
187          INTEGER nonlinFreeSurf
188          INTEGER select_rStar
189          INTEGER tempAdvScheme, tempVertAdvScheme
190          INTEGER saltAdvScheme, saltVertAdvScheme
191          INTEGER selectKEscheme
192          INTEGER debugLevel
193    
194    C
195          INTEGER debLevZero
196          PARAMETER(debLevZero=0)
197          INTEGER debLevA
198          PARAMETER(debLevA=1)
199          INTEGER debLevB
200          PARAMETER(debLevB=2)
201    
202  C--   COMMON /PARM_L/ Logical valued parameters used by the model.  C--   COMMON /PARM_L/ Logical valued parameters used by the model.
203  C     usingCartesianGrid - If TRUE grid generation will be in a cartesian  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian
204  C                          coordinate frame.  C                          coordinate frame.
205  C     usingSphericalPolarGrid - If TRUE grid generation will be in a  C     usingSphericalPolarGrid :: If TRUE grid generation will be in a
206  C                               spherical polar frame.  C                               spherical polar frame.
207  C     momViscosity  - Flag which turns momentum friction terms on and off.  C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical
208  C     momAdvection  - Flag which turns advection of momentum on and off.  C     no_slip_sides :: Impose "no-slip" at lateral boundaries.
209  C     momForcing    - Flag which turns external forcing of momentum on  C     no_slip_bottom :: Impose "no-slip" at bottom boundary.
210    C     staggerTimeStep :: enable a Stagger time stepping T,S Rho then U,V
211    C     momViscosity  :: Flag which turns momentum friction terms on and off.
212    C     momAdvection  :: Flag which turns advection of momentum on and off.
213    C     momForcing    :: Flag which turns external forcing of momentum on
214  C                     and off.  C                     and off.
215  C     useCoriolis   - Flag which turns the coriolis terms on and off.  C     momPressureForcing :: Flag which turns pressure term in momentum equation
216  C     tempDiffusion - Flag which turns diffusion of temperature on  C                          on and off.
217    C     metricTerms   :: Flag which turns metric terms on or off.
218    C     usingSphericalPolarMTerms :: If TRUE use spherical polar metric terms.
219    C     useNHMTerms :: If TRUE use non-hydrostatic metric terms.
220    C     useCoriolis   :: Flag which turns the coriolis terms on and off.
221    C     tempAdvection :: Flag which turns advection of temperature on
222  C                     and off.  C                     and off.
223  C     tempAdvection - Flag which turns advection of temperature on  C     tempForcing   :: Flag which turns external forcing of temperature on
224  C                     and off.  C                     and off.
225  C     tempForcing   - Flag which turns external forcing of temperature on  C     saltAdvection :: Flag which turns advection of salinity on
226  C                     and off.  C                     and off.
227  C     saltDiffusion - Flag which turns diffusion of salinit on  C     saltForcing   :: Flag which turns external forcing of salinity on
 C                     and off.  
 C     saltAdvection - Flag which turns advection of salinit on  
 C                     and off.  
 C     saltForcing   - Flag which turns external forcing of salinit on  
228  C                     and off.  C                     and off.
229    C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux
230    C                         as a real Fresh Water (=> changes the Sea Level)
231    C                         if F, converts P+R-E to salt flux (no SL effect)
232    C     useFullLeith   :: Set to true to use full Leith viscosity(may be unstable
233    C                       on irregular grids)
234    C     useAreaViscLength :: Set to true to use old scaling for viscous
235    C              lengths, e.g., L2=Raz.  May be preferable for cube sphere.
236    C     useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms
237    C     rigidLid            :: Set to true to use rigid lid
238    C     implicitFreeSurface :: Set to true to use implcit free surface
239    C     exactConserv        :: Set to true to conserve exactly the total Volume
240    C     uniformLin_PhiSurf  :: Set to true to use a uniform Bo_surf in the
241    C                           linear relation Phi_surf = Bo_surf*eta
242    C     momStepping   :: Turns momentum equation time-stepping off
243    C     tempStepping  :: Turns temperature equation time-stepping off
244    C     saltStepping  :: Turns salinity equation time-stepping off
245    C     useConstantF  :: Coriolis parameter set to f0
246    C     useBetaPlaneF :: Coriolis parameter set to f0 + beta.y
247    C     useSphereF    :: Coriolis parameter set to 2.omega.sin(phi)
248    C     useCDscheme   :: use CD-scheme to calculate Coriolis terms.
249    C     useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986)
250    C     useJamartMomAdv :: Use wet-point method for V.I. non-linear term
251    C     SadournyCoriolis :: use the enstrophy conserving scheme by Sadourny
252    C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term
253    C     highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection)
254    C     upwindShear        :: use 1rst order upwind interp. (V.I., vertical advection)
255    C     useAbsVorticity :: work with f+zeta in Coriolis terms
256    C     implicitDiffusion :: Turns implicit vertical diffusion on
257    C     implicitViscosity :: Turns implicit vertical viscosity on
258    C     tempImplVertAdv :: Turns on implicit vertical advection for Temperature
259    C     saltImplVertAdv :: Turns on implicit vertical advection for Salinity
260    C     momImplVertAdv  :: Turns on implicit vertical advection for Momentum
261    C     multiDimAdvection :: Flag that enable multi-dimension advection
262    C     useMultiDimAdvec  :: True if multi-dim advection is used at least once
263    C     forcing_In_AB     :: if False, put forcing (Temp,Salt,Tracers) contribution
264    C                          out off Adams-Bashforth time stepping.
265    C     doAB_onGtGs       :: if the Adams-Bashforth time stepping is used, always
266    C                          apply AB on tracer tendencies (rather than on Tracer)
267    C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup
268    C     doThetaClimRelax :: Set true if relaxation to temperature
269    C                        climatology is required.
270    C     doSaltClimRelax  :: Set true if relaxation to salinity
271    C                        climatology is required.
272    C     periodicExternalForcing :: Set true if forcing is time-dependant
273    C     usingPCoords     :: Set to indicate that we are working in a pressure
274    C                        type coordinate (p or p*).
275    C     usingZCoords     :: Set to indicate that we are working in a height
276    C                        type coordinate (z or z*)
277    C     fluidIsAir       :: Set to indicate that the fluid major constituent
278    C                        is air
279    C     fluidIsWater     :: Set to indicate that the fluid major constituent
280    C                        is water
281    C     useDynP_inEos_Zc :: use the dynamical pressure in EOS (with Z-coord.)
282    C                         this requires specific code for restart & exchange
283    C     setCenterDr    :: set cell Center depth and put Interface at the middle
284    C     nonHydrostatic :: Using non-hydrostatic terms
285    C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm
286    C     globalFiles    :: Selects between "global" and "tiled" files
287    C     useSingleCpuIO :: On SGI platforms, option globalFiles is either
288    C                       slow (f77) or does not work (f90).  When
289    C                       useSingleCpuIO is set, mdsio_writefield.F
290    C                       outputs from master mpi process only.
291    C     allowFreezing  :: Allows surface water to freeze and form ice
292    C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)
293    C     pickup_write_mdsio :: use mdsio to write pickups
294    C     pickup_read_mdsio  :: use mdsio to read  pickups
295    C     pickup_write_immed :: echo the pickup immediately (for conversion)
296    C     timeave_mdsio      :: use mdsio for timeave output
297    C     snapshot_mdsio     :: use mdsio for "snapshot" (dumpfreq/diagfreq) output
298    C     monitor_stdio      :: use stdio for monitor output
299    C     calendarDumps :: When set, approximate months (30-31 days) and years (360-372 days)
300    C                      for parameters chkPtFreq, pChkPtFreq, taveFreq, SEAICE_taveFreq,
301    C                      KPP_taveFreq, and freq in pkg/diagnostics are converted to exact
302    C                      calendar months and years.  Requires pkg/cal.
303    C     dumpInitAndLast :: dumps model state to files at Initial (nIter0)
304    C                        & Last iteration, in addition multiple of dumpFreq iter.
305        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
306       & momViscosity, momAdvection, momForcing, useCoriolis,       & usingCurvilinearGrid, usingCylindricalGrid,
307       & tempDiffusion, tempAdvection, tempForcing,       & no_slip_sides,no_slip_bottom,
308       & saltDiffusion, saltAdvection, saltForcing       & staggerTimeStep,
309         & momViscosity, momAdvection, momForcing, useCoriolis,
310         & momPressureForcing, vectorInvariantMomentum,
311         & tempAdvection, tempForcing,
312         & saltAdvection, saltForcing,
313         & useRealFreshWaterFlux,
314         & useFullLeith, useStrainTensionVisc,
315         & useAreaViscLength,
316         & rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf,
317         & momStepping, tempStepping, saltStepping,
318         & metricTerms, usingSphericalPolarMTerms, useNHMTerms,
319         & useConstantF, useBetaPlaneF, useSphereF,
320         & useCDscheme,
321         & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,
322         & SadournyCoriolis, upwindVorticity, highOrderVorticity,
323         & useAbsVorticity, upwindShear,
324         & implicitDiffusion, implicitViscosity,
325         & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,
326         & multiDimAdvection, useMultiDimAdvec,
327         & forcing_In_AB, doAB_onGtGs,
328         & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,
329         & periodicExternalForcing,
330         & fluidIsAir, fluidIsWater,
331         & usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,
332         & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,
333         & allowFreezing, useOldFreezing,
334         & usePickupBeforeC35, usePickupBeforeC54, startFromPickupAB2,
335         & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,
336         & timeave_mdsio, snapshot_mdsio, monitor_stdio,
337         & outputTypesInclusive, dumpInitAndLast, debugMode,
338         & inAdMode, inAdTrue, inAdFalse, inAdExact,
339         & calendarDumps
340    
341        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
342        LOGICAL usingSphericalPolarGrid        LOGICAL usingSphericalPolarGrid
343          LOGICAL usingCylindricalGrid
344          LOGICAL usingCurvilinearGrid
345          LOGICAL usingSphericalPolarMTerms
346          LOGICAL useNHMTerms
347          LOGICAL no_slip_sides
348          LOGICAL no_slip_bottom
349          LOGICAL staggerTimeStep
350        LOGICAL momViscosity        LOGICAL momViscosity
351        LOGICAL momAdvection        LOGICAL momAdvection
352        LOGICAL momForcing        LOGICAL momForcing
353          LOGICAL momPressureForcing
354        LOGICAL useCoriolis        LOGICAL useCoriolis
355        LOGICAL tempDiffusion        LOGICAL vectorInvariantMomentum
356        LOGICAL tempAdvection        LOGICAL tempAdvection
357        LOGICAL tempForcing        LOGICAL tempForcing
       LOGICAL saltDiffusion  
358        LOGICAL saltAdvection        LOGICAL saltAdvection
359        LOGICAL saltForcing        LOGICAL saltForcing
360          LOGICAL useRealFreshWaterFlux
361          LOGICAL useFullLeith
362          LOGICAL useStrainTensionVisc
363          LOGICAL useAreaViscLength
364          LOGICAL rigidLid
365          LOGICAL implicitFreeSurface
366          LOGICAL exactConserv
367          LOGICAL uniformLin_PhiSurf
368          LOGICAL momStepping
369          LOGICAL tempStepping
370          LOGICAL saltStepping
371          LOGICAL metricTerms
372          LOGICAL useConstantF
373          LOGICAL useBetaPlaneF
374          LOGICAL useSphereF
375          LOGICAL useCDscheme
376          LOGICAL useEnergyConservingCoriolis
377          LOGICAL useJamartWetPoints
378          LOGICAL useJamartMomAdv
379          LOGICAL SadournyCoriolis
380          LOGICAL upwindVorticity
381          LOGICAL highOrderVorticity
382          LOGICAL useAbsVorticity
383          LOGICAL upwindShear
384          LOGICAL implicitDiffusion
385          LOGICAL implicitViscosity
386          LOGICAL tempImplVertAdv
387          LOGICAL saltImplVertAdv
388          LOGICAL momImplVertAdv
389          LOGICAL multiDimAdvection
390          LOGICAL useMultiDimAdvec
391          LOGICAL forcing_In_AB
392          LOGICAL doAB_onGtGs
393          LOGICAL doThetaClimRelax
394          LOGICAL doSaltClimRelax
395          LOGICAL doTr1ClimRelax
396          LOGICAL periodicExternalForcing
397          LOGICAL fluidIsAir
398          LOGICAL fluidIsWater
399          LOGICAL usingPCoords
400          LOGICAL usingZCoords
401          LOGICAL useDynP_inEos_Zc
402          LOGICAL setCenterDr
403          LOGICAL nonHydrostatic
404          LOGICAL quasiHydrostatic
405          LOGICAL globalFiles
406          LOGICAL useSingleCpuIO
407          LOGICAL allowFreezing
408          LOGICAL useOldFreezing
409          LOGICAL usePickupBeforeC35
410          LOGICAL usePickupBeforeC54
411          LOGICAL startFromPickupAB2
412          LOGICAL dumpInitAndLast
413          LOGICAL debugMode
414          LOGICAL pickup_read_mdsio, pickup_write_mdsio
415          LOGICAL pickup_write_immed
416          LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio
417          LOGICAL outputTypesInclusive
418          LOGICAL inAdMode, inAdTrue, inAdFalse, inAdExact
419          LOGICAL calendarDumps
420    
421  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
422  C     cg2dTargetResidual  C     cg2dTargetResidual
423  C               - Target residual for cg2d solver.  C          :: Target residual for cg2d solver; no unit (RHS normalisation)
424  C     delZ      - Vertical grid spacing ( m ) - delZ is the distance  C     cg2dTargetResWunit
425  C                 between "w" surfaces.  C          :: Target residual for cg2d solver; W unit (No RHS normalisation)
426  C     delX      - Separation between cell faces (m) or (deg), depending  C     cg3dTargetResidual
427    C               :: Target residual for cg3d solver.
428    C     cg2dpcOffDFac :: Averaging weight for preconditioner off-diagonal.
429    C     Note. 20th May 1998
430    C           I made a weird discovery! In the model paper we argue
431    C           for the form of the preconditioner used here ( see
432    C           A Finite-volume, Incompressible Navier-Stokes Model
433    C           ...., Marshall et. al ). The algebra gives a simple
434    C           0.5 factor for the averaging of ac and aCw to get a
435    C           symmettric pre-conditioner. By using a factor of 0.51
436    C           i.e. scaling the off-diagonal terms in the
437    C           preconditioner down slightly I managed to get the
438    C           number of iterations for convergence in a test case to
439    C           drop form 192 -> 134! Need to investigate this further!
440    C           For now I have introduced a parameter cg2dpcOffDFac which
441    C           defaults to 0.51 but can be set at runtime.
442    C     delR      :: Vertical grid spacing ( units of r ).
443    C     delRc     :: Vertical grid spacing between cell centers (r unit).
444    C     delX      :: Separation between cell faces (m) or (deg), depending
445  C     delY        on input flags.  C     delY        on input flags.
446  C     gravity   - Accel. due to gravity ( m/s^2 )  C     gravity   :: Accel. due to gravity ( m/s^2 )
447  C     ronil     - Reference density  C     recip_gravity and its inverse
448  C     startTime - Start time for model ( s )  C     gBaro     :: Accel. due to gravity used in barotropic equation ( m/s^2 )
449  C     phiMin    - Latitude of southern most cell face.  C     rhoNil    :: Reference density for the linear equation of state
450  C     thetaMin  - Longitude of western most cell face (this  C     rhoConst  :: Vertically constant reference density
451    C     rhoConstFresh :: Constant reference density for fresh water (rain)
452    C     tRef      :: reference vertical profile for potential temperature
453    C     sRef      :: reference vertical profile for salinity/specific humidity
454    C     phiMin    :: Latitude of southern most cell face.
455    C     thetaMin  :: Longitude of western most cell face (this
456  C                 is an "inert" parameter but it is included  C                 is an "inert" parameter but it is included
457  C                 to make geographical references simple.)  C                 to make geographical references simple.)
458  C     rSphere   - Radius of sphere for a spherical polar grid ( m ).  C     rSphere   :: Radius of sphere for a spherical polar grid ( m ).
459  C     f0        - Reference coriolis parameter ( 1/s )  C     recip_RSphere  :: Reciprocal radius of sphere ( m ).
460    C     f0        :: Reference coriolis parameter ( 1/s )
461  C                 ( Southern edge f for beta plane )  C                 ( Southern edge f for beta plane )
462  C     beta      - df/dy ( s^-1.m^-1 )  C     beta      :: df/dy ( s^-1.m^-1 )
463  C     viscAh    - Eddy viscosity coeff. for mixing of  C     omega     :: Angular velocity ( rad/s )
464    C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)
465    C     viscAh    :: Eddy viscosity coeff. for mixing of
466  C                 momentum laterally ( m^2/s )  C                 momentum laterally ( m^2/s )
467  C     viscAz    - Eddy viscosity coeff. for mixing of  C     viscAhW   :: Eddy viscosity coeff. for mixing of vertical
468  C                 momentum vertically ( m^2/s )  C                 momentum laterally, no effect for hydrostatic
469  C     viscA4    - Biharmonic viscosity coeff. for mixing of  C                 model, defaults to viscAh if unset ( m^2/s )
470    C     viscAr    :: Eddy viscosity coeff. for mixing of
471    C                 momentum vertically ( units of r^2/s )
472    C     viscA4    :: Biharmonic viscosity coeff. for mixing of
473  C                 momentum laterally ( m^4/s )  C                 momentum laterally ( m^4/s )
474  C     diffKhT   - Laplacian diffusion coeff. for mixing of  C     viscA4W   :: Biharmonic viscosity coeff. for mixing of vertical
475    C                 momentum laterally, no effect for hydrostatic
476    C                 model, defaults to viscA4 if unset ( m^2/s )
477    C     viscAhD   :: Eddy viscosity coeff. for mixing of momentum laterally
478    C                  (act on Divergence part) ( m^2/s )
479    C     viscAhZ   :: Eddy viscosity coeff. for mixing of momentum laterally
480    C                  (act on Vorticity  part) ( m^2/s )
481    C     viscA4D   :: Biharmonic viscosity coeff. for mixing of momentum laterally
482    C                  (act on Divergence part) ( m^4/s )
483    C     viscA4Z   :: Biharmonic viscosity coeff. for mixing of momentum laterally
484    C                  (act on Vorticity  part) ( m^4/s )
485    C     viscC2leith :: Leith non-dimensional viscosity factor (grad(vort))
486    C     viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div))
487    C     viscC2smag  :: Smagorinsky non-dimensional viscosity factor (harmonic)
488    C     viscC4smag  :: Smagorinsky non-dimensional viscosity factor (biharmonic)
489    C     viscAhMax :: Maximum eddy viscosity coeff. for mixing of
490    C                   momentum laterally ( m^2/s )
491    C     viscAhReMax :: Maximum gridscale Reynolds number for eddy viscosity
492    C                   coeff. for mixing of momentum laterally (non-dim)
493    C     viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ...
494    C     viscAhGridMin::  in terms of non-dimensional grid-size dependent visc.
495    C     viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of
496    C                 momentum laterally ( m^4/s )
497    C     viscA4ReMax :: Maximum Gridscale Reynolds number for
498    C                  biharmonic viscosity coeff. momentum laterally (non-dim)
499    C     viscAhGrid:: non-dimensional grid-size dependent viscosity
500    C     viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity
501    C     viscA4GridMax:: maximum and minimum biharmonic viscosity coefficients ...
502    C     viscA4GridMin::  in terms of non-dimensional grid-size dependent viscosity
503    C     viscC4leith :: Leith non-dimensional viscosity factor (grad(vort))
504    C     viscC4leithD :: Modified Leith non-dimensional viscosity factor (grad(div))
505    C     diffKhT   :: Laplacian diffusion coeff. for mixing of
506  C                 heat laterally ( m^2/s )  C                 heat laterally ( m^2/s )
507  C     diffKzT   - Laplacian diffusion coeff. for mixing of  C     diffKrNrT :: vertical profile of Laplacian diffusion coeff.
508  C                 heat vertically ( m^2/s )  C                 for mixing of heat vertically ( units of r^2/s )
509  C     diffK4T   - Biharmonic diffusion coeff. for mixing of  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of
510  C                 heat laterally ( m^4/s )  C                 heat laterally ( m^4/s )
511  C     diffKhS  -  Laplacian diffusion coeff. for mixing of  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of
512  C                 salt laterally ( m^2/s )  C                 salt laterally ( m^2/s )
513  C     diffKzS   - Laplacian diffusion coeff. for mixing of  C     diffKrNrS :: vertical profile of Laplacian diffusion coeff.
514  C                 salt vertically ( m^2/s )  C                 for mixing of salt vertically ( units of r^2/s ),
515  C     diffK4S   - Biharmonic diffusion coeff. for mixing of  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of
516  C                 salt laterally ( m^4/s )  C                 salt laterally ( m^4/s )
517  C     tAlpha    - Linear EOS thermal expansion coefficient ( 1/degree ).  C     diffKrBL79surf :: T/S surface diffusivity (m^2/s) Bryan and Lewis, 1979
518  C     sBeta     - Linear EOS haline contraction coefficient.  C     diffKrBL79deep :: T/S deep diffusivity (m^2/s) Bryan and Lewis, 1979
519  C     deltaT    - Default timestep ( s )  C     diffKrBL79scl  :: depth scale for arctan fn (m) Bryan and Lewis, 1979
520  C     deltaTMom    - Timestep for momemtum equations ( s )  C     diffKrBL79Ho   :: depth offset for arctan fn (m) Bryan and Lewis, 1979
521  C     deltaTtracer - Timestep for tracer equations ( s )  C     deltaT    :: Default timestep ( s )
522  C     tauCD     - CD scheme coupling timescale ( 1/s )  C     deltaTClock  :: Timestep used as model "clock". This determines the
523  C     rCD       - CD scheme normalised coupling parameter ( 0-1 )  C                    IO frequencies and is used in tagging output. It can
524  C     startTime - Starting time for this integration ( s ).  C                    be totally different to the dynamical time. Typically
525  C     endTime   - Ending time for this integration ( s ).  C                    it will be the deep-water timestep for accelerated runs.
526  C     chkPtFreq - Frequency of check pointing ( s ).  C                    Frequency of checkpointing and dumping of the model state
527  C     dumpFreq  - Frequency with which model state is written to  C                    are referenced to this clock. ( s )
528  C                 post-processing files ( s ).  C     deltaTMom    :: Timestep for momemtum equations ( s )
529        COMMON /PARM_R/ cg2dTargetResidual, delZ, delX, delY, deltaT,  C     dTtracerLev  :: Timestep for tracer equations ( s ), function of level k
530       & deltaTmom, deltaTtracer, abeps, startTime, phiMin, thetaMin,  C     deltaTfreesurf :: Timestep for free-surface equation ( s )
531       & rSphere, f0, fCori, beta, viscAh, viscAz, viscA4, diffKhT, diffKzT,  C     freesurfFac  :: Parameter to turn implicit free surface term on or off
532       & diffK4T, diffKhS, diffKzS, diffK4S, delT, tauCD, rCD,  C                    freesurfac = 1. uses implicit free surface
533       & tAlpha, sBeta, gravity, rhonil, tRef, sRef,  C                    freesurfac = 0. uses rigid lid
534       & endTime, chkPtFreq, dumpFreq  C     abEps        :: Adams-Bashforth-2 stabilizing weight
535        REAL cg2dTargetResidual  C     alph_AB      :: Adams-Bashforth-3 primary factor
536        REAL delZ(Nz)  C     beta_AB      :: Adams-Bashforth-3 secondary factor
537        REAL delX(Nx)  C     implicSurfPress :: parameter of the Crank-Nickelson time stepping :
538        REAL delY(Ny)  C                     Implicit part of Surface Pressure Gradient ( 0-1 )
539        REAL deltaT  C     implicDiv2Dflow :: parameter of the Crank-Nickelson time stepping :
540        REAL deltaTmom  C                     Implicit part of barotropic flow Divergence ( 0-1 )
541        REAL deltaTtracer  C     hFacMin      :: Minimum fraction size of a cell (affects hFacC etc...)
542        REAL abeps  C     hFacMinDz    :: Minimum dimesional size of a cell (affects hFacC etc..., m)
543        REAL phiMin  C     hFacMinDp    :: Minimum dimesional size of a cell (affects hFacC etc..., Pa)
544        REAL thetaMin  C     hFacMinDr    :: Minimum dimesional size of a cell (affects hFacC etc..., units of r)
545        REAL rSphere  C     hFacInf      :: Threshold (inf and sup) for fraction size of surface cell
546        REAL f0  C     hFacSup        that control vanishing and creating levels
547        _RL  beta  C     tauCD        :: CD scheme coupling timescale ( 1/s )
548        REAL viscAh  C     rCD          :: CD scheme normalised coupling parameter ( 0-1 )
549        REAL viscAz  C     baseTime      :: model base time (time origin) = time @ iteration zero
550        REAL viscA4  C     startTime     :: Starting time for this integration ( s ).
551        REAL diffKhT  C     endTime       :: Ending time for this integration ( s ).
552        REAL diffKzT  C     chkPtFreq     :: Frequency of rolling check pointing ( s ).
553        REAL diffK4T  C     pChkPtFreq    :: Frequency of permanent check pointing ( s ).
554        REAL diffKhS  C     dumpFreq      :: Frequency with which model state is written to
555        REAL diffKzS  C                     post-processing files ( s ).
556        REAL diffK4S  C     diagFreq      :: Frequency with which model writes diagnostic output
557        REAL delt  C                     of intermediate quantities.
558        REAL tauCD  C     afFacMom      :: Advection of momentum term tracer parameter
559        REAL rCD  C     vfFacMom      :: Momentum viscosity tracer parameter
560        REAL tAlpha  C     pfFacMom      :: Momentum pressure forcing tracer parameter
561        REAL sBeta  C     cfFacMom      :: Coriolis term tracer parameter
562        REAL gravity  C     foFacMom      :: Momentum forcing tracer parameter
563        REAL rhonil  C     mtFacMom      :: Metric terms tracer parameter
564        REAL tRef(Nz)  C     cosPower      :: Power of cosine of latitude to multiply viscosity
565        REAL sRef(Nz)  C     cAdjFreq      :: Frequency of convective adjustment
566        real Fcori(1-OLx:sNx+OLx,1-OLy:sNy+OLy,nSx,nSy)  C
567        REAL startTime  C     taveFreq      :: Frequency with which time-averaged model state
568        REAL endTime  C                      is written to post-processing files ( s ).
569        REAL chkPtFreq  C     tave_lastIter :: (for state variable only) fraction of the last time
570        REAL dumpFreq  C                      step (of each taveFreq period) put in the time average.
571        COMMON /PARM_A/ HeatCapacity_Cp,  C                      (fraction for 1rst iter = 1 - tave_lastIter)
572       &                Lamba_theta  C     tauThetaClimRelax :: Relaxation to climatology time scale ( s ).
573        REAL HeatCapacity_Cp  C     tauSaltClimRelax :: Relaxation to climatology time scale ( s ).
574        REAL Lamba_theta  C     latBandClimRelax :: latitude band where Relaxation to Clim. is applied,
575    C                         i.e. where |yC| <= latBandClimRelax
576    C     externForcingPeriod :: Is the period of which forcing varies (eg. 1 month)
577    C     externForcingCycle :: Is the repeat time of the forcing (eg. 1 year)
578    C                          (note: externForcingCycle must be an integer
579    C                           number times externForcingPeriod)
580    C     convertFW2Salt :: salinity, used to convert Fresh-Water Flux to Salt Flux
581    C                       (use model surface (local) value if set to -1)
582    C     temp_EvPrRn :: temperature of Rain & Evap.
583    C     salt_EvPrRn :: salinity of Rain & Evap.
584    C        (notes: a) tracer content of Rain/Evap only used if both
585    C                     NonLin_FrSurf & useRealFreshWater are set.
586    C                b) use model surface (local) value if set to UNSET_RL)
587    C     horiVertRatio      :: Ratio on units in vertical to units in horizontal.
588    C     recip_horiVertRatio  ( 1 if horiz in m and vertical in m ).
589    C                          ( g*rho if horiz in m and vertical in Pa ).
590    C     Ro_SeaLevel        :: standard position of Sea-Level in "R" coordinate, used as
591    C                          starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)
592    C     sideDragFactor     :: side-drag scaling factor (used only if no_slip_sides)
593    C                           (default=2: full drag ; =1: gives half-slip BC)
594    C     bottomDragLinear   :: Drag coefficient built in to core dynamics
595    C     bottomDragQuadratic   ( linear: 1/s, quadratic: 1/m )
596    C     nh_Am2        :: scales the non-hydrostatic terms and changes internal scales
597    C                      (i.e. allows convection at different Rayleigh numbers)
598          COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
599         & cg2dpcOffDFac, cg3dTargetResidual,
600         & delR, delRc, delX, delY,
601         & deltaT, deltaTmom, dTtracerLev, deltaTfreesurf, deltaTClock,
602         & abEps, alph_AB, beta_AB,
603         & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,
604         & viscAh, viscAhW, viscAhMax,
605         & viscAhGrid, viscAhGridMax, viscAhGridMin,
606         & viscC2leith, viscC2leithD,
607         & viscC2smag, viscC4smag,
608         & viscAhD, viscAhZ, viscA4D, viscA4Z,
609         & viscA4, viscA4W,
610         & viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin,
611         & viscAhRemax, viscA4Remax,
612         & viscC4leith, viscC4leithD, viscAr,
613         & diffKhT, diffK4T, diffKrNrT,
614         & diffKhS, diffK4S, diffKrNrS,
615         & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,
616         & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,
617         & hFacMin, hFacMinDz, hFacInf, hFacSup,
618         & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
619         & recip_rhoConst, rhoConst,
620         & rhoConstFresh, convertEmP2rUnit, tRef, sRef,
621         & baseTime, startTime, endTime,
622         & chkPtFreq, pchkPtFreq, dumpFreq, adjDumpFreq,
623         & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq,
624         & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
625         & cosPower, cAdjFreq, omega, rotationPeriod,
626         & tauThetaClimRelax,
627         & tauSaltClimRelax,
628         & tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax,
629         & externForcingCycle, externForcingPeriod,
630         & convertFW2Salt, temp_EvPrRn, salt_EvPrRn,
631         & hFacMinDr, hFacMinDp,
632         & horiVertRatio, recip_horiVertRatio,
633         & ivdc_kappa, Ro_SeaLevel,
634         & sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2,
635         & tCylIn, tCylOut
636    
637          _RL cg2dTargetResidual
638          _RL cg2dTargetResWunit
639          _RL cg3dTargetResidual
640          _RL cg2dpcOffDFac
641          _RL delR(Nr)
642          _RL delRc(Nr+1)
643          _RL delX(Nx)
644          _RL delY(Ny)
645          _RL deltaT
646          _RL deltaTClock
647          _RL deltaTmom
648          _RL dTtracerLev(Nr)
649          _RL deltaTfreesurf
650          _RL abEps, alph_AB, beta_AB
651          _RL phiMin
652          _RL thetaMin
653          _RL rSphere
654          _RL recip_RSphere
655          _RL f0
656          _RL freeSurfFac
657          _RL implicSurfPress
658          _RL implicDiv2Dflow
659          _RL hFacMin
660          _RL hFacMinDz
661          _RL hFacMinDp
662          _RL hFacMinDr
663          _RL hFacInf
664          _RL hFacSup
665          _RL beta
666          _RL viscAh
667          _RL viscAhW
668          _RL viscAhD
669          _RL viscAhZ
670          _RL viscAhMax
671          _RL viscAhReMax
672          _RL viscAhGrid
673          _RL viscAhGridMax
674          _RL viscAhGridMin
675          _RL viscC2leith
676          _RL viscC2leithD
677          _RL viscC2smag
678          _RL viscC4smag
679          _RL viscAr
680          _RL viscA4
681          _RL viscA4W
682          _RL viscA4D
683          _RL viscA4Z
684          _RL viscA4Max
685          _RL viscA4ReMax
686          _RL viscA4Grid, viscA4GridMax, viscA4GridMin
687          _RL viscC4leith
688          _RL viscC4leithD
689          _RL diffKhT
690          _RL diffKrNrT(Nr)
691          _RL diffK4T
692          _RL diffKhS
693          _RL diffKrNrS(Nr)
694          _RL diffK4S
695          _RL diffKrBL79surf
696          _RL diffKrBL79deep
697          _RL diffKrBL79scl
698          _RL diffKrBL79Ho
699          _RL delt
700          _RL tauCD
701          _RL rCD
702          _RL gravity
703          _RL recip_gravity
704          _RL gBaro
705          _RL rhonil
706          _RL recip_rhonil
707          _RL rhoConst
708          _RL recip_rhoConst
709          _RL rhoConstFresh
710          _RL convertEmP2rUnit
711          _RL tRef(Nr)
712          _RL sRef(Nr)
713          _RL baseTime
714          _RL startTime
715          _RL endTime
716          _RL chkPtFreq
717          _RL pChkPtFreq
718          _RL dumpFreq
719          _RL adjDumpFreq
720          _RL diagFreq
721          _RL taveFreq
722          _RL tave_lastIter
723          _RL monitorFreq
724          _RL adjMonitorFreq
725          _RL afFacMom
726          _RL vfFacMom
727          _RL pfFacMom
728          _RL cfFacMom
729          _RL foFacMom
730          _RL mTFacMom
731          _RL cosPower
732          _RL cAdjFreq
733          _RL omega
734          _RL rotationPeriod
735          _RL tauThetaClimRelax
736          _RL tauSaltClimRelax
737          _RL tauTr1ClimRelax
738          _RL lambdaTr1ClimRelax
739          _RL latBandClimRelax
740          _RL externForcingCycle
741          _RL externForcingPeriod
742          _RL convertFW2Salt
743          _RL temp_EvPrRn
744          _RL salt_EvPrRn
745          _RL horiVertRatio
746          _RL recip_horiVertRatio
747          _RL ivdc_kappa
748          _RL Ro_SeaLevel
749          _RL sideDragFactor
750          _RL bottomDragLinear
751          _RL bottomDragQuadratic
752          _RL nh_Am2
753          _RL tCylIn
754          _RL tCylOut
755    
756    C--   COMMON /PARM_A/ Thermodynamics constants ?
757          COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp
758          _RL HeatCapacity_Cp
759          _RL recip_Cp
760    
761    C--   COMMON /PARM_ATM/ Atmospheric physical parameters (Ideal Gas EOS, ...)
762    C     celsius2K :: convert centigrade (Celsius) degree to Kelvin
763    C     atm_Po    :: standard reference pressure
764    C     atm_Cp    :: specific heat (Cp) of the (dry) air at constant pressure
765    C     atm_Rd    :: gas constant for dry air
766    C     atm_kappa :: kappa = R/Cp (R: constant of Ideal Gas EOS)
767    C     atm_Rq    :: water vapour specific volume anomaly relative to dry air
768    C                  (e.g. typical value = (29/18 -1) 10^-3 with q [g/kg])
769    C     integr_GeoPot :: option to select the way we integrate the geopotential
770    C                     (still a subject of discussions ...)
771    C     selectFindRoSurf :: select the way surf. ref. pressure (=Ro_surf) is
772    C             derived from the orography. Implemented: 0,1 (see INI_P_GROUND)
773          COMMON /PARM_ATM/
774         &            celsius2K,
775         &            atm_Cp, atm_Rd, atm_kappa, atm_Rq, atm_Po,
776         &            integr_GeoPot, selectFindRoSurf
777          _RL celsius2K
778          _RL atm_Po, atm_Cp, atm_Rd, atm_kappa, atm_Rq
779          INTEGER integr_GeoPot, selectFindRoSurf
780    
781    C Logical flags for selecting packages
782          LOGICAL useOPPS
783          LOGICAL usePP81
784          LOGICAL useMY82
785          LOGICAL useGGL90
786          LOGICAL useKPP
787          LOGICAL useGMRedi
788          LOGICAL useOBCS
789          LOGICAL useAIM
790          LOGICAL useLand
791          LOGICAL useGrdchk
792          LOGICAL useECCO
793          LOGICAL useSHAP_FILT
794          LOGICAL useZONAL_FILT
795          LOGICAL useFLT
796          LOGICAL usePTRACERS
797          LOGICAL useGCHEM
798          LOGICAL useSBO
799          LOGICAL useSEAICE
800          LOGICAL useBulkForce
801          LOGICAL useThSIce
802          LOGICAL usefizhi
803          LOGICAL usegridalt
804          LOGICAL usediagnostics
805          LOGICAL useEBM
806          LOGICAL useMNC
807          LOGICAL useMATRIX
808          LOGICAL useRunClock
809          COMMON /PARM_PACKAGES/
810         &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,
811         &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
812         &        usePTRACERS, useGCHEM,
813         &        useSBO, useSEAICE, useThSIce, useBulkForce,
814         &        usefizhi, usegridalt, usediagnostics, useEBM, useMNC,
815         &        usePP81, useMY82, useOPPS, useGGL90, useMATRIX,
816         &        useRunClock
817    
818    CEH3 ;;; Local Variables: ***
819    CEH3 ;;; mode:fortran ***
820    CEH3 ;;; End: ***
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.170
  ViewVC Help
Powered by ViewVC 1.1.22