MOM6: +Add BACKSCATTER_UNDERBOUND

- NOAA-GFDL/MOM6@36dda7e4e +Add BACKSCATTER_UNDERBOUND
- NOAA-GFDL/MOM6@31afce04b (*)Avoid diagnostic seg fault in a 10km deep ocean
- NOAA-GFDL/MOM6@215b9609c Refactor PressureForce_FV before ice shelf fixes

coupled_AM2_LM3_SIS2/AM2_SIS2_MOM6i_1deg/MOM_parameter_doc.all

@@ -1157,6 +1157,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

coupled_AM2_LM3_SIS2/Concurrent_ice_1deg/MOM_parameter_doc.all

@@ -1157,6 +1157,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

coupled_AM2_LM3_SIS2/Intersperse_ice_1deg/MOM_parameter_doc.all

@@ -1157,6 +1157,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ice_ocean_SIS2/Baltic/MOM_parameter_doc.all

@@ -1157,6 +1157,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ice_ocean_SIS2/Baltic_ALE_z_offline_tracers/MOM_parameter_doc.all

@@ -1368,6 +1368,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ice_ocean_SIS2/Baltic_OM4_025/MOM_parameter_doc.all

@@ -1392,6 +1392,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ice_ocean_SIS2/Baltic_OM4_05/MOM_parameter_doc.all

@@ -1411,6 +1411,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ice_ocean_SIS2/OM4_025/MOM_parameter_doc.all

@@ -1392,6 +1392,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ice_ocean_SIS2/OM4_05/MOM_parameter_doc.all

@@ -1411,6 +1411,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ice_ocean_SIS2/OM_1deg/MOM_parameter_doc.all

@@ -1385,6 +1385,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ice_ocean_SIS2/SIS2/MOM_parameter_doc.all

@@ -1123,6 +1123,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ice_ocean_SIS2/SIS2_bergs_cgrid/MOM_parameter_doc.all

@@ -1123,6 +1123,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ice_ocean_SIS2/SIS2_cgrid/MOM_parameter_doc.all

@@ -1123,6 +1123,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ocean_only/circle_obcs/MOM_parameter_doc.all

@@ -991,6 +991,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ocean_only/global_ALE/hycom/MOM_parameter_doc.all

@@ -1418,6 +1418,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True

ocean_only/global_ALE/layer/MOM_parameter_doc.all

@@ -1313,6 +1313,13 @@ BETTER_BOUND_AH = True ! [Boolean] default = True
 RE_AH = 0.0 ! [nondim] default = 0.0
  ! If nonzero, the biharmonic coefficient is scaled so that the biharmonic
  ! Reynolds number is equal to this.
+BACKSCATTER_UNDERBOUND = True ! [Boolean] default = True
+ ! If true, the bounds on the biharmonic viscosity are allowed to increase where
+ ! the Laplacian viscosity is negative (due to backscatter parameterizations)
+ ! beyond the largest timestep-dependent stable values of biharmonic viscosity
+ ! when no Laplacian viscosity is applied. The default is true for historical
+ ! reasons, but this option probably should not be used because it can contribute
+ ! to numerical instabilities.
 SMAG_BI_CONST = 0.06 ! [nondim] default = 0.0
  ! The nondimensional biharmonic Smagorinsky constant, typically 0.015 - 0.06.
 BOUND_CORIOLIS_BIHARM = True ! [Boolean] default = True