/[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.155 by edhill, Sun Apr 17 00:25:45 2005 UTC revision 1.178 by jmc, Fri Feb 10 22:54:11 2006 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)
59  C                   (mainly used for the atmosphere = ground height).  C                   (mainly used for the atmosphere = ground height).
60    C     shelfIceFile:: File containing the topography of the shelfice draught
61    C                    (unit=m)
62  C     hydrogThetaFile :: File containing initial hydrographic data for potential  C     hydrogThetaFile :: File containing initial hydrographic data for potential
63  C                       temperature.  C                       temperature.
64  C     hydrogSaltFile  :: File containing initial hydrographic data for salinity.  C     hydrogSaltFile  :: File containing initial hydrographic data for salinity.
# Line 78  C     eosType         :: choose the equa Line 84  C     eosType         :: choose the equa
84  C                        LINEAR, POLY3, UNESCO, JMD95Z, JMD95P, MDJWF, IDEALGAS  C                        LINEAR, POLY3, UNESCO, JMD95Z, JMD95P, MDJWF, IDEALGAS
85  C     the_run_name    :: string identifying the name of the model "run"  C     the_run_name    :: string identifying the name of the model "run"
86        COMMON /PARM_C/ checkPtSuff,        COMMON /PARM_C/ checkPtSuff,
87       &                bathyFile, topoFile,       &                delXFile, delYFile, horizGridFile,
88         &                bathyFile, topoFile, shelfIceFile,
89       &                hydrogThetaFile, hydrogSaltFile,       &                hydrogThetaFile, hydrogSaltFile,
90       &                zonalWindFile, meridWindFile, thetaClimFile,       &                zonalWindFile, meridWindFile, thetaClimFile,
91       &                saltClimFile, buoyancyRelation,       &                saltClimFile, buoyancyRelation,
# Line 92  C     the_run_name    :: string identify Line 99  C     the_run_name    :: string identify
99       &                mdsioLocalDir,       &                mdsioLocalDir,
100       &                the_run_name       &                the_run_name
101        CHARACTER*(5) checkPtSuff(maxNoChkptLev)        CHARACTER*(5) checkPtSuff(maxNoChkptLev)
102        CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile        CHARACTER*(MAX_LEN_FNAM) delXFile
103          CHARACTER*(MAX_LEN_FNAM) delYFile
104          CHARACTER*(MAX_LEN_FNAM) horizGridFile
105          CHARACTER*(MAX_LEN_FNAM) bathyFile, topoFile, shelfIceFile
106        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile        CHARACTER*(MAX_LEN_FNAM) hydrogThetaFile
107        CHARACTER*(MAX_LEN_FNAM) hydrogSaltFile        CHARACTER*(MAX_LEN_FNAM) hydrogSaltFile
108        CHARACTER*(MAX_LEN_FNAM) zonalWindFile        CHARACTER*(MAX_LEN_FNAM) zonalWindFile
# Line 143  C                           =0 Linear fr Line 153  C                           =0 Linear fr
153  C     select_rStar        :: option related to r* vertical coordinate  C     select_rStar        :: option related to r* vertical coordinate
154  C                           =0 (default) use r coord. ; > 0 use r*  C                           =0 (default) use r coord. ; > 0 use r*
155  C     tempAdvScheme       :: Temp. Horiz.Advection scheme selector  C     tempAdvScheme       :: Temp. Horiz.Advection scheme selector
156  C     tempVertAdvScheme   :: Temp. Vert. Advection scheme selector  C     tempVertAdvScheme   :: Temp. Vert. Advection scheme selector
157  C     saltAdvScheme       :: Salt. Horiz.advection scheme selector  C     saltAdvScheme       :: Salt. Horiz.advection scheme selector
158  C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector  C     saltVertAdvScheme   :: Salt. Vert. Advection scheme selector
159    C     selectKEscheme      :: Kinetic Energy scheme selector (Vector Inv.)
160  C     debugLevel          :: debug level selector: higher -> more writing  C     debugLevel          :: debug level selector: higher -> more writing
161    
162        COMMON /PARM_I/        COMMON /PARM_I/
# Line 160  C     debugLevel          :: debug level Line 171  C     debugLevel          :: debug level
171       &        nonlinFreeSurf, select_rStar,       &        nonlinFreeSurf, select_rStar,
172       &        tempAdvScheme, tempVertAdvScheme,       &        tempAdvScheme, tempVertAdvScheme,
173       &        saltAdvScheme, saltVertAdvScheme,       &        saltAdvScheme, saltVertAdvScheme,
174         &        selectKEscheme,
175       &        debugLevel       &        debugLevel
176        INTEGER cg2dMaxIters        INTEGER cg2dMaxIters
177        INTEGER cg2dChkResFreq        INTEGER cg2dChkResFreq
# Line 178  C     debugLevel          :: debug level Line 190  C     debugLevel          :: debug level
190        INTEGER select_rStar        INTEGER select_rStar
191        INTEGER tempAdvScheme, tempVertAdvScheme        INTEGER tempAdvScheme, tempVertAdvScheme
192        INTEGER saltAdvScheme, saltVertAdvScheme        INTEGER saltAdvScheme, saltVertAdvScheme
193          INTEGER selectKEscheme
194        INTEGER debugLevel        INTEGER debugLevel
195    
196  C  C
# Line 196  C                               spherica Line 209  C                               spherica
209  C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical  C     usingCylindricalGrid :: If TRUE grid generation will be Cylindrical
210  C     no_slip_sides :: Impose "no-slip" at lateral boundaries.  C     no_slip_sides :: Impose "no-slip" at lateral boundaries.
211  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  
212  C     momViscosity  :: Flag which turns momentum friction terms on and off.  C     momViscosity  :: Flag which turns momentum friction terms on and off.
213  C     momAdvection  :: Flag which turns advection of momentum on and off.  C     momAdvection  :: Flag which turns advection of momentum on and off.
214  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 230  C                     and off.
230  C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux  C     useRealFreshWaterFlux :: if True (=Natural BCS), treats P+R-E flux
231  C                         as a real Fresh Water (=> changes the Sea Level)  C                         as a real Fresh Water (=> changes the Sea Level)
232  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)
233  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
234  C                       on irregular grids)  C                       on irregular grids)
235  C     useAnisotropicViscAGridMax :: Set to true to use Alistair's latest  C     useAreaViscLength :: Set to true to use old scaling for viscous
236  C              anisotropic length scale.  It is used only for maximum viscosity  C              lengths, e.g., L2=Raz.  May be preferable for cube sphere.
237  C              calculations.  Alistair recommends a value of viscA*GridMax=.25  C     useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms
238  C     rigidLid            :: Set to true to use rigid lid  C     rigidLid            :: Set to true to use rigid lid
239  C     implicitFreeSurface :: Set to true to use implcit free surface  C     implicitFreeSurface :: Set to true to use implicit free surface
240  C     exactConserv        :: Set to true to conserve exactly the total Volume  C     exactConserv        :: Set to true to conserve exactly the total Volume
241  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
242  C                           linear relation Phi_surf = Bo_surf*eta  C                           linear relation Phi_surf = Bo_surf*eta
243    C     use3Dsolver   :: set to true to use 3-D pressure solver
244    C     implicitIntGravWave :: treat Internal Gravity Wave implicitly
245    C     staggerTimeStep :: enable a Stagger time stepping T,S Rho then U,V
246  C     momStepping   :: Turns momentum equation time-stepping off  C     momStepping   :: Turns momentum equation time-stepping off
247  C     tempStepping  :: Turns temperature equation time-stepping off  C     tempStepping  :: Turns temperature equation time-stepping off
248  C     saltStepping  :: Turns salinity equation time-stepping off  C     saltStepping  :: Turns salinity equation time-stepping off
# Line 239  C     useJamartWetPoints :: Use wet-poin Line 254  C     useJamartWetPoints :: Use wet-poin
254  C     useJamartMomAdv :: Use wet-point method for V.I. non-linear term  C     useJamartMomAdv :: Use wet-point method for V.I. non-linear term
255  C     SadournyCoriolis :: use the enstrophy conserving scheme by Sadourny  C     SadournyCoriolis :: use the enstrophy conserving scheme by Sadourny
256  C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term  C     upwindVorticity :: bias interpolation of vorticity in the Coriolis term
257  C     highOrderVorticity :: use 3rd/4th order interp. of vorticity in Coriolis  C     highOrderVorticity :: use 3rd/4th order interp. of vorticity (V.I., advection)
258    C     upwindShear        :: use 1rst order upwind interp. (V.I., vertical advection)
259  C     useAbsVorticity :: work with f+zeta in Coriolis terms  C     useAbsVorticity :: work with f+zeta in Coriolis terms
260  C     implicitDiffusion :: Turns implicit vertical diffusion on  C     implicitDiffusion :: Turns implicit vertical diffusion on
261  C     implicitViscosity :: Turns implicit vertical viscosity on  C     implicitViscosity :: Turns implicit vertical viscosity on
# Line 248  C     saltImplVertAdv :: Turns on implic Line 264  C     saltImplVertAdv :: Turns on implic
264  C     momImplVertAdv  :: Turns on implicit vertical advection for Momentum  C     momImplVertAdv  :: Turns on implicit vertical advection for Momentum
265  C     multiDimAdvection :: Flag that enable multi-dimension advection  C     multiDimAdvection :: Flag that enable multi-dimension advection
266  C     useMultiDimAdvec  :: True if multi-dim advection is used at least once  C     useMultiDimAdvec  :: True if multi-dim advection is used at least once
267  C     forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution  C     forcing_In_AB     :: if False, put forcing (Temp,Salt,Tracers) contribution
268  C                      out off Adams-Bashforth time stepping.  C                          out off Adams-Bashforth time stepping.
269    C     doAB_onGtGs       :: if the Adams-Bashforth time stepping is used, always
270    C                          apply AB on tracer tendencies (rather than on Tracer)
271  C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup  C     startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup
272  C     doThetaClimRelax :: Set true if relaxation to temperature  C     doThetaClimRelax :: Set true if relaxation to temperature
273  C                        climatology is required.  C                        climatology is required.
# Line 276  C                       useSingleCpuIO i Line 294  C                       useSingleCpuIO i
294  C                       outputs from master mpi process only.  C                       outputs from master mpi process only.
295  C     allowFreezing  :: Allows surface water to freeze and form ice  C     allowFreezing  :: Allows surface water to freeze and form ice
296  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)  
297  C     pickup_write_mdsio :: use mdsio to write pickups  C     pickup_write_mdsio :: use mdsio to write pickups
298  C     pickup_read_mdsio  :: use mdsio to read  pickups  C     pickup_read_mdsio  :: use mdsio to read  pickups
299  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 304  C     calendarDumps :: When set, approxi
304  C                      for parameters chkPtFreq, pChkPtFreq, taveFreq, SEAICE_taveFreq,  C                      for parameters chkPtFreq, pChkPtFreq, taveFreq, SEAICE_taveFreq,
305  C                      KPP_taveFreq, and freq in pkg/diagnostics are converted to exact  C                      KPP_taveFreq, and freq in pkg/diagnostics are converted to exact
306  C                      calendar months and years.  Requires pkg/cal.  C                      calendar months and years.  Requires pkg/cal.
307    C     dumpInitAndLast :: dumps model state to files at Initial (nIter0)
308    C                        & Last iteration, in addition multiple of dumpFreq iter.
309    C     balanceEmPmR    :: substract global mean of EmPmR at every time step
310    C     balanceQnet     :: substract global mean of Qnet at every time step
311    C     balancePrintMean:: print substracted global means to STDOUT
312    
313        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,        COMMON /PARM_L/ usingCartesianGrid, usingSphericalPolarGrid,
314       & usingCurvilinearGrid, usingCylindricalGrid,       & usingCurvilinearGrid, usingCylindricalGrid,
315       & no_slip_sides,no_slip_bottom,       & no_slip_sides,no_slip_bottom,
      & staggerTimeStep,  
