--- MITgcm/verification/MLAdjust/input.1.leith/data	2005/09/19 15:14:53	1.2
+++ MITgcm/verification/MLAdjust/input.1.leith/data	2005/09/27 14:35:30	1.3
@@ -4,28 +4,46 @@
 #
 # Continuous equation parameters
  &PARM01
- viscAr=5.E-3,
-# This is the minimum viscosity used total=viscah+viscleith
-# Often U*deltaX/2 for grid noise (Griffies and Hallberg,2000)
- viscAh=3.E1,
+ viscAr=5.E-5,
+# This is the background viscosity used total=viscah+viscleith
+# viscAh=1.E2,
+# viscA4=1.E2,
+# This is the background viscosity used on vorticity points
+# viscAhZ=1.E2,
+# viscA4Z=1.E2,
+# This is the background viscosity used on divergence points
+# viscAhD=1.E2,
+# viscA4D=1.E2,
+# viscAhd and viscahz are set to viscah if they are unset initially
+#
+# This is a maximum gridscale Reynolds number limiter
+# ViscAhRemax=10,
+# ViscA4Remax=10,
+# This is a minimum autoscaling viscosity
+# ViscAhGridMin=.01,
+# ViscA4GridMin=.01,
 # This is the standard Leith Viscosity, prop to grad(curl v_h)
- viscC2leith=5.,
+# viscC2leith=0.2,
+ viscC4leith=0.2,
 # This is an additional Leith Viscosity, prop to grad(div.v_h)
 # It traps instabilities that have no signal in grad(curl v_h)
-# viscC2leithD=5.,
-# These are the minimum viscosities used with Smagorinsky
-# total=viscastrain+viscsmag, etc.
-# viscAstrain=3.E1,
-# viscAtension=3.E1,
+# viscC2leithD=0.2,
+ viscC4leithD=0.2,
 # This is harmonic Smagorinsky Coefficient
-# 0.15=value from 3-d turbulence theory; 
-# 3-4=stable oceanic values (Griffies and Hallberg, 2000)
-# viscC2Smag=4.,
+# 0.2-0.9=value from 3-d turbulence theory; Smagorinsky (1993)
+# 2.2-4=stable oceanic values (Griffies and Hallberg, 2000)
+# 0.7 c2smag is approximately the same as 0.12 c2leith+c2leithd
+# 1.0 c2smag is approximately the same as 0.25 c2leith+c2leithd
+# viscC2Smag=1.0,
+ viscC4Smag=1.0,
 # This is the max Viscosity used, be it Smag or Leith, and it scales 
 # with grid size and timestep
-# Often 0.5 for stability (Griffies and Hallberg,2000)
- viscAhGridMax=0.5,
+# viscAhGridMax=1.0,
+ viscA4GridMax=1.0,
+# This uses the full grad(vort) and grad(div)
  useFullLeith=.TRUE.,
+# This uses the Strain-Tension form for calculating viscous terms
+# useStrainTensionVisc=.TRUE.,
  no_slip_sides=.FALSE.,
  no_slip_bottom=.FALSE.,
  ivdc_kappa=10.,
@@ -64,13 +82,16 @@
 # Time stepping parameters
  &PARM03
  niter0=36,
- endTime=57600.,
+#endTime=86400.,
+#- for testing, only run for 12 iterations:
+ nTimeSteps=12,
  deltaT=1200.,
  abEps=0.01,
- pchkptFreq=7200.,
- taveFreq=7200.,
- dumpFreq=7200.,
- diagFreq=7200.,
+#- for testing, reduce amount of output by commenting out output frequencies:
+# pchkptFreq=43200.,
+#taveFreq=7200.,
+#dumpFreq=7200.,
+#diagFreq=7200.,
  monitorFreq=1200.,
  cadjfreq=0.,
  /