/[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.223 by jmc, Fri Sep 5 19:55:56 2008 UTC revision 1.278 by jmc, Sat Jan 3 23:53:50 2015 UTC
# Line 14  C     here. Notes describing the paramet Line 14  C     here. Notes describing the paramet
14    
15  CEOP  CEOP
16    
 C     Macros for special grid options  
 #include "PARAMS_MACROS.h"  
   
17  C--   Contants  C--   Contants
18  C     Useful physical values  C     Useful physical values
19        Real*8 PI        Real*8 PI
# Line 33  C     pickupSuff      :: force to start Line 30  C     pickupSuff      :: force to start
30  C                        and read pickup files with this suffix (max 10 Char.)  C                        and read pickup files with this suffix (max 10 Char.)
31  C     mdsioLocalDir   :: read-write tiled file from/to this directory name  C     mdsioLocalDir   :: read-write tiled file from/to this directory name
32  C                        (+ 4 digits Processor-Rank) instead of current dir.  C                        (+ 4 digits Processor-Rank) instead of current dir.
33    C     adTapeDir       :: read-write checkpointing tape files from/to this
34    C                        directory name instead of current dir. Conflicts
35    C                        mdsioLocalDir, so only one of the two can be set.
36    C                        In contrast to mdsioLocalDir, if specified adTapeDir
37    C                        must exist before the model starts.
38  C     tRefFile      :: File containing reference Potential Temperat.  tRef (1.D)  C     tRefFile      :: File containing reference Potential Temperat.  tRef (1.D)
39  C     sRefFile      :: File containing reference salinity/spec.humid. sRef (1.D)  C     sRefFile      :: File containing reference salinity/spec.humid. sRef (1.D)
40  C     rhoRefFile    :: File containing reference density profile rhoRef (1.D)  C     rhoRefFile    :: File containing reference density profile rhoRef (1.D)
41  C     delRFile      :: File containing vertical grid spacing delR  (1.D array)  C     delRFile      :: File containing vertical grid spacing delR  (1.D array)
42  C     delRcFile     :: File containing vertical grid spacing delRc (1.D array)  C     delRcFile     :: File containing vertical grid spacing delRc (1.D array)
43    C     hybSigmFile   :: File containing hybrid-sigma vertical coord. coeff. (2x 1.D)
44  C     delXFile      :: File containing X-spacing grid definition (1.D array)  C     delXFile      :: File containing X-spacing grid definition (1.D array)
45  C     delYFile      :: File containing Y-spacing grid definition (1.D array)  C     delYFile      :: File containing Y-spacing grid definition (1.D array)
46  C     horizGridFile :: File containing horizontal-grid definition  C     horizGridFile :: File containing horizontal-grid definition
# Line 46  C     bathyFile       :: File containing Line 49  C     bathyFile       :: File containing
49  C                        is taken from inline function.  C                        is taken from inline function.
50  C     topoFile        :: File containing the topography of the surface (unit=m)  C     topoFile        :: File containing the topography of the surface (unit=m)
51  C                        (mainly used for the atmosphere = ground height).  C                        (mainly used for the atmosphere = ground height).
52  C     shelfIceFile    :: File containing the topography of the shelfice draught  C     addWwallFile    :: File containing 2-D additional Western  cell-edge wall
53  C                        (unit=m)  C     addSwallFile    :: File containing 2-D additional Southern cell-edge wall
54    C                        (e.g., to add "thin-wall" where it is =1)
55  C     hydrogThetaFile :: File containing initial hydrographic data (3-D)  C     hydrogThetaFile :: File containing initial hydrographic data (3-D)
56  C                        for potential temperature.  C                        for potential temperature.
57  C     hydrogSaltFile  :: File containing initial hydrographic data (3-D)  C     hydrogSaltFile  :: File containing initial hydrographic data (3-D)
58  C                        for salinity.  C                        for salinity.
59  C     diffKrFile      :: File containing 3D specification of vertical diffusivity  C     diffKrFile      :: File containing 3D specification of vertical diffusivity
60    C     viscAhDfile     :: File containing 3D specification of horizontal viscosity
61    C     viscAhZfile     :: File containing 3D specification of horizontal viscosity
62    C     viscA4Dfile     :: File containing 3D specification of horizontal viscosity
63    C     viscA4Zfile     :: File containing 3D specification of horizontal viscosity
64  C     zonalWindFile   :: File containing zonal wind data  C     zonalWindFile   :: File containing zonal wind data
65  C     meridWindFile   :: File containing meridional wind data  C     meridWindFile   :: File containing meridional wind data
66  C     thetaClimFile   :: File containing surface theta climataology used  C     thetaClimFile   :: File containing surface theta climataology used
67  C                       in relaxation term -lambda(theta-theta*)  C                        in relaxation term -lambda(theta-theta*)
68  C     saltClimFile    :: File containing surface salt climataology used  C     saltClimFile    :: File containing surface salt climataology used
69  C                       in relaxation term -lambda(salt-salt*)  C                        in relaxation term -lambda(salt-salt*)
70  C     surfQfile       :: File containing surface heat flux, excluding SW  C     surfQfile       :: File containing surface heat flux, excluding SW
71  C                        (old version, kept for backward compatibility)  C                        (old version, kept for backward compatibility)
72  C     surfQnetFile    :: File containing surface net heat flux  C     surfQnetFile    :: File containing surface net heat flux
73  C     surfQswFile     :: File containing surface shortwave radiation  C     surfQswFile     :: File containing surface shortwave radiation
 C     dQdTfile        :: File containing thermal relaxation coefficient  
74  C     EmPmRfile       :: File containing surface fresh water flux  C     EmPmRfile       :: File containing surface fresh water flux
75  C           NOTE: for backward compatibility EmPmRfile is specified in  C           NOTE: for backward compatibility EmPmRfile is specified in
76  C                 m/s when using external_fields_load.F.  It is converted  C                 m/s when using external_fields_load.F.  It is converted
77  C                 to kg/m2/s by multiplying by rhoConstFresh.  C                 to kg/m2/s by multiplying by rhoConstFresh.
78  C     saltFluxFile    :: File containing surface salt flux  C     saltFluxFile    :: File containing surface salt flux
79  C     pLoadFile       :: File containing pressure loading  C     pLoadFile       :: File containing pressure loading
80    C     addMassFile     :: File containing source/sink of fluid in the interior
81  C     eddyPsiXFile    :: File containing zonal Eddy streamfunction data  C     eddyPsiXFile    :: File containing zonal Eddy streamfunction data
82  C     eddyPsiYFile    :: File containing meridional Eddy streamfunction data  C     eddyPsiYFile    :: File containing meridional Eddy streamfunction data
83  C     the_run_name    :: string identifying the name of the model "run"  C     the_run_name    :: string identifying the name of the model "run"
84        COMMON /PARM_C/        COMMON /PARM_C/
85       &                buoyancyRelation, eosType,       &                buoyancyRelation, eosType,
86       &                pickupSuff, mdsioLocalDir,       &                pickupSuff, mdsioLocalDir, adTapeDir,
87       &                tRefFile, sRefFile, rhoRefFile,       &                tRefFile, sRefFile, rhoRefFile,
88       &                delRFile, delRcFile,       &                delRFile, delRcFile, hybSigmFile,
89       &                delXFile, delYFile, horizGridFile,       &                delXFile, delYFile, horizGridFile,
90       &                bathyFile, topoFile, shelfIceFile,       &                bathyFile, topoFile, addWwallFile, addSwallFile,
91         &                viscAhDfile, viscAhZfile,
92         &                viscA4Dfile, viscA4Zfile,
93       &                hydrogThetaFile, hydrogSaltFile, diffKrFile,       &                hydrogThetaFile, hydrogSaltFile, diffKrFile,
94       &                zonalWindFile, meridWindFile, thetaClimFile,       &                zonalWindFile, meridWindFile, thetaClimFile,
95       &                saltClimFile,       &                saltClimFile,
# Line 87  C     the_run_name    :: string identify Line 97  C     the_run_name    :: string identify
97       &                surfQfile, surfQnetFile, surfQswFile,       &                surfQfile, surfQnetFile, surfQswFile,
98       &                lambdaThetaFile, lambdaSaltFile,       &                lambdaThetaFile, lambdaSaltFile,
99       &                uVelInitFile, vVelInitFile, pSurfInitFile,       &                uVelInitFile, vVelInitFile, pSurfInitFile,
100       &                dQdTfile, ploadFile,       &                pLoadFile, addMassFile,
101       &                eddyPsiXFile, eddyPsiYFile,       &                eddyPsiXFile, eddyPsiYFile, geothermalFile,
102       &                the_run_name       &                the_run_name
103        CHARACTER*(MAX_LEN_FNAM) buoyancyRelation        CHARACTER*(MAX_LEN_FNAM) buoyancyRelation
104        CHARACTER*(6)  eosType        CHARACTER*(6)  eosType
105        CHARACTER*(10) pickupSuff        CHARACTER*(10) pickupSuff
106        CHARACTER*(MAX_LEN_FNAM) mdsioLocalDir        CHARACTER*(MAX_LEN_FNAM) mdsioLocalDir
107          CHARACTER*(MAX_LEN_FNAM) adTapeDir
108        CHARACTER*(MAX_LEN_FNAM) tRefFile        CHARACTER*(MAX_LEN_FNAM) tRefFile
109        CHARACTER*(MAX_LEN_FNAM) sRefFile        CHARACTER*(MAX_LEN_FNAM) sRefFile
110        CHARACTER*(MAX_LEN_FNAM) rhoRefFile        CHARACTER*(MAX_LEN_FNAM) rhoRefFile
111        CHARACTER*(MAX_LEN_FNAM) delRFile        CHARACTER*(MAX_LEN_FNAM) delRFile
112        CHARACTER*(MAX_LEN_FNAM) delRcFile        CHARACTER*(MAX_LEN_FNAM) delRcFile
113          CHARACTER*(MAX_LEN_FNAM) hybSigmFile
114        CHARACTER*(MAX_LEN_FNAM) delXFile        CHARACTER*(MAX_LEN_FNAM) delXFile
115        CHARACTER*(MAX_LEN_FNAM) delYFile        CHARACTER*(MAX_LEN_FNAM) delYFile
116        CHARACTER*(MAX_LEN_FNAM) horizGridFile        CHARACTER*(MAX_LEN_FNAM) horizGridFile
117        CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile, shelfIceFile        CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile
118          CHARACTER*(MAX_LEN_FNAM) addWwallFile, addSwallFile
119        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile, hydrogSaltFile        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile, hydrogSaltFile
120        CHARACTER*(MAX_LEN_FNAM) diffKrFile        CHARACTER*(MAX_LEN_FNAM) diffKrFile
121          CHARACTER*(MAX_LEN_FNAM) viscAhDfile
122          CHARACTER*(MAX_LEN_FNAM) viscAhZfile
123          CHARACTER*(MAX_LEN_FNAM) viscA4Dfile
124          CHARACTER*(MAX_LEN_FNAM) viscA4Zfile
125        CHARACTER*(MAX_LEN_FNAM) zonalWindFile        CHARACTER*(MAX_LEN_FNAM) zonalWindFile
126        CHARACTER*(MAX_LEN_FNAM) meridWindFile        CHARACTER*(MAX_LEN_FNAM) meridWindFile
127        CHARACTER*(MAX_LEN_FNAM) thetaClimFile        CHARACTER*(MAX_LEN_FNAM) thetaClimFile
# Line 117  C     the_run_name    :: string identify Line 134  C     the_run_name    :: string identify
134        CHARACTER*(MAX_LEN_FNAM) uVelInitFile        CHARACTER*(MAX_LEN_FNAM) uVelInitFile
135        CHARACTER*(MAX_LEN_FNAM) vVelInitFile        CHARACTER*(MAX_LEN_FNAM) vVelInitFile
136        CHARACTER*(MAX_LEN_FNAM) pSurfInitFile        CHARACTER*(MAX_LEN_FNAM) pSurfInitFile
137        CHARACTER*(MAX_LEN_FNAM) dQdTfile        CHARACTER*(MAX_LEN_FNAM) pLoadFile
138        CHARACTER*(MAX_LEN_FNAM) ploadFile        CHARACTER*(MAX_LEN_FNAM) addMassFile
139        CHARACTER*(MAX_LEN_FNAM) eddyPsiXFile        CHARACTER*(MAX_LEN_FNAM) eddyPsiXFile
140        CHARACTER*(MAX_LEN_FNAM) eddyPsiYFile        CHARACTER*(MAX_LEN_FNAM) eddyPsiYFile
141          CHARACTER*(MAX_LEN_FNAM) geothermalFile
142        CHARACTER*(MAX_LEN_FNAM) lambdaThetaFile        CHARACTER*(MAX_LEN_FNAM) lambdaThetaFile
143        CHARACTER*(MAX_LEN_FNAM) lambdaSaltFile        CHARACTER*(MAX_LEN_FNAM) lambdaSaltFile
144        CHARACTER*(MAX_LEN_PREC/2) the_run_name        CHARACTER*(MAX_LEN_PREC/2) the_run_name
145    
146  C--   COMMON /PARM_I/ Integer valued parameters used by the model.  C--   COMMON /PARM_I/ Integer valued parameters used by the model.
147  C     cg2dMaxIters        :: Maximum number of iterations in the  C     cg2dMaxIters        :: Maximum number of iterations in the
148  C                           two-dimensional con. grad solver.  C                            two-dimensional con. grad solver.
149  C     cg2dChkResFreq      :: Frequency with which to check residual  C     cg2dChkResFreq      :: Frequency with which to check residual
150  C                           in con. grad solver.  C                            in con. grad solver.
151  C     cg2dPreCondFreq     :: Frequency for updating cg2d preconditioner  C     cg2dPreCondFreq     :: Frequency for updating cg2d preconditioner
152  C                            (non-linear free-surf.)  C                            (non-linear free-surf.)
153    C     cg2dUseMinResSol    :: =0 : use last-iteration/converged solution
154    C                            =1 : use solver minimum-residual solution
155  C     cg3dMaxIters        :: Maximum number of iterations in the  C     cg3dMaxIters        :: Maximum number of iterations in the
156  C                           three-dimensional con. grad solver.  C                            three-dimensional con. grad solver.
157  C     cg3dChkResFreq      :: Frequency with which to check residual  C     cg3dChkResFreq      :: Frequency with which to check residual
158  C                           in con. grad solver.  C                            in con. grad solver.
159    C     printResidualFreq   :: Frequency for printing residual in CG iterations
160  C     nIter0              :: Start time-step number of for this run  C     nIter0              :: Start time-step number of for this run
161  C     nTimeSteps          :: Number of timesteps to execute  C     nTimeSteps          :: Number of timesteps to execute
162  C     writeStatePrec      :: Precision used for writing model state.  C     writeStatePrec      :: Precision used for writing model state.
163  C     writeBinaryPrec     :: Precision used for writing binary files  C     writeBinaryPrec     :: Precision used for writing binary files
164  C     readBinaryPrec      :: Precision used for reading binary files  C     readBinaryPrec      :: Precision used for reading binary files
165    C     selectCoriMap       :: select setting of Coriolis parameter map:
166    C                           =0 f-Plane (Constant Coriolis, = f0)
167    C                           =1 Beta-Plane Coriolis (= f0 + beta.y)
168    C                           =2 Spherical Coriolis (= 2.omega.sin(phi))
169    C                           =3 Read Coriolis 2-d fields from files.
170    C     selectSigmaCoord    :: option related to sigma vertical coordinate
171  C     nonlinFreeSurf      :: option related to non-linear free surface  C     nonlinFreeSurf      :: option related to non-linear free surface
172  C                           =0 Linear free surface ; >0 Non-linear  C                           =0 Linear free surface ; >0 Non-linear
173  C     select_rStar        :: option related to r* vertical coordinate  C     select_rStar        :: option related to r* vertical coordinate
174  C                           =0 (default) use r coord. ; > 0 use r*  C                           =0 (default) use r coord. ; > 0 use r*
175    C     selectNHfreeSurf    :: option for Non-Hydrostatic (free-)Surface formulation:
176    C                           =0 (default) hydrostatic surf. ; > 0 add NH effects.
177  C     selectAddFluid      :: option to add mass source/sink of fluid in the interior  C     selectAddFluid      :: option to add mass source/sink of fluid in the interior
178  C                            (3-D generalisation of oceanic real-fresh water flux)  C                            (3-D generalisation of oceanic real-fresh water flux)
179  C                           =0 off ; =1 add fluid ; =-1 virtual flux (no mass added)  C                           =0 off ; =1 add fluid ; =-1 virtual flux (no mass added)
180  C     momForcingOutAB     :: =1: take momentum forcing contribution  C     momForcingOutAB     :: =1: take momentum forcing contribution
181  C                           out of (=0: in) Adams-Bashforth time stepping.  C                            out of (=0: in) Adams-Bashforth time stepping.
182  C     tracForcingOutAB    :: =1: take tracer (Temp,Salt,pTracers) forcing contribution  C     tracForcingOutAB    :: =1: take tracer (Temp,Salt,pTracers) forcing contribution
183  C                           out of (=0: in) Adams-Bashforth time stepping.  C                            out of (=0: in) Adams-Bashforth time stepping.
184  C     tempAdvScheme       :: Temp. Horiz.Advection scheme selector  C     tempAdvScheme       :: Temp. Horiz.Advection scheme selector
185  C     tempVertAdvScheme   :: Temp. Vert. Advection scheme selector  C     tempVertAdvScheme   :: Temp. Vert. Advection scheme selector
186  C     saltAdvScheme       :: Salt. Horiz.advection scheme selector  C     saltAdvScheme       :: Salt. Horiz.advection scheme selector
187  C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector  C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector
188  C     selectKEscheme      :: Kinetic Energy scheme selector (Vector Inv.)  C     selectKEscheme      :: Kinetic Energy scheme selector (Vector Inv.)
189  C     selectVortScheme    :: Scheme selector for Vorticity term (Vector Inv.)  C     selectVortScheme    :: Scheme selector for Vorticity term (Vector Inv.)
190    C     selectBotDragQuadr  :: quadratic bottom drag discretisation option:
191    C                           =0: average KE from grid center to U & V location
192    C                           =1: use local velocity norm @ U & V location
193    C                           =2: same with wet-point averaging of other component
194  C     monitorSelect       :: select group of variables to monitor  C     monitorSelect       :: select group of variables to monitor
195  C                            =1 : dynvars ; =2 : + vort ; =3 : + surface  C                            =1 : dynvars ; =2 : + vort ; =3 : + surface
196  C     debugLevel          :: debug level selector: higher -> more writing  C-    debugLevel          :: controls printing of algorithm intermediate results
197    C                            and statistics ; higher -> more writing
198    
199        COMMON /PARM_I/        COMMON /PARM_I/
200       &        cg2dMaxIters,       &        cg2dMaxIters, cg2dChkResFreq,
201       &        cg2dChkResFreq, cg2dPreCondFreq,       &        cg2dPreCondFreq, cg2dUseMinResSol,
202       &        cg3dMaxIters,       &        cg3dMaxIters, cg3dChkResFreq,
203       &        cg3dChkResFreq,       &        printResidualFreq,
204       &        nIter0, nTimeSteps, nEndIter,       &        nIter0, nTimeSteps, nEndIter,
205       &        writeStatePrec,       &        writeStatePrec,
206       &        writeBinaryPrec, readBinaryPrec,       &        writeBinaryPrec, readBinaryPrec,
207         &        selectCoriMap,
208         &        selectSigmaCoord,
209       &        nonlinFreeSurf, select_rStar,       &        nonlinFreeSurf, select_rStar,
210         &        selectNHfreeSurf,
211       &        selectAddFluid,       &        selectAddFluid,
212       &        momForcingOutAB, tracForcingOutAB,       &        momForcingOutAB, tracForcingOutAB,
213       &        tempAdvScheme, tempVertAdvScheme,       &        tempAdvScheme, tempVertAdvScheme,
214       &        saltAdvScheme, saltVertAdvScheme,       &        saltAdvScheme, saltVertAdvScheme,
215       &        selectKEscheme, selectVortScheme,       &        selectKEscheme, selectVortScheme,
216         &        selectBotDragQuadr,
217       &        monitorSelect, debugLevel       &        monitorSelect, debugLevel
218        INTEGER cg2dMaxIters        INTEGER cg2dMaxIters
219        INTEGER cg2dChkResFreq        INTEGER cg2dChkResFreq
220        INTEGER cg2dPreCondFreq        INTEGER cg2dPreCondFreq
221          INTEGER cg2dUseMinResSol
222        INTEGER cg3dMaxIters        INTEGER cg3dMaxIters
223        INTEGER cg3dChkResFreq        INTEGER cg3dChkResFreq
224          INTEGER printResidualFreq
225        INTEGER nIter0        INTEGER nIter0
226        INTEGER nTimeSteps        INTEGER nTimeSteps
227        INTEGER nEndIter        INTEGER nEndIter
228        INTEGER writeStatePrec        INTEGER writeStatePrec
229        INTEGER writeBinaryPrec        INTEGER writeBinaryPrec
230        INTEGER readBinaryPrec        INTEGER readBinaryPrec
231          INTEGER selectCoriMap
232          INTEGER selectSigmaCoord
233        INTEGER nonlinFreeSurf        INTEGER nonlinFreeSurf
234        INTEGER select_rStar        INTEGER select_rStar
235          INTEGER selectNHfreeSurf
236        INTEGER selectAddFluid        INTEGER selectAddFluid
237        INTEGER momForcingOutAB, tracForcingOutAB        INTEGER momForcingOutAB, tracForcingOutAB
238        INTEGER tempAdvScheme, tempVertAdvScheme        INTEGER tempAdvScheme, tempVertAdvScheme
239        INTEGER saltAdvScheme, saltVertAdvScheme        INTEGER saltAdvScheme, saltVertAdvScheme
240        INTEGER selectKEscheme        INTEGER selectKEscheme
241        INTEGER selectVortScheme        INTEGER selectVortScheme
242          INTEGER selectBotDragQuadr
243        INTEGER monitorSelect        INTEGER monitorSelect
244        INTEGER debugLevel        INTEGER debugLevel
245    
246  C--   COMMON /PARM_L/ Logical valued parameters used by the model.  C--   COMMON /PARM_L/ Logical valued parameters used by the model.
247  C- Coordinate + Grid params:  C- Coordinate + Grid params:
248  C     fluidIsAir       :: Set to indicate that the fluid major constituent  C     fluidIsAir       :: Set to indicate that the fluid major constituent
249  C                        is air  C                         is air
250  C     fluidIsWater     :: Set to indicate that the fluid major constituent  C     fluidIsWater     :: Set to indicate that the fluid major constituent
251  C                        is water  C                         is water
252  C     usingPCoords     :: Set to indicate that we are working in a pressure  C     usingPCoords     :: Set to indicate that we are working in a pressure
253  C                        type coordinate (p or p*).  C                         type coordinate (p or p*).
254  C     usingZCoords     :: Set to indicate that we are working in a height  C     usingZCoords     :: Set to indicate that we are working in a height
255  C                        type coordinate (z or z*)  C                         type coordinate (z or z*)
256  C     useDynP_inEos_Zc :: use the dynamical pressure in EOS (with Z-coord.)  C     useDynP_inEos_Zc :: use the dynamical pressure in EOS (with Z-coord.)
257  C                         this requires specific code for restart & exchange  C                         this requires specific code for restart & exchange
258  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian
259  C                          coordinate frame.  C                           coordinate frame.
260  C     usingSphericalPolarGrid :: If TRUE grid generation will be in a  C     usingSphericalPolarGrid :: If TRUE grid generation will be in a
261  C                               spherical polar frame.  C                                spherical polar frame.
262  C     rotateGrid      :: rotate grid coordinates to geographical coordinates  C     rotateGrid      :: rotate grid coordinates to geographical coordinates
263  C                        according to Euler angles phiEuler, thetaEuler, psiEuler  C                        according to Euler angles phiEuler, thetaEuler, psiEuler
 C     usingCurvilinearGrid :: If TRUE, use a curvilinear grid (to be provided)  
264  C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical  C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical
265    C     usingCurvilinearGrid :: If TRUE, use a curvilinear grid (to be provided)
266    C     hasWetCSCorners :: domain contains CS-type corners where dynamics is solved
267  C     deepAtmosphere :: deep model (drop the shallow-atmosphere approximation)  C     deepAtmosphere :: deep model (drop the shallow-atmosphere approximation)
268  C     setInterFDr    :: set Interface depth (put cell-Center at the middle)  C     setInterFDr    :: set Interface depth (put cell-Center at the middle)
269  C     setCenterDr    :: set cell-Center depth (put Interface at the middle)  C     setCenterDr    :: set cell-Center depth (put Interface at the middle)
270  C- Momentum params:  C- Momentum params:
271  C     no_slip_sides  :: Impose "no-slip" at lateral boundaries.  C     no_slip_sides  :: Impose "no-slip" at lateral boundaries.
272  C     no_slip_bottom :: Impose "no-slip" at bottom boundary.  C     no_slip_bottom :: Impose "no-slip" at bottom boundary.
273    C     bottomVisc_pCell :: account for partial-cell in bottom visc. (no-slip BC)
274    C     useSmag3D      :: Use isotropic 3-D Smagorinsky
275  C     useFullLeith   :: Set to true to use full Leith viscosity(may be unstable  C     useFullLeith   :: Set to true to use full Leith viscosity(may be unstable
276  C                       on irregular grids)  C                       on irregular grids)
277  C     useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms  C     useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms
# Line 240  C     metricTerms   :: Flag which turns Line 287  C     metricTerms   :: Flag which turns
287  C     useNHMTerms   :: If TRUE use non-hydrostatic metric terms.  C     useNHMTerms   :: If TRUE use non-hydrostatic metric terms.
288  C     useCoriolis   :: Flag which turns the coriolis terms on and off.  C     useCoriolis   :: Flag which turns the coriolis terms on and off.
289  C     use3dCoriolis :: Turns the 3-D coriolis terms (in Omega.cos Phi) on - off  C     use3dCoriolis :: Turns the 3-D coriolis terms (in Omega.cos Phi) on - off
 C     useConstantF  :: Coriolis parameter set to f0  
 C     useBetaPlaneF :: Coriolis parameter set to f0 + beta.y  
 C     useSphereF    :: Coriolis parameter set to 2.omega.sin(phi)  
290  C     useCDscheme   :: use CD-scheme to calculate Coriolis terms.  C     useCDscheme   :: use CD-scheme to calculate Coriolis terms.
291  C     vectorInvariantMomentum :: use Vector-Invariant form (mom_vecinv package)  C     vectorInvariantMomentum :: use Vector-Invariant form (mom_vecinv package)
292  C                                (default = F = use mom_fluxform package)  C                                (default = F = use mom_fluxform package)
# Line 253  C     highOrderVorticity :: use 3rd/4th Line 297  C     highOrderVorticity :: use 3rd/4th
297  C     useAbsVorticity :: work with f+zeta in Coriolis terms  C     useAbsVorticity :: work with f+zeta in Coriolis terms
298  C     upwindShear        :: use 1rst order upwind interp. (V.I., vertical advection)  C     upwindShear        :: use 1rst order upwind interp. (V.I., vertical advection)
299  C     momStepping    :: Turns momentum equation time-stepping off  C     momStepping    :: Turns momentum equation time-stepping off
300    C     calc_wVelocity :: Turns of vertical velocity calculation off
301  C- Temp. & Salt params:  C- Temp. & Salt params:
302  C     tempStepping   :: Turns temperature equation time-stepping off  C     tempStepping   :: Turns temperature equation time-stepping on/off
303  C     saltStepping   :: Turns salinity equation time-stepping off  C     saltStepping   :: Turns salinity equation time-stepping on/off
304    C     addFrictionHeating :: account for frictional heating
305  C     tempAdvection  :: Flag which turns advection of temperature on and off.  C     tempAdvection  :: Flag which turns advection of temperature on and off.
306    C     tempVertDiff4  :: use vertical bi-harmonic diffusion for temperature
307  C     tempIsActiveTr :: Pot.Temp. is a dynamically active tracer  C     tempIsActiveTr :: Pot.Temp. is a dynamically active tracer
308  C     tempForcing    :: Flag which turns external forcing of temperature on  C     tempForcing    :: Flag which turns external forcing of temperature on/off
 C                       and off.  
309  C     saltAdvection  :: Flag which turns advection of salinity on and off.  C     saltAdvection  :: Flag which turns advection of salinity on and off.
310    C     saltVertDiff4  :: use vertical bi-harmonic diffusion for salinity
311  C     saltIsActiveTr :: Salinity  is a dynamically active tracer  C     saltIsActiveTr :: Salinity  is a dynamically active tracer
312  C     saltForcing    :: Flag which turns external forcing of salinity on  C     saltForcing    :: Flag which turns external forcing of salinity on/off
313  C                       and off.  C     maskIniTemp    :: apply mask to initial Pot.Temp.
314  C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux  C     maskIniSalt    :: apply mask to initial salinity
315  C                         as a real Fresh Water (=> changes the Sea Level)  C     checkIniTemp   :: check for points with identically zero initial Pot.Temp.
316  C                         if F, converts P+R-E to salt flux (no SL effect)  C     checkIniSalt   :: check for points with identically zero initial salinity
317  C- Time-stepping params:  C- Pressure solver related parameters (PARM02)
318    C     useSRCGSolver  :: Set to true to use conjugate gradient
319    C                       solver with single reduction (only one call of
320    C                       s/r mpi_allreduce), default is false
321    C- Time-stepping & free-surface params:
322  C     rigidLid            :: Set to true to use rigid lid  C     rigidLid            :: Set to true to use rigid lid
323  C     implicitFreeSurface :: Set to true to use implicit free surface  C     implicitFreeSurface :: Set to true to use implicit free surface
324    C     uniformLin_PhiSurf  :: Set to true to use a uniform Bo_surf in the
325    C                            linear relation Phi_surf = Bo_surf*eta
326    C     uniformFreeSurfLev  :: TRUE if free-surface level-index is uniform (=1)
327  C     exactConserv        :: Set to true to conserve exactly the total Volume  C     exactConserv        :: Set to true to conserve exactly the total Volume
328  C     linFSConserveTr     :: Set to true to correct source/sink of tracer  C     linFSConserveTr     :: Set to true to correct source/sink of tracer
329  C                            at the surface due to Linear Free Surface  C                            at the surface due to Linear Free Surface
330  C     uniformLin_PhiSurf  :: Set to true to use a uniform Bo_surf in the  C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux
331  C                            linear relation Phi_surf = Bo_surf*eta  C                         as a real Fresh Water (=> changes the Sea Level)
332    C                         if F, converts P+R-E to salt flux (no SL effect)
333  C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm  C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm
334  C     nonHydrostatic   :: Using non-hydrostatic algorithm  C     nonHydrostatic   :: Using non-hydrostatic algorithm
335  C     use3Dsolver      :: set to true to use 3-D pressure solver  C     use3Dsolver      :: set to true to use 3-D pressure solver
336  C     implicitIntGravWave :: treat Internal Gravity Wave implicitly  C     implicitIntGravWave :: treat Internal Gravity Wave implicitly
337  C     staggerTimeStep   :: enable a Stagger time stepping U,V (& W) then T,S  C     staggerTimeStep   :: enable a Stagger time stepping U,V (& W) then T,S
338    C     doResetHFactors   :: Do reset thickness factors @ beginning of each time-step
339  C     implicitDiffusion :: Turns implicit vertical diffusion on  C     implicitDiffusion :: Turns implicit vertical diffusion on
340  C     implicitViscosity :: Turns implicit vertical viscosity on  C     implicitViscosity :: Turns implicit vertical viscosity on
341  C     tempImplVertAdv :: Turns on implicit vertical advection for Temperature  C     implBottomFriction :: Turns on implicit bottom friction (drag & no-slip BC)
342  C     saltImplVertAdv :: Turns on implicit vertical advection for Salinity  C     tempImplVertAdv   :: Turns on implicit vertical advection for Temperature
343  C     momImplVertAdv  :: Turns on implicit vertical advection for Momentum  C     saltImplVertAdv   :: Turns on implicit vertical advection for Salinity
344    C     momImplVertAdv    :: Turns on implicit vertical advection for Momentum
345  C     multiDimAdvection :: Flag that enable multi-dimension advection  C     multiDimAdvection :: Flag that enable multi-dimension advection
346  C     useMultiDimAdvec  :: True if multi-dim advection is used at least once  C     useMultiDimAdvec  :: True if multi-dim advection is used at least once
347  C     momDissip_In_AB   :: if False, put Dissipation tendency contribution  C     momDissip_In_AB   :: if False, put Dissipation tendency contribution
# Line 299  C     doThetaClimRelax :: Set true if re Line 356  C     doThetaClimRelax :: Set true if re
356  C                        climatology is required.  C                        climatology is required.
357  C     doSaltClimRelax  :: Set true if relaxation to salinity  C     doSaltClimRelax  :: Set true if relaxation to salinity
358  C                        climatology is required.  C                        climatology is required.
359    C     balanceThetaClimRelax :: substract global mean effect at every time step
360    C     balanceSaltClimRelax :: substract global mean effect at every time step
361  C     allowFreezing  :: Allows surface water to freeze and form ice  C     allowFreezing  :: Allows surface water to freeze and form ice
 C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)  
