/[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.157 by adcroft, Tue May 31 14:49:38 2005 UTC revision 1.173 by jmc, Tue Dec 13 16:33:43 2005 UTC
# Line 49  C     Checkpoint data Line 49  C     Checkpoint data
49    
50  C--   COMMON /PARM_C/ Character valued parameters used by the model.  C--   COMMON /PARM_C/ Character valued parameters used by the model.
51  C     checkPtSuff :: List of checkpoint file suffices  C     checkPtSuff :: List of checkpoint file suffices
52    C     delXFile      :: File containing X-spacing grid definition (1.D array)
53    C     delYFile      :: File containing Y-spacing grid definition (1.D array)
54    C     horizGridFile :: File containing horizontal-grid definition
55    C                        (only when using curvilinear_grid)
56  C     bathyFile   :: File containing bathymetry. If not defined bathymetry  C     bathyFile   :: File containing bathymetry. If not defined bathymetry
57  C                   is taken from inline function.  C                   is taken from inline function.
58  C     topoFile    :: File containing the topography of the surface (unit=m)  C     topoFile    :: File containing the topography of the surface (unit=m)
# Line 78  C     eosType         :: choose the equa Line 82  C     eosType         :: choose the equa
82  C                        LINEAR, POLY3, UNESCO, JMD95Z, JMD95P, MDJWF, IDEALGAS  C                        LINEAR, POLY3, UNESCO, JMD95Z, JMD95P, MDJWF, IDEALGAS
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/ checkPtSuff,        COMMON /PARM_C/ checkPtSuff,
85         &                delXFile, delYFile, horizGridFile,
86       &                bathyFile, topoFile,       &                bathyFile, topoFile,
87       &                hydrogThetaFile, hydrogSaltFile,       &                hydrogThetaFile, hydrogSaltFile,
88       &                zonalWindFile, meridWindFile, thetaClimFile,       &                zonalWindFile, meridWindFile, thetaClimFile,
# Line 92  C     the_run_name    :: string identify Line 97  C     the_run_name    :: string identify
97       &                mdsioLocalDir,       &                mdsioLocalDir,
98       &                the_run_name       &                the_run_name
99        CHARACTER*(5) checkPtSuff(maxNoChkptLev)        CHARACTER*(5) checkPtSuff(maxNoChkptLev)
100          CHARACTER*(MAX_LEN_FNAM) delXFile
101          CHARACTER*(MAX_LEN_FNAM) delYFile
102          CHARACTER*(MAX_LEN_FNAM) horizGridFile
103        CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile        CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile
104        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile
105        CHARACTER*(MAX_LEN_FNAM) hydrogSaltFile        CHARACTER*(MAX_LEN_FNAM) hydrogSaltFile
# Line 143  C                           =0 Linear fr Line 151  C                           =0 Linear fr
151  C     select_rStar        :: option related to r* vertical coordinate  C     select_rStar        :: option related to r* vertical coordinate
152  C                           =0 (default) use r coord. ; > 0 use r*  C                           =0 (default) use r coord. ; > 0 use r*
153  C     tempAdvScheme       :: Temp. Horiz.Advection scheme selector  C     tempAdvScheme       :: Temp. Horiz.Advection scheme selector
154  C     tempVertAdvScheme   :: Temp. Vert. Advection scheme selector  C     tempVertAdvScheme   :: Temp. Vert. Advection scheme selector
155  C     saltAdvScheme       :: Salt. Horiz.advection scheme selector  C     saltAdvScheme       :: Salt. Horiz.advection scheme selector
156  C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector  C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector
157    C     selectKEscheme      :: Kinetic Energy scheme selector (Vector Inv.)
158  C     debugLevel          :: debug level selector: higher -> more writing  C     debugLevel          :: debug level selector: higher -> more writing
159    
160        COMMON /PARM_I/        COMMON /PARM_I/
# Line 160  C     debugLevel          :: debug level Line 169  C     debugLevel          :: debug level
169       &        nonlinFreeSurf, select_rStar,       &        nonlinFreeSurf, select_rStar,
170       &        tempAdvScheme, tempVertAdvScheme,       &        tempAdvScheme, tempVertAdvScheme,
171       &        saltAdvScheme, saltVertAdvScheme,       &        saltAdvScheme, saltVertAdvScheme,
172         &        selectKEscheme,
173       &        debugLevel       &        debugLevel
174        INTEGER cg2dMaxIters        INTEGER cg2dMaxIters
175        INTEGER cg2dChkResFreq        INTEGER cg2dChkResFreq
# Line 178  C     debugLevel          :: debug level Line 188  C     debugLevel          :: debug level
188        INTEGER select_rStar        INTEGER select_rStar
189        INTEGER tempAdvScheme, tempVertAdvScheme        INTEGER tempAdvScheme, tempVertAdvScheme
190        INTEGER saltAdvScheme, saltVertAdvScheme        INTEGER saltAdvScheme, saltVertAdvScheme
191          INTEGER selectKEscheme
192        INTEGER debugLevel        INTEGER debugLevel
193    
194  C  C
# Line 196  C                               spherica Line 207  C                               spherica
207  C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical  C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical
208  C     no_slip_sides :: Impose "no-slip" at lateral boundaries.  C     no_slip_sides :: Impose "no-slip" at lateral boundaries.
209  C     no_slip_bottom :: Impose "no-slip" at bottom boundary.  C     no_slip_bottom :: Impose "no-slip" at bottom boundary.
 C     staggerTimeStep :: enable a Stagger time stepping T,S Rho then U,V  
210  C     momViscosity  :: Flag which turns momentum friction terms on and off.  C     momViscosity  :: Flag which turns momentum friction terms on and off.
211  C     momAdvection  :: Flag which turns advection of momentum on and off.  C     momAdvection  :: Flag which turns advection of momentum on and off.
212  C     momForcing    :: Flag which turns external forcing of momentum on  C     momForcing    :: Flag which turns external forcing of momentum on
# Line 218  C                     and off. Line 228  C                     and off.
228  C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux  C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux
229  C                         as a real Fresh Water (=> changes the Sea Level)  C                         as a real Fresh Water (=> changes the Sea Level)
230  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)
231  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
232  C                       on irregular grids)  C                       on irregular grids)
233  C     useAnisotropicViscAGridMax :: Set to true to use Alistair's latest  C     useAreaViscLength :: Set to true to use old scaling for viscous
234  C              anisotropic length scale.  It is used only for maximum viscosity  C              lengths, e.g., L2=Raz.  May be preferable for cube sphere.
235  C              calculations.  Alistair recommends a value of viscA*GridMax=.25  C     useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms
236  C     rigidLid            :: Set to true to use rigid lid  C     rigidLid            :: Set to true to use rigid lid
237  C     implicitFreeSurface :: Set to true to use implcit free surface  C     implicitFreeSurface :: Set to true to use implicit free surface
238  C     exactConserv        :: Set to true to conserve exactly the total Volume  C     exactConserv        :: Set to true to conserve exactly the total Volume
239  C     uniformLin_PhiSurf  :: Set to true to use a uniform Bo_surf in the  C     uniformLin_PhiSurf  :: Set to true to use a uniform Bo_surf in the
240  C                           linear relation Phi_surf = Bo_surf*eta  C                           linear relation Phi_surf = Bo_surf*eta
241    C     implicitIntGravWave :: treat Internal Gravity Wave implicitly
242    C     staggerTimeStep :: enable a Stagger time stepping T,S Rho then U,V
243  C     momStepping   :: Turns momentum equation time-stepping off  C     momStepping   :: Turns momentum equation time-stepping off
244  C     tempStepping  :: Turns temperature equation time-stepping off  C     tempStepping  :: Turns temperature equation time-stepping off
245  C     saltStepping  :: Turns salinity equation time-stepping off  C     saltStepping  :: Turns salinity equation time-stepping off
# Line 239  C     useJamartWetPoints :: Use wet-poin Line 251  C     useJamartWetPoints :: Use wet-poin
251  C     useJamartMomAdv :: Use wet-point method for V.I. non-linear term  C     useJamartMomAdv :: Use wet-point method for V.I. non-linear term
252  C     SadournyCoriolis :: use the enstrophy conserving scheme by Sadourny  C     SadournyCoriolis :: use the enstrophy conserving scheme by Sadourny
253  C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term  C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term
254  C     highOrderVorticity :: use 3rd/4th order interp. of vorticity in Coriolis  C     highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection)
255    C     upwindShear        :: use 1rst order upwind interp. (V.I., vertical advection)
256  C     useAbsVorticity :: work with f+zeta in Coriolis terms  C     useAbsVorticity :: work with f+zeta in Coriolis terms
257  C     implicitDiffusion :: Turns implicit vertical diffusion on  C     implicitDiffusion :: Turns implicit vertical diffusion on
258  C     implicitViscosity :: Turns implicit vertical viscosity on  C     implicitViscosity :: Turns implicit vertical viscosity on
# Line 248  C     saltImplVertAdv :: Turns on implic Line 261  C     saltImplVertAdv :: Turns on implic
261  C     momImplVertAdv  :: Turns on implicit vertical advection for Momentum  C     momImplVertAdv  :: Turns on implicit vertical advection for Momentum
262  C     multiDimAdvection :: Flag that enable multi-dimension advection  C     multiDimAdvection :: Flag that enable multi-dimension advection
263  C     useMultiDimAdvec  :: True if multi-dim advection is used at least once  C     useMultiDimAdvec  :: True if multi-dim advection is used at least once
264  C     forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution  C     forcing_In_AB     :: if False, put forcing (Temp,Salt,Tracers) contribution
265  C                      out off Adams-Bashforth time stepping.  C                          out off Adams-Bashforth time stepping.
266    C     doAB_onGtGs       :: if the Adams-Bashforth time stepping is used, always
267    C                          apply AB on tracer tendencies (rather than on Tracer)
268  C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup  C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup
269  C     doThetaClimRelax :: Set true if relaxation to temperature  C     doThetaClimRelax :: Set true if relaxation to temperature
270  C                        climatology is required.  C                        climatology is required.
# Line 276  C                       useSingleCpuIO i Line 291  C                       useSingleCpuIO i
291  C                       outputs from master mpi process only.  C                       outputs from master mpi process only.
292  C     allowFreezing  :: Allows surface water to freeze and form ice  C     allowFreezing  :: Allows surface water to freeze and form ice
293  C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)  C     useOldFreezing :: use the old version (before checkpoint52a_pre, 2003-11-12)
 C     groundAtK1  :: put the surface(k=1) at the Lower Boundary (=ground)  
