MOM6: +Rotationally symmetric epipycnal diffusion option

- NOAA-GFDL/MOM6@a3fd1f3a7 +Rotationally symmetric epipycnal diffusion option
- NOAA-GFDL/MOM6@f9372f3d6 +Add get_netcdf_filename for a get_field_nc error
- NOAA-GFDL/MOM6@5e34f486b Set flags properly for rotated tripolar grids
- NOAA-GFDL/MOM6@2b59089ea Set up auxiliary domain for unrotated grid

coupled_AM2_LM3_SIS2/AM2_SIS2_MOM6i_1deg/MOM_parameter_doc.all

@@ -1948,6 +1948,15 @@ RECALC_NEUTRAL_SURF = False ! [Boolean] default = False
  ! If true, then recalculate the neutral surfaces if the
  ! diffusive CFL is exceeded. If false, assume that the
  ! positions of the surfaces do not change
+HOR_DIFF_ANSWER_DATE = 20240101 ! default = 20240101
+ ! The vintage of the order of arithmetic to use for the tracer diffusion.
+ ! Values of 20240330 or below recover the answers from the original form of the
+ ! along-isopycnal mixed layer to interior mixing code, while higher values use
+ ! mathematically equivalent expressions that recover rotational symmetry when
+ ! DIFFUSE_ML_TO_INTERIOR is true.
+HOR_DIFF_LIMIT_BUG = True ! [Boolean] default = True
+ ! If true and the answer date is 20240330 or below, use a rotational symmetry
+ ! breaking bug when limiting the tracer properties in tracer_epipycnal_ML_diff.
 ML_KHTR_SCALE = 0.0 ! [nondim] default = 1.0
  ! With Diffuse_ML_interior, the ratio of the truly horizontal diffusivity in the
  ! mixed layer to the epipycnal diffusivity. The valid range is 0 to 1.

coupled_AM2_LM3_SIS2/Concurrent_ice_1deg/MOM_parameter_doc.all

@@ -1948,6 +1948,15 @@ RECALC_NEUTRAL_SURF = False ! [Boolean] default = False
  ! If true, then recalculate the neutral surfaces if the
  ! diffusive CFL is exceeded. If false, assume that the
  ! positions of the surfaces do not change
+HOR_DIFF_ANSWER_DATE = 20240101 ! default = 20240101
+ ! The vintage of the order of arithmetic to use for the tracer diffusion.
+ ! Values of 20240330 or below recover the answers from the original form of the
+ ! along-isopycnal mixed layer to interior mixing code, while higher values use
+ ! mathematically equivalent expressions that recover rotational symmetry when
+ ! DIFFUSE_ML_TO_INTERIOR is true.
+HOR_DIFF_LIMIT_BUG = True ! [Boolean] default = True
+ ! If true and the answer date is 20240330 or below, use a rotational symmetry
+ ! breaking bug when limiting the tracer properties in tracer_epipycnal_ML_diff.
 ML_KHTR_SCALE = 0.0 ! [nondim] default = 1.0
  ! With Diffuse_ML_interior, the ratio of the truly horizontal diffusivity in the
  ! mixed layer to the epipycnal diffusivity. The valid range is 0 to 1.

coupled_AM2_LM3_SIS2/Intersperse_ice_1deg/MOM_parameter_doc.all

@@ -1948,6 +1948,15 @@ RECALC_NEUTRAL_SURF = False ! [Boolean] default = False
  ! If true, then recalculate the neutral surfaces if the
  ! diffusive CFL is exceeded. If false, assume that the
  ! positions of the surfaces do not change
+HOR_DIFF_ANSWER_DATE = 20240101 ! default = 20240101
+ ! The vintage of the order of arithmetic to use for the tracer diffusion.
+ ! Values of 20240330 or below recover the answers from the original form of the
+ ! along-isopycnal mixed layer to interior mixing code, while higher values use
+ ! mathematically equivalent expressions that recover rotational symmetry when
+ ! DIFFUSE_ML_TO_INTERIOR is true.
+HOR_DIFF_LIMIT_BUG = True ! [Boolean] default = True
+ ! If true and the answer date is 20240330 or below, use a rotational symmetry
+ ! breaking bug when limiting the tracer properties in tracer_epipycnal_ML_diff.
 ML_KHTR_SCALE = 0.0 ! [nondim] default = 1.0
  ! With Diffuse_ML_interior, the ratio of the truly horizontal diffusivity in the
  ! mixed layer to the epipycnal diffusivity. The valid range is 0 to 1.

