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.