294  C     pickup_write_mdsio :: use mdsio to write pickups  C     pickup_write_mdsio :: use mdsio to write pickups
295  C     pickup_read_mdsio  :: use mdsio to read  pickups  C     pickup_read_mdsio  :: use mdsio to read  pickups
296  C     pickup_write_immed :: echo the pickup immediately (for conversion)  C     pickup_write_immed :: echo the pickup immediately (for conversion)
# Line 287  C     calendarDumps :: When set, approxi Line 301  C     calendarDumps :: When set, approxi
301  C                      for parameters chkPtFreq, pChkPtFreq, taveFreq, SEAICE_taveFreq,  C                      for parameters chkPtFreq, pChkPtFreq, taveFreq, SEAICE_taveFreq,
302  C                      KPP_taveFreq, and freq in pkg/diagnostics are converted to exact  C                      KPP_taveFreq, and freq in pkg/diagnostics are converted to exact
303  C                      calendar months and years.  Requires pkg/cal.  C                      calendar months and years.  Requires pkg/cal.
304    C     dumpInitAndLast :: dumps model state to files at Initial (nIter0)
305    C                        & Last iteration, in addition multiple of dumpFreq iter.
306        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
307       & usingCurvilinearGrid, usingCylindricalGrid,       & usingCurvilinearGrid, usingCylindricalGrid,
308       & no_slip_sides,no_slip_bottom,       & no_slip_sides,no_slip_bottom,
      & staggerTimeStep,  