ice_ocean_SIS2/Baltic/MOM_parameter_doc.all

@@ -1948,6 +1948,15 @@ RECALC_NEUTRAL_SURF = False ! [Boolean] default = False
  ! If true, then recalculate the neutral surfaces if the
  ! diffusive CFL is exceeded. If false, assume that the
  ! positions of the surfaces do not change
+HOR_DIFF_ANSWER_DATE = 20240101 ! default = 20240101
+ ! The vintage of the order of arithmetic to use for the tracer diffusion.
+ ! Values of 20240330 or below recover the answers from the original form of the
+ ! along-isopycnal mixed layer to interior mixing code, while higher values use
+ ! mathematically equivalent expressions that recover rotational symmetry when
+ ! DIFFUSE_ML_TO_INTERIOR is true.
+HOR_DIFF_LIMIT_BUG = True ! [Boolean] default = True
+ ! If true and the answer date is 20240330 or below, use a rotational symmetry
+ ! breaking bug when limiting the tracer properties in tracer_epipycnal_ML_diff.
 ML_KHTR_SCALE = 0.0 ! [nondim] default = 1.0
  ! With Diffuse_ML_interior, the ratio of the truly horizontal diffusivity in the
  ! mixed layer to the epipycnal diffusivity. The valid range is 0 to 1.

ice_ocean_SIS2/SIS2/MOM_parameter_doc.all

@@ -1920,6 +1920,15 @@ RECALC_NEUTRAL_SURF = False ! [Boolean] default = False
  ! If true, then recalculate the neutral surfaces if the
  ! diffusive CFL is exceeded. If false, assume that the
  ! positions of the surfaces do not change
+HOR_DIFF_ANSWER_DATE = 20240101 ! default = 20240101
+ ! The vintage of the order of arithmetic to use for the tracer diffusion.
+ ! Values of 20240330 or below recover the answers from the original form of the
+ ! along-isopycnal mixed layer to interior mixing code, while higher values use
+ ! mathematically equivalent expressions that recover rotational symmetry when
+ ! DIFFUSE_ML_TO_INTERIOR is true.
+HOR_DIFF_LIMIT_BUG = True ! [Boolean] default = True
+ ! If true and the answer date is 20240330 or below, use a rotational symmetry
+ ! breaking bug when limiting the tracer properties in tracer_epipycnal_ML_diff.
 ML_KHTR_SCALE = 0.0 ! [nondim] default = 1.0
  ! With Diffuse_ML_interior, the ratio of the truly horizontal diffusivity in the
  ! mixed layer to the epipycnal diffusivity. The valid range is 0 to 1.

ice_ocean_SIS2/SIS2_bergs_cgrid/MOM_parameter_doc.all

@@ -1924,6 +1924,15 @@ RECALC_NEUTRAL_SURF = False ! [Boolean] default = False
  ! If true, then recalculate the neutral surfaces if the
  ! diffusive CFL is exceeded. If false, assume that the
  ! positions of the surfaces do not change
+HOR_DIFF_ANSWER_DATE = 20240101 ! default = 20240101
+ ! The vintage of the order of arithmetic to use for the tracer diffusion.
+ ! Values of 20240330 or below recover the answers from the original form of the
+ ! along-isopycnal mixed layer to interior mixing code, while higher values use
+ ! mathematically equivalent expressions that recover rotational symmetry when
+ ! DIFFUSE_ML_TO_INTERIOR is true.
+HOR_DIFF_LIMIT_BUG = True ! [Boolean] default = True
+ ! If true and the answer date is 20240330 or below, use a rotational symmetry
+ ! breaking bug when limiting the tracer properties in tracer_epipycnal_ML_diff.
 ML_KHTR_SCALE = 0.0 ! [nondim] default = 1.0
  ! With Diffuse_ML_interior, the ratio of the truly horizontal diffusivity in the
  ! mixed layer to the epipycnal diffusivity. The valid range is 0 to 1.

