/[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.276 by heimbach, Thu Aug 7 18:43:33 2014 UTC revision 1.288 by jmc, Wed Oct 4 20:30:04 2017 UTC
# Line 33  C                        (+ 4 digits Pro Line 33  C                        (+ 4 digits Pro
33  C     adTapeDir       :: read-write checkpointing tape files from/to this  C     adTapeDir       :: read-write checkpointing tape files from/to this
34  C                        directory name instead of current dir. Conflicts  C                        directory name instead of current dir. Conflicts
35  C                        mdsioLocalDir, so only one of the two can be set.  C                        mdsioLocalDir, so only one of the two can be set.
 C                        In contrast to mdsioLocalDir, if specified adTapeDir  
 C                        must exist before the model starts.  
36  C     tRefFile      :: File containing reference Potential Temperat.  tRef (1.D)  C     tRefFile      :: File containing reference Potential Temperat.  tRef (1.D)
37  C     sRefFile      :: File containing reference salinity/spec.humid. sRef (1.D)  C     sRefFile      :: File containing reference salinity/spec.humid. sRef (1.D)
38  C     rhoRefFile    :: File containing reference density profile rhoRef (1.D)  C     rhoRefFile    :: File containing reference density profile rhoRef (1.D)
39    C     gravityFile   :: File containing gravity vertical profile (1.D)
40  C     delRFile      :: File containing vertical grid spacing delR  (1.D array)  C     delRFile      :: File containing vertical grid spacing delR  (1.D array)
41  C     delRcFile     :: File containing vertical grid spacing delRc (1.D array)  C     delRcFile     :: File containing vertical grid spacing delRc (1.D array)
42  C     hybSigmFile   :: File containing hybrid-sigma vertical coord. coeff. (2x 1.D)  C     hybSigmFile   :: File containing hybrid-sigma vertical coord. coeff. (2x 1.D)
# Line 84  C     the_run_name    :: string identify Line 83  C     the_run_name    :: string identify
83        COMMON /PARM_C/        COMMON /PARM_C/
84       &                buoyancyRelation, eosType,       &                buoyancyRelation, eosType,
85       &                pickupSuff, mdsioLocalDir, adTapeDir,       &                pickupSuff, mdsioLocalDir, adTapeDir,
86       &                tRefFile, sRefFile, rhoRefFile,       &                tRefFile, sRefFile, rhoRefFile, gravityFile,
87       &                delRFile, delRcFile, hybSigmFile,       &                delRFile, delRcFile, hybSigmFile,
88       &                delXFile, delYFile, horizGridFile,       &                delXFile, delYFile, horizGridFile,
89       &                bathyFile, topoFile, addWwallFile, addSwallFile,       &                bathyFile, topoFile, addWwallFile, addSwallFile,
# Line 108  C     the_run_name    :: string identify Line 107  C     the_run_name    :: string identify
107        CHARACTER*(MAX_LEN_FNAM) tRefFile        CHARACTER*(MAX_LEN_FNAM) tRefFile
108        CHARACTER*(MAX_LEN_FNAM) sRefFile        CHARACTER*(MAX_LEN_FNAM) sRefFile
109        CHARACTER*(MAX_LEN_FNAM) rhoRefFile        CHARACTER*(MAX_LEN_FNAM) rhoRefFile
110          CHARACTER*(MAX_LEN_FNAM) gravityFile
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        CHARACTER*(MAX_LEN_FNAM) hybSigmFile
# Line 159  C                            in con. gra Line 159  C                            in con. gra
159  C     printResidualFreq   :: Frequency for printing residual in CG iterations  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     nTimeSteps_l2       :: Number of inner timesteps to execute per timestep
 C     writeBinaryPrec     :: Precision used for writing binary files  
 C     readBinaryPrec      :: Precision used for reading binary files  
163  C     selectCoriMap       :: select setting of Coriolis parameter map:  C     selectCoriMap       :: select setting of Coriolis parameter map:
164  C                           =0 f-Plane (Constant Coriolis, = f0)  C                           =0 f-Plane (Constant Coriolis, = f0)
165  C                           =1 Beta-Plane Coriolis (= f0 + beta.y)  C                           =1 Beta-Plane Coriolis (= f0 + beta.y)
# Line 174  C     select_rStar        :: option rela Line 172  C     select_rStar        :: option rela
172  C                           =0 (default) use r coord. ; > 0 use r*  C                           =0 (default) use r coord. ; > 0 use r*
173  C     selectNHfreeSurf    :: option for Non-Hydrostatic (free-)Surface formulation:  C     selectNHfreeSurf    :: option for Non-Hydrostatic (free-)Surface formulation:
174  C                           =0 (default) hydrostatic surf. ; > 0 add NH effects.  C                           =0 (default) hydrostatic surf. ; > 0 add NH effects.
175    C     selectP_inEOS_Zc    :: select which pressure to use in EOS (for z-coords)
176    C                           =0: simply: -g*rhoConst*z
177    C                           =1: use pRef = integral{-g*rho(Tref,Sref,pRef)*dz}
178    C                           =2: use hydrostatic dynamical pressure
179    C                           =3: use full (Hyd+NH) dynamical pressure
180  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
181  C                            (3-D generalisation of oceanic real-fresh water flux)  C                            (3-D generalisation of oceanic real-fresh water flux)
182  C                           =0 off ; =1 add fluid ; =-1 virtual flux (no mass added)  C                           =0 off ; =1 add fluid ; =-1 virtual flux (no mass added)
183    C     selectImplicitDrag  :: select Implicit treatment of bottom/top drag
184    C                           = 0: fully explicit
185    C                           = 1: implicit on provisional velocity
186    C                                (i.e., before grad.Eta increment)
187    C                           = 2: fully implicit (combined with Impl Surf.Press)
188  C     momForcingOutAB     :: =1: take momentum forcing contribution  C     momForcingOutAB     :: =1: take momentum forcing contribution
189  C                            out of (=0: in) Adams-Bashforth time stepping.  C                            out of (=0: in) Adams-Bashforth time stepping.
190  C     tracForcingOutAB    :: =1: take tracer (Temp,Salt,pTracers) forcing contribution  C     tracForcingOutAB    :: =1: take tracer (Temp,Salt,pTracers) forcing contribution
# Line 187  C     saltAdvScheme       :: Salt. Horiz Line 195  C     saltAdvScheme       :: Salt. Horiz
195  C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector  C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector
196  C     selectKEscheme      :: Kinetic Energy scheme selector (Vector Inv.)  C     selectKEscheme      :: Kinetic Energy scheme selector (Vector Inv.)
197  C     selectVortScheme    :: Scheme selector for Vorticity term (Vector Inv.)  C     selectVortScheme    :: Scheme selector for Vorticity term (Vector Inv.)
198    C     selectBotDragQuadr  :: quadratic bottom drag discretisation option:
199    C                           =0: average KE from grid center to U & V location
200    C                           =1: use local velocity norm @ U & V location
201    C                           =2: same with wet-point averaging of other component
202    C     readBinaryPrec      :: Precision used for reading binary files
203    C     writeStatePrec      :: Precision used for writing model state.
204    C     writeBinaryPrec     :: Precision used for writing binary files
205    C     rwSuffixType        :: controls the format of the mds file suffix.
206    C                          =0 (default): use iteration number (myIter, I10.10);
207    C                          =1: 100*myTime (100th sec); =2: myTime (seconds);
208    C                          =3: myTime/360 (10th of hr); =4: myTime/3600 (hours).
209  C     monitorSelect       :: select group of variables to monitor  C     monitorSelect       :: select group of variables to monitor
210  C                            =1 : dynvars ; =2 : + vort ; =3 : + surface  C                            =1 : dynvars ; =2 : + vort ; =3 : + surface
211  C-    debugLevel          :: controls printing of algorithm intermediate results  C-    debugLevel          :: controls printing of algorithm intermediate results
212  C                            and statistics ; higher -> more writing  C                            and statistics ; higher -> more writing
213    C-    plotLevel           :: controls printing of field maps ; higher -> more flds
214    
215        COMMON /PARM_I/        COMMON /PARM_I/
216       &        cg2dMaxIters, cg2dChkResFreq,       &        cg2dMaxIters, cg2dChkResFreq,
217       &        cg2dPreCondFreq, cg2dUseMinResSol,       &        cg2dPreCondFreq, cg2dUseMinResSol,
218       &        cg3dMaxIters, cg3dChkResFreq,       &        cg3dMaxIters, cg3dChkResFreq,
219       &        printResidualFreq,       &        printResidualFreq,
220       &        nIter0, nTimeSteps, nEndIter,       &        nIter0, nTimeSteps, nTimeSteps_l2, nEndIter,
      &        writeStatePrec,  
      &        writeBinaryPrec, readBinaryPrec,  
221       &        selectCoriMap,       &        selectCoriMap,
222       &        selectSigmaCoord,       &        selectSigmaCoord,
223       &        nonlinFreeSurf, select_rStar,       &        nonlinFreeSurf, select_rStar,
224       &        selectNHfreeSurf,       &        selectNHfreeSurf, selectP_inEOS_Zc,
225       &        selectAddFluid,       &        selectAddFluid, selectImplicitDrag,
226       &        momForcingOutAB, tracForcingOutAB,       &        momForcingOutAB, tracForcingOutAB,
227       &        tempAdvScheme, tempVertAdvScheme,       &        tempAdvScheme, tempVertAdvScheme,
228       &        saltAdvScheme, saltVertAdvScheme,       &        saltAdvScheme, saltVertAdvScheme,
229       &        selectKEscheme, selectVortScheme,       &        selectKEscheme, selectVortScheme,
230       &        monitorSelect, debugLevel       &        selectBotDragQuadr,
231         &        readBinaryPrec, writeBinaryPrec, writeStatePrec,
232         &        rwSuffixType, monitorSelect, debugLevel, plotLevel
233        INTEGER cg2dMaxIters        INTEGER cg2dMaxIters
234        INTEGER cg2dChkResFreq        INTEGER cg2dChkResFreq
235        INTEGER cg2dPreCondFreq        INTEGER cg2dPreCondFreq
# Line 219  C                            and statist Line 239  C                            and statist
239        INTEGER printResidualFreq        INTEGER printResidualFreq
240        INTEGER nIter0        INTEGER nIter0
241        INTEGER nTimeSteps        INTEGER nTimeSteps
242          INTEGER nTimeSteps_l2
243        INTEGER nEndIter        INTEGER nEndIter
       INTEGER writeStatePrec  
       INTEGER writeBinaryPrec  
       INTEGER readBinaryPrec  
244        INTEGER selectCoriMap        INTEGER selectCoriMap
245        INTEGER selectSigmaCoord        INTEGER selectSigmaCoord
246        INTEGER nonlinFreeSurf        INTEGER nonlinFreeSurf
247        INTEGER select_rStar        INTEGER select_rStar
248        INTEGER selectNHfreeSurf        INTEGER selectNHfreeSurf
249          INTEGER selectP_inEOS_Zc
250        INTEGER selectAddFluid        INTEGER selectAddFluid
251          INTEGER selectImplicitDrag
252        INTEGER momForcingOutAB, tracForcingOutAB        INTEGER momForcingOutAB, tracForcingOutAB
253        INTEGER tempAdvScheme, tempVertAdvScheme        INTEGER tempAdvScheme, tempVertAdvScheme
254        INTEGER saltAdvScheme, saltVertAdvScheme        INTEGER saltAdvScheme, saltVertAdvScheme
255        INTEGER selectKEscheme        INTEGER selectKEscheme
256        INTEGER selectVortScheme        INTEGER selectVortScheme
257          INTEGER selectBotDragQuadr
258          INTEGER readBinaryPrec
259          INTEGER writeStatePrec
260          INTEGER writeBinaryPrec
261          INTEGER rwSuffixType
262        INTEGER monitorSelect        INTEGER monitorSelect
263        INTEGER debugLevel        INTEGER debugLevel
264          INTEGER plotLevel
265    
266  C--   COMMON /PARM_L/ Logical valued parameters used by the model.  C--   COMMON /PARM_L/ Logical valued parameters used by the model.
267  C- Coordinate + Grid params:  C- Coordinate + Grid params:
# Line 247  C     usingPCoords     :: Set to indicat Line 273  C     usingPCoords     :: Set to indicat
273  C                         type coordinate (p or p*).  C                         type coordinate (p or p*).
274  C     usingZCoords     :: Set to indicate that we are working in a height  C     usingZCoords     :: Set to indicate that we are working in a height
275  C                         type coordinate (z or z*)  C                         type coordinate (z or z*)
 C     useDynP_inEos_Zc :: use the dynamical pressure in EOS (with Z-coord.)  
 C                         this requires specific code for restart & exchange  
276  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian
277  C                           coordinate frame.  C                           coordinate frame.
278  C     usingSphericalPolarGrid :: If TRUE grid generation will be in a  C     usingSphericalPolarGrid :: If TRUE grid generation will be in a
# Line 261  C     hasWetCSCorners :: domain contains Line 285  C     hasWetCSCorners :: domain contains
285  C     deepAtmosphere :: deep model (drop the shallow-atmosphere approximation)  C     deepAtmosphere :: deep model (drop the shallow-atmosphere approximation)
286  C     setInterFDr    :: set Interface depth (put cell-Center at the middle)  C     setInterFDr    :: set Interface depth (put cell-Center at the middle)
287  C     setCenterDr    :: set cell-Center depth (put Interface at the middle)  C     setCenterDr    :: set cell-Center depth (put Interface at the middle)
288    C     useMin4hFacEdges :: set hFacW,hFacS as minimum of adjacent hFacC factor
289  C- Momentum params:  C- Momentum params:
290  C     no_slip_sides  :: Impose "no-slip" at lateral boundaries.  C     no_slip_sides  :: Impose "no-slip" at lateral boundaries.
291  C     no_slip_bottom :: Impose "no-slip" at bottom boundary.  C     no_slip_bottom :: Impose "no-slip" at bottom boundary.
292    C     bottomVisc_pCell :: account for partial-cell in bottom visc. (no-slip BC)
293  C     useSmag3D      :: Use isotropic 3-D Smagorinsky  C     useSmag3D      :: Use isotropic 3-D Smagorinsky
294  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
295  C                       on irregular grids)  C                       on irregular grids)
# Line 272  C     useAreaViscLength :: Set to true t Line 298  C     useAreaViscLength :: Set to true t
298  C                          e.g., L2=Raz.  May be preferable for cube sphere.  C                          e.g., L2=Raz.  May be preferable for cube sphere.
299  C     momViscosity  :: Flag which turns momentum friction terms on and off.  C     momViscosity  :: Flag which turns momentum friction terms on and off.
300  C     momAdvection  :: Flag which turns advection of momentum on and off.  C     momAdvection  :: Flag which turns advection of momentum on and off.
301  C     momForcing    :: Flag which turns external forcing of momentum on  C     momForcing    :: Flag which turns external forcing of momentum on and off.
302  C                      and off.  C     momTidalForcing    :: Flag which turns tidal forcing on and off.
303  C     momPressureForcing :: Flag which turns pressure term in momentum equation  C     momPressureForcing :: Flag which turns pressure term in momentum equation
304  C                          on and off.  C                          on and off.
305  C     metricTerms   :: Flag which turns metric terms on or off.  C     metricTerms   :: Flag which turns metric terms on or off.
# Line 288  C     useJamartMomAdv :: Use wet-point m Line 314  C     useJamartMomAdv :: Use wet-point m
314  C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term  C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term
315  C     highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection)  C     highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection)
316  C     useAbsVorticity :: work with f+zeta in Coriolis terms  C     useAbsVorticity :: work with f+zeta in Coriolis terms
317  C     upwindShear        :: use 1rst order upwind interp. (V.I., vertical advection)  C     upwindShear     :: use 1rst order upwind interp. (V.I., vertical advection)
318  C     momStepping    :: Turns momentum equation time-stepping off  C     momStepping    :: Turns momentum equation time-stepping off
319  C     calc_wVelocity :: Turns of vertical velocity calculation off  C     calc_wVelocity :: Turns vertical velocity calculation off
320  C- Temp. & Salt params:  C- Temp. & Salt params:
321  C     tempStepping   :: Turns temperature equation time-stepping on/off  C     tempStepping   :: Turns temperature equation time-stepping on/off
322  C     saltStepping   :: Turns salinity equation time-stepping on/off  C     saltStepping   :: Turns salinity equation time-stepping on/off
# Line 323  C                            at the surf Line 349  C                            at the surf
349  C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux  C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux
350  C                         as a real Fresh Water (=> changes the Sea Level)  C                         as a real Fresh Water (=> changes the Sea Level)
351  C                         if F, converts P+R-E to salt flux (no SL effect)  C                         if F, converts P+R-E to salt flux (no SL effect)
352    C     storePhiHyd4Phys :: store hydrostatic potential for use in Physics/EOS
353    C                         this requires specific code for restart & exchange
354  C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm  C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm
355  C     nonHydrostatic   :: Using non-hydrostatic algorithm  C     nonHydrostatic   :: Using non-hydrostatic algorithm
356  C     use3Dsolver      :: set to true to use 3-D pressure solver  C     use3Dsolver      :: set to true to use 3-D pressure solver
357  C     implicitIntGravWave :: treat Internal Gravity Wave implicitly  C     implicitIntGravWave :: treat Internal Gravity Wave implicitly
358  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
359    C     applyExchUV_early :: Apply EXCH to U,V earlier, just before integr_continuity
360  C     doResetHFactors   :: Do reset thickness factors @ beginning of each time-step  C     doResetHFactors   :: Do reset thickness factors @ beginning of each time-step
361  C     implicitDiffusion :: Turns implicit vertical diffusion on  C     implicitDiffusion :: Turns implicit vertical diffusion on
362  C     implicitViscosity :: Turns implicit vertical viscosity on  C     implicitViscosity :: Turns implicit vertical viscosity on
363  C     tempImplVertAdv :: Turns on implicit vertical advection for Temperature  C     tempImplVertAdv   :: Turns on implicit vertical advection for Temperature
364  C     saltImplVertAdv :: Turns on implicit vertical advection for Salinity  C     saltImplVertAdv   :: Turns on implicit vertical advection for Salinity
365  C     momImplVertAdv  :: Turns on implicit vertical advection for Momentum  C     momImplVertAdv    :: Turns on implicit vertical advection for Momentum
366  C     multiDimAdvection :: Flag that enable multi-dimension advection  C     multiDimAdvection :: Flag that enable multi-dimension advection
367  C     useMultiDimAdvec  :: True if multi-dim advection is used at least once  C     useMultiDimAdvec  :: True if multi-dim advection is used at least once
368  C     momDissip_In_AB   :: if False, put Dissipation tendency contribution  C     momDissip_In_AB   :: if False, put Dissipation tendency contribution
# Line 370  C     snapshot_mdsio     :: use mdsio fo Line 399  C     snapshot_mdsio     :: use mdsio fo
399  C     monitor_stdio      :: use stdio for monitor output  C     monitor_stdio      :: use stdio for monitor output
400  C     dumpInitAndLast :: dumps model state to files at Initial (nIter0)  C     dumpInitAndLast :: dumps model state to files at Initial (nIter0)
401  C                        & Last iteration, in addition multiple of dumpFreq iter.  C                        & Last iteration, in addition multiple of dumpFreq iter.
 C     printDomain     :: controls printing of domain fields (bathy, hFac ...).  