309       & momViscosity, momAdvection, momForcing, useCoriolis,       & momViscosity, momAdvection, momForcing, useCoriolis,
310       & momPressureForcing, vectorInvariantMomentum,       & momPressureForcing, vectorInvariantMomentum,
311       & tempAdvection, tempForcing,       & tempAdvection, tempForcing,
312       & saltAdvection, saltForcing,       & saltAdvection, saltForcing,
313       & useRealFreshWaterFlux,       & useRealFreshWaterFlux,
314       & useFullLeith, useAnisotropicViscAGridMax,       & useFullLeith, useStrainTensionVisc,
315         & useAreaViscLength,
316       & rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf,       & rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf,
317         & implicitIntGravWave, staggerTimeStep,
318       & momStepping, tempStepping, saltStepping,       & momStepping, tempStepping, saltStepping,
319       & metricTerms, usingSphericalPolarMTerms, useNHMTerms,       & metricTerms, usingSphericalPolarMTerms, useNHMTerms,
320       & useConstantF, useBetaPlaneF, useSphereF,       & useConstantF, useBetaPlaneF, useSphereF,
321       & useCDscheme,       & useCDscheme,
322       & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,       & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,
323       & SadournyCoriolis, upwindVorticity, highOrderVorticity,       & SadournyCoriolis, upwindVorticity, highOrderVorticity,
324       & useAbsVorticity,       & useAbsVorticity, upwindShear,
325       & implicitDiffusion, implicitViscosity,       & implicitDiffusion, implicitViscosity,
326       & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,       & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,
327       & multiDimAdvection, useMultiDimAdvec, forcing_In_AB,       & multiDimAdvection, useMultiDimAdvec,
328         & forcing_In_AB, doAB_onGtGs,
329       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,
330       & periodicExternalForcing,       & periodicExternalForcing,
331       & fluidIsAir, fluidIsWater,       & fluidIsAir, fluidIsWater,
332       & usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,       & usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,
333       & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,       & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,
334       & allowFreezing, useOldFreezing, groundAtK1,       & allowFreezing, useOldFreezing,
335       & usePickupBeforeC35, usePickupBeforeC54, startFromPickupAB2,       & usePickupBeforeC35, usePickupBeforeC54, startFromPickupAB2,
336       & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,       & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,
337       & timeave_mdsio, snapshot_mdsio, monitor_stdio,       & timeave_mdsio, snapshot_mdsio, monitor_stdio,
338       & outputTypesInclusive, debugMode,       & outputTypesInclusive, dumpInitAndLast, debugMode,
339       & inAdMode, inAdTrue, inAdFalse, inAdExact,       & inAdMode, inAdTrue, inAdFalse, inAdExact,
340       & calendarDumps       & calendarDumps
341    
# Line 329  C                      calendar months a Line 347  C                      calendar months a
347        LOGICAL useNHMTerms        LOGICAL useNHMTerms
348        LOGICAL no_slip_sides        LOGICAL no_slip_sides
349        LOGICAL no_slip_bottom        LOGICAL no_slip_bottom
       LOGICAL staggerTimeStep  
