MOM6 parameter changes for 0.4.0

[dougiesquire]

With the updated version of MOM6 in the recently released ACCESS-OM3 0.4.0, there are a number of new bug fixes and recommended settings. These default to values that preserve answers and so are not used unless explicitly turned on in MOM_input. However, the CICE and CDEPS/CMEPS updates in ACCESS-OM3 0.4.0 broke reproducibility with 0.3.1 so there's no real value in using MOM6 settings that preserve MOM6 answers. Instead, we should move to these recommended settings:

• EQN_OF_STATE = "WRIGHT_REDUCED" ! default = "WRIGHT"
                                  ! EQN_OF_STATE determines which ocean equation of state should be used.
                                  ! Currently, the valid choices are "LINEAR", "UNESCO", "JACKETT_MCD", "WRIGHT",
                                  ! "WRIGHT_REDUCED", "WRIGHT_FULL", "NEMO", "ROQUET_RHO", "ROQUET_SPV" and
                                  ! "TEOS10".  This is only used if USE_EOS is true.
  

  Justification: GFDL to main (2024-05-31) mom-ocean/MOM6#1631 (comment)

• VISC_REM_BUG = False            !   [Boolean] default = True
                                  ! If true, visc_rem_[uv] in split mode is incorrectly calculated or accounted
                                  ! for in two places. This parameter controls the defaults of two individual
                                  ! flags, VISC_REM_TIMESTEP_BUG in MOM_dynamics_split_RK2(b) and
                                  ! VISC_REM_BT_WEIGHT_BUG in MOM_barotropic.
  

  Justification: see parameter description

• BACKSCATTER_UNDERBOUND = False  !   [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.
  

  Justification: see parameter description

• FRICTWORK_BUG = False           !   [Boolean] default = True
                                  ! If true, retain an answer-changing bug in calculating the FrictWork, which
                                  ! cancels the h in thickness flux and the h at velocity point. This isnot
                                  ! recommended.
  

  Justification: see parameter description

• REMAPPING_USE_OM4_SUBCELLS = False !   [Boolean] default = True
                                  ! This selects the remapping algorithm used in OM4 that does not use the full
                                  ! reconstruction for the top- and lower-most sub-layers, but instead assumes
                                  ! they are always vanished (untrue) and so just uses their edge values. We
                                  ! recommend setting this option to false.
  
  Z_INIT_REMAPPING_USE_OM4_SUBCELLS = False !   [Boolean] default = True
                                  ! If true, use the OM4 remapping-via-subcells algorithm for initialization. See
                                  ! REMAPPING_USE_OM4_SUBCELLS for more details. We recommend setting this option
                                  ! to false.
  
  DIAG_REMAPPING_USE_OM4_SUBCELLS = False !   [Boolean] default = True
                                  ! If true, use the OM4 remapping-via-subcells algorithm for diagnostics. See
                                  ! REMAPPING_USE_OM4_SUBCELLS for details. We recommend setting this option to
                                  ! false.
  
  EBT_REMAPPING_USE_OM4_SUBCELLS = False !   [Boolean] default = True
                                  ! If true, use the OM4 remapping-via-subcells algorithm for calculating EBT
                                  ! structure. See REMAPPING_USE_OM4_SUBCELLS for details. We recommend setting
                                  ! this option to false.
  
  INTWAVE_REMAPPING_USE_OM4_SUBCELLS = False !   [Boolean] default = True
                                  ! If true, use the OM4 remapping-via-subcells algorithm for calculating EBT
                                  ! structure. See REMAPPING_USE_OM4_SUBCELLS for details. We recommend setting
                                  ! this option to false.
  
  HBD_REMAPPING_USE_OM4_SUBCELLS = False !   [Boolean] default = True
                                  ! If true, use the OM4 remapping-via-subcells algorithm for horizontal boundary
                                  ! diffusion. See REMAPPING_USE_OM4_SUBCELLS for details. We recommend setting
                                  ! this option to false.
  

  Justification: see parameter description

I'll start with a PR to the dev-1deg_jra55do_ryf branch, then will cherry-pick into the open PRs for dev-1deg_jra55do_iaf, dev-1deg_jra55do_ryf_wombatlite and dev-025deg_jra55do_ryf

[dougiesquire]

Note, we have also been defaulting to:

TIDE_USE_SAL_SCALAR = True      !   [Boolean] default = True
                                ! If true and TIDES is true, use the scalar approximation when calculating
                                ! self-attraction and loading.

and setting

TIDE_SAL_SCALAR_VALUE = 0.094   !   [m m-1]
                                ! The constant of proportionality between sea surface height (really it should
                                ! be bottom pressure) anomalies and bottom geopotential anomalies. This is only
                                ! used if TIDES and TIDE_USE_SAL_SCALAR are true.

However, these have now been deprecated (see NOAA-GFDL/MOM6#382).

The replacement parameters take default values equal to what we were using, so we can just use the defaults:

SAL_SCALAR_APPROX = True        !   [Boolean] default = True
                                ! If true, use the scalar approximation to calculate self-attraction and
                                ! loading.
SAL_SCALAR_VALUE = 0.094        !   [m m-1] default = 0.094
                                ! The constant of proportionality between sea surface height (really it should
                                ! be bottom pressure) anomalies and bottom geopotential anomalies. This is only
                                ! used if USE_SAL_SCALAR is true or USE_PREVIOUS_TIDES is true.

[dougiesquire]

@minghangli-uni also suggested:

AUTO_MASKTABLE = True           !   [Boolean] default = False
                                ! Turn on automatic mask table generation to eliminate land blocks.

This is currently set in our 0.25deg config, but not the 1deg configs