402    
403        COMMON /PARM_L/        COMMON /PARM_L/
404       & fluidIsAir, fluidIsWater,       & fluidIsAir, fluidIsWater,
405       & usingPCoords, usingZCoords, useDynP_inEos_Zc,       & usingPCoords, usingZCoords,
406       & usingCartesianGrid, usingSphericalPolarGrid, rotateGrid,       & usingCartesianGrid, usingSphericalPolarGrid, rotateGrid,
407       & usingCylindricalGrid, usingCurvilinearGrid, hasWetCSCorners,       & usingCylindricalGrid, usingCurvilinearGrid, hasWetCSCorners,
408       & deepAtmosphere, setInterFDr, setCenterDr,       & deepAtmosphere, setInterFDr, setCenterDr, useMin4hFacEdges,
409       & no_slip_sides, no_slip_bottom, useSmag3D,       & no_slip_sides, no_slip_bottom, bottomVisc_pCell, useSmag3D,
410       & useFullLeith, useStrainTensionVisc, useAreaViscLength,       & useFullLeith, useStrainTensionVisc, useAreaViscLength,
411       & momViscosity, momAdvection, momForcing,       & momViscosity, momAdvection, momForcing, momTidalForcing,
412       & momPressureForcing, metricTerms, useNHMTerms,       & momPressureForcing, metricTerms, useNHMTerms,
413       & useCoriolis, use3dCoriolis,       & useCoriolis, use3dCoriolis,
414       & useCDscheme, vectorInvariantMomentum,       & useCDscheme, vectorInvariantMomentum,
# Line 396  C     printDomain     :: controls printi Line 424  C     printDomain     :: controls printi
424       & rigidLid, implicitFreeSurface,       & rigidLid, implicitFreeSurface,
425       & uniformLin_PhiSurf, uniformFreeSurfLev,       & uniformLin_PhiSurf, uniformFreeSurfLev,
426       & exactConserv, linFSConserveTr, useRealFreshWaterFlux,       & exactConserv, linFSConserveTr, useRealFreshWaterFlux,
427       & quasiHydrostatic, nonHydrostatic, use3Dsolver,       & storePhiHyd4Phys, quasiHydrostatic, nonHydrostatic,
428       & implicitIntGravWave, staggerTimeStep, doResetHFactors,       & use3Dsolver, implicitIntGravWave, staggerTimeStep,
429         & applyExchUV_early, doResetHFactors,
430       & implicitDiffusion, implicitViscosity,       & implicitDiffusion, implicitViscosity,
431       & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,       & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,
432       & multiDimAdvection, useMultiDimAdvec,       & multiDimAdvection, useMultiDimAdvec,
# Line 412  C     printDomain     :: controls printi Line 441  C     printDomain     :: controls printi
441       & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,       & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,
442       & writePickupAtEnd,       & writePickupAtEnd,
443       & timeave_mdsio, snapshot_mdsio, monitor_stdio,       & timeave_mdsio, snapshot_mdsio, monitor_stdio,
444       & outputTypesInclusive, dumpInitAndLast,       & outputTypesInclusive, dumpInitAndLast
      & printDomain  