ice_ocean_SIS2/SIS2_cgrid/MOM_parameter_doc.all

@@ -1920,6 +1920,15 @@ RECALC_NEUTRAL_SURF = False ! [Boolean] default = False
  ! If true, then recalculate the neutral surfaces if the
  ! diffusive CFL is exceeded. If false, assume that the
  ! positions of the surfaces do not change
+HOR_DIFF_ANSWER_DATE = 20240101 ! default = 20240101
+ ! The vintage of the order of arithmetic to use for the tracer diffusion.
+ ! Values of 20240330 or below recover the answers from the original form of the
+ ! along-isopycnal mixed layer to interior mixing code, while higher values use
+ ! mathematically equivalent expressions that recover rotational symmetry when
+ ! DIFFUSE_ML_TO_INTERIOR is true.
+HOR_DIFF_LIMIT_BUG = True ! [Boolean] default = True
+ ! If true and the answer date is 20240330 or below, use a rotational symmetry
+ ! breaking bug when limiting the tracer properties in tracer_epipycnal_ML_diff.
 ML_KHTR_SCALE = 0.0 ! [nondim] default = 1.0
  ! With Diffuse_ML_interior, the ratio of the truly horizontal diffusivity in the
  ! mixed layer to the epipycnal diffusivity. The valid range is 0 to 1.

ocean_only/benchmark/MOM_parameter_doc.all

@@ -1717,6 +1717,15 @@ RECALC_NEUTRAL_SURF = False ! [Boolean] default = False
  ! If true, then recalculate the neutral surfaces if the
  ! diffusive CFL is exceeded. If false, assume that the
  ! positions of the surfaces do not change
+HOR_DIFF_ANSWER_DATE = 20240101 ! default = 20240101
+ ! The vintage of the order of arithmetic to use for the tracer diffusion.
+ ! Values of 20240330 or below recover the answers from the original form of the
+ ! along-isopycnal mixed layer to interior mixing code, while higher values use
+ ! mathematically equivalent expressions that recover rotational symmetry when
+ ! DIFFUSE_ML_TO_INTERIOR is true.
+HOR_DIFF_LIMIT_BUG = True ! [Boolean] default = True
+ ! If true and the answer date is 20240330 or below, use a rotational symmetry
+ ! breaking bug when limiting the tracer properties in tracer_epipycnal_ML_diff.
 ML_KHTR_SCALE = 0.0 ! [nondim] default = 1.0
  ! With Diffuse_ML_interior, the ratio of the truly horizontal diffusivity in the
  ! mixed layer to the epipycnal diffusivity. The valid range is 0 to 1.

ocean_only/nonBous_global/MOM_parameter_doc.all

@@ -1899,6 +1899,15 @@ RECALC_NEUTRAL_SURF = False ! [Boolean] default = False
  ! If true, then recalculate the neutral surfaces if the
  ! diffusive CFL is exceeded. If false, assume that the
  ! positions of the surfaces do not change
+HOR_DIFF_ANSWER_DATE = 20240101 ! default = 20240101
+ ! The vintage of the order of arithmetic to use for the tracer diffusion.
+ ! Values of 20240330 or below recover the answers from the original form of the
+ ! along-isopycnal mixed layer to interior mixing code, while higher values use
+ ! mathematically equivalent expressions that recover rotational symmetry when
+ ! DIFFUSE_ML_TO_INTERIOR is true.
+HOR_DIFF_LIMIT_BUG = True ! [Boolean] default = True
+ ! If true and the answer date is 20240330 or below, use a rotational symmetry
+ ! breaking bug when limiting the tracer properties in tracer_epipycnal_ML_diff.
 ML_KHTR_SCALE = 0.0 ! [nondim] default = 1.0
  ! With Diffuse_ML_interior, the ratio of the truly horizontal diffusivity in the
  ! mixed layer to the epipycnal diffusivity. The valid range is 0 to 1.