MOM6: +Add optional argument du_cor to continuity_PPM

- NOAA-GFDL/MOM6@78772923d +Add optional argument du_cor to continuity_PPM
- NOAA-GFDL/MOM6@6b6f8dcde Move continuity_fluxes to MOM_continuity_PPM.F90
- NOAA-GFDL/MOM6@a70482cfc +Add continuity_zonal_convergence
- NOAA-GFDL/MOM6@bbda30f12 +Add zonal_BT_mass_flux & meridional_BT_mass_flux
- NOAA-GFDL/MOM6@d3496fee3 +Call zonal_edge_thickness outside of zonal_mass_flux
- NOAA-GFDL/MOM6@6b86ab84d +Create zonal_edge_thickness in MOM_continuity_PPM
- NOAA-GFDL/MOM6@fc0bef1c2 +Pass-through from continuity_PPM via continuity
- NOAA-GFDL/MOM6@b898297e0 +(*)Add continuity_fluxes & continuity_adjust_vel
- NOAA-GFDL/MOM6@11759d6a2 Fix unpaired parentheses in a formatting string
- NOAA-GFDL/MOM6@34b880faa Add units to comments in MOM_neutral_diffusion
- NOAA-GFDL/MOM6@7cc191a6f Add rule and workflow to compare results of micro-timing tests
- NOAA-GFDL/MOM6@d8579d292 Use a class to deliver elemental procedures for EOS
- NOAA-GFDL/MOM6@cce4b3da5 Fix a bug that left OBC%debug uninitialized

coupled_AM2_LM3_SIS2/AM2_SIS2_MOM6i_1deg/MOM_parameter_doc.all

@@ -952,13 +952,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

coupled_AM2_LM3_SIS2/AM2_SIS2_MOM6i_1deg/MOM_parameter_doc.short

@@ -271,8 +271,6 @@ PORBAR_ANSWER_DATE = 99991231 ! default = 20181231
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.15E-09
  ! The tolerance for the differences between the barotropic and baroclinic

coupled_AM2_LM3_SIS2/Concurrent_ice_1deg/MOM_parameter_doc.all

@@ -952,13 +952,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

coupled_AM2_LM3_SIS2/Concurrent_ice_1deg/MOM_parameter_doc.short

@@ -251,8 +251,6 @@ KHTH_MAX = 900.0 ! [m2 s-1] default = 0.0
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.15E-09
  ! The tolerance for the differences between the barotropic and baroclinic

coupled_AM2_LM3_SIS2/Intersperse_ice_1deg/MOM_parameter_doc.all

@@ -952,13 +952,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

coupled_AM2_LM3_SIS2/Intersperse_ice_1deg/MOM_parameter_doc.short

@@ -254,8 +254,6 @@ KHTH_MAX = 900.0 ! [m2 s-1] default = 0.0
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.15E-09
  ! The tolerance for the differences between the barotropic and baroclinic

ice_ocean_SIS2/Baltic/MOM_parameter_doc.all

@@ -952,13 +952,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ice_ocean_SIS2/Baltic/MOM_parameter_doc.short

@@ -262,8 +262,6 @@ KHTH_MAX = 900.0 ! [m2 s-1] default = 0.0
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.15E-09
  ! The tolerance for the differences between the barotropic and baroclinic

ice_ocean_SIS2/Baltic_ALE_z_offline_tracers/MOM_parameter_doc.all

@@ -1093,13 +1093,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ice_ocean_SIS2/Baltic_ALE_z_offline_tracers/MOM_parameter_doc.short

@@ -296,8 +296,6 @@ KHTH_MAX = 900.0 ! [m2 s-1] default = 0.0
 TIDES = True ! [Boolean] default = False
  ! If true, apply tidal momentum forcing.
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 2.5E-09
  ! The tolerance for the differences between the barotropic and baroclinic

ice_ocean_SIS2/Baltic_OM4_025/MOM_parameter_doc.all

@@ -1187,13 +1187,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ice_ocean_SIS2/Baltic_OM4_025/MOM_parameter_doc.short

@@ -392,8 +392,6 @@ PORBAR_ANSWER_DATE = 99991231 ! default = 20181231
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-09
  ! The tolerance for the differences between the barotropic and baroclinic

ice_ocean_SIS2/Baltic_OM4_05/MOM_parameter_doc.all

@@ -1201,13 +1201,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ice_ocean_SIS2/Baltic_OM4_05/MOM_parameter_doc.short

@@ -409,8 +409,6 @@ PORBAR_ANSWER_DATE = 99991231 ! default = 20181231
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-09
  ! The tolerance for the differences between the barotropic and baroclinic

ice_ocean_SIS2/OM4_025/MOM_parameter_doc.all

@@ -1187,13 +1187,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ice_ocean_SIS2/OM4_025/MOM_parameter_doc.short

@@ -399,8 +399,6 @@ PORBAR_ANSWER_DATE = 99991231 ! default = 20181231
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-09
  ! The tolerance for the differences between the barotropic and baroclinic

ice_ocean_SIS2/OM4_05/MOM_parameter_doc.all