445    
446        LOGICAL fluidIsAir        LOGICAL fluidIsAir
447        LOGICAL fluidIsWater        LOGICAL fluidIsWater
448        LOGICAL usingPCoords        LOGICAL usingPCoords
449        LOGICAL usingZCoords        LOGICAL usingZCoords
       LOGICAL useDynP_inEos_Zc  
450        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
451        LOGICAL usingSphericalPolarGrid, rotateGrid        LOGICAL usingSphericalPolarGrid, rotateGrid
452        LOGICAL usingCylindricalGrid        LOGICAL usingCylindricalGrid
# Line 427  C     printDomain     :: controls printi Line 454  C     printDomain     :: controls printi
454        LOGICAL deepAtmosphere        LOGICAL deepAtmosphere
455        LOGICAL setInterFDr        LOGICAL setInterFDr
456        LOGICAL setCenterDr        LOGICAL setCenterDr
457          LOGICAL useMin4hFacEdges
458    
459        LOGICAL no_slip_sides        LOGICAL no_slip_sides
460        LOGICAL no_slip_bottom        LOGICAL no_slip_bottom
461          LOGICAL bottomVisc_pCell
462        LOGICAL useSmag3D        LOGICAL useSmag3D
463        LOGICAL useFullLeith        LOGICAL useFullLeith
464        LOGICAL useStrainTensionVisc        LOGICAL useStrainTensionVisc
# Line 437  C     printDomain     :: controls printi Line 466  C     printDomain     :: controls printi
466        LOGICAL momViscosity        LOGICAL momViscosity
467        LOGICAL momAdvection        LOGICAL momAdvection
468        LOGICAL momForcing        LOGICAL momForcing
469          LOGICAL momTidalForcing
470        LOGICAL momPressureForcing        LOGICAL momPressureForcing
471        LOGICAL metricTerms        LOGICAL metricTerms
472        LOGICAL useNHMTerms        LOGICAL useNHMTerms
# Line 477  C     printDomain     :: controls printi Line 507  C     printDomain     :: controls printi
507        LOGICAL exactConserv        LOGICAL exactConserv
508        LOGICAL linFSConserveTr        LOGICAL linFSConserveTr
509        LOGICAL useRealFreshWaterFlux        LOGICAL useRealFreshWaterFlux
510          LOGICAL storePhiHyd4Phys
511        LOGICAL quasiHydrostatic        LOGICAL quasiHydrostatic
512        LOGICAL nonHydrostatic        LOGICAL nonHydrostatic
513        LOGICAL use3Dsolver        LOGICAL use3Dsolver
514        LOGICAL implicitIntGravWave        LOGICAL implicitIntGravWave
515        LOGICAL staggerTimeStep        LOGICAL staggerTimeStep
516          LOGICAL applyExchUV_early
517        LOGICAL doResetHFactors        LOGICAL doResetHFactors
518        LOGICAL implicitDiffusion        LOGICAL implicitDiffusion
519        LOGICAL implicitViscosity        LOGICAL implicitViscosity
# Line 510  C     printDomain     :: controls printi Line 542  C     printDomain     :: controls printi
542        LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio        LOGICAL timeave_mdsio, snapshot_mdsio, monitor_stdio
543        LOGICAL outputTypesInclusive        LOGICAL outputTypesInclusive
544        LOGICAL dumpInitAndLast        LOGICAL dumpInitAndLast
       LOGICAL printDomain  