316       & momViscosity, momAdvection, momForcing, useCoriolis,       & momViscosity, momAdvection, momForcing, useCoriolis,
317       & momPressureForcing, vectorInvariantMomentum,       & momPressureForcing, vectorInvariantMomentum,
318       & tempAdvection, tempForcing,       & tempAdvection, tempForcing,
319       & saltAdvection, saltForcing,       & saltAdvection, saltForcing,
320       & useRealFreshWaterFlux,       & useRealFreshWaterFlux,
321       & useFullLeith, useAnisotropicViscAGridMax,       & useFullLeith, useStrainTensionVisc,
322         & useAreaViscLength,
323       & rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf,       & rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf,
324         & use3Dsolver, implicitIntGravWave, staggerTimeStep,
325       & momStepping, tempStepping, saltStepping,       & momStepping, tempStepping, saltStepping,
326       & metricTerms, usingSphericalPolarMTerms, useNHMTerms,       & metricTerms, usingSphericalPolarMTerms, useNHMTerms,
327       & useConstantF, useBetaPlaneF, useSphereF,       & useConstantF, useBetaPlaneF, useSphereF,
328       & useCDscheme,       & useCDscheme,
329       & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,       & useEnergyConservingCoriolis, useJamartWetPoints, useJamartMomAdv,
330       & SadournyCoriolis, upwindVorticity, highOrderVorticity,       & SadournyCoriolis, upwindVorticity, highOrderVorticity,
331       & useAbsVorticity,       & useAbsVorticity, upwindShear,
332       & implicitDiffusion, implicitViscosity,       & implicitDiffusion, implicitViscosity,
333       & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,       & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,
334       & multiDimAdvection, useMultiDimAdvec, forcing_In_AB,       & multiDimAdvection, useMultiDimAdvec,
335         & forcing_In_AB, doAB_onGtGs,
336       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,       & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax,
337       & periodicExternalForcing,       & periodicExternalForcing,
338       & fluidIsAir, fluidIsWater,       & fluidIsAir, fluidIsWater,
339       & usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,       & usingPCoords, usingZCoords, useDynP_inEos_Zc, setCenterDr,
340       & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,       & nonHydrostatic, quasiHydrostatic, globalFiles, useSingleCpuIO,
341       & allowFreezing, useOldFreezing, groundAtK1,       & allowFreezing, useOldFreezing,
342       & usePickupBeforeC35, usePickupBeforeC54, startFromPickupAB2,       & usePickupBeforeC35, usePickupBeforeC54, startFromPickupAB2,
343       & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,       & pickup_read_mdsio, pickup_write_mdsio, pickup_write_immed,
344       & timeave_mdsio, snapshot_mdsio, monitor_stdio,       & timeave_mdsio, snapshot_mdsio, monitor_stdio,
345       & outputTypesInclusive, debugMode,       & outputTypesInclusive, dumpInitAndLast, debugMode,
346       & inAdMode, inAdTrue, inAdFalse, inAdExact,       & inAdMode, inAdTrue, inAdFalse, inAdExact,
347       & calendarDumps       & calendarDumps,
348         & balanceEmPmR, balanceQnet, balancePrintMean
349    
350        LOGICAL usingCartesianGrid        LOGICAL usingCartesianGrid
351        LOGICAL usingSphericalPolarGrid        LOGICAL usingSphericalPolarGrid
# Line 329  C                      calendar months a Line 355  C                      calendar months a
355        LOGICAL useNHMTerms        LOGICAL useNHMTerms
356        LOGICAL no_slip_sides        LOGICAL no_slip_sides
357        LOGICAL no_slip_bottom        LOGICAL no_slip_bottom
       LOGICAL staggerTimeStep  
