--- MITgcm/model/inc/PARAMS.h	2003/12/15 23:00:43	1.99
+++ MITgcm/model/inc/PARAMS.h	2004/02/07 16:27:13	1.102
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/inc/PARAMS.h,v 1.99 2003/12/15 23:00:43 molod Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/model/inc/PARAMS.h,v 1.102 2004/02/07 16:27:13 adcroft Exp $
 C $Name:  $
 C
 CBOP
@@ -215,6 +215,9 @@
 C     useJamartWetPoints :: Use wet-point method for Coriolis (Jamart and Ozer, 1986)
 C     implicitDiffusion :: Turns implicit vertical diffusion on
 C     implicitViscosity :: Turns implicit vertical viscosity on
+C     tempImplVertAdv :: Turns on implicit vertical advection for Temperature
+C     saltImplVertAdv :: Turns on implicit vertical advection for Salinity
+C     momImplVertAdv  :: Turns on implicit vertical advection for Momentum
 C     multiDimAdvection :: Flag that enable multi-dimension advection
 C     forcing_In_AB :: if False, put forcing (Temp,Salt,Tracers) contribution
 C                      out off Adams-Bashforth time stepping.
@@ -256,6 +259,7 @@
      & useCDscheme,
      & useEnergyConservingCoriolis, useJamartWetPoints,
      & implicitDiffusion, implicitViscosity,
+     & tempImplVertAdv, saltImplVertAdv, momImplVertAdv,
      & multiDimAdvection, forcing_In_AB,
      & doThetaClimRelax, doSaltClimRelax, doTr1ClimRelax, 
      & periodicExternalForcing, 
@@ -302,6 +306,9 @@
       LOGICAL useJamartWetPoints
       LOGICAL implicitDiffusion
       LOGICAL implicitViscosity
+      LOGICAL tempImplVertAdv 
+      LOGICAL saltImplVertAdv 
+      LOGICAL momImplVertAdv
       LOGICAL multiDimAdvection
       LOGICAL forcing_In_AB
       LOGICAL doThetaClimRelax
@@ -373,6 +380,12 @@
 C                 momentum vertically ( units of r^2/s )
 C     viscA4    :: Biharmonic viscosity coeff. for mixing of
 C                 momentum laterally ( m^4/s )
+C     viscAhMax :: Maximum eddy viscosity coeff. for mixing of
+C                 momentum laterally ( m^2/s )
+C     viscA4Max :: Maximum biharmonic viscosity coeff. for mixing of
+C                 momentum laterally ( m^4/s )
+C     viscAhGrid:: non-dimensional grid-size dependent viscosity
+C     viscA4Grid:: non-dimensional grid-size dependent bi-harmonic viscosity
 C     diffKhT   :: Laplacian diffusion coeff. for mixing of
 C                 heat laterally ( m^2/s )
 C     diffKrT   :: Laplacian diffusion coeff. for mixing of
@@ -464,7 +477,9 @@
      & abeps, startTime, 
      & phiMin, thetaMin, rSphere, recip_RSphere, f0, beta,
      & fCori, fCoriG, fCoriCos,
-     & viscAh,  viscA4,  viscAr, viscAstrain, viscAtension,
+     & viscAh,  viscAhMax, viscAhGrid,
+     & viscA4,  viscA4Max, viscA4Grid,
+     & viscAr, viscAstrain, viscAtension,
      & diffKhT, diffK4T, diffKrT,
      & diffKhS, diffK4S, diffKrS,
      & delT, tauCD, rCD, freeSurfFac, implicSurfPress, implicDiv2Dflow,
@@ -516,10 +531,14 @@
       _RL hFacSup
       _RL beta
       _RL viscAh
+      _RL viscAhMax
+      _RL viscAhGrid
       _RL viscAstrain
       _RL viscAtension
       _RL viscAr
       _RL viscA4 
+      _RL viscA4Max
+      _RL viscA4Grid
       _RL diffKhT 
       _RL diffKrT
       _RL diffK4T 
@@ -625,8 +644,9 @@
       LOGICAL useThSIce
       LOGICAL usefizhi
       LOGICAL usegridalt
+      LOGICAL usediagnostics
       COMMON /PARM_PACKAGES/
      &        useKPP, useGMRedi, useOBCS, useAIM, useLand, useECCO, 
      &        useSHAP_FILT, useZONAL_FILT, useGrdchk, useFLT,
-     &        usePTRACERS,  useSBO, useSEAICE, 
-     &        useThSIce, useBulkForce, usefizhi,  usegridalt
+     &        usePTRACERS,  useSBO, useSEAICE, useThSIce, useBulkForce, 
+     &        usefizhi,  usegridalt, usediagnostics