/[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.257 by jmc, Sat Jan 21 18:46:27 2012 UTC revision 1.266 by jmc, Fri Nov 9 22:29:32 2012 UTC
# Line 64  C     viscA4Zfile     :: File containing Line 64  C     viscA4Zfile     :: File containing
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
# Line 78  C                 m/s when using externa Line 78  C                 m/s when using externa
78  C                 to kg/m2/s by multiplying by rhoConstFresh.  C                 to kg/m2/s by multiplying by rhoConstFresh.
79  C     saltFluxFile    :: File containing surface salt flux  C     saltFluxFile    :: File containing surface salt flux
80  C     pLoadFile       :: File containing pressure loading  C     pLoadFile       :: File containing pressure loading
81    C     addMassFile     :: File containing source/sink of fluid in the interior
82  C     eddyPsiXFile    :: File containing zonal Eddy streamfunction data  C     eddyPsiXFile    :: File containing zonal Eddy streamfunction data
83  C     eddyPsiYFile    :: File containing meridional Eddy streamfunction data  C     eddyPsiYFile    :: File containing meridional Eddy streamfunction data
84  C     the_run_name    :: string identifying the name of the model "run"  C     the_run_name    :: string identifying the name of the model "run"
# Line 97  C     the_run_name    :: string identify Line 98  C     the_run_name    :: string identify
98       &                surfQfile, surfQnetFile, surfQswFile,       &                surfQfile, surfQnetFile, surfQswFile,
99       &                lambdaThetaFile, lambdaSaltFile,       &                lambdaThetaFile, lambdaSaltFile,
100       &                uVelInitFile, vVelInitFile, pSurfInitFile,       &                uVelInitFile, vVelInitFile, pSurfInitFile,
101       &                dQdTfile, ploadFile,       &                dQdTfile, pLoadFile, addMassFile,
102       &                eddyPsiXFile, eddyPsiYFile,       &                eddyPsiXFile, eddyPsiYFile,
103       &                the_run_name       &                the_run_name
104        CHARACTER*(MAX_LEN_FNAM) buoyancyRelation        CHARACTER*(MAX_LEN_FNAM) buoyancyRelation
# Line 134  C     the_run_name    :: string identify Line 135  C     the_run_name    :: string identify
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) dQdTfile
138        CHARACTER*(MAX_LEN_FNAM) ploadFile        CHARACTER*(MAX_LEN_FNAM) pLoadFile
139          CHARACTER*(MAX_LEN_FNAM) addMassFile
140        CHARACTER*(MAX_LEN_FNAM) eddyPsiXFile        CHARACTER*(MAX_LEN_FNAM) eddyPsiXFile
141        CHARACTER*(MAX_LEN_FNAM) eddyPsiYFile        CHARACTER*(MAX_LEN_FNAM) eddyPsiYFile
142        CHARACTER*(MAX_LEN_FNAM) lambdaThetaFile        CHARACTER*(MAX_LEN_FNAM) lambdaThetaFile
# Line 143  C     the_run_name    :: string identify Line 145  C     the_run_name    :: string identify
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     cg3dMaxIters        :: Maximum number of iterations in the  C     cg3dMaxIters        :: Maximum number of iterations in the
154  C                           three-dimensional con. grad solver.  C                            three-dimensional con. grad solver.
155  C     cg3dChkResFreq      :: Frequency with which to check residual  C     cg3dChkResFreq      :: Frequency with which to check residual
156  C                           in con. grad solver.  C                            in con. grad solver.
157  C     printResidualFreq   :: Frequency for printing residual in CG iterations  C     printResidualFreq   :: Frequency for printing residual in CG iterations
158  C     nIter0              :: Start time-step number of for this run  C     nIter0              :: Start time-step number of for this run
159  C     nTimeSteps          :: Number of timesteps to execute  C     nTimeSteps          :: Number of timesteps to execute
# Line 174  C     selectAddFluid      :: option to a Line 176  C     selectAddFluid      :: option to a
176  C                            (3-D generalisation of oceanic real-fresh water flux)  C                            (3-D generalisation of oceanic real-fresh water flux)
177  C                           =0 off ; =1 add fluid ; =-1 virtual flux (no mass added)  C                           =0 off ; =1 add fluid ; =-1 virtual flux (no mass added)
178  C     momForcingOutAB     :: =1: take momentum forcing contribution  C     momForcingOutAB     :: =1: take momentum forcing contribution
179  C                           out of (=0: in) Adams-Bashforth time stepping.  C                            out of (=0: in) Adams-Bashforth time stepping.
180  C     tracForcingOutAB    :: =1: take tracer (Temp,Salt,pTracers) forcing contribution  C     tracForcingOutAB    :: =1: take tracer (Temp,Salt,pTracers) forcing contribution
181  C                           out of (=0: in) Adams-Bashforth time stepping.  C                            out of (=0: in) Adams-Bashforth time stepping.
182  C     tempAdvScheme       :: Temp. Horiz.Advection scheme selector  C     tempAdvScheme       :: Temp. Horiz.Advection scheme selector
183  C     tempVertAdvScheme   :: Temp. Vert. Advection scheme selector  C     tempVertAdvScheme   :: Temp. Vert. Advection scheme selector
184  C     saltAdvScheme       :: Salt. Horiz.advection scheme selector  C     saltAdvScheme       :: Salt. Horiz.advection scheme selector
# Line 234  C                            and statist Line 236  C                            and statist
236  C--   COMMON /PARM_L/ Logical valued parameters used by the model.  C--   COMMON /PARM_L/ Logical valued parameters used by the model.
237  C- Coordinate + Grid params:  C- Coordinate + Grid params:
238  C     fluidIsAir       :: Set to indicate that the fluid major constituent  C     fluidIsAir       :: Set to indicate that the fluid major constituent
239  C                        is air  C                         is air
240  C     fluidIsWater     :: Set to indicate that the fluid major constituent  C     fluidIsWater     :: Set to indicate that the fluid major constituent
241  C                        is water  C                         is water
242  C     usingPCoords     :: Set to indicate that we are working in a pressure  C     usingPCoords     :: Set to indicate that we are working in a pressure
243  C                        type coordinate (p or p*).  C                         type coordinate (p or p*).
244  C     usingZCoords     :: Set to indicate that we are working in a height  C     usingZCoords     :: Set to indicate that we are working in a height
245  C                        type coordinate (z or z*)  C                         type coordinate (z or z*)
246  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.)
247  C                         this requires specific code for restart & exchange  C                         this requires specific code for restart & exchange
248  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian  C     usingCartesianGrid :: If TRUE grid generation will be in a cartesian
249  C                          coordinate frame.  C                           coordinate frame.
250  C     usingSphericalPolarGrid :: If TRUE grid generation will be in a  C     usingSphericalPolarGrid :: If TRUE grid generation will be in a
251  C                               spherical polar frame.  C                                spherical polar frame.
252  C     rotateGrid      :: rotate grid coordinates to geographical coordinates  C     rotateGrid      :: rotate grid coordinates to geographical coordinates
253  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)  
254  C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical  C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical
255    C     usingCurvilinearGrid :: If TRUE, use a curvilinear grid (to be provided)
256    C     hasWetCSCorners :: domain contains CS-type corners where dynamics is solved
257  C     deepAtmosphere :: deep model (drop the shallow-atmosphere approximation)  C     deepAtmosphere :: deep model (drop the shallow-atmosphere approximation)
258  C     setInterFDr    :: set Interface depth (put cell-Center at the middle)  C     setInterFDr    :: set Interface depth (put cell-Center at the middle)
259  C     setCenterDr    :: set cell-Center depth (put Interface at the middle)  C     setCenterDr    :: set cell-Center depth (put Interface at the middle)
# Line 284  C     upwindShear        :: use 1rst ord Line 287  C     upwindShear        :: use 1rst ord
287  C     momStepping    :: Turns momentum equation time-stepping off  C     momStepping    :: Turns momentum equation time-stepping off
288  C     calc_wVelocity :: Turns of vertical velocity calculation off  C     calc_wVelocity :: Turns of vertical velocity calculation off
289  C- Temp. & Salt params:  C- Temp. & Salt params:
290  C     tempStepping   :: Turns temperature equation time-stepping off  C     tempStepping   :: Turns temperature equation time-stepping on/off
291  C     saltStepping   :: Turns salinity equation time-stepping off  C     saltStepping   :: Turns salinity equation time-stepping on/off
292    C     addFrictionHeating :: account for frictional heating
293  C     tempAdvection  :: Flag which turns advection of temperature on and off.  C     tempAdvection  :: Flag which turns advection of temperature on and off.
294    C     tempVertDiff4  :: use vertical bi-harmonic diffusion for temperature
295  C     tempIsActiveTr :: Pot.Temp. is a dynamically active tracer  C     tempIsActiveTr :: Pot.Temp. is a dynamically active tracer
296  C     tempForcing    :: Flag which turns external forcing of temperature on  C     tempForcing    :: Flag which turns external forcing of temperature on/off
 C                       and off.  