545    
546  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
547  C     cg2dTargetResidual  C     cg2dTargetResidual
# Line 542  C                :: most cell face (Lat- Line 573  C                :: most cell face (Lat-
573  C                :: parameter but it makes geographical references simple.)  C                :: parameter but it makes geographical references simple.)
574  C     ygOrigin   :: Origin of the Y-axis (Cartesian Grid) / Latitude of Southern  C     ygOrigin   :: Origin of the Y-axis (Cartesian Grid) / Latitude of Southern
575  C                :: most face (Lat-Lon grid).  C                :: most face (Lat-Lon grid).
576  C     gravity    :: Accel. due to gravity ( m/s^2 )  C     rSphere    :: Radius of sphere for a spherical polar grid ( m ).
577  C     recip_gravity and its inverse  C     recip_rSphere :: Reciprocal radius of sphere ( m^-1 ).
578    C     radius_fromHorizGrid :: sphere Radius of input horiz. grid (Curvilinear Grid)
579    C     seaLev_Z   :: the reference height of sea-level (usually zero)
580    C     top_Pres   :: pressure (P-Coords) or reference pressure (Z-Coords) at the top
581    C     rSigmaBnd  :: vertical position (in r-unit) of r/sigma transition (Hybrid-Sigma)
582    C     gravity    :: Acceleration due to constant gravity ( m/s^2 )
583    C     recip_gravity :: Reciprocal gravity acceleration ( s^2/m )
584  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 )
585    C     gravFacC   :: gravity factor (vs surf. gravity) vert. profile at cell-Center
586    C     gravFacF   :: gravity factor (vs surf. gravity) vert. profile at cell-interF
587  C     rhoNil     :: Reference density for the linear equation of state  C     rhoNil     :: Reference density for the linear equation of state
588  C     rhoConst   :: Vertically constant reference density (Boussinesq)  C     rhoConst   :: Vertically constant reference density (Boussinesq)
589  C     thetaConst :: Constant reference for potential temperature  C     rho1Ref    :: reference vertical profile for density (anelastic)
590  C     rhoFacC    :: normalized (by rhoConst) reference density at cell-Center  C     rhoFacC    :: normalized (by rhoConst) reference density at cell-Center
591  C     rhoFacF    :: normalized (by rhoConst) reference density at cell-interFace  C     rhoFacF    :: normalized (by rhoConst) reference density at cell-interFace
592  C     rhoConstFresh :: Constant reference density for fresh water (rain)  C     rhoConstFresh :: Constant reference density for fresh water (rain)
593  C     rho1Ref    :: reference vertical profile for density  C     thetaConst :: Constant reference for potential temperature
594  C     tRef       :: reference vertical profile for potential temperature  C     tRef       :: reference vertical profile for potential temperature
595  C     sRef       :: reference vertical profile for salinity/specific humidity  C     sRef       :: reference vertical profile for salinity/specific humidity
596  C     phiRef     :: reference potential (pressure/rho, geopotential) profile  C     pRef4EOS   :: reference pressure used in EOS (case selectP_inEOS_Zc=1)
597    C     phiRef     :: reference potential (press/rho, geopot) profile (m^2/s^2)
598  C     dBdrRef    :: vertical gradient of reference buoyancy  [(m/s/r)^2]:  C     dBdrRef    :: vertical gradient of reference buoyancy  [(m/s/r)^2]:
599  C                :: z-coord: = N^2_ref = Brunt-Vaissala frequency [s^-2]  C                :: z-coord: = N^2_ref = Brunt-Vaissala frequency [s^-2]
600  C                :: p-coord: = -(d.alpha/dp)_ref          [(m^2.s/kg)^2]  C                :: p-coord: = -(d.alpha/dp)_ref          [(m^2.s/kg)^2]
# Line 572  C     rUnit2mass :: units conversion fac Line 612  C     rUnit2mass :: units conversion fac
612  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].
613  C                :: z-coord: = rhoConst  ( [m] * rho = [kg/m2] ) ;  C                :: z-coord: = rhoConst  ( [m] * rho = [kg/m2] ) ;
614  C                :: p-coord: = 1/gravity ( [Pa] /  g = [kg/m2] ) ;  C                :: p-coord: = 1/gravity ( [Pa] /  g = [kg/m2] ) ;
 C     rSphere    :: Radius of sphere for a spherical polar grid ( m ).  
 C     recip_rSphere  :: Reciprocal radius of sphere ( m ).  
 C     radius_fromHorizGrid :: sphere Radius of input horiz. grid (Curvilinear Grid)  