362  C     periodicExternalForcing :: Set true if forcing is time-dependant  C     periodicExternalForcing :: Set true if forcing is time-dependant
363  C- I/O parameters -  C- I/O parameters -
364  C     globalFiles    :: Selects between "global" and "tiled" files  C     globalFiles    :: Selects between "global" and "tiled" files.
365  C     useSingleCpuIO :: On SGI platforms, option globalFiles is either  C                       On some platforms with MPI, option globalFiles is either
366  C                       slow (f77) or does not work (f90).  When  C                       slow or does not work. Use useSingleCpuIO instead.
367  C                       useSingleCpuIO is set, mdsio_writefield.F  C     useSingleCpuIO :: moved to EEPARAMS.h
 C                       outputs from master mpi process only.  
368  C     pickupStrictlyMatch :: check and stop if pickup-file do not stricly match  C     pickupStrictlyMatch :: check and stop if pickup-file do not stricly match
369  C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup  C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup
370  C     usePickupBeforeC54 :: start from old-pickup files, generated with code from  C     usePickupBeforeC54 :: start from old-pickup files, generated with code from
# Line 321  C     snapshot_mdsio     :: use mdsio fo Line 378  C     snapshot_mdsio     :: use mdsio fo
378  C     monitor_stdio      :: use stdio for monitor output  C     monitor_stdio      :: use stdio for monitor output
379  C     dumpInitAndLast :: dumps model state to files at Initial (nIter0)  C     dumpInitAndLast :: dumps model state to files at Initial (nIter0)
380  C                        & Last iteration, in addition multiple of dumpFreq iter.  C                        & Last iteration, in addition multiple of dumpFreq iter.
381    C     printDomain     :: controls printing of domain fields (bathy, hFac ...).
382    
383        COMMON /PARM_L/        COMMON /PARM_L/
384       & fluidIsAir, fluidIsWater,       & fluidIsAir, fluidIsWater,
385       & usingPCoords, usingZCoords, useDynP_inEos_Zc,       & usingPCoords, usingZCoords, useDynP_inEos_Zc,
386       & usingCartesianGrid, usingSphericalPolarGrid, rotateGrid,       & usingCartesianGrid, usingSphericalPolarGrid, rotateGrid,
387       & usingCurvilinearGrid, usingCylindricalGrid,       & usingCylindricalGrid, usingCurvilinearGrid, hasWetCSCorners,
388       & deepAtmosphere, setInterFDr, setCenterDr,       & deepAtmosphere, setInterFDr, setCenterDr,
389       & no_slip_sides, no_slip_bottom,       & no_slip_sides, no_slip_bottom, bottomVisc_pCell, useSmag3D,
390       & useFullLeith, useStrainTensionVisc, useAreaViscLength,       & useFullLeith, useStrainTensionVisc, useAreaViscLength,
391       & momViscosity, momAdvection, momForcing,       & momViscosity, momAdvection, momForcing,
392       & momPressureForcing, metricTerms, useNHMTerms,       & momPressureForcing, metricTerms, useNHMTerms,
393       & useCoriolis, use3dCoriolis,       & useCoriolis, use3dCoriolis,
      & useConstantF, useBetaPlaneF, useSphereF,  