297  C     saltAdvection  :: Flag which turns advection of salinity on and off.  C     saltAdvection  :: Flag which turns advection of salinity on and off.
298    C     saltVertDiff4  :: use vertical bi-harmonic diffusion for salinity
299  C     saltIsActiveTr :: Salinity  is a dynamically active tracer  C     saltIsActiveTr :: Salinity  is a dynamically active tracer
300  C     saltForcing    :: Flag which turns external forcing of salinity on  C     saltForcing    :: Flag which turns external forcing of salinity on/off
 C                       and off.  
301  C     maskIniTemp    :: apply mask to initial Pot.Temp.  C     maskIniTemp    :: apply mask to initial Pot.Temp.
302  C     maskIniSalt    :: apply mask to initial salinity  C     maskIniSalt    :: apply mask to initial salinity
303  C     checkIniTemp   :: check for points with identically zero initial Pot.Temp.  C     checkIniTemp   :: check for points with identically zero initial Pot.Temp.
304  C     checkIniSalt   :: check for points with identically zero initial salinity  C     checkIniSalt   :: check for points with identically zero initial salinity
 C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux  
 C                         as a real Fresh Water (=> changes the Sea Level)  
 C                         if F, converts P+R-E to salt flux (no SL effect)  
305  C- Pressure solver related parameters (PARM02)  C- Pressure solver related parameters (PARM02)
306  C     useSRCGSolver  :: Set to true to use conjugate gradient  C     useSRCGSolver  :: Set to true to use conjugate gradient
307  C                       solver with single reduction (only one call of  C                       solver with single reduction (only one call of
308  C                       s/r mpi_allreduce), default is false  C                       s/r mpi_allreduce), default is false
309  C- Time-stepping params:  C- Time-stepping & free-surface params:
310  C     rigidLid            :: Set to true to use rigid lid  C     rigidLid            :: Set to true to use rigid lid
311  C     implicitFreeSurface :: Set to true to use implicit free surface  C     implicitFreeSurface :: Set to true to use implicit free surface
312    C     uniformLin_PhiSurf  :: Set to true to use a uniform Bo_surf in the
313    C                            linear relation Phi_surf = Bo_surf*eta
314    C     uniformFreeSurfLev  :: TRUE if free-surface level-index is uniform (=1)
315  C     exactConserv        :: Set to true to conserve exactly the total Volume  C     exactConserv        :: Set to true to conserve exactly the total Volume
316  C     linFSConserveTr     :: Set to true to correct source/sink of tracer  C     linFSConserveTr     :: Set to true to correct source/sink of tracer
317  C                            at the surface due to Linear Free Surface  C                            at the surface due to Linear Free Surface
318  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
319  C                            linear relation Phi_surf = Bo_surf*eta  C                         as a real Fresh Water (=> changes the Sea Level)
320    C                         if F, converts P+R-E to salt flux (no SL effect)
321  C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm  C     quasiHydrostatic :: Using non-hydrostatic terms in hydrostatic algorithm
322  C     nonHydrostatic   :: Using non-hydrostatic algorithm  C     nonHydrostatic   :: Using non-hydrostatic algorithm
323  C     use3Dsolver      :: set to true to use 3-D pressure solver  C     use3Dsolver      :: set to true to use 3-D pressure solver
# Line 338  C     doThetaClimRelax :: Set true if re Line 343  C     doThetaClimRelax :: Set true if re
343  C                        climatology is required.  C                        climatology is required.
344  C     doSaltClimRelax  :: Set true if relaxation to salinity  C     doSaltClimRelax  :: Set true if relaxation to salinity
345  C                        climatology is required.  C                        climatology is required.
346    C     balanceThetaClimRelax :: substract global mean effect at every time step
347    C     balanceSaltClimRelax :: substract global mean effect at every time step
348  C     allowFreezing  :: Allows surface water to freeze and form ice  C     allowFreezing  :: Allows surface water to freeze and form ice
 C     allowInteriorFreezing :: Allow water at depth to freeze and rise to the surface  