615  C     f0         :: Reference coriolis parameter ( 1/s )  C     f0         :: Reference coriolis parameter ( 1/s )
616  C                   ( Southern edge f for beta plane )  C                   ( Southern edge f for beta plane )
617  C     beta       :: df/dy ( s^-1.m^-1 )  C     beta       :: df/dy ( s^-1.m^-1 )
# Line 669  C     alph_AB      :: Adams-Bashforth-3 Line 706  C     alph_AB      :: Adams-Bashforth-3
706  C     beta_AB      :: Adams-Bashforth-3 secondary factor  C     beta_AB      :: Adams-Bashforth-3 secondary factor
707  C     implicSurfPress :: parameter of the Crank-Nickelson time stepping :  C     implicSurfPress :: parameter of the Crank-Nickelson time stepping :
708  C                     Implicit part of Surface Pressure Gradient ( 0-1 )  C                     Implicit part of Surface Pressure Gradient ( 0-1 )
709  C     implicDiv2Dflow :: parameter of the Crank-Nickelson time stepping :  C     implicDiv2DFlow :: parameter of the Crank-Nickelson time stepping :
710  C                     Implicit part of barotropic flow Divergence ( 0-1 )  C                     Implicit part of barotropic flow Divergence ( 0-1 )
711  C     implicitNHPress :: parameter of the Crank-Nickelson time stepping :  C     implicitNHPress :: parameter of the Crank-Nickelson time stepping :
712  C                     Implicit part of Non-Hydrostatic Pressure Gradient ( 0-1 )  C                     Implicit part of Non-Hydrostatic Pressure Gradient ( 0-1 )
# Line 726  C     hMixCriteria:: criteria for mixed- Line 763  C     hMixCriteria:: criteria for mixed-
763  C     dRhoSmall   :: parameter for mixed-layer diagnostic  C     dRhoSmall   :: parameter for mixed-layer diagnostic
764  C     hMixSmooth  :: Smoothing parameter for mixed-layer diag (default=0=no smoothing)  C     hMixSmooth  :: Smoothing parameter for mixed-layer diag (default=0=no smoothing)
765  C     ivdc_kappa  :: implicit vertical diffusivity for convection [m^2/s]  C     ivdc_kappa  :: implicit vertical diffusivity for convection [m^2/s]
 C     Ro_SeaLevel :: standard position of Sea-Level in "R" coordinate, used as  
 C                    starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)  
 C     rSigmaBnd   :: vertical position (in r-unit) of r/sigma transition (Hybrid-Sigma)  
