--- MITgcm/model/inc/PARAMS.h 2005/09/16 19:33:05 1.165 +++ MITgcm/model/inc/PARAMS.h 2005/10/24 22:38:49 1.170 @@ -1,4 +1,4 @@ -C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/inc/PARAMS.h,v 1.165 2005/09/16 19:33:05 baylor Exp $ +C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/inc/PARAMS.h,v 1.170 2005/10/24 22:38:49 jmc Exp $ C $Name: $ C @@ -151,9 +151,10 @@ C select_rStar :: option related to r* vertical coordinate C =0 (default) use r coord. ; > 0 use r* C tempAdvScheme :: Temp. Horiz.Advection scheme selector -C tempVertAdvScheme :: Temp. Vert. Advection scheme selector +C tempVertAdvScheme :: Temp. Vert. Advection scheme selector C saltAdvScheme :: Salt. Horiz.advection scheme selector C saltVertAdvScheme :: Salt. Vert. Advection scheme selector +C selectKEscheme :: Kinetic Energy scheme selector (Vector Inv.) C debugLevel :: debug level selector: higher -> more writing COMMON /PARM_I/ @@ -168,6 +169,7 @@ & nonlinFreeSurf, select_rStar, & tempAdvScheme, tempVertAdvScheme, & saltAdvScheme, saltVertAdvScheme, + & selectKEscheme, & debugLevel INTEGER cg2dMaxIters INTEGER cg2dChkResFreq @@ -186,6 +188,7 @@ INTEGER select_rStar INTEGER tempAdvScheme, tempVertAdvScheme INTEGER saltAdvScheme, saltVertAdvScheme + INTEGER selectKEscheme INTEGER debugLevel C @@ -226,11 +229,10 @@ 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) -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 C on irregular grids) -C useAnisotropicViscAGridMax :: Set to true to use Alistair's latest -C anisotropic length scale. It is used only for maximum viscosity -C calculations. Alistair recommends a value of viscA*GridMax=.25 +C useAreaViscLength :: Set to true to use old scaling for viscous +C lengths, e.g., L2=Raz. May be preferable for cube sphere. C useStrainTensionVisc:: Set to true to use Strain-Tension viscous terms C rigidLid :: Set to true to use rigid lid C implicitFreeSurface :: Set to true to use implcit free surface @@ -258,8 +260,10 @@ C momImplVertAdv :: Turns on implicit vertical advection for Momentum C multiDimAdvection :: Flag that enable multi-dimension advection C useMultiDimAdvec :: True if multi-dim advection is used at least once -C forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution -C out off Adams-Bashforth time stepping. +C forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution +C out off Adams-Bashforth time stepping. +C doAB_onGtGs :: if the Adams-Bashforth time stepping is used, always +C apply AB on tracer tendencies (rather than on Tracer) C startFromPickupAB2 :: with AB-3 code, start from an AB-2 pickup C doThetaClimRelax :: Set true if relaxation to temperature C climatology is required. @@ -307,7 +311,8 @@ & tempAdvection, tempForcing, & saltAdvection, saltForcing, & useRealFreshWaterFlux, - & useFullLeith, useAnisotropicViscAGridMax, useStrainTensionVisc, + & useFullLeith, useStrainTensionVisc, + & useAreaViscLength, & rigidLid, implicitFreeSurface, exactConserv, uniformLin_PhiSurf, & momStepping, tempStepping, saltStepping, & metricTerms, usingSphericalPolarMTerms, useNHMTerms, @@ -318,7 +323,8 @@ & useAbsVorticity, upwindShear, & implicitDiffusion, implicitViscosity, & tempImplVertAdv, saltImplVertAdv, momImplVertAdv, - & multiDimAdvection, useMultiDimAdvec, forcing_In_AB, + & multiDimAdvection, useMultiDimAdvec, + & forcing_In_AB, doAB_onGtGs, & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax, & periodicExternalForcing, & fluidIsAir, fluidIsWater, @@ -353,8 +359,8 @@ LOGICAL saltForcing LOGICAL useRealFreshWaterFlux LOGICAL useFullLeith - LOGICAL useAnisotropicViscAGridMax LOGICAL useStrainTensionVisc + LOGICAL useAreaViscLength LOGICAL rigidLid LOGICAL implicitFreeSurface LOGICAL exactConserv @@ -383,6 +389,7 @@ LOGICAL multiDimAdvection LOGICAL useMultiDimAdvec LOGICAL forcing_In_AB + LOGICAL doAB_onGtGs LOGICAL doThetaClimRelax LOGICAL doSaltClimRelax LOGICAL doTr1ClimRelax @@ -476,15 +483,19 @@ C viscA4Z :: Biharmonic viscosity coeff. for mixing of momentum laterally C (act on Vorticity part) ( m^4/s ) C viscC2leith :: Leith non-dimensional viscosity factor (grad(vort)) -C viscC2leithD :: Modified Leith non-dimensional viscosity factor (grad(div)) +C viscC2leithD :: Modified Leith non-dimensional visc. factor (grad(div)) C viscC2smag :: Smagorinsky non-dimensional viscosity factor (harmonic) C viscC4smag :: Smagorinsky non-dimensional viscosity factor (biharmonic) C viscAhMax :: Maximum eddy viscosity coeff. for mixing of -C momentum laterally ( m^2/s ) +C momentum laterally ( m^2/s ) +C viscAhReMax :: Maximum gridscale Reynolds number for eddy viscosity +C coeff. for mixing of momentum laterally (non-dim) C viscAhGridMax:: maximum and minimum harmonic viscosity coefficients ... -C viscAhGridMin:: in terms of non-dimensional grid-size dependent viscosity +C viscAhGridMin:: in terms of non-dimensional grid-size dependent visc. C viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of C momentum laterally ( m^4/s ) +C viscA4ReMax :: Maximum Gridscale Reynolds number for +C biharmonic viscosity coeff. momentum laterally (non-dim) C viscAhGrid:: non-dimensional grid-size dependent viscosity C viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity C viscA4GridMax:: maximum and minimum biharmonic viscosity coefficients ... @@ -578,8 +589,12 @@ C ( g*rho if horiz in m and vertical in Pa ). C Ro_SeaLevel :: standard position of Sea-Level in "R" coordinate, used as C starting value (k=1) for vertical coordinate (rf(1)=Ro_SeaLevel) +C sideDragFactor :: side-drag scaling factor (used only if no_slip_sides) +C (default=2: full drag ; =1: gives half-slip BC) C bottomDragLinear :: Drag coefficient built in to core dynamics -C --"-"-- Quadratic ( linear: 1/s, quadratic: 1/m ) +C bottomDragQuadratic ( linear: 1/s, quadratic: 1/m ) +C nh_Am2 :: scales the non-hydrostatic terms and changes internal scales +C (i.e. allows convection at different Rayleigh numbers) COMMON /PARM_R/ cg2dTargetResidual, cg2dTargetResWunit, & cg2dpcOffDFac, cg3dTargetResidual, & delR, delRc, delX, delY, @@ -593,6 +608,7 @@ & viscAhD, viscAhZ, viscA4D, viscA4Z, & viscA4, viscA4W, & viscA4Max, viscA4Grid, viscA4GridMax, viscA4GridMin, + & viscAhRemax, viscA4Remax, & viscC4leith, viscC4leithD, viscAr, & diffKhT, diffK4T, diffKrNrT, & diffKhS, diffK4S, diffKrNrS, @@ -615,10 +631,9 @@ & hFacMinDr, hFacMinDp, & horiVertRatio, recip_horiVertRatio, & ivdc_kappa, Ro_SeaLevel, - & bottomDragLinear,bottomDragQuadratic,nh_Am2, + & sideDragFactor, bottomDragLinear, bottomDragQuadratic, nh_Am2, & tCylIn, tCylOut - _RL nh_Am2 _RL cg2dTargetResidual _RL cg2dTargetResWunit _RL cg3dTargetResidual @@ -653,6 +668,7 @@ _RL viscAhD _RL viscAhZ _RL viscAhMax + _RL viscAhReMax _RL viscAhGrid _RL viscAhGridMax _RL viscAhGridMin @@ -666,6 +682,7 @@ _RL viscA4D _RL viscA4Z _RL viscA4Max + _RL viscA4ReMax _RL viscA4Grid, viscA4GridMax, viscA4GridMin _RL viscC4leith _RL viscC4leithD @@ -729,8 +746,10 @@ _RL recip_horiVertRatio _RL ivdc_kappa _RL Ro_SeaLevel + _RL sideDragFactor _RL bottomDragLinear _RL bottomDragQuadratic + _RL nh_Am2 _RL tCylIn _RL tCylOut