349  C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)  C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)
350  C     periodicExternalForcing :: Set true if forcing is time-dependant  C     periodicExternalForcing :: Set true if forcing is time-dependant
351  C- I/O parameters -  C- I/O parameters -
# Line 366  C     printDomain     :: controls printi Line 372  C     printDomain     :: controls printi
372       & fluidIsAir, fluidIsWater,       & fluidIsAir, fluidIsWater,
373       & usingPCoords, usingZCoords, useDynP_inEos_Zc,       & usingPCoords, usingZCoords, useDynP_inEos_Zc,
374       & usingCartesianGrid, usingSphericalPolarGrid, rotateGrid,       & usingCartesianGrid, usingSphericalPolarGrid, rotateGrid,
375       & usingCurvilinearGrid, usingCylindricalGrid,       & usingCylindricalGrid, usingCurvilinearGrid, hasWetCSCorners,
376       & deepAtmosphere, setInterFDr, setCenterDr,       & deepAtmosphere, setInterFDr, setCenterDr,
377       & no_slip_sides, no_slip_bottom,       & no_slip_sides, no_slip_bottom,
378       & useFullLeith, useStrainTensionVisc, useAreaViscLength,       & useFullLeith, useStrainTensionVisc, useAreaViscLength,
# Line 378  C     printDomain     :: controls printi Line 384  C     printDomain     :: controls printi
384       & upwindVorticity, highOrderVorticity,       & upwindVorticity, highOrderVorticity,
385       & useAbsVorticity, upwindShear,       & useAbsVorticity, upwindShear,
386       & momStepping, calc_wVelocity, tempStepping, saltStepping,       & momStepping, calc_wVelocity, tempStepping, saltStepping,
387       & tempAdvection, tempIsActiveTr, tempForcing,       & addFrictionHeating,
388       & saltAdvection, saltIsActiveTr, saltForcing,       & tempAdvection, tempVertDiff4, tempIsActiveTr, tempForcing,
389         & saltAdvection, saltVertDiff4, saltIsActiveTr, saltForcing,
390       & maskIniTemp, maskIniSalt, checkIniTemp, checkIniSalt,       & maskIniTemp, maskIniSalt, checkIniTemp, checkIniSalt,
391       & useRealFreshWaterFlux, useSRCGSolver,       & useSRCGSolver,
392       & rigidLid, implicitFreeSurface, exactConserv, linFSConserveTr,       & rigidLid, implicitFreeSurface,
393       & uniformLin_PhiSurf,       & uniformLin_PhiSurf, uniformFreeSurfLev,
394         & exactConserv, linFSConserveTr, useRealFreshWaterFlux,
395       & quasiHydrostatic, nonHydrostatic, use3Dsolver,       & quasiHydrostatic, nonHydrostatic, use3Dsolver,
396       & implicitIntGravWave, staggerTimeStep, doResetHFactors,       & implicitIntGravWave, staggerTimeStep, doResetHFactors,
397       & implicitDiffusion, implicitViscosity,       & implicitDiffusion, implicitViscosity,
# Line 391  C     printDomain     :: controls printi Line 399  C     printDomain     :: controls printi
399       & multiDimAdvection, useMultiDimAdvec,       & multiDimAdvection, useMultiDimAdvec,
400       & momDissip_In_AB, doAB_onGtGs,       & momDissip_In_AB, doAB_onGtGs,
401       & balanceEmPmR, balanceQnet, balancePrintMean,       & balanceEmPmR, balanceQnet, balancePrintMean,
402         & balanceThetaClimRelax, balanceSaltClimRelax,
403       & doThetaClimRelax, doSaltClimRelax,       & doThetaClimRelax, doSaltClimRelax,
404       & allowFreezing, allowInteriorFreezing, useOldFreezing,       & allowFreezing, useOldFreezing,
405       & periodicExternalForcing,       & periodicExternalForcing,
406       & globalFiles,       & globalFiles,
407       & pickupStrictlyMatch, usePickupBeforeC54, startFromPickupAB2,       & pickupStrictlyMatch, usePickupBeforeC54, startFromPickupAB2,
# Line 400  C     printDomain     :: controls printi Line 409  C     printDomain     :: controls printi
409       & writePickupAtEnd,       & writePickupAtEnd,
410       & timeave_mdsio, snapshot_mdsio, monitor_stdio,       & timeave_mdsio, snapshot_mdsio, monitor_stdio,
411       & outputTypesInclusive, dumpInitAndLast,       & outputTypesInclusive, dumpInitAndLast,
412       & printDomain,       & printDomain
      & inAdMode, inAdTrue, inAdFalse, inAdExact  