358        LOGICAL momViscosity        LOGICAL momViscosity
359        LOGICAL momAdvection        LOGICAL momAdvection
360        LOGICAL momForcing        LOGICAL momForcing
# Line 341  C                      calendar months a Line 366  C                      calendar months a
366        LOGICAL saltAdvection        LOGICAL saltAdvection
367        LOGICAL saltForcing        LOGICAL saltForcing
368        LOGICAL useRealFreshWaterFlux        LOGICAL useRealFreshWaterFlux
369        LOGICAL useFullLeith, useAnisotropicViscAGridMax        LOGICAL useFullLeith
370          LOGICAL useStrainTensionVisc
371          LOGICAL useAreaViscLength
372        LOGICAL rigidLid        LOGICAL rigidLid
373        LOGICAL implicitFreeSurface        LOGICAL implicitFreeSurface
374        LOGICAL exactConserv        LOGICAL exactConserv
375        LOGICAL uniformLin_PhiSurf        LOGICAL uniformLin_PhiSurf
376          LOGICAL use3Dsolver
377          LOGICAL implicitIntGravWave
378          LOGICAL staggerTimeStep
379        LOGICAL momStepping        LOGICAL momStepping
380        LOGICAL tempStepping        LOGICAL tempStepping
381        LOGICAL saltStepping        LOGICAL saltStepping
# Line 361  C                      calendar months a Line 391  C                      calendar months a
391        LOGICAL upwindVorticity        LOGICAL upwindVorticity
392        LOGICAL highOrderVorticity        LOGICAL highOrderVorticity
393        LOGICAL useAbsVorticity        LOGICAL useAbsVorticity
394          LOGICAL upwindShear
395        LOGICAL implicitDiffusion        LOGICAL implicitDiffusion
396        LOGICAL implicitViscosity        LOGICAL implicitViscosity
397        LOGICAL tempImplVertAdv        LOGICAL tempImplVertAdv
# Line 369  C                      calendar months a Line 400  C                      calendar months a
400        LOGICAL multiDimAdvection        LOGICAL multiDimAdvection
401        LOGICAL useMultiDimAdvec        LOGICAL useMultiDimAdvec
402        LOGICAL forcing_In_AB        LOGICAL forcing_In_AB
403          LOGICAL doAB_onGtGs
404        LOGICAL doThetaClimRelax        LOGICAL doThetaClimRelax
405        LOGICAL doSaltClimRelax        LOGICAL doSaltClimRelax
406        LOGICAL doTr1ClimRelax        LOGICAL doTr1ClimRelax
# Line 385  C                      calendar months a Line 417  C                      calendar months a
417        LOGICAL useSingleCpuIO        LOGICAL useSingleCpuIO
418        LOGICAL allowFreezing        LOGICAL allowFreezing
419        LOGICAL useOldFreezing        LOGICAL useOldFreezing
       LOGICAL groundAtK1  