394       & useCDscheme, vectorInvariantMomentum,       & useCDscheme, vectorInvariantMomentum,
395       & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,       & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,
396       & upwindVorticity, highOrderVorticity,       & upwindVorticity, highOrderVorticity,
397       & useAbsVorticity, upwindShear,       & useAbsVorticity, upwindShear,
398       & momStepping, tempStepping, saltStepping,       & momStepping, calc_wVelocity, tempStepping, saltStepping,
399       & tempAdvection, tempIsActiveTr, tempForcing,       & addFrictionHeating,
400       & saltAdvection, saltIsActiveTr, saltForcing,       & tempAdvection, tempVertDiff4, tempIsActiveTr, tempForcing,
401       & useRealFreshWaterFlux,       & saltAdvection, saltVertDiff4, saltIsActiveTr, saltForcing,
402       & rigidLid, implicitFreeSurface, exactConserv, linFSConserveTr,       & maskIniTemp, maskIniSalt, checkIniTemp, checkIniSalt,
403       & uniformLin_PhiSurf,       & useSRCGSolver,
404       & quasiHydrostatic, nonHydrostatic,       & rigidLid, implicitFreeSurface,
405       & use3Dsolver, implicitIntGravWave, staggerTimeStep,       & uniformLin_PhiSurf, uniformFreeSurfLev,
406       & implicitDiffusion, implicitViscosity,       & exactConserv, linFSConserveTr, useRealFreshWaterFlux,
407         & quasiHydrostatic, nonHydrostatic, use3Dsolver,
408         & implicitIntGravWave, staggerTimeStep, doResetHFactors,
409         & implicitDiffusion, implicitViscosity, implBottomFriction,
410       & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,       & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,
411       & multiDimAdvection, useMultiDimAdvec,       & multiDimAdvection, useMultiDimAdvec,
412       & momDissip_In_AB, doAB_onGtGs,       & momDissip_In_AB, doAB_onGtGs,
413       & balanceEmPmR, balanceQnet, balancePrintMean,       & balanceEmPmR, balanceQnet, balancePrintMean,
414         & balanceThetaClimRelax, balanceSaltClimRelax,
415       & doThetaClimRelax, doSaltClimRelax,       & doThetaClimRelax, doSaltClimRelax,
416       & allowFreezing, useOldFreezing,       & allowFreezing,
417       & periodicExternalForcing,       & periodicExternalForcing,
418       & globalFiles, useSingleCpuIO,       & globalFiles,
419       & pickupStrictlyMatch, usePickupBeforeC54, startFromPickupAB2,       & pickupStrictlyMatch, usePickupBeforeC54, startFromPickupAB2,
420       & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,       & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,
421       & writePickupAtEnd,       & writePickupAtEnd,
422       & timeave_mdsio, snapshot_mdsio, monitor_stdio,       & timeave_mdsio, snapshot_mdsio, monitor_stdio,
423       & outputTypesInclusive, dumpInitAndLast, debugMode,       & outputTypesInclusive, dumpInitAndLast,
424       & inAdMode, inAdTrue, inAdFalse, inAdExact       & printDomain
425    
426        LOGICAL fluidIsAir        LOGICAL fluidIsAir
427        LOGICAL fluidIsWater        LOGICAL fluidIsWater
# Line 370  C                        & Last iteratio Line 431  C                        & Last iteratio
431        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
432        LOGICAL usingSphericalPolarGrid, rotateGrid        LOGICAL usingSphericalPolarGrid, rotateGrid
433        LOGICAL usingCylindricalGrid        LOGICAL usingCylindricalGrid
434        LOGICAL usingCurvilinearGrid        LOGICAL usingCurvilinearGrid, hasWetCSCorners
435        LOGICAL deepAtmosphere        LOGICAL deepAtmosphere
436        LOGICAL setInterFDr        LOGICAL setInterFDr
437        LOGICAL setCenterDr        LOGICAL setCenterDr
438        LOGICAL useNHMTerms  
439        LOGICAL no_slip_sides        LOGICAL no_slip_sides
440        LOGICAL no_slip_bottom        LOGICAL no_slip_bottom
441          LOGICAL bottomVisc_pCell
442          LOGICAL useSmag3D
443          LOGICAL useFullLeith
444          LOGICAL useStrainTensionVisc
445          LOGICAL useAreaViscLength
446        LOGICAL momViscosity        LOGICAL momViscosity
447        LOGICAL momAdvection        LOGICAL momAdvection
448        LOGICAL momForcing        LOGICAL momForcing
449        LOGICAL momPressureForcing        LOGICAL momPressureForcing
450          LOGICAL metricTerms
451          LOGICAL useNHMTerms
452    
453        LOGICAL useCoriolis        LOGICAL useCoriolis
454          LOGICAL use3dCoriolis
455          LOGICAL useCDscheme
456        LOGICAL vectorInvariantMomentum        LOGICAL vectorInvariantMomentum
457          LOGICAL useEnergyConservingCoriolis
458          LOGICAL useJamartWetPoints
459          LOGICAL useJamartMomAdv
460          LOGICAL upwindVorticity
461          LOGICAL highOrderVorticity
462          LOGICAL useAbsVorticity
463          LOGICAL upwindShear
464          LOGICAL momStepping
465          LOGICAL calc_wVelocity
466          LOGICAL tempStepping
467          LOGICAL saltStepping
468          LOGICAL addFrictionHeating
469        LOGICAL tempAdvection        LOGICAL tempAdvection
470          LOGICAL tempVertDiff4
471        LOGICAL tempIsActiveTr        LOGICAL tempIsActiveTr
472        LOGICAL tempForcing        LOGICAL tempForcing
473        LOGICAL saltAdvection        LOGICAL saltAdvection
474          LOGICAL saltVertDiff4
475        LOGICAL saltIsActiveTr        LOGICAL saltIsActiveTr
476        LOGICAL saltForcing        LOGICAL saltForcing
477        LOGICAL useRealFreshWaterFlux        LOGICAL maskIniTemp
478        LOGICAL useFullLeith        LOGICAL maskIniSalt
479        LOGICAL useStrainTensionVisc        LOGICAL checkIniTemp
480        LOGICAL useAreaViscLength        LOGICAL checkIniSalt
481          LOGICAL useSRCGSolver
482        LOGICAL rigidLid        LOGICAL rigidLid
483        LOGICAL implicitFreeSurface        LOGICAL implicitFreeSurface
484          LOGICAL uniformLin_PhiSurf
485          LOGICAL uniformFreeSurfLev
486        LOGICAL exactConserv        LOGICAL exactConserv
487        LOGICAL linFSConserveTr        LOGICAL linFSConserveTr
488        LOGICAL uniformLin_PhiSurf        LOGICAL useRealFreshWaterFlux
489        LOGICAL quasiHydrostatic        LOGICAL quasiHydrostatic
490        LOGICAL nonHydrostatic        LOGICAL nonHydrostatic
491        LOGICAL use3Dsolver        LOGICAL use3Dsolver
492        LOGICAL implicitIntGravWave        LOGICAL implicitIntGravWave
493        LOGICAL staggerTimeStep        LOGICAL staggerTimeStep
494        LOGICAL momStepping        LOGICAL doResetHFactors
       LOGICAL tempStepping  
       LOGICAL saltStepping  
       LOGICAL metricTerms  
       LOGICAL useConstantF  
       LOGICAL useBetaPlaneF  
       LOGICAL useSphereF  
       LOGICAL use3dCoriolis  
       LOGICAL useCDscheme  
       LOGICAL useEnergyConservingCoriolis  
       LOGICAL useJamartWetPoints  
       LOGICAL useJamartMomAdv  
       LOGICAL upwindVorticity  
       LOGICAL highOrderVorticity  
       LOGICAL useAbsVorticity  
       LOGICAL upwindShear  