350        LOGICAL momViscosity        LOGICAL momViscosity
351        LOGICAL momAdvection        LOGICAL momAdvection
352        LOGICAL momForcing        LOGICAL momForcing
# Line 341  C                      calendar months a Line 358  C                      calendar months a
358        LOGICAL saltAdvection        LOGICAL saltAdvection
359        LOGICAL saltForcing        LOGICAL saltForcing
360        LOGICAL useRealFreshWaterFlux        LOGICAL useRealFreshWaterFlux
361        LOGICAL useFullLeith, useAnisotropicViscAGridMax        LOGICAL useFullLeith
362          LOGICAL useStrainTensionVisc
363          LOGICAL useAreaViscLength
364        LOGICAL rigidLid        LOGICAL rigidLid
365        LOGICAL implicitFreeSurface        LOGICAL implicitFreeSurface
366        LOGICAL exactConserv        LOGICAL exactConserv
367        LOGICAL uniformLin_PhiSurf        LOGICAL uniformLin_PhiSurf
368          LOGICAL implicitIntGravWave
369          LOGICAL staggerTimeStep
370        LOGICAL momStepping        LOGICAL momStepping
371        LOGICAL tempStepping        LOGICAL tempStepping
372        LOGICAL saltStepping        LOGICAL saltStepping
# Line 361  C                      calendar months a Line 382  C                      calendar months a
382        LOGICAL upwindVorticity        LOGICAL upwindVorticity
383        LOGICAL highOrderVorticity        LOGICAL highOrderVorticity
384        LOGICAL useAbsVorticity        LOGICAL useAbsVorticity
385          LOGICAL upwindShear
386        LOGICAL implicitDiffusion        LOGICAL implicitDiffusion
387        LOGICAL implicitViscosity        LOGICAL implicitViscosity
388        LOGICAL tempImplVertAdv        LOGICAL tempImplVertAdv
# Line 369  C                      calendar months a Line 391  C                      calendar months a
391        LOGICAL multiDimAdvection        LOGICAL multiDimAdvection
392        LOGICAL useMultiDimAdvec        LOGICAL useMultiDimAdvec
393        LOGICAL forcing_In_AB        LOGICAL forcing_In_AB
394          LOGICAL doAB_onGtGs
395        LOGICAL doThetaClimRelax        LOGICAL doThetaClimRelax
396        LOGICAL doSaltClimRelax        LOGICAL doSaltClimRelax
397        LOGICAL doTr1ClimRelax        LOGICAL doTr1ClimRelax
# Line 385  C                      calendar months a Line 408  C                      calendar months a
408        LOGICAL useSingleCpuIO        LOGICAL useSingleCpuIO
409        LOGICAL allowFreezing        LOGICAL allowFreezing
410        LOGICAL useOldFreezing        LOGICAL useOldFreezing
       LOGICAL groundAtK1  