420        LOGICAL usePickupBeforeC35        LOGICAL usePickupBeforeC35
421        LOGICAL usePickupBeforeC54        LOGICAL usePickupBeforeC54
422        LOGICAL startFromPickupAB2        LOGICAL startFromPickupAB2
423          LOGICAL dumpInitAndLast
424        LOGICAL debugMode        LOGICAL debugMode
425        LOGICAL pickup_read_mdsio, pickup_write_mdsio        LOGICAL pickup_read_mdsio, pickup_write_mdsio
426        LOGICAL pickup_write_immed        LOGICAL pickup_write_immed
# Line 397  C                      calendar months a Line 429  C                      calendar months a
429        LOGICAL inAdMode, inAdTrue, inAdFalse, inAdExact        LOGICAL inAdMode, inAdTrue, inAdFalse, inAdExact
430        LOGICAL calendarDumps        LOGICAL calendarDumps
431    
432          LOGICAL balanceEmPmR
433          LOGICAL balanceQnet
434          LOGICAL balancePrintMean
435    
436  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.  C--   COMMON /PARM_R/ "Real" valued parameters used by the model.
437  C     cg2dTargetResidual  C     cg2dTargetResidual
438  C          :: Target residual for cg2d solver; no unit (RHS normalisation)  C          :: Target residual for cg2d solver; no unit (RHS normalisation)
# Line 430  C     rhoConst  :: Vertically constant r Line 466  C     rhoConst  :: Vertically constant r
466  C     rhoConstFresh :: Constant reference density for fresh water (rain)  C     rhoConstFresh :: Constant reference density for fresh water (rain)
467  C     tRef      :: reference vertical profile for potential temperature  C     tRef      :: reference vertical profile for potential temperature
468  C     sRef      :: reference vertical profile for salinity/specific humidity  C     sRef      :: reference vertical profile for salinity/specific humidity
469    C     phiRef    :: reference potential (pressure/rho, geopotential) profile
470    C     dBdrRef   :: vertical gradient of reference boyancy  [(m/s/r)^2)]:
471    C               :: z-coord: = N^2_ref = Brunt-Vaissala frequency [s^-2]
472    C               :: p-coord: = -(d.alpha/dp)_ref          [(m^2.s/kg)^2]
473  C     phiMin    :: Latitude of southern most cell face.  C     phiMin    :: Latitude of southern most cell face.
474  C     thetaMin  :: Longitude of western most cell face (this  C     thetaMin  :: Longitude of western most cell face (this
475  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 502  C                  (act on Divergence pa
502  C     viscA4Z   :: Biharmonic viscosity coeff. for mixing of momentum laterally  C     viscA4Z   :: Biharmonic viscosity coeff. for mixing of momentum laterally
503  C                  (act on Vorticity  part) ( m^4/s )  C                  (act on Vorticity  part) ( m^4/s )
504  C     viscC2leith :: Leith non-dimensional viscosity factor (grad(vort))  C     viscC2leith :: Leith non-dimensional viscosity factor (grad(vort))
505  C     viscC2leithD :: Modified Leith non-dimensional viscosity factor (grad(div))  C     viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div))
506  C     viscC2smag  :: Smagorinsky non-dimensional viscosity factor  C     viscC2smag  :: Smagorinsky non-dimensional viscosity factor (harmonic)
507    C     viscC4smag  :: Smagorinsky non-dimensional viscosity factor (biharmonic)
508  C     viscAhMax :: Maximum eddy viscosity coeff. for mixing of  C     viscAhMax :: Maximum eddy viscosity coeff. for mixing of
509  C                 momentum laterally ( m^2/s )  C                   momentum laterally ( m^2/s )
510    C     viscAhReMax :: Maximum gridscale Reynolds number for eddy viscosity
511    C                   coeff. for mixing of momentum laterally (non-dim)
512  C     viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ...  C     viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ...
513  C     viscAhGridMin::  in terms of non-dimensional grid-size dependent viscosity  C     viscAhGridMin::  in terms of non-dimensional grid-size dependent visc.
514  C     viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of  C     viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of
515  C                 momentum laterally ( m^4/s )  C                 momentum laterally ( m^4/s )
516    C     viscA4ReMax :: Maximum Gridscale Reynolds number for
517    C                  biharmonic viscosity coeff. momentum laterally (non-dim)
518  C     viscAhGrid:: non-dimensional grid-size dependent viscosity  C     viscAhGrid:: non-dimensional grid-size dependent viscosity
519  C     viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity  C     viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity
520  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 608  C     recip_horiVertRatio  ( 1 if horiz
608  C                          ( g*rho if horiz in m and vertical in Pa ).  C                          ( g*rho if horiz in m and vertical in Pa ).
609  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
610  C                          starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)  C                          starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel)
611    C     sideDragFactor     :: side-drag scaling factor (used only if no_slip_sides)
612    C                           (default=2: full drag ; =1: gives half-slip BC)
613  C     bottomDragLinear   :: Drag coefficient built in to core dynamics  C     bottomDragLinear   :: Drag coefficient built in to core dynamics
614  C      --"-"--  Quadratic  ( linear: 1/s, quadratic: 1/m )  C     bottomDragQuadratic   ( linear: 1/s, quadratic: 1/m )
615    C     nh_Am2        :: scales the non-hydrostatic terms and changes internal scales
616    C                      (i.e. allows convection at different Rayleigh numbers)
617        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,        COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit,
618       & cg2dpcOffDFac, cg3dTargetResidual,       & cg2dpcOffDFac, cg3dTargetResidual,
619       & delR, delRc, delX, delY,       & delR, delRc, delX, delY,
# Line 574  C      --"-"--  Quadratic  ( linear: 1/s Line 623  C      --"-"--  Quadratic  ( linear: 1/s
623       & viscAh, viscAhW, viscAhMax,       & viscAh, viscAhW, viscAhMax,
624       & viscAhGrid, viscAhGridMax, viscAhGridMin,       & viscAhGrid, viscAhGridMax, viscAhGridMin,
625       & viscC2leith, viscC2leithD,       & viscC2leith, viscC2leithD,
626       & viscC2smag,       & viscC2smag, viscC4smag,
627       & viscAhD, viscAhZ, viscA4D, viscA4Z,       & viscAhD, viscAhZ, viscA4D, viscA4Z,
628       & viscA4, viscA4W,       & viscA4, viscA4W,
629       & viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin,       & viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin,
630       & viscC4leith, viscC4leithD, viscAr, viscAstrain, viscAtension,       & viscAhRemax, viscA4Remax,
631         & viscC4leith, viscC4leithD, viscAr,
632       & diffKhT, diffK4T, diffKrNrT,       & diffKhT, diffK4T, diffKrNrT,
633       & diffKhS, diffK4S, diffKrNrS,       & diffKhS, diffK4S, diffKrNrS,
634       & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,       & diffKrBL79surf, diffKrBL79deep, diffKrBL79scl, diffKrBL79Ho,
# Line 586  C      --"-"--  Quadratic  ( linear: 1/s Line 636  C      --"-"--  Quadratic  ( linear: 1/s
636       & hFacMin, hFacMinDz, hFacInf, hFacSup,       & hFacMin, hFacMinDz, hFacInf, hFacSup,
637       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,       & gravity, recip_Gravity, gBaro, rhonil, recip_rhonil,
638       & recip_rhoConst, rhoConst,       & recip_rhoConst, rhoConst,
639       & rhoConstFresh, convertEmP2rUnit, tRef, sRef,       & rhoConstFresh, convertEmP2rUnit, tRef, sRef, phiRef, dBdrRef,
640       & baseTime, startTime, endTime,       & baseTime, startTime, endTime,
641       & chkPtFreq, pchkPtFreq, dumpFreq, adjDumpFreq,       & chkPtFreq, pchkPtFreq, dumpFreq, adjDumpFreq,
642       & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq,       & diagFreq, taveFreq, tave_lastIter, monitorFreq, adjMonitorFreq,
643       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,       & afFacMom, vfFacMom, pfFacMom, cfFacMom, foFacMom, mtFacMom,
644       & cosPower, cAdjFreq, omega, rotationPeriod,       & cosPower, cAdjFreq, omega, rotationPeriod,
645       & tauThetaClimRelax,       & tauThetaClimRelax, tauThetaClimRelax3Dim,
646       & tauSaltClimRelax,       & tauSaltClimRelax, tauSaltClimRelax3Dim,
647       & tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax,       & tauTr1ClimRelax, lambdaTr1ClimRelax, latBandClimRelax,
648       & externForcingCycle, externForcingPeriod,       & externForcingCycle, externForcingPeriod,
649       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn,       & convertFW2Salt, temp_EvPrRn, salt_EvPrRn,
650       & hFacMinDr, hFacMinDp,       & hFacMinDr, hFacMinDp,
651       & horiVertRatio, recip_horiVertRatio,       & horiVertRatio, recip_horiVertRatio,
652       & ivdc_kappa, Ro_SeaLevel,       & ivdc_kappa, Ro_SeaLevel,
653       & bottomDragLinear,bottomDragQuadratic,       & sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2,
654       & tCylIn, tCylOut       & tCylIn, tCylOut
655    
656        _RL cg2dTargetResidual        _RL cg2dTargetResidual
# Line 637  C      --"-"--  Quadratic  ( linear: 1/s Line 687  C      --"-"--  Quadratic  ( linear: 1/s
687        _RL viscAhD        _RL viscAhD
688        _RL viscAhZ        _RL viscAhZ
689        _RL viscAhMax        _RL viscAhMax
690          _RL viscAhReMax
691        _RL viscAhGrid        _RL viscAhGrid
692        _RL viscAhGridMax        _RL viscAhGridMax
693        _RL viscAhGridMin        _RL viscAhGridMin
694        _RL viscC2leith        _RL viscC2leith
695        _RL viscC2leithD        _RL viscC2leithD
696        _RL viscC2smag        _RL viscC2smag
697        _RL viscAstrain        _RL viscC4smag
       _RL viscAtension  