495        LOGICAL implicitDiffusion        LOGICAL implicitDiffusion
496        LOGICAL implicitViscosity        LOGICAL implicitViscosity
497          LOGICAL implBottomFriction
498        LOGICAL tempImplVertAdv        LOGICAL tempImplVertAdv
499        LOGICAL saltImplVertAdv        LOGICAL saltImplVertAdv
500        LOGICAL momImplVertAdv        LOGICAL momImplVertAdv
# Line 433  C                        & Last iteratio Line 507  C                        & Last iteratio
507        LOGICAL balancePrintMean        LOGICAL balancePrintMean
508        LOGICAL doThetaClimRelax        LOGICAL doThetaClimRelax
509        LOGICAL doSaltClimRelax        LOGICAL doSaltClimRelax
510          LOGICAL balanceThetaClimRelax
511          LOGICAL balanceSaltClimRelax
512        LOGICAL allowFreezing        LOGICAL allowFreezing
       LOGICAL useOldFreezing  
513        LOGICAL periodicExternalForcing        LOGICAL periodicExternalForcing
514        LOGICAL globalFiles        LOGICAL globalFiles
       LOGICAL useSingleCpuIO  
515        LOGICAL pickupStrictlyMatch        LOGICAL pickupStrictlyMatch
516        LOGICAL usePickupBeforeC54        LOGICAL usePickupBeforeC54
517        LOGICAL startFromPickupAB2        LOGICAL startFromPickupAB2
# Line 446  C                        & Last iteratio Line 520  C                        & Last iteratio
520        LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio        LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio
521        LOGICAL outputTypesInclusive        LOGICAL outputTypesInclusive
522        LOGICAL dumpInitAndLast        LOGICAL dumpInitAndLast
523        LOGICAL debugMode        LOGICAL printDomain
       LOGICAL inAdMode, inAdTrue, inAdFalse, inAdExact  
