--- MITgcm_contrib/submesoscale/code/gmredi_readparms.F	2008/05/30 21:51:24	1.1
+++ MITgcm_contrib/submesoscale/code/gmredi_readparms.F	2008/05/30 22:13:42	1.2
@@ -1,4 +1,4 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/submesoscale/code/gmredi_readparms.F,v 1.1 2008/05/30 21:51:24 dimitri Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm_contrib/submesoscale/code/gmredi_readparms.F,v 1.2 2008/05/30 22:13:42 dimitri Exp $
 C $Name:  $
 
 #include "GMREDI_OPTIONS.h"
@@ -45,6 +45,7 @@
 C     GM_slopeSqCutoff :: slope^2 cut-off value
       NAMELIST /GM_PARM01/
      &          GM_AdvForm, GM_AdvSeparate,
+     &          GM_ExtraDiag,
      &          GM_isopycK,
      &          GM_background_K,
      &          GM_taper_scheme,
@@ -57,7 +58,10 @@
      &          GM_Visbeck_maxval_K,
      &          GM_Scrit,
      &          GM_Sd,
-     &          GM_MNC
+     &          GM_MNC,
+     &          GM_SM_Ce,
+     &          GM_SM_Lf,
+     &          GM_SM_tau
 c     &          GM_dumpFreq,
 c     &          GM_taveFreq
 
@@ -69,6 +73,7 @@
 
 C--   Default values GM/Redi
       GM_AdvForm          = .FALSE.
+      GM_ExtraDiag        = .FALSE.
       GM_AdvSeparate      = .FALSE.
       GM_isopycK          = -999.
       GM_background_K     = 0. _d 0
@@ -80,6 +85,18 @@
       GM_Scrit            = 0.004 _d 0
       GM_Sd               = 0.001 _d 0
       GM_MNC              = useMNC
+CBFK Submesoscale Paramters
+C     GM_SM_Ce is the efficiency factor.  Nondimensional, typically 0.06
+C     Will be converted to a dimensional factor later
+C     Set to zero to kill Submeso at runtime.
+      GM_SM_Ce            = 0.06  _d 0
+C     GM_SM_Lf is the minimum Front Width scale.  Typically 1-10km
+      GM_SM_Lf            = 5.00  _d 3
+C     GM_SM_Lmax is the maximum upscale.  Typically 100km 
+      GM_SM_Lmax            = 1.00  _d 5
+C     GM_SM_tau is the momentum mixing timescale: 1/d to 1/week
+      GM_SM_tau           = 8.64  _d 4
+C     Will be converted to an inverse squared timescale later
 
 C--   Default values GM/Redi I/O control
 c     GM_dumpFreq         = -1.
@@ -121,7 +138,8 @@
       ELSE
         GM_skewflx = 1.
         GM_advect  = 0.
-        GM_ExtraDiag = GM_isopycK.NE.GM_background_K
+        GM_ExtraDiag = (GM_isopycK.NE.GM_background_K)
+     &      .or.(GM_SM_Ce.gt.0 _d 0)
       ENDIF
 
 C     Make sure that we locally honor the global MNC on/off flag