698        _RL viscAr        _RL viscAr
699        _RL viscA4        _RL viscA4
700        _RL viscA4W        _RL viscA4W
701        _RL viscA4D        _RL viscA4D
702        _RL viscA4Z        _RL viscA4Z
703        _RL viscA4Max        _RL viscA4Max
704          _RL viscA4ReMax
705        _RL viscA4Grid, viscA4GridMax, viscA4GridMin        _RL viscA4Grid, viscA4GridMax, viscA4GridMin
706        _RL viscC4leith        _RL viscC4leith
707        _RL viscC4leithD        _RL viscC4leithD
# Line 678  C      --"-"--  Quadratic  ( linear: 1/s Line 729  C      --"-"--  Quadratic  ( linear: 1/s
729        _RL convertEmP2rUnit        _RL convertEmP2rUnit
730        _RL tRef(Nr)        _RL tRef(Nr)
731        _RL sRef(Nr)        _RL sRef(Nr)
732          _RL phiRef(2*Nr+1)
733          _RL dBdrRef(Nr)
734        _RL baseTime        _RL baseTime
735        _RL startTime        _RL startTime
736        _RL endTime        _RL endTime
# Line 702  C      --"-"--  Quadratic  ( linear: 1/s Line 755  C      --"-"--  Quadratic  ( linear: 1/s
755        _RL rotationPeriod        _RL rotationPeriod
756        _RL tauThetaClimRelax        _RL tauThetaClimRelax
757        _RL tauSaltClimRelax        _RL tauSaltClimRelax
758          _RL tauThetaClimRelax3Dim
759          _RL tauSaltClimRelax3Dim
760        _RL tauTr1ClimRelax        _RL tauTr1ClimRelax
761        _RL lambdaTr1ClimRelax        _RL lambdaTr1ClimRelax
762        _RL latBandClimRelax        _RL latBandClimRelax
# Line 714  C      --"-"--  Quadratic  ( linear: 1/s Line 769  C      --"-"--  Quadratic  ( linear: 1/s
769        _RL recip_horiVertRatio        _RL recip_horiVertRatio
770        _RL ivdc_kappa        _RL ivdc_kappa
771        _RL Ro_SeaLevel        _RL Ro_SeaLevel
772          _RL sideDragFactor
773        _RL bottomDragLinear        _RL bottomDragLinear
774        _RL bottomDragQuadratic        _RL bottomDragQuadratic
775          _RL nh_Am2
776        _RL tCylIn        _RL tCylIn
777        _RL tCylOut        _RL tCylOut
778    
# Line 761  C Logical flags for selecting packages Line 818  C Logical flags for selecting packages
818        LOGICAL useFLT        LOGICAL useFLT
819        LOGICAL usePTRACERS        LOGICAL usePTRACERS
820        LOGICAL useGCHEM        LOGICAL useGCHEM
821          LOGICAL useRBCS
822        LOGICAL useSBO        LOGICAL useSBO
823        LOGICAL useSEAICE        LOGICAL useSEAICE
824          LOGICAL useShelfIce
825        LOGICAL useBulkForce        LOGICAL useBulkForce
826        LOGICAL useThSIce        LOGICAL useThSIce
827        LOGICAL usefizhi        LOGICAL usefizhi
# Line 770  C Logical flags for selecting packages Line 829  C Logical flags for selecting packages
829        LOGICAL usediagnostics        LOGICAL usediagnostics
830        LOGICAL useEBM        LOGICAL useEBM
831        LOGICAL useMNC        LOGICAL useMNC
832          LOGICAL useMATRIX
833          LOGICAL useRunClock
834          LOGICAL useEMBED_FILES
835        COMMON /PARM_PACKAGES/        COMMON /PARM_PACKAGES/
836       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,       &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO,
837       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,       &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
838       &        usePTRACERS, useGCHEM,       &        usePTRACERS, useGCHEM, useRBCS,
839       &        useSBO, useSEAICE, useThSIce, useBulkForce,       &        useSBO, useSEAICE, useShelfIce, useThSIce, useBulkForce,
840       &        usefizhi, usegridalt, usediagnostics, useEBM, useMNC,       &        usefizhi, usegridalt, usediagnostics, useEBM, useMNC,
841       &        usePP81, useMY82, useOPPS, useGGL90       &        usePP81, useMY82, useOPPS, useGGL90, useMATRIX,
842         &        useRunClock, useEMBED_FILES
843    
844  CEH3 ;;; Local Variables: ***  CEH3 ;;; Local Variables: ***
845  CEH3 ;;; mode:fortran ***  CEH3 ;;; mode:fortran ***
Legend:
Removed from v.1.155  
changed lines
  Added in v.1.178
  ViewVC Help
Powered by ViewVC 1.1.22