524    
525  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
526  C     cg2dTargetResidual  C     cg2dTargetResidual
# Line 473  C           defaults to 0.51 but can be Line 546  C           defaults to 0.51 but can be
546  C     delR      :: Vertical grid spacing ( units of r ).  C     delR      :: Vertical grid spacing ( units of r ).
547  C     delRc     :: Vertical grid spacing between cell centers (r unit).  C     delRc     :: Vertical grid spacing between cell centers (r unit).
548  C     delX      :: Separation between cell faces (m) or (deg), depending  C     delX      :: Separation between cell faces (m) or (deg), depending
549  C     delY        on input flags.  C     delY         on input flags. Note: moved to header file SET_GRID.h
550  C     gravity   :: Accel. due to gravity ( m/s^2 )  C     xgOrigin   :: Origin of the X-axis (Cartesian Grid) / Longitude of Western
551    C                :: most cell face (Lat-Lon grid) (Note: this is an "inert"
552    C                :: parameter but it makes geographical references simple.)
553    C     ygOrigin   :: Origin of the Y-axis (Cartesian Grid) / Latitude of Southern
554    C                :: most face (Lat-Lon grid).
555    C     gravity    :: Accel. due to gravity ( m/s^2 )
556  C     recip_gravity and its inverse  C     recip_gravity and its inverse
557  C     gBaro     :: Accel. due to gravity used in barotropic equation ( m/s^2 )  C     gBaro      :: Accel. due to gravity used in barotropic equation ( m/s^2 )
558  C     rhoNil    :: Reference density for the linear equation of state  C     rhoNil     :: Reference density for the linear equation of state
559  C     rhoConst  :: Vertically constant reference density  C     rhoConst   :: Vertically constant reference density (Boussinesq)
560  C     rhoFacC   :: normalized (by rhoConst) reference density at cell-Center  C     thetaConst :: Constant reference for potential temperature
561  C     rhoFacF   :: normalized (by rhoConst) reference density at cell-interFace  C     rhoFacC    :: normalized (by rhoConst) reference density at cell-Center
562    C     rhoFacF    :: normalized (by rhoConst) reference density at cell-interFace
563  C     rhoConstFresh :: Constant reference density for fresh water (rain)  C     rhoConstFresh :: Constant reference density for fresh water (rain)
564  C     rho1Ref   :: reference vertical profile for density  C     rho1Ref    :: reference vertical profile for density
565  C     tRef      :: reference vertical profile for potential temperature  C     tRef       :: reference vertical profile for potential temperature
566  C     sRef      :: reference vertical profile for salinity/specific humidity  C     sRef       :: reference vertical profile for salinity/specific humidity
567  C     phiRef    :: reference potential (pressure/rho, geopotential) profile  C     phiRef     :: reference potential (pressure/rho, geopotential) profile
568  C     dBdrRef   :: vertical gradient of reference boyancy  [(m/s/r)^2)]:  C     dBdrRef    :: vertical gradient of reference buoyancy  [(m/s/r)^2]:
569  C               :: z-coord: = N^2_ref = Brunt-Vaissala frequency [s^-2]  C                :: z-coord: = N^2_ref = Brunt-Vaissala frequency [s^-2]
570  C               :: p-coord: = -(d.alpha/dp)_ref          [(m^2.s/kg)^2]  C                :: p-coord: = -(d.alpha/dp)_ref          [(m^2.s/kg)^2]
571  C     rVel2wUnit :: units conversion factor (Non-Hydrostatic code),  C     rVel2wUnit :: units conversion factor (Non-Hydrostatic code),
572  C                :: from r-coordinate vertical velocity to vertical velocity [m/s].  C                :: from r-coordinate vertical velocity to vertical velocity [m/s].
573  C                :: z-coord: = 1 ; p-coord: wSpeed [m/s] = rVel [Pa/s] * rVel2wUnit  C                :: z-coord: = 1 ; p-coord: wSpeed [m/s] = rVel [Pa/s] * rVel2wUnit
# Line 503  C     rUnit2mass :: units conversion fac Line 582  C     rUnit2mass :: units conversion fac
582  C                :: from vertical r-coordinate unit to mass per unit area [kg/m2].  C                :: from vertical r-coordinate unit to mass per unit area [kg/m2].
583  C                :: z-coord: = rhoConst  ( [m] * rho = [kg/m2] ) ;  C                :: z-coord: = rhoConst  ( [m] * rho = [kg/m2] ) ;
584  C                :: p-coord: = 1/gravity ( [Pa] /  g = [kg/m2] ) ;  C                :: p-coord: = 1/gravity ( [Pa] /  g = [kg/m2] ) ;
585  C     phiMin    :: Latitude of southern most cell face.  C     rSphere    :: Radius of sphere for a spherical polar grid ( m ).
 C     thetaMin  :: Longitude of western most cell face (this  
 C                 is an "inert" parameter but it is included  
 C                 to make geographical references simple.)  
 C     rSphere   :: Radius of sphere for a spherical polar grid ( m ).  