411        LOGICAL usePickupBeforeC35        LOGICAL usePickupBeforeC35
412        LOGICAL usePickupBeforeC54        LOGICAL usePickupBeforeC54
413        LOGICAL startFromPickupAB2        LOGICAL startFromPickupAB2
414          LOGICAL dumpInitAndLast
415        LOGICAL debugMode        LOGICAL debugMode
416        LOGICAL pickup_read_mdsio, pickup_write_mdsio        LOGICAL pickup_read_mdsio, pickup_write_mdsio
417        LOGICAL pickup_write_immed        LOGICAL pickup_write_immed
# Line 430  C     rhoConst  :: Vertically constant r Line 453  C     rhoConst  :: Vertically constant r
453  C     rhoConstFresh :: Constant reference density for fresh water (rain)  C     rhoConstFresh :: Constant reference density for fresh water (rain)
454  C     tRef      :: reference vertical profile for potential temperature  C     tRef      :: reference vertical profile for potential temperature
455  C     sRef      :: reference vertical profile for salinity/specific humidity  C     sRef      :: reference vertical profile for salinity/specific humidity
456    C     phiRef    :: reference potential (pressure/rho, geopotential) profile
457    C     dBdrRef   :: vertical gradient of reference boyancy  [(m/s/r)^2)]:
458    C               :: z-coord: = N^2_ref = Brunt-Vaissala frequency [s^-2]
459    C               :: p-coord: = -(d.alpha/dp)_ref          [(m^2.s/kg)^2]
460  C     phiMin    :: Latitude of southern most cell face.  C     phiMin    :: Latitude of southern most cell face.
461  C     thetaMin  :: Longitude of western most cell face (this  C     thetaMin  :: Longitude of western most cell face (this
462  C                 is an "inert" parameter but it is included  C                 is an "inert" parameter but it is included
# Line 462  C                  (act on Divergence pa Line 489  C                  (act on Divergence pa
489  C     viscA4Z   :: Biharmonic viscosity coeff. for mixing of momentum laterally  C     viscA4Z   :: Biharmonic viscosity coeff. for mixing of momentum laterally
490  C                  (act on Vorticity  part) ( m^4/s )  C                  (act on Vorticity  part) ( m^4/s )
491  C     viscC2leith :: Leith non-dimensional viscosity factor (grad(vort))  C     viscC2leith :: Leith non-dimensional viscosity factor (grad(vort))
492  C     viscC2leithD :: Modified Leith non-dimensional viscosity factor (grad(div))  C     viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div))
493  C     viscC2smag  :: Smagorinsky non-dimensional viscosity factor  C     viscC2smag  :: Smagorinsky non-dimensional viscosity factor (harmonic)
494    C     viscC4smag  :: Smagorinsky non-dimensional viscosity factor (biharmonic)
495  C     viscAhMax :: Maximum eddy viscosity coeff. for mixing of  C     viscAhMax :: Maximum eddy viscosity coeff. for mixing of
496  C                 momentum laterally ( m^2/s )  C                   momentum laterally ( m^2/s )
497    C     viscAhReMax :: Maximum gridscale Reynolds number for eddy viscosity
498    C                   coeff. for mixing of momentum laterally (non-dim)
499  C     viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ...  C     viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ...
500  C     viscAhGridMin::  in terms of non-dimensional grid-size dependent viscosity  C     viscAhGridMin::  in terms of non-dimensional grid-size dependent visc.
501  C     viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of  C     viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of
502  C                 momentum laterally ( m^4/s )  C                 momentum laterally ( m^4/s )
503    C     viscA4ReMax :: Maximum Gridscale Reynolds number for
504    C                  biharmonic viscosity coeff. momentum laterally (non-dim)
505  C     viscAhGrid:: non-dimensional grid-size dependent viscosity  C     viscAhGrid:: non-dimensional grid-size dependent viscosity
506  C     viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity  C     viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity
507  C     viscA4GridMax:: maximum and minimum biharmonic viscosity coefficients ...  C     viscA4GridMax:: maximum and minimum biharmonic viscosity coefficients ...
# Line 563  C     recip_horiVertRatio  ( 1 if horiz Line 595  C     recip_horiVertRatio  ( 1 if horiz
595  C                          ( g*rho if horiz in m and vertical in Pa ).  C                          ( g*rho if horiz in m and vertical in Pa ).
596  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
597  C                          starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)  C                          starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)
598    C     sideDragFactor     :: side-drag scaling factor (used only if no_slip_sides)
599    C                           (default=2: full drag ; =1: gives half-slip BC)
600  C     bottomDragLinear   :: Drag coefficient built in to core dynamics  C     bottomDragLinear   :: Drag coefficient built in to core dynamics
601  C      --"-"--  Quadratic  ( linear: 1/s, quadratic: 1/m )  C     bottomDragQuadratic   ( linear: 1/s, quadratic: 1/m )
602    C     nh_Am2        :: scales the non-hydrostatic terms and changes internal scales
603    C                      (i.e. allows convection at different Rayleigh numbers)
604        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
605       & cg2dpcOffDFac, cg3dTargetResidual,       & cg2dpcOffDFac, cg3dTargetResidual,
606       & delR, delRc, delX, delY,       & delR, delRc, delX, delY,
# Line 574  C      --"-"--  Quadratic  ( linear: 1/s Line 610  C      --"-"--  Quadratic  ( linear: 1/s
610       & viscAh, viscAhW, viscAhMax,       & viscAh, viscAhW, viscAhMax,
611       & viscAhGrid, viscAhGridMax, viscAhGridMin,       & viscAhGrid, viscAhGridMax, viscAhGridMin,
612       & viscC2leith, viscC2leithD,       & viscC2leith, viscC2leithD,
613       & viscC2smag,       & viscC2smag, viscC4smag,
614       & viscAhD, viscAhZ, viscA4D, viscA4Z,       & viscAhD, viscAhZ, viscA4D, viscA4Z,
615       & viscA4, viscA4W,       & viscA4, viscA4W,
616       & viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin,       & viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin,
617       & viscC4leith, viscC4leithD, viscAr, viscAstrain, viscAtension,       & viscAhRemax, viscA4Remax,
618         & viscC4leith, viscC4leithD, viscAr,
619       & diffKhT, diffK4T, diffKrNrT,       & diffKhT, diffK4T, diffKrNrT,
620       & diffKhS, diffK4S, diffKrNrS,       & diffKhS, diffK4S, diffKrNrS,
621       & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,       & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,
# Line 586  C      --"-"--  Quadratic  ( linear: 1/s Line 623  C      --"-"--  Quadratic  ( linear: 1/s
623       & hFacMin, hFacMinDz, hFacInf, hFacSup,       & hFacMin, hFacMinDz, hFacInf, hFacSup,
624       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
625       & recip_rhoConst, rhoConst,       & recip_rhoConst, rhoConst,
626       & rhoConstFresh, convertEmP2rUnit, tRef, sRef,       & rhoConstFresh, convertEmP2rUnit, tRef, sRef, phiRef, dBdrRef,
627       & baseTime, startTime, endTime,       & baseTime, startTime, endTime,
628       & chkPtFreq, pchkPtFreq, dumpFreq, adjDumpFreq,       & chkPtFreq, pchkPtFreq, dumpFreq, adjDumpFreq,
629       & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq,       & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq,
# Line 600  C      --"-"--  Quadratic  ( linear: 1/s Line 637  C      --"-"--  Quadratic  ( linear: 1/s
637       & hFacMinDr, hFacMinDp,       & hFacMinDr, hFacMinDp,
638       & horiVertRatio, recip_horiVertRatio,       & horiVertRatio, recip_horiVertRatio,
639       & ivdc_kappa, Ro_SeaLevel,       & ivdc_kappa, Ro_SeaLevel,
640       & bottomDragLinear,bottomDragQuadratic,nh_Am2,       & sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2,
641       & tCylIn, tCylOut       & tCylIn, tCylOut
642    
       _RL nh_Am2  