766  C     sideDragFactor     :: side-drag scaling factor (used only if no_slip_sides)  C     sideDragFactor     :: side-drag scaling factor (used only if no_slip_sides)
767  C                           (default=2: full drag ; =1: gives half-slip BC)  C                           (default=2: full drag ; =1: gives half-slip BC)
768  C     bottomDragLinear    :: Linear    bottom-drag coefficient (units of [r]/s)  C     bottomDragLinear    :: Linear    bottom-drag coefficient (units of [r]/s)
# Line 746  C     thetaEuler    :: Euler angle, rota Line 780  C     thetaEuler    :: Euler angle, rota
780  C     psiEuler      :: Euler angle, rotation about new z-axis  C     psiEuler      :: Euler angle, rotation about new z-axis
781        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
782       & cg2dpcOffDFac, cg3dTargetResidual,       & cg2dpcOffDFac, cg3dTargetResidual,
783       & delR, delRc, xgOrigin, ygOrigin,       & delR, delRc, xgOrigin, ygOrigin, rSphere, recip_rSphere,
784         & radius_fromHorizGrid, seaLev_Z, top_Pres, rSigmaBnd,
785       & deltaT, deltaTMom, dTtracerLev, deltaTFreeSurf, deltaTClock,       & deltaT, deltaTMom, dTtracerLev, deltaTFreeSurf, deltaTClock,
786       & abEps, alph_AB, beta_AB,       & abEps, alph_AB, beta_AB,
      & rSphere, recip_rSphere, radius_fromHorizGrid,  