586  C     recip_rSphere  :: Reciprocal radius of sphere ( m ).  C     recip_rSphere  :: Reciprocal radius of sphere ( m ).
587  C     f0        :: Reference coriolis parameter ( 1/s )  C     radius_fromHorizGrid :: sphere Radius of input horiz. grid (Curvilinear Grid)
588  C                 ( Southern edge f for beta plane )  C     f0         :: Reference coriolis parameter ( 1/s )
589  C     beta      :: df/dy ( s^-1.m^-1 )  C                   ( Southern edge f for beta plane )
590  C     omega     :: Angular velocity ( rad/s )  C     beta       :: df/dy ( s^-1.m^-1 )
591    C     fPrime     :: Second Coriolis parameter ( 1/s ), related to Y-component
592    C                   of rotation (reference value = 2.Omega.Cos(Phi))
593    C     omega      :: Angular velocity ( rad/s )
594  C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)  C     rotationPeriod :: Rotation period (s) (= 2.pi/omega)
595  C     viscAr    :: Eddy viscosity coeff. for mixing of  C     viscArNr   :: vertical profile of Eddy viscosity coeff.
596  C                 momentum vertically ( units of r^2/s )  C                   for vertical mixing of momentum ( units of r^2/s )
597  C     viscAh    :: Eddy viscosity coeff. for mixing of  C     viscAh     :: Eddy viscosity coeff. for mixing of
598  C                 momentum laterally ( m^2/s )  C                   momentum laterally ( m^2/s )
599  C     viscAhW   :: Eddy viscosity coeff. for mixing of vertical  C     viscAhW    :: Eddy viscosity coeff. for mixing of vertical
600  C                 momentum laterally, no effect for hydrostatic  C                   momentum laterally, no effect for hydrostatic
601  C                 model, defaults to viscAh if unset ( m^2/s )  C                   model, defaults to viscAhD if unset ( m^2/s )
602  C                 Not used if variable horiz. viscosity is used.  C                   Not used if variable horiz. viscosity is used.
603  C     viscA4    :: Biharmonic viscosity coeff. for mixing of  C     viscA4     :: Biharmonic viscosity coeff. for mixing of
604  C                 momentum laterally ( m^4/s )  C                   momentum laterally ( m^4/s )
605  C     viscA4W   :: Biharmonic viscosity coeff. for mixing of vertical  C     viscA4W    :: Biharmonic viscosity coeff. for mixing of vertical
606  C                 momentum laterally, no effect for hydrostatic  C                   momentum laterally, no effect for hydrostatic
607  C                 model, defaults to viscA4 if unset ( m^2/s )  C                   model, defaults to viscA4D if unset ( m^2/s )
608  C                 Not used if variable horiz. viscosity is used.  C                   Not used if variable horiz. viscosity is used.
609  C     viscAhD   :: Eddy viscosity coeff. for mixing of momentum laterally  C     viscAhD    :: Eddy viscosity coeff. for mixing of momentum laterally
610  C                  (act on Divergence part) ( m^2/s )  C                   (act on Divergence part) ( m^2/s )
611  C     viscAhZ   :: Eddy viscosity coeff. for mixing of momentum laterally  C     viscAhZ    :: Eddy viscosity coeff. for mixing of momentum laterally
612  C                  (act on Vorticity  part) ( m^2/s )  C                   (act on Vorticity  part) ( m^2/s )
613  C     viscA4D   :: Biharmonic viscosity coeff. for mixing of momentum laterally  C     viscA4D    :: Biharmonic viscosity coeff. for mixing of momentum laterally
614  C                  (act on Divergence part) ( m^4/s )  C                   (act on Divergence part) ( m^4/s )
615  C     viscA4Z   :: Biharmonic viscosity coeff. for mixing of momentum laterally  C     viscA4Z    :: Biharmonic viscosity coeff. for mixing of momentum laterally
616  C                  (act on Vorticity  part) ( m^4/s )  C                   (act on Vorticity  part) ( m^4/s )
617    C     smag3D_coeff :: Isotropic 3-D Smagorinsky coefficient (-)
618  C     viscC2leith  :: Leith non-dimensional viscosity factor (grad(vort))  C     viscC2leith  :: Leith non-dimensional viscosity factor (grad(vort))
619  C     viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div))  C     viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div))
620  C     viscC4leith  :: Leith non-dimensional viscosity factor (grad(vort))  C     viscC4leith  :: Leith non-dimensional viscosity factor (grad(vort))
# Line 558  C     viscA4GridMax:: maximum and minimu Line 637  C     viscA4GridMax:: maximum and minimu
637  C     viscA4GridMin::  in terms of non-dimensional grid-size dependent viscosity  C     viscA4GridMin::  in terms of non-dimensional grid-size dependent viscosity
638  C     diffKhT   :: Laplacian diffusion coeff. for mixing of  C     diffKhT   :: Laplacian diffusion coeff. for mixing of
639  C                 heat laterally ( m^2/s )  C                 heat laterally ( m^2/s )
 C     diffKrNrT :: vertical profile of Laplacian diffusion coeff.  
 C                 for mixing of heat vertically ( units of r^2/s )  
640  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of
641  C                 heat laterally ( m^4/s )  C                 heat laterally ( m^4/s )
642    C     diffKrNrT :: vertical profile of Laplacian diffusion coeff.
643    C                 for mixing of heat vertically ( units of r^2/s )
644    C     diffKr4T  :: vertical profile of Biharmonic diffusion coeff.
645    C                 for mixing of heat vertically ( units of r^4/s )
646  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of
647  C                 salt laterally ( m^2/s )  C                 salt laterally ( m^2/s )
 C     diffKrNrS :: vertical profile of Laplacian diffusion coeff.  
 C                 for mixing of salt vertically ( units of r^2/s ),  