@@ -1201,13 +1201,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ice_ocean_SIS2/OM4_05/MOM_parameter_doc.short

@@ -414,8 +414,6 @@ PORBAR_ANSWER_DATE = 99991231 ! default = 20181231
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-09
  ! The tolerance for the differences between the barotropic and baroclinic

ice_ocean_SIS2/OM_1deg/MOM_parameter_doc.all

@@ -1178,13 +1178,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ice_ocean_SIS2/OM_1deg/MOM_parameter_doc.short

@@ -426,8 +426,6 @@ PORBAR_ANSWER_DATE = 99991231 ! default = 20181231
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.75E-09
  ! The tolerance for the differences between the barotropic and baroclinic

ice_ocean_SIS2/SIS2/MOM_parameter_doc.all

@@ -924,13 +924,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ice_ocean_SIS2/SIS2/MOM_parameter_doc.short

@@ -289,8 +289,6 @@ KHTH_MAX = 900.0 ! [m2 s-1] default = 0.0
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.15E-09
  ! The tolerance for the differences between the barotropic and baroclinic

ice_ocean_SIS2/SIS2_bergs_cgrid/MOM_parameter_doc.all

@@ -924,13 +924,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ice_ocean_SIS2/SIS2_bergs_cgrid/MOM_parameter_doc.short

@@ -289,8 +289,6 @@ KHTH_MAX = 900.0 ! [m2 s-1] default = 0.0
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.15E-09
  ! The tolerance for the differences between the barotropic and baroclinic

ice_ocean_SIS2/SIS2_cgrid/MOM_parameter_doc.all

@@ -924,13 +924,6 @@ STORE_CORIOLIS_ACCEL = True ! [Boolean] default = True
 FPMIX = False ! [Boolean] default = False
  ! If true, apply profiles of momentum flux magnitude and direction
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ice_ocean_SIS2/SIS2_cgrid/MOM_parameter_doc.short

@@ -289,8 +289,6 @@ KHTH_MAX = 900.0 ! [m2 s-1] default = 0.0
 
 ! === module MOM_dynamics_split_RK2 ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 ETA_TOLERANCE = 1.0E-06 ! [m] default = 3.15E-09
  ! The tolerance for the differences between the barotropic and baroclinic

ocean_only/CVmix_SCM_tests/cooling_only/BML/MOM_parameter_doc.all

@@ -753,13 +753,6 @@ TIDES = False ! [Boolean] default = False
 CALCULATE_SAL = False ! [Boolean] default = False
  ! If true, calculate self-attraction and loading.
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ocean_only/CVmix_SCM_tests/cooling_only/BML/MOM_parameter_doc.short

@@ -196,8 +196,6 @@ KV = 0.0 ! [m2 s-1]
 
 ! === module MOM_dynamics_unsplit ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 UPWIND_1ST_CONTINUITY = True ! [Boolean] default = False
  ! If true, CONTINUITY_PPM becomes a 1st-order upwind continuity solver. This

ocean_only/CVmix_SCM_tests/cooling_only/EPBL/MOM_parameter_doc.all

@@ -862,13 +862,6 @@ TIDES = False ! [Boolean] default = False
 CALCULATE_SAL = False ! [Boolean] default = False
  ! If true, calculate self-attraction and loading.
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ocean_only/CVmix_SCM_tests/cooling_only/EPBL/MOM_parameter_doc.short

@@ -225,8 +225,6 @@ KV = 0.0 ! [m2 s-1]
 
 ! === module MOM_dynamics_unsplit ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 UPWIND_1ST_CONTINUITY = True ! [Boolean] default = False
  ! If true, CONTINUITY_PPM becomes a 1st-order upwind continuity solver. This

ocean_only/CVmix_SCM_tests/cooling_only/KPP/MOM_parameter_doc.all

@@ -862,13 +862,6 @@ TIDES = False ! [Boolean] default = False
 CALCULATE_SAL = False ! [Boolean] default = False
  ! If true, calculate self-attraction and loading.
 
-! === module MOM_continuity ===
-CONTINUITY_SCHEME = "PPM" ! default = "PPM"
- ! CONTINUITY_SCHEME selects the discretization for the continuity solver. The
- ! only valid value currently is:
- ! PPM - use a positive-definite (or monotonic)
- ! piecewise parabolic reconstruction solver.
-
 ! === module MOM_continuity_PPM ===
 MONOTONIC_CONTINUITY = False ! [Boolean] default = False
  ! If true, CONTINUITY_PPM uses the Colella and Woodward monotonic limiter. The

ocean_only/CVmix_SCM_tests/cooling_only/KPP/MOM_parameter_doc.short

@@ -225,8 +225,6 @@ KV = 0.0 ! [m2 s-1]
 
 ! === module MOM_dynamics_unsplit ===
 
-! === module MOM_continuity ===
-
 ! === module MOM_continuity_PPM ===
 UPWIND_1ST_CONTINUITY = True ! [Boolean] default = False
  ! If true, CONTINUITY_PPM becomes a 1st-order upwind continuity solver. This