787       & f0, beta, fPrime, omega, rotationPeriod,       & f0, beta, fPrime, omega, rotationPeriod,
788       & viscFacAdj, viscAh, viscAhW, smag3D_coeff,       & viscFacAdj, viscAh, viscAhW, smag3D_coeff,
789       & viscAhMax, viscAhGrid, viscAhGridMax, viscAhGridMin,       & viscAhMax, viscAhGrid, viscAhGridMax, viscAhGridMin,
# Line 766  C     psiEuler      :: Euler angle, rota Line 800  C     psiEuler      :: Euler angle, rota
800       & BL79LatVary,       & BL79LatVary,
801       & diffKrBLEQsurf, diffKrBLEQdeep, diffKrBLEQscl, diffKrBLEQHo,       & diffKrBLEQsurf, diffKrBLEQdeep, diffKrBLEQscl, diffKrBLEQHo,
802       & tauCD, rCD, epsAB_CD,       & tauCD, rCD, epsAB_CD,
803       & freeSurfFac, implicSurfPress, implicDiv2Dflow, implicitNHPress,       & freeSurfFac, implicSurfPress, implicDiv2DFlow, implicitNHPress,
804       & hFacMin, hFacMinDz, hFacInf, hFacSup,       & hFacMin, hFacMinDz, hFacInf, hFacSup,
805       & gravity, recip_gravity, gBaro,       & gravity, recip_gravity, gBaro,
806       & rhoNil, rhoConst, recip_rhoConst, thetaConst,       & gravFacC, recip_gravFacC, gravFacF, recip_gravFacF,
807       & rhoFacC, recip_rhoFacC, rhoFacF, recip_rhoFacF,       & rhoNil, rhoConst, recip_rhoConst, rho1Ref,
808       & rhoConstFresh, rho1Ref, tRef, sRef, phiRef, dBdrRef,       & rhoFacC, recip_rhoFacC, rhoFacF, recip_rhoFacF, rhoConstFresh,
809         & thetaConst, tRef, sRef, pRef4EOS, phiRef, dBdrRef,
810       & rVel2wUnit, wUnit2rVel, mass2rUnit, rUnit2mass,       & rVel2wUnit, wUnit2rVel, mass2rUnit, rUnit2mass,
811       & baseTime, startTime, endTime,       & baseTime, startTime, endTime,
812       & chkPtFreq, pChkPtFreq, dumpFreq, adjDumpFreq,       & chkPtFreq, pChkPtFreq, dumpFreq, adjDumpFreq,
# Line 783  C     psiEuler      :: Euler angle, rota Line 818  C     psiEuler      :: Euler angle, rota
818       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn,       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn,
819       & temp_addMass, salt_addMass, hFacMinDr, hFacMinDp,       & temp_addMass, salt_addMass, hFacMinDr, hFacMinDp,
820       & ivdc_kappa, hMixCriteria, dRhoSmall, hMixSmooth,       & ivdc_kappa, hMixCriteria, dRhoSmall, hMixSmooth,
      & Ro_SeaLevel, rSigmaBnd,  
821       & sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2,       & sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2,
822       & smoothAbsFuncRange,       & smoothAbsFuncRange,
823       & tCylIn, tCylOut,       & tCylIn, tCylOut,
# Line 797  C     psiEuler      :: Euler angle, rota Line 831  C     psiEuler      :: Euler angle, rota
831        _RL delRc(Nr+1)        _RL delRc(Nr+1)
832        _RL xgOrigin        _RL xgOrigin
833        _RL ygOrigin        _RL ygOrigin
834          _RL rSphere
835          _RL recip_rSphere
836          _RL radius_fromHorizGrid
837          _RL seaLev_Z
838          _RL top_Pres
839          _RL rSigmaBnd
840        _RL deltaT        _RL deltaT
841        _RL deltaTClock        _RL deltaTClock
842        _RL deltaTMom        _RL deltaTMom
843        _RL dTtracerLev(Nr)        _RL dTtracerLev(Nr)
844        _RL deltaTFreeSurf        _RL deltaTFreeSurf
845        _RL abEps, alph_AB, beta_AB        _RL abEps, alph_AB, beta_AB
       _RL rSphere  
       _RL recip_rSphere  
       _RL radius_fromHorizGrid  