648  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of
649  C                 salt laterally ( m^4/s )  C                 salt laterally ( m^4/s )
650    C     diffKrNrS :: vertical profile of Laplacian diffusion coeff.
651    C                 for mixing of salt vertically ( units of r^2/s ),
652    C     diffKr4S  :: vertical profile of Biharmonic diffusion coeff.
653    C                 for mixing of salt vertically ( units of r^4/s )
654  C     diffKrBL79surf :: T/S surface diffusivity (m^2/s) Bryan and Lewis, 1979  C     diffKrBL79surf :: T/S surface diffusivity (m^2/s) Bryan and Lewis, 1979
655  C     diffKrBL79deep :: T/S deep diffusivity (m^2/s) Bryan and Lewis, 1979  C     diffKrBL79deep :: T/S deep diffusivity (m^2/s) Bryan and Lewis, 1979
656  C     diffKrBL79scl  :: depth scale for arctan fn (m) Bryan and Lewis, 1979  C     diffKrBL79scl  :: depth scale for arctan fn (m) Bryan and Lewis, 1979
# Line 587  C                    Frequency of checkp Line 670  C                    Frequency of checkp
670  C                    are referenced to this clock. ( s )  C                    are referenced to this clock. ( s )
671  C     deltaTMom    :: Timestep for momemtum equations ( s )  C     deltaTMom    :: Timestep for momemtum equations ( s )
672  C     dTtracerLev  :: Timestep for tracer equations ( s ), function of level k  C     dTtracerLev  :: Timestep for tracer equations ( s ), function of level k
673  C     deltaTfreesurf :: Timestep for free-surface equation ( s )  C     deltaTFreeSurf :: Timestep for free-surface equation ( s )
674  C     freesurfFac  :: Parameter to turn implicit free surface term on or off  C     freeSurfFac  :: Parameter to turn implicit free surface term on or off
675  C                    freesurfac = 1. uses implicit free surface  C                     freeSurFac = 1. uses implicit free surface
676  C                    freesurfac = 0. uses rigid lid  C                     freeSurFac = 0. uses rigid lid
677  C     abEps        :: Adams-Bashforth-2 stabilizing weight  C     abEps        :: Adams-Bashforth-2 stabilizing weight
678  C     alph_AB      :: Adams-Bashforth-3 primary factor  C     alph_AB      :: Adams-Bashforth-3 primary factor
679  C     beta_AB      :: Adams-Bashforth-3 secondary factor  C     beta_AB      :: Adams-Bashforth-3 secondary factor
# Line 598  C     implicSurfPress :: parameter of th Line 681  C     implicSurfPress :: parameter of th
681  C                     Implicit part of Surface Pressure Gradient ( 0-1 )  C                     Implicit part of Surface Pressure Gradient ( 0-1 )
682  C     implicDiv2Dflow :: parameter of the Crank-Nickelson time stepping :  C     implicDiv2Dflow :: parameter of the Crank-Nickelson time stepping :
683  C                     Implicit part of barotropic flow Divergence ( 0-1 )  C                     Implicit part of barotropic flow Divergence ( 0-1 )
684    C     implicitNHPress :: parameter of the Crank-Nickelson time stepping :
685    C                     Implicit part of Non-Hydrostatic Pressure Gradient ( 0-1 )
686  C     hFacMin      :: Minimum fraction size of a cell (affects hFacC etc...)  C     hFacMin      :: Minimum fraction size of a cell (affects hFacC etc...)
687  C     hFacMinDz    :: Minimum dimesional size of a cell (affects hFacC etc..., m)  C     hFacMinDz    :: Minimum dimensional size of a cell (affects hFacC etc..., m)
688  C     hFacMinDp    :: Minimum dimesional size of a cell (affects hFacC etc..., Pa)  C     hFacMinDp    :: Minimum dimensional size of a cell (affects hFacC etc..., Pa)
689  C     hFacMinDr    :: Minimum dimesional size of a cell (affects hFacC etc..., units of r)  C     hFacMinDr    :: Minimum dimensional size of a cell (-> hFacC etc..., r units)
690  C     hFacInf      :: Threshold (inf and sup) for fraction size of surface cell  C     hFacInf      :: Threshold (inf and sup) for fraction size of surface cell
691  C     hFacSup        that control vanishing and creating levels  C     hFacSup          that control vanishing and creating levels
692  C     tauCD        :: CD scheme coupling timescale ( 1/s )  C     tauCD         :: CD scheme coupling timescale ( s )
693  C     rCD          :: CD scheme normalised coupling parameter ( 0-1 )  C     rCD           :: CD scheme normalised coupling parameter (= 1 - deltaT/tauCD)
694    C     epsAB_CD      :: Adams-Bashforth-2 stabilizing weight used in CD scheme
695  C     baseTime      :: model base time (time origin) = time @ iteration zero  C     baseTime      :: model base time (time origin) = time @ iteration zero
696  C     startTime     :: Starting time for this integration ( s ).  C     startTime     :: Starting time for this integration ( s ).
697  C     endTime       :: Ending time for this integration ( s ).  C     endTime       :: Ending time for this integration ( s ).
698  C     chkPtFreq     :: Frequency of rolling check pointing ( s ).  C     chkPtFreq     :: Frequency of rolling check pointing ( s ).
699  C     pChkPtFreq    :: Frequency of permanent check pointing ( s ).  C     pChkPtFreq    :: Frequency of permanent check pointing ( s ).
700  C     dumpFreq      :: Frequency with which model state is written to  C     dumpFreq      :: Frequency with which model state is written to
701  C                     post-processing files ( s ).  C                      post-processing files ( s ).
702  C     diagFreq      :: Frequency with which model writes diagnostic output  C     diagFreq      :: Frequency with which model writes diagnostic output
703  C                     of intermediate quantities.  C                      of intermediate quantities.
704  C     afFacMom      :: Advection of momentum term tracer parameter  C     afFacMom      :: Advection of momentum term tracer parameter
705  C     vfFacMom      :: Momentum viscosity tracer parameter  C     vfFacMom      :: Momentum viscosity tracer parameter
706  C     pfFacMom      :: Momentum pressure forcing tracer parameter  C     pfFacMom      :: Momentum pressure forcing tracer parameter
# Line 641  C     convertFW2Salt :: salinity, used t Line 727  C     convertFW2Salt :: salinity, used t
727  C                       (use model surface (local) value if set to -1)  C                       (use model surface (local) value if set to -1)
728  C     temp_EvPrRn :: temperature of Rain & Evap.  C     temp_EvPrRn :: temperature of Rain & Evap.
729  C     salt_EvPrRn :: salinity of Rain & Evap.  C     salt_EvPrRn :: salinity of Rain & Evap.
730    C     temp_addMass :: temperature of addMass array
731    C     salt_addMass :: salinity of addMass array
732  C        (notes: a) tracer content of Rain/Evap only used if both  C        (notes: a) tracer content of Rain/Evap only used if both
733  C                     NonLin_FrSurf & useRealFreshWater are set.  C                     NonLin_FrSurf & useRealFreshWater are set.
734  C                b) use model surface (local) value if set to UNSET_RL)  C                b) use model surface (local) value if set to UNSET_RL)
735  C     hMixCrit    :: criteria for mixed-layer diagnostic  C     hMixCriteria:: criteria for mixed-layer diagnostic
736    C     dRhoSmall   :: parameter for mixed-layer diagnostic
737    C     hMixSmooth  :: Smoothing parameter for mixed-layer diag (default=0=no smoothing)
738  C     ivdc_kappa  :: implicit vertical diffusivity for convection [m^2/s]  C     ivdc_kappa  :: implicit vertical diffusivity for convection [m^2/s]
739  C     Ro_SeaLevel :: standard position of Sea-Level in "R" coordinate, used as  C     Ro_SeaLevel :: standard position of Sea-Level in "R" coordinate, used as
740  C                    starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)  C                    starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)
741    C     rSigmaBnd   :: vertical position (in r-unit) of r/sigma transition (Hybrid-Sigma)
742  C     sideDragFactor     :: side-drag scaling factor (used only if no_slip_sides)  C     sideDragFactor     :: side-drag scaling factor (used only if no_slip_sides)
743  C                           (default=2: full drag ; =1: gives half-slip BC)  C                           (default=2: full drag ; =1: gives half-slip BC)
744  C     bottomDragLinear    :: Linear    bottom-drag coefficient (units of [r]/s)  C     bottomDragLinear    :: Linear    bottom-drag coefficient (units of [r]/s)
# Line 658  C                           is to be rep Line 749  C                           is to be rep
749  C                           (affects myabs, mymin, mymax)  C                           (affects myabs, mymin, mymax)
750  C     nh_Am2        :: scales the non-hydrostatic terms and changes internal scales  C     nh_Am2        :: scales the non-hydrostatic terms and changes internal scales
751  C                      (i.e. allows convection at different Rayleigh numbers)  C                      (i.e. allows convection at different Rayleigh numbers)
752    C     tCylIn        :: Temperature of the cylinder inner boundary
753    C     tCylOut       :: Temperature of the cylinder outer boundary
754  C     phiEuler      :: Euler angle, rotation about original z-axis  C     phiEuler      :: Euler angle, rotation about original z-axis
755  C     thetaEuler    :: Euler angle, rotation about new x-axis  C     thetaEuler    :: Euler angle, rotation about new x-axis
756  C     psiEuler      :: Euler angle, rotation about new z-axis  C     psiEuler      :: Euler angle, rotation about new z-axis
757        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
758       & cg2dpcOffDFac, cg3dTargetResidual,       & cg2dpcOffDFac, cg3dTargetResidual,
759       & delR, delRc, delX, delY,       & delR, delRc, xgOrigin, ygOrigin,
760       & deltaT, deltaTmom, dTtracerLev, deltaTfreesurf, deltaTClock,       & deltaT, deltaTMom, dTtracerLev, deltaTFreeSurf, deltaTClock,
761       & abEps, alph_AB, beta_AB,       & abEps, alph_AB, beta_AB,
762       & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,       & rSphere, recip_rSphere, radius_fromHorizGrid,
763       & viscAh, viscAhW, viscAhMax,       & f0, beta, fPrime, omega, rotationPeriod,
764       & viscAhGrid, viscAhGridMax, viscAhGridMin,       & viscFacAdj, viscAh, viscAhW, smag3D_coeff,
765         & viscAhMax, viscAhGrid, viscAhGridMax, viscAhGridMin,
766       & viscC2leith, viscC2leithD,       & viscC2leith, viscC2leithD,
767       & viscC2smag, viscC4smag,       & viscC2smag, viscC4smag,
768       & viscAhD, viscAhZ, viscA4D, viscA4Z,       & viscAhD, viscAhZ, viscA4D, viscA4Z,
769       & viscA4, viscA4W, viscA4Max,       & viscA4, viscA4W, viscA4Max,
770       & viscA4Grid, viscA4GridMax, viscA4GridMin,       & viscA4Grid, viscA4GridMax, viscA4GridMin,
771       & viscAhReMax, viscA4ReMax,       & viscAhReMax, viscA4ReMax,
772       & viscC4leith, viscC4leithD, viscAr,       & viscC4leith, viscC4leithD, viscArNr,
773       & diffKhT, diffK4T, diffKrNrT,       & diffKhT, diffK4T, diffKrNrT, diffKr4T,
774       & diffKhS, diffK4S, diffKrNrS,       & diffKhS, diffK4S, diffKrNrS, diffKr4S,
775       & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,       & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,
776       & BL79LatVary,       & BL79LatVary,
777       & diffKrBLEQsurf, diffKrBLEQdeep, diffKrBLEQscl, diffKrBLEQHo,       & diffKrBLEQsurf, diffKrBLEQdeep, diffKrBLEQscl, diffKrBLEQHo,
778       & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,       & tauCD, rCD, epsAB_CD,
779         & freeSurfFac, implicSurfPress, implicDiv2Dflow, implicitNHPress,
780       & hFacMin, hFacMinDz, hFacInf, hFacSup,       & hFacMin, hFacMinDz, hFacInf, hFacSup,
781       & gravity, recip_gravity, gBaro,       & gravity, recip_gravity, gBaro,
782       & rhonil, recip_rhonil, rhoConst, recip_rhoConst,       & rhoNil, rhoConst, recip_rhoConst, thetaConst,
783       & rhoFacC, recip_rhoFacC, rhoFacF, recip_rhoFacF,       & rhoFacC, recip_rhoFacC, rhoFacF, recip_rhoFacF,
784       & rhoConstFresh, rho1Ref, tRef, sRef, phiRef, dBdrRef,       & rhoConstFresh, rho1Ref, tRef, sRef, phiRef, dBdrRef,
785       & rVel2wUnit, wUnit2rVel, mass2rUnit, rUnit2mass,       & rVel2wUnit, wUnit2rVel, mass2rUnit, rUnit2mass,
# Line 692  C     psiEuler      :: Euler angle, rota Line 787  C     psiEuler      :: Euler angle, rota
787       & chkPtFreq, pChkPtFreq, dumpFreq, adjDumpFreq,       & chkPtFreq, pChkPtFreq, dumpFreq, adjDumpFreq,
788       & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq,       & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq,
789       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
790       & cosPower, cAdjFreq, omega, rotationPeriod,       & cosPower, cAdjFreq,
791       & tauThetaClimRelax, tauSaltClimRelax, latBandClimRelax,       & tauThetaClimRelax, tauSaltClimRelax, latBandClimRelax,
792       & externForcingCycle, externForcingPeriod,       & externForcingCycle, externForcingPeriod,
793       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn,       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn,
794       & hFacMinDr, hFacMinDp,       & temp_addMass, salt_addMass, hFacMinDr, hFacMinDp,
795       & ivdc_kappa, hMixCriteria, Ro_SeaLevel,       & ivdc_kappa, hMixCriteria, dRhoSmall, hMixSmooth,
796         & Ro_SeaLevel, rSigmaBnd,
797       & sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2,       & sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2,
798       & smoothAbsFuncRange,       & smoothAbsFuncRange,
799       & tCylIn, tCylOut,       & tCylIn, tCylOut,
# Line 709  C     psiEuler      :: Euler angle, rota Line 805  C     psiEuler      :: Euler angle, rota
805        _RL cg2dpcOffDFac        _RL cg2dpcOffDFac
806        _RL delR(Nr)        _RL delR(Nr)
807        _RL delRc(Nr+1)        _RL delRc(Nr+1)
808        _RL delX(Nx)        _RL xgOrigin
809        _RL delY(Ny)        _RL ygOrigin
810        _RL deltaT        _RL deltaT
811        _RL deltaTClock        _RL deltaTClock
812        _RL deltaTmom        _RL deltaTMom
813        _RL dTtracerLev(Nr)        _RL dTtracerLev(Nr)
814        _RL deltaTfreesurf        _RL deltaTFreeSurf
815        _RL abEps, alph_AB, beta_AB        _RL abEps, alph_AB, beta_AB
       _RL phiMin  
       _RL thetaMin  