643        _RL cg2dTargetResidual        _RL cg2dTargetResidual
644        _RL cg2dTargetResWunit        _RL cg2dTargetResWunit
645        _RL cg3dTargetResidual        _RL cg3dTargetResidual
# Line 638  C      --"-"--  Quadratic  ( linear: 1/s Line 674  C      --"-"--  Quadratic  ( linear: 1/s
674        _RL viscAhD        _RL viscAhD
675        _RL viscAhZ        _RL viscAhZ
676        _RL viscAhMax        _RL viscAhMax
677          _RL viscAhReMax
678        _RL viscAhGrid        _RL viscAhGrid
679        _RL viscAhGridMax        _RL viscAhGridMax
680        _RL viscAhGridMin        _RL viscAhGridMin
681        _RL viscC2leith        _RL viscC2leith
682        _RL viscC2leithD        _RL viscC2leithD
683        _RL viscC2smag        _RL viscC2smag
684        _RL viscAstrain        _RL viscC4smag
       _RL viscAtension  
685        _RL viscAr        _RL viscAr
686        _RL viscA4        _RL viscA4
687        _RL viscA4W        _RL viscA4W
688        _RL viscA4D        _RL viscA4D
689        _RL viscA4Z        _RL viscA4Z
690        _RL viscA4Max        _RL viscA4Max
691          _RL viscA4ReMax
692        _RL viscA4Grid, viscA4GridMax, viscA4GridMin        _RL viscA4Grid, viscA4GridMax, viscA4GridMin
693        _RL viscC4leith        _RL viscC4leith
694        _RL viscC4leithD        _RL viscC4leithD
# Line 679  C      --"-"--  Quadratic  ( linear: 1/s Line 716  C      --"-"--  Quadratic  ( linear: 1/s
716        _RL convertEmP2rUnit        _RL convertEmP2rUnit
717        _RL tRef(Nr)        _RL tRef(Nr)
718        _RL sRef(Nr)        _RL sRef(Nr)
719          _RL phiRef(2*Nr+1)
720          _RL dBdrRef(Nr)
721        _RL baseTime        _RL baseTime
722        _RL startTime        _RL startTime
723        _RL endTime        _RL endTime
# Line 715  C      --"-"--  Quadratic  ( linear: 1/s Line 754  C      --"-"--  Quadratic  ( linear: 1/s
754        _RL recip_horiVertRatio        _RL recip_horiVertRatio
755        _RL ivdc_kappa        _RL ivdc_kappa
756        _RL Ro_SeaLevel        _RL Ro_SeaLevel
757          _RL sideDragFactor
758        _RL bottomDragLinear        _RL bottomDragLinear
759        _RL bottomDragQuadratic        _RL bottomDragQuadratic
760          _RL nh_Am2
761        _RL tCylIn        _RL tCylIn
762        _RL tCylOut        _RL tCylOut
763    
# Line 762  C Logical flags for selecting packages Line 803  C Logical flags for selecting packages
803        LOGICAL useFLT        LOGICAL useFLT
804        LOGICAL usePTRACERS        LOGICAL usePTRACERS
805        LOGICAL useGCHEM        LOGICAL useGCHEM
806          LOGICAL useRBCS
807        LOGICAL useSBO        LOGICAL useSBO
808        LOGICAL useSEAICE        LOGICAL useSEAICE
809        LOGICAL useBulkForce        LOGICAL useBulkForce
# Line 772  C Logical flags for selecting packages Line 814  C Logical flags for selecting packages
814        LOGICAL useEBM        LOGICAL useEBM
815        LOGICAL useMNC        LOGICAL useMNC
816        LOGICAL useMATRIX        LOGICAL useMATRIX
817          LOGICAL useRunClock
818        COMMON /PARM_PACKAGES/        COMMON /PARM_PACKAGES/
819       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,
820       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
821       &        usePTRACERS, useGCHEM,       &        usePTRACERS, useGCHEM, useRBCS,
822       &        useSBO, useSEAICE, useThSIce, useBulkForce,       &        useSBO, useSEAICE, useThSIce, useBulkForce,
823       &        usefizhi, usegridalt, usediagnostics, useEBM, useMNC,       &        usefizhi, usegridalt, usediagnostics, useEBM, useMNC,
824       &        usePP81, useMY82, useOPPS, useGGL90, useMATRIX       &        usePP81, useMY82, useOPPS, useGGL90, useMATRIX,
825         &        useRunClock
826    
827  CEH3 ;;; Local Variables: ***  CEH3 ;;; Local Variables: ***
828  CEH3 ;;; mode:fortran ***  CEH3 ;;; mode:fortran ***
Legend:
Removed from v.1.157  
changed lines
  Added in v.1.173
  ViewVC Help
Powered by ViewVC 1.1.22