846        _RL f0        _RL f0
847        _RL beta        _RL beta
848        _RL fPrime        _RL fPrime
# Line 813  C     psiEuler      :: Euler angle, rota Line 850  C     psiEuler      :: Euler angle, rota
850        _RL rotationPeriod        _RL rotationPeriod
851        _RL freeSurfFac        _RL freeSurfFac
852        _RL implicSurfPress        _RL implicSurfPress
853        _RL implicDiv2Dflow        _RL implicDiv2DFlow
854        _RL implicitNHPress        _RL implicitNHPress
855        _RL hFacMin        _RL hFacMin
856        _RL hFacMinDz        _RL hFacMinDz
# Line 862  C     psiEuler      :: Euler angle, rota Line 899  C     psiEuler      :: Euler angle, rota
899        _RL diffKrBLEQscl        _RL diffKrBLEQscl
900        _RL diffKrBLEQHo        _RL diffKrBLEQHo
901        _RL tauCD, rCD, epsAB_CD        _RL tauCD, rCD, epsAB_CD
902        _RL gravity        _RL gravity,       recip_gravity
       _RL recip_gravity  
903        _RL gBaro        _RL gBaro
904          _RL gravFacC(Nr),   recip_gravFacC(Nr)
905          _RL gravFacF(Nr+1), recip_gravFacF(Nr+1)
906        _RL rhoNil        _RL rhoNil
907        _RL rhoConst,      recip_rhoConst        _RL rhoConst,      recip_rhoConst
908        _RL thetaConst        _RL rho1Ref(Nr)
909        _RL rhoFacC(Nr),   recip_rhoFacC(Nr)        _RL rhoFacC(Nr),   recip_rhoFacC(Nr)
910        _RL rhoFacF(Nr+1), recip_rhoFacF(Nr+1)        _RL rhoFacF(Nr+1), recip_rhoFacF(Nr+1)
911        _RL rhoConstFresh        _RL rhoConstFresh
912        _RL rho1Ref(Nr)        _RL thetaConst
913        _RL tRef(Nr)        _RL tRef(Nr)
914        _RL sRef(Nr)        _RL sRef(Nr)
915          _RL pRef4EOS(Nr)
916        _RL phiRef(2*Nr+1)        _RL phiRef(2*Nr+1)
917        _RL dBdrRef(Nr)        _RL dBdrRef(Nr)
918        _RL rVel2wUnit(Nr+1), wUnit2rVel(Nr+1)        _RL rVel2wUnit(Nr+1), wUnit2rVel(Nr+1)
# Line 912  C     psiEuler      :: Euler angle, rota Line 951  C     psiEuler      :: Euler angle, rota
951        _RL hMixCriteria        _RL hMixCriteria
952        _RL dRhoSmall        _RL dRhoSmall
953        _RL hMixSmooth        _RL hMixSmooth
       _RL Ro_SeaLevel  
       _RL rSigmaBnd  
954        _RL sideDragFactor        _RL sideDragFactor
955        _RL bottomDragLinear        _RL bottomDragLinear
956        _RL bottomDragQuadratic        _RL bottomDragQuadratic
# Line 965  C Logical flags for selecting packages Line 1002  C Logical flags for selecting packages
1002        LOGICAL useBulkForce        LOGICAL useBulkForce
1003        LOGICAL useEBM        LOGICAL useEBM
1004        LOGICAL useCheapAML        LOGICAL useCheapAML
1005          LOGICAL useAUTODIFF
1006        LOGICAL useGrdchk        LOGICAL useGrdchk
1007        LOGICAL useSMOOTH        LOGICAL useSMOOTH
1008        LOGICAL usePROFILES        LOGICAL usePROFILES
1009        LOGICAL useECCO        LOGICAL useECCO
1010          LOGICAL useCTRL
1011        LOGICAL useSBO        LOGICAL useSBO
1012        LOGICAL useFLT        LOGICAL useFLT
1013        LOGICAL usePTRACERS        LOGICAL usePTRACERS
# Line 1001  C Logical flags for selecting packages Line 1040  C Logical flags for selecting packages
1040       &        useOPPS, usePP81, useKL10, useMY82, useGGL90, useKPP,       &        useOPPS, usePP81, useKL10, useMY82, useGGL90, useKPP,
1041       &        useGMRedi, useBBL, useDOWN_SLOPE,       &        useGMRedi, useBBL, useDOWN_SLOPE,
1042       &        useCAL, useEXF, useBulkForce, useEBM, useCheapAML,       &        useCAL, useEXF, useBulkForce, useEBM, useCheapAML,
1043       &        useGrdchk,useSMOOTH,usePROFILES,useECCO,useSBO, useFLT,       &        useGrdchk, useSMOOTH, usePROFILES, useECCO, useCTRL,
1044         &        useSBO, useFLT, useAUTODIFF,
1045       &        usePTRACERS, useGCHEM, useRBCS, useOffLine, useMATRIX,       &        usePTRACERS, useGCHEM, useRBCS, useOffLine, useMATRIX,
1046       &        useFRAZIL, useSEAICE, useSALT_PLUME, useShelfIce,       &        useFRAZIL, useSEAICE, useSALT_PLUME, useShelfIce,
1047       &        useStreamIce, useICEFRONT, useThSIce, useLand,       &        useStreamIce, useICEFRONT, useThSIce, useLand,
Legend:
Removed from v.1.276  
changed lines
  Added in v.1.288
  ViewVC Help
Powered by ViewVC 1.1.22