816        _RL rSphere        _RL rSphere
817        _RL recip_rSphere        _RL recip_rSphere
818          _RL radius_fromHorizGrid
819        _RL f0        _RL f0
820          _RL beta
821          _RL fPrime
822          _RL omega
823          _RL rotationPeriod
824        _RL freeSurfFac        _RL freeSurfFac
825        _RL implicSurfPress        _RL implicSurfPress
826        _RL implicDiv2Dflow        _RL implicDiv2Dflow
827          _RL implicitNHPress
828        _RL hFacMin        _RL hFacMin
829        _RL hFacMinDz        _RL hFacMinDz
830        _RL hFacMinDp        _RL hFacMinDp
831        _RL hFacMinDr        _RL hFacMinDr
832        _RL hFacInf        _RL hFacInf
833        _RL hFacSup        _RL hFacSup
834        _RL beta        _RL viscArNr(Nr)
835        _RL viscAr        _RL viscFacAdj
836        _RL viscAh        _RL viscAh
837        _RL viscAhW        _RL viscAhW
838        _RL viscAhD        _RL viscAhD
839        _RL viscAhZ        _RL viscAhZ
840          _RL smag3D_coeff
841        _RL viscAhMax        _RL viscAhMax
842        _RL viscAhReMax        _RL viscAhReMax
843        _RL viscAhGrid, viscAhGridMax, viscAhGridMin        _RL viscAhGrid, viscAhGridMax, viscAhGridMin
# Line 754  C     psiEuler      :: Euler angle, rota Line 855  C     psiEuler      :: Euler angle, rota
855        _RL viscC4leithD        _RL viscC4leithD
856        _RL viscC4smag        _RL viscC4smag
857        _RL diffKhT        _RL diffKhT
       _RL diffKrNrT(Nr)  
858        _RL diffK4T        _RL diffK4T
859          _RL diffKrNrT(Nr)
860          _RL diffKr4T(Nr)
861        _RL diffKhS        _RL diffKhS
       _RL diffKrNrS(Nr)  
862        _RL diffK4S        _RL diffK4S
863          _RL diffKrNrS(Nr)
864          _RL diffKr4S(Nr)
865        _RL diffKrBL79surf        _RL diffKrBL79surf
866        _RL diffKrBL79deep        _RL diffKrBL79deep
867        _RL diffKrBL79scl        _RL diffKrBL79scl
# Line 768  C     psiEuler      :: Euler angle, rota Line 871  C     psiEuler      :: Euler angle, rota
871        _RL diffKrBLEQdeep        _RL diffKrBLEQdeep
872        _RL diffKrBLEQscl        _RL diffKrBLEQscl
873        _RL diffKrBLEQHo        _RL diffKrBLEQHo
874        _RL delt        _RL tauCD, rCD, epsAB_CD
       _RL tauCD  
       _RL rCD  
875        _RL gravity        _RL gravity
876        _RL recip_gravity        _RL recip_gravity
877        _RL gBaro        _RL gBaro
878        _RL rhonil,        recip_rhonil        _RL rhoNil
879        _RL rhoConst,      recip_rhoConst        _RL rhoConst,      recip_rhoConst
880          _RL thetaConst
881        _RL rhoFacC(Nr),   recip_rhoFacC(Nr)        _RL rhoFacC(Nr),   recip_rhoFacC(Nr)
882        _RL rhoFacF(Nr+1), recip_rhoFacF(Nr+1)        _RL rhoFacF(Nr+1), recip_rhoFacF(Nr+1)
883        _RL rhoConstFresh        _RL rhoConstFresh
# Line 806  C     psiEuler      :: Euler angle, rota Line 908  C     psiEuler      :: Euler angle, rota
908        _RL mtFacMom        _RL mtFacMom
909        _RL cosPower        _RL cosPower
910        _RL cAdjFreq        _RL cAdjFreq
       _RL omega  
       _RL rotationPeriod  
911        _RL tauThetaClimRelax        _RL tauThetaClimRelax
912        _RL tauSaltClimRelax        _RL tauSaltClimRelax
913        _RL latBandClimRelax        _RL latBandClimRelax
# Line 816  C     psiEuler      :: Euler angle, rota Line 916  C     psiEuler      :: Euler angle, rota
916        _RL convertFW2Salt        _RL convertFW2Salt
917        _RL temp_EvPrRn        _RL temp_EvPrRn
918        _RL salt_EvPrRn        _RL salt_EvPrRn
919          _RL temp_addMass
920          _RL salt_addMass
921        _RL ivdc_kappa        _RL ivdc_kappa
922        _RL hMixCriteria        _RL hMixCriteria
923          _RL dRhoSmall
924          _RL hMixSmooth
925        _RL Ro_SeaLevel        _RL Ro_SeaLevel
926          _RL rSigmaBnd
927        _RL sideDragFactor        _RL sideDragFactor
928        _RL bottomDragLinear        _RL bottomDragLinear
929        _RL bottomDragQuadratic        _RL bottomDragQuadratic
930        _RL smoothAbsFuncRange        _RL smoothAbsFuncRange
931        _RL nh_Am2        _RL nh_Am2
932        _RL tCylIn        _RL tCylIn, tCylOut
       _RL tCylOut  
933        _RL phiEuler, thetaEuler, psiEuler        _RL phiEuler, thetaEuler, psiEuler
934    
935  C--   COMMON /PARM_A/ Thermodynamics constants ?  C--   COMMON /PARM_A/ Thermodynamics constants ?
936        COMMON /PARM_A/ HeatCapacity_Cp,recip_Cp        COMMON /PARM_A/ HeatCapacity_Cp
937        _RL HeatCapacity_Cp        _RL HeatCapacity_Cp
       _RL recip_Cp  
938    
939  C--   COMMON /PARM_ATM/ Atmospheric physical parameters (Ideal Gas EOS, ...)  C--   COMMON /PARM_ATM/ Atmospheric physical parameters (Ideal Gas EOS, ...)
940  C     celsius2K :: convert centigrade (Celsius) degree to Kelvin  C     celsius2K :: convert centigrade (Celsius) degree to Kelvin
# Line 860  C Logical flags for selecting packages Line 963  C Logical flags for selecting packages
963        LOGICAL useZONAL_FILT        LOGICAL useZONAL_FILT
964        LOGICAL useOPPS        LOGICAL useOPPS
965        LOGICAL usePP81        LOGICAL usePP81
966          LOGICAL useKL10
967        LOGICAL useMY82        LOGICAL useMY82
968        LOGICAL useGGL90        LOGICAL useGGL90
969        LOGICAL useKPP        LOGICAL useKPP
970        LOGICAL useGMRedi        LOGICAL useGMRedi
971        LOGICAL useDOWN_SLOPE        LOGICAL useDOWN_SLOPE
972          LOGICAL useBBL
973        LOGICAL useCAL        LOGICAL useCAL
974        LOGICAL useEXF        LOGICAL useEXF
975        LOGICAL useBulkForce        LOGICAL useBulkForce
976        LOGICAL useEBM        LOGICAL useEBM
977        LOGICAL useCheapAML        LOGICAL useCheapAML
978        LOGICAL useGrdchk        LOGICAL useGrdchk
979          LOGICAL useSMOOTH
980          LOGICAL usePROFILES
981        LOGICAL useECCO        LOGICAL useECCO
982          LOGICAL useCTRL
983        LOGICAL useSBO        LOGICAL useSBO
984        LOGICAL useFLT        LOGICAL useFLT
985        LOGICAL usePTRACERS        LOGICAL usePTRACERS
# Line 879  C Logical flags for selecting packages Line 987  C Logical flags for selecting packages
987        LOGICAL useRBCS        LOGICAL useRBCS
988        LOGICAL useOffLine        LOGICAL useOffLine
989        LOGICAL useMATRIX        LOGICAL useMATRIX
990          LOGICAL useFRAZIL
991        LOGICAL useSEAICE        LOGICAL useSEAICE
992        LOGICAL useSALT_PLUME        LOGICAL useSALT_PLUME
993        LOGICAL useShelfIce        LOGICAL useShelfIce
994          LOGICAL useStreamIce
995          LOGICAL useICEFRONT
996        LOGICAL useThSIce        LOGICAL useThSIce
997          LOGICAL useLand
998        LOGICAL useATM2d        LOGICAL useATM2d
999        LOGICAL useAIM        LOGICAL useAIM
1000        LOGICAL useLand        LOGICAL useAtm_Phys
1001        LOGICAL useFizhi        LOGICAL useFizhi
1002        LOGICAL useGridAlt        LOGICAL useGridAlt
1003        LOGICAL useDiagnostics        LOGICAL useDiagnostics
1004        LOGICAL useREGRID        LOGICAL useREGRID
1005          LOGICAL useLayers
1006        LOGICAL useMNC        LOGICAL useMNC
1007        LOGICAL useRunClock        LOGICAL useRunClock
1008        LOGICAL useEMBED_FILES        LOGICAL useEMBED_FILES
1009        LOGICAL useMYPACKAGE        LOGICAL useMYPACKAGE
1010        COMMON /PARM_PACKAGES/        COMMON /PARM_PACKAGES/
1011       &        useGAD, useOBCS, useSHAP_FILT, useZONAL_FILT,       &        useGAD, useOBCS, useSHAP_FILT, useZONAL_FILT,
1012       &        useOPPS, usePP81, useMY82, useGGL90, useKPP,       &        useOPPS, usePP81, useKL10, useMY82, useGGL90, useKPP,
1013       &        useGMRedi, useDOWN_SLOPE,       &        useGMRedi, useBBL, useDOWN_SLOPE,
1014       &        useCAL, useEXF, useBulkForce, useEBM, useCheapAML,       &        useCAL, useEXF, useBulkForce, useEBM, useCheapAML,
1015       &        useGrdchk, useECCO, useSBO, useFLT,       &        useGrdchk, useSMOOTH, usePROFILES, useECCO, useCTRL,
1016         &        useSBO, useFLT,
1017       &        usePTRACERS, useGCHEM, useRBCS, useOffLine, useMATRIX,       &        usePTRACERS, useGCHEM, useRBCS, useOffLine, useMATRIX,
1018       &        useSEAICE, useSALT_PLUME, useShelfIce, useThSIce,       &        useFRAZIL, useSEAICE, useSALT_PLUME, useShelfIce,
1019       &        useATM2D, useAIM, useLand, useFizhi, useGridAlt,       &        useStreamIce, useICEFRONT, useThSIce, useLand,
1020       &        useDiagnostics, useREGRID, useMNC,       &        useATM2D, useAIM, useAtm_Phys, useFizhi, useGridAlt,
1021         &        useDiagnostics, useREGRID, useLayers, useMNC,
1022       &        useRunClock, useEMBED_FILES,       &        useRunClock, useEMBED_FILES,
1023       &        useMYPACKAGE       &        useMYPACKAGE
1024    
1025  CEH3 ;;; Local Variables: ***  CEH3 ;;; Local Variables: ***
1026  CEH3 ;;; mode:fortran ***  CEH3 ;;; mode:fortran ***
1027  CEH3 ;;; End: ***  CEH3 ;;; End: ***
Legend:
Removed from v.1.223  
changed lines
  Added in v.1.278
  ViewVC Help
Powered by ViewVC 1.1.22