413    
414        LOGICAL fluidIsAir        LOGICAL fluidIsAir
415        LOGICAL fluidIsWater        LOGICAL fluidIsWater
# Line 411  C     printDomain     :: controls printi Line 419  C     printDomain     :: controls printi
419        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
420        LOGICAL usingSphericalPolarGrid, rotateGrid        LOGICAL usingSphericalPolarGrid, rotateGrid
421        LOGICAL usingCylindricalGrid        LOGICAL usingCylindricalGrid
422        LOGICAL usingCurvilinearGrid        LOGICAL usingCurvilinearGrid, hasWetCSCorners
423        LOGICAL deepAtmosphere        LOGICAL deepAtmosphere
424        LOGICAL setInterFDr        LOGICAL setInterFDr
425        LOGICAL setCenterDr        LOGICAL setCenterDr
426        LOGICAL useNHMTerms  
427        LOGICAL no_slip_sides        LOGICAL no_slip_sides
428        LOGICAL no_slip_bottom        LOGICAL no_slip_bottom
429        LOGICAL useFullLeith        LOGICAL useFullLeith
# Line 425  C     printDomain     :: controls printi Line 433  C     printDomain     :: controls printi
433        LOGICAL momAdvection        LOGICAL momAdvection
434        LOGICAL momForcing        LOGICAL momForcing
435        LOGICAL momPressureForcing        LOGICAL momPressureForcing
436          LOGICAL metricTerms
437          LOGICAL useNHMTerms
438    
439        LOGICAL useCoriolis        LOGICAL useCoriolis
440          LOGICAL use3dCoriolis
441          LOGICAL useCDscheme
442        LOGICAL vectorInvariantMomentum        LOGICAL vectorInvariantMomentum
443          LOGICAL useEnergyConservingCoriolis
444          LOGICAL useJamartWetPoints
445          LOGICAL useJamartMomAdv
446          LOGICAL upwindVorticity
447          LOGICAL highOrderVorticity
448          LOGICAL useAbsVorticity
449          LOGICAL upwindShear
450          LOGICAL momStepping
451          LOGICAL calc_wVelocity
452          LOGICAL tempStepping
453          LOGICAL saltStepping
454          LOGICAL addFrictionHeating
455        LOGICAL tempAdvection        LOGICAL tempAdvection
456          LOGICAL tempVertDiff4
457        LOGICAL tempIsActiveTr        LOGICAL tempIsActiveTr
458        LOGICAL tempForcing        LOGICAL tempForcing
459        LOGICAL saltAdvection        LOGICAL saltAdvection
460          LOGICAL saltVertDiff4
461        LOGICAL saltIsActiveTr        LOGICAL saltIsActiveTr
462        LOGICAL saltForcing        LOGICAL saltForcing
463        LOGICAL maskIniTemp        LOGICAL maskIniTemp
464        LOGICAL maskIniSalt        LOGICAL maskIniSalt
465        LOGICAL checkIniTemp        LOGICAL checkIniTemp
466        LOGICAL checkIniSalt        LOGICAL checkIniSalt
       LOGICAL useRealFreshWaterFlux  
467        LOGICAL useSRCGSolver        LOGICAL useSRCGSolver
468        LOGICAL rigidLid        LOGICAL rigidLid
469        LOGICAL implicitFreeSurface        LOGICAL implicitFreeSurface
470          LOGICAL uniformLin_PhiSurf
471          LOGICAL uniformFreeSurfLev
472        LOGICAL exactConserv        LOGICAL exactConserv
473        LOGICAL linFSConserveTr        LOGICAL linFSConserveTr
474        LOGICAL uniformLin_PhiSurf        LOGICAL useRealFreshWaterFlux
475        LOGICAL quasiHydrostatic        LOGICAL quasiHydrostatic
476        LOGICAL nonHydrostatic        LOGICAL nonHydrostatic
477        LOGICAL use3Dsolver        LOGICAL use3Dsolver
478        LOGICAL implicitIntGravWave        LOGICAL implicitIntGravWave
479        LOGICAL staggerTimeStep        LOGICAL staggerTimeStep
480        LOGICAL doResetHFactors        LOGICAL doResetHFactors
       LOGICAL momStepping  
       LOGICAL calc_wVelocity  
       LOGICAL tempStepping  
       LOGICAL saltStepping  
       LOGICAL metricTerms  
       LOGICAL use3dCoriolis  
       LOGICAL useCDscheme  
       LOGICAL useEnergyConservingCoriolis  
       LOGICAL useJamartWetPoints  
       LOGICAL useJamartMomAdv  
       LOGICAL upwindVorticity  
       LOGICAL highOrderVorticity  
       LOGICAL useAbsVorticity  
       LOGICAL upwindShear  
481        LOGICAL implicitDiffusion        LOGICAL implicitDiffusion
482        LOGICAL implicitViscosity        LOGICAL implicitViscosity
483        LOGICAL tempImplVertAdv        LOGICAL tempImplVertAdv
# Line 478  C     printDomain     :: controls printi Line 492  C     printDomain     :: controls printi
492        LOGICAL balancePrintMean        LOGICAL balancePrintMean
493        LOGICAL doThetaClimRelax        LOGICAL doThetaClimRelax
494        LOGICAL doSaltClimRelax        LOGICAL doSaltClimRelax
495          LOGICAL balanceThetaClimRelax
496          LOGICAL balanceSaltClimRelax
497        LOGICAL allowFreezing        LOGICAL allowFreezing
       LOGICAL allowInteriorFreezing  
498        LOGICAL useOldFreezing        LOGICAL useOldFreezing
499        LOGICAL periodicExternalForcing        LOGICAL periodicExternalForcing
500        LOGICAL globalFiles        LOGICAL globalFiles
# Line 492  C     printDomain     :: controls printi Line 507  C     printDomain     :: controls printi
507        LOGICAL outputTypesInclusive        LOGICAL outputTypesInclusive
508        LOGICAL dumpInitAndLast        LOGICAL dumpInitAndLast
509        LOGICAL printDomain        LOGICAL printDomain
       LOGICAL inAdMode, inAdTrue, inAdFalse, inAdExact  
510    
511  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
512  C     cg2dTargetResidual  C     cg2dTargetResidual
# Line 607  C     viscA4GridMax:: maximum and minimu Line 621  C     viscA4GridMax:: maximum and minimu
621  C     viscA4GridMin::  in terms of non-dimensional grid-size dependent viscosity  C     viscA4GridMin::  in terms of non-dimensional grid-size dependent viscosity
622  C     diffKhT   :: Laplacian diffusion coeff. for mixing of  C     diffKhT   :: Laplacian diffusion coeff. for mixing of
623  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 )  
624  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of  C     diffK4T   :: Biharmonic diffusion coeff. for mixing of
625  C                 heat laterally ( m^4/s )  C                 heat laterally ( m^4/s )
626    C     diffKrNrT :: vertical profile of Laplacian diffusion coeff.
627    C                 for mixing of heat vertically ( units of r^2/s )
628    C     diffKr4T  :: vertical profile of Biharmonic diffusion coeff.
629    C                 for mixing of heat vertically ( units of r^4/s )
630  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of  C     diffKhS  ::  Laplacian diffusion coeff. for mixing of
631  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 ),  
632  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of  C     diffK4S   :: Biharmonic diffusion coeff. for mixing of
633  C                 salt laterally ( m^4/s )  C                 salt laterally ( m^4/s )
634    C     diffKrNrS :: vertical profile of Laplacian diffusion coeff.
635    C                 for mixing of salt vertically ( units of r^2/s ),
636    C     diffKr4S  :: vertical profile of Biharmonic diffusion coeff.
637    C                 for mixing of salt vertically ( units of r^4/s )
638  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
639  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
640  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 636  C                    Frequency of checkp Line 654  C                    Frequency of checkp
654  C                    are referenced to this clock. ( s )  C                    are referenced to this clock. ( s )
655  C     deltaTMom    :: Timestep for momemtum equations ( s )  C     deltaTMom    :: Timestep for momemtum equations ( s )
656  C     dTtracerLev  :: Timestep for tracer equations ( s ), function of level k  C     dTtracerLev  :: Timestep for tracer equations ( s ), function of level k
657  C     deltaTfreesurf :: Timestep for free-surface equation ( s )  C     deltaTFreeSurf :: Timestep for free-surface equation ( s )
658  C     freesurfFac  :: Parameter to turn implicit free surface term on or off  C     freeSurfFac  :: Parameter to turn implicit free surface term on or off
659  C                    freesurfac = 1. uses implicit free surface  C                     freeSurFac = 1. uses implicit free surface
660  C                    freesurfac = 0. uses rigid lid  C                     freeSurFac = 0. uses rigid lid
661  C     abEps        :: Adams-Bashforth-2 stabilizing weight  C     abEps        :: Adams-Bashforth-2 stabilizing weight
662  C     alph_AB      :: Adams-Bashforth-3 primary factor  C     alph_AB      :: Adams-Bashforth-3 primary factor
663  C     beta_AB      :: Adams-Bashforth-3 secondary factor  C     beta_AB      :: Adams-Bashforth-3 secondary factor
# Line 664  C     endTime       :: Ending time for t Line 682  C     endTime       :: Ending time for t
682  C     chkPtFreq     :: Frequency of rolling check pointing ( s ).  C     chkPtFreq     :: Frequency of rolling check pointing ( s ).
683  C     pChkPtFreq    :: Frequency of permanent check pointing ( s ).  C     pChkPtFreq    :: Frequency of permanent check pointing ( s ).
684  C     dumpFreq      :: Frequency with which model state is written to  C     dumpFreq      :: Frequency with which model state is written to
685  C                     post-processing files ( s ).  C                      post-processing files ( s ).
686  C     diagFreq      :: Frequency with which model writes diagnostic output  C     diagFreq      :: Frequency with which model writes diagnostic output
687  C                     of intermediate quantities.  C                      of intermediate quantities.
688  C     afFacMom      :: Advection of momentum term tracer parameter  C     afFacMom      :: Advection of momentum term tracer parameter
689  C     vfFacMom      :: Momentum viscosity tracer parameter  C     vfFacMom      :: Momentum viscosity tracer parameter
690  C     pfFacMom      :: Momentum pressure forcing tracer parameter  C     pfFacMom      :: Momentum pressure forcing tracer parameter
# Line 715  C                           is to be rep Line 733  C                           is to be rep
733  C                           (affects myabs, mymin, mymax)  C                           (affects myabs, mymin, mymax)
734  C     nh_Am2        :: scales the non-hydrostatic terms and changes internal scales  C     nh_Am2        :: scales the non-hydrostatic terms and changes internal scales
735  C                      (i.e. allows convection at different Rayleigh numbers)  C                      (i.e. allows convection at different Rayleigh numbers)
736    C     tCylIn        :: Temperature of the cylinder inner boundary
737    C     tCylOut       :: Temperature of the cylinder outer boundary
738  C     phiEuler      :: Euler angle, rotation about original z-axis  C     phiEuler      :: Euler angle, rotation about original z-axis
739  C     thetaEuler    :: Euler angle, rotation about new x-axis  C     thetaEuler    :: Euler angle, rotation about new x-axis
740  C     psiEuler      :: Euler angle, rotation about new z-axis  C     psiEuler      :: Euler angle, rotation about new z-axis
741        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
742       & cg2dpcOffDFac, cg3dTargetResidual,       & cg2dpcOffDFac, cg3dTargetResidual,
743       & delR, delRc, xgOrigin, ygOrigin,       & delR, delRc, xgOrigin, ygOrigin,
744       & deltaT, deltaTmom, dTtracerLev, deltaTfreesurf, deltaTClock,       & deltaT, deltaTMom, dTtracerLev, deltaTFreeSurf, deltaTClock,
745       & abEps, alph_AB, beta_AB,       & abEps, alph_AB, beta_AB,
746       & rSphere, recip_rSphere, radius_fromHorizGrid,       & rSphere, recip_rSphere, radius_fromHorizGrid,
747       & f0, beta, fPrime, omega, rotationPeriod,       & f0, beta, fPrime, omega, rotationPeriod,
# Line 734  C     psiEuler      :: Euler angle, rota Line 754  C     psiEuler      :: Euler angle, rota
754       & viscA4Grid, viscA4GridMax, viscA4GridMin,       & viscA4Grid, viscA4GridMax, viscA4GridMin,
755       & viscAhReMax, viscA4ReMax,       & viscAhReMax, viscA4ReMax,
756       & viscC4leith, viscC4leithD, viscArNr,       & viscC4leith, viscC4leithD, viscArNr,
757       & diffKhT, diffK4T, diffKrNrT,       & diffKhT, diffK4T, diffKrNrT, diffKr4T,
758       & diffKhS, diffK4S, diffKrNrS,       & diffKhS, diffK4S, diffKrNrS, diffKr4S,
759       & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,       & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,
760       & BL79LatVary,       & BL79LatVary,
761       & diffKrBLEQsurf, diffKrBLEQdeep, diffKrBLEQscl, diffKrBLEQHo,       & diffKrBLEQsurf, diffKrBLEQdeep, diffKrBLEQscl, diffKrBLEQHo,
# Line 773  C     psiEuler      :: Euler angle, rota Line 793  C     psiEuler      :: Euler angle, rota
793        _RL ygOrigin        _RL ygOrigin
794        _RL deltaT        _RL deltaT
795        _RL deltaTClock        _RL deltaTClock
796        _RL deltaTmom        _RL deltaTMom
797        _RL dTtracerLev(Nr)        _RL dTtracerLev(Nr)
798        _RL deltaTfreesurf        _RL deltaTFreeSurf
799        _RL abEps, alph_AB, beta_AB        _RL abEps, alph_AB, beta_AB
800        _RL rSphere        _RL rSphere
801        _RL recip_rSphere        _RL recip_rSphere
# Line 818  C     psiEuler      :: Euler angle, rota Line 838  C     psiEuler      :: Euler angle, rota
838        _RL viscC4leithD        _RL viscC4leithD
839        _RL viscC4smag        _RL viscC4smag
840        _RL diffKhT        _RL diffKhT
       _RL diffKrNrT(Nr)  
841        _RL diffK4T        _RL diffK4T
842          _RL diffKrNrT(Nr)
843          _RL diffKr4T(Nr)
844        _RL diffKhS        _RL diffKhS
       _RL diffKrNrS(Nr)  
845        _RL diffK4S        _RL diffK4S
846          _RL diffKrNrS(Nr)
847          _RL diffKr4S(Nr)
848        _RL diffKrBL79surf        _RL diffKrBL79surf
849        _RL diffKrBL79deep        _RL diffKrBL79deep
850        _RL diffKrBL79scl        _RL diffKrBL79scl
# Line 889  C     psiEuler      :: Euler angle, rota Line 911  C     psiEuler      :: Euler angle, rota
911        _RL bottomDragQuadratic        _RL bottomDragQuadratic
912        _RL smoothAbsFuncRange        _RL smoothAbsFuncRange
913        _RL nh_Am2        _RL nh_Am2
914        _RL tCylIn        _RL tCylIn, tCylOut
       _RL tCylOut  
915        _RL phiEuler, thetaEuler, psiEuler        _RL phiEuler, thetaEuler, psiEuler
916    
917  C--   COMMON /PARM_A/ Thermodynamics constants ?  C--   COMMON /PARM_A/ Thermodynamics constants ?
# Line 938  C Logical flags for selecting packages Line 959  C Logical flags for selecting packages
959        LOGICAL useCheapAML        LOGICAL useCheapAML
960        LOGICAL useGrdchk        LOGICAL useGrdchk
961        LOGICAL useSMOOTH        LOGICAL useSMOOTH
962          LOGICAL usePROFILES
963        LOGICAL useECCO        LOGICAL useECCO
964        LOGICAL useSBO        LOGICAL useSBO
965        LOGICAL useFLT        LOGICAL useFLT
# Line 946  C Logical flags for selecting packages Line 968  C Logical flags for selecting packages
968        LOGICAL useRBCS        LOGICAL useRBCS
969        LOGICAL useOffLine        LOGICAL useOffLine
970        LOGICAL useMATRIX        LOGICAL useMATRIX
971          LOGICAL useFRAZIL
972        LOGICAL useSEAICE        LOGICAL useSEAICE
973        LOGICAL useSALT_PLUME        LOGICAL useSALT_PLUME
974        LOGICAL useShelfIce        LOGICAL useShelfIce
975          LOGICAL useStreamIce
976        LOGICAL useICEFRONT        LOGICAL useICEFRONT
977        LOGICAL useThSIce        LOGICAL useThSIce
978        LOGICAL useATM2d        LOGICAL useATM2d
# Line 968  C Logical flags for selecting packages Line 992  C Logical flags for selecting packages
992       &        useOPPS, usePP81, useMY82, useGGL90, useKPP,       &        useOPPS, usePP81, useMY82, useGGL90, useKPP,
993       &        useGMRedi, useBBL, useDOWN_SLOPE,       &        useGMRedi, useBBL, useDOWN_SLOPE,
994       &        useCAL, useEXF, useBulkForce, useEBM, useCheapAML,       &        useCAL, useEXF, useBulkForce, useEBM, useCheapAML,
995       &        useGrdchk, useSMOOTH, useECCO, useSBO, useFLT,       &        useGrdchk,useSMOOTH,usePROFILES,useECCO,useSBO, useFLT,
996       &        usePTRACERS, useGCHEM, useRBCS, useOffLine, useMATRIX,       &        usePTRACERS, useGCHEM, useRBCS, useOffLine, useMATRIX,
997       &        useSEAICE, useSALT_PLUME, useShelfIce,       &        useFRAZIL, useSEAICE, useSALT_PLUME, useShelfIce,
998       &        useICEFRONT, useThSIce,       &        useStreamIce, useICEFRONT, useThSIce,
999       &        useATM2D, useAIM, useLand, useFizhi, useGridAlt,       &        useATM2D, useAIM, useLand, useFizhi, useGridAlt,
1000       &        useDiagnostics, useREGRID, useLayers, useMNC,       &        useDiagnostics, useREGRID, useLayers, useMNC,
1001       &        useRunClock, useEMBED_FILES,       &        useRunClock, useEMBED_FILES,
1002       &        useMYPACKAGE       &        useMYPACKAGE
1003    
 C     Logical flags for turning off parts of the code in adjoint mode  
       LOGICAL useKPPinAdMode, useKPPinFwdMode  
       LOGICAL useGMrediInAdMode, useGMrediInFwdMode  
       LOGICAL useSEAICEinAdMode, useSEAICEinFwdMode  
       COMMON /PARM_PACKAGES_ADJ/  
      &       useKPPinAdMode, useKPPinFwdMode,  
      &       useGMrediInAdMode, useGMrediInFwdMode,  
      &       useSEAICEinAdMode, useSEAICEinFwdMode  
   
1004  CEH3 ;;; Local Variables: ***  CEH3 ;;; Local Variables: ***
1005  CEH3 ;;; mode:fortran ***  CEH3 ;;; mode:fortran ***
1006  CEH3 ;;; End: ***  CEH3 ;;; End: ***
Legend:
Removed from v.1.257  
changed lines
  Added in v.1.266
  ViewVC Help
Powered by ViewVC 1.1.22