diff --git a/param_templates/MOM_input.yaml b/param_templates/MOM_input.yaml
index ed05736..0db23cb 100644
--- a/param_templates/MOM_input.yaml
+++ b/param_templates/MOM_input.yaml
@@ -1587,27 +1587,6 @@ Global:
         units: kg/m3
         value:
             $OCN_GRID == "tx2_3v2": 0.02
-    FOX_KEMPER_ML_RESTRAT_COEF:
-        description: |
-            "[nondim] default = 0.0
-            A nondimensional coefficient that is proportional to
-            the ratio of the deformation radius to the dominant
-            lengthscale of the submesoscale mixed layer
-            instabilities, times the minimum of the ratio of the
-            mesoscale eddy kinetic energy to the large-scale
-            geostrophic kinetic energy or 1 plus the square of the
-            grid spacing over the deformation radius, as detailed
-            by Fox-Kemper et al. (2010)
-            The file that specifies the vertical grid for
-            depth-space diagnostics, or blank to disable
-            depth-space output.
-            The number of depth-space levels.  This is determined
-            from the size of the variable zw in the output grid file."
-        datatype: real
-        units: nondim
-        value:
-            $OCN_GRID == "tx2_3v2": 1.0
-            $OCN_GRID == "tx0.25v1": 1.0
     MLE_FRONT_LENGTH:
         description: |
             "[m] default = 0.0
@@ -1620,14 +1599,6 @@ Global:
         value:
             $OCN_GRID == "tx2_3v2": 1000.0
             $OCN_GRID == "tx0.25v1": 500.0
-    MLE_FRONT_LENGTH_FROM_FILE:
-        description: |
-            "[Boolean] default = False
-            If true, the MLE front-length scale is read from a file."
-        datatype: logical
-        units: Boolean
-        value:
-            $OCN_GRID == "tx2_3v2": False
     MLE_MLD_DECAY_TIME:
         description: |
             "[s] default = 0.0
@@ -2639,14 +2610,6 @@ Global:
         units: Boolean
         value:
             $OCN_GRID == "tx2_3v2": False
-    MEKE_RESTORING_TIMESCALE:
-        description: |
-            "[s] default = 1.0E+06
-            The timescale used to nudge MEKE toward its equilibrium value."
-        datatype: real
-        units: nondim
-        value:
-            $OCN_GRID == "tx2_3v2": 1.0E+07
     MEKE_ADVECTION_FACTOR:
         description: |
             "[nondim] default = 0.0
diff --git a/param_templates/json/MOM_input.json b/param_templates/json/MOM_input.json
index 66099cc..d9d6bfc 100644
--- a/param_templates/json/MOM_input.json
+++ b/param_templates/json/MOM_input.json
@@ -1213,15 +1213,6 @@
             "$OCN_GRID == \"tx2_3v2\"": 0.02
          }
       },
-      "FOX_KEMPER_ML_RESTRAT_COEF": {
-         "description": "\"[nondim] default = 0.0\nA nondimensional coefficient that is proportional to\nthe ratio of the deformation radius to the dominant\nlengthscale of the submesoscale mixed layer\ninstabilities, times the minimum of the ratio of the\nmesoscale eddy kinetic energy to the large-scale\ngeostrophic kinetic energy or 1 plus the square of the\ngrid spacing over the deformation radius, as detailed\nby Fox-Kemper et al. (2010)\nThe file that specifies the vertical grid for\ndepth-space diagnostics, or blank to disable\ndepth-space output.\nThe number of depth-space levels.  This is determined\nfrom the size of the variable zw in the output grid file.\"\n",
-         "datatype": "real",
-         "units": "nondim",
-         "value": {
-            "$OCN_GRID == \"tx2_3v2\"": 1.0,
-            "$OCN_GRID == \"tx0.25v1\"": 1.0
-         }
-      },
       "MLE_FRONT_LENGTH": {
          "description": "\"[m] default = 0.0\nIf non-zero, is the frontal-length scale used to calculate the\nupscaling of buoyancy gradients that is otherwise represented\nby the parameter FOX_KEMPER_ML_RESTRAT_COEF. If MLE_FRONT_LENGTH is\nnon-zero, it is recommended to set FOX_KEMPER_ML_RESTRAT_COEF=1.\"\n",
          "datatype": "real",
@@ -1231,14 +1222,6 @@
             "$OCN_GRID == \"tx0.25v1\"": 500.0
          }
       },
-      "MLE_FRONT_LENGTH_FROM_FILE": {
-         "description": "\"[Boolean] default = False\nIf true, the MLE front-length scale is read from a file.\"\n",
-         "datatype": "logical",
-         "units": "Boolean",
-         "value": {
-            "$OCN_GRID == \"tx2_3v2\"": false
-         }
-      },
       "MLE_MLD_DECAY_TIME": {
          "description": "\"[s] default = 0.0\nThe time-scale for a running-mean filter applied to the mixed-layer\ndepth used in the MLE restratification parameterization. When\nthe MLD deepens below the current running-mean the running-mean\nis instantaneously set to the current MLD.\"\n",
          "datatype": "real",
@@ -2110,14 +2093,6 @@
             "$OCN_GRID == \"tx2_3v2\"": false
          }
       },
-      "MEKE_RESTORING_TIMESCALE": {
-         "description": "\"[s] default = 1.0E+06\nThe timescale used to nudge MEKE toward its equilibrium value.\"\n",
-         "datatype": "real",
-         "units": "nondim",
-         "value": {
-            "$OCN_GRID == \"tx2_3v2\"": 10000000.0
-         }
-      },
       "MEKE_ADVECTION_FACTOR": {
          "description": "\"[nondim] default = 0.0\nA scale factor in front of advection of eddy energy. Zero turns advection off.\nUsing unity would be normal but other values could accommodate a mismatch\nbetween the advecting barotropic flow and the vertical structure of MEKE.\"\n",
          "datatype": "real",