Skip to content

Commit 70a0af7

Browse files
gustavo-marquesgustavo-marques
authored
More miscellaneous runtime parameter changes for beta05 (#219)
* Update tx2_3v2 defaults for beta05 - Added new topographic edits and updated channel list. - Enabled COMPUTE_MONIN_OBUKHOV = True to improve the OBL seasonal cycle by limiting the OBL depth to no deeper than the Monin-Obukhov depth. - Updated background mixing parameters: - KV = 1.0E-6 - KD = 1.0E-7 - PRANDTL_BKGND = 10 - Configured MLE backup settings: - Set MLE_FRONT_LENGTH_FROM_FILE = True - Updated MLE_FL_FILE to "mle-lf-clim-tx2_3v2_20250108.nc" - Enabled Leithy parameterization with: - USE_LEITHY = True - LEITH_BI_CONST = 90.0 - LEITHY_CK = 1.0 (full backscatter) - LEITH_AH = False * Add input files to the data list * Add missing input files to the data list * Remove duplication (TOPO_EDITS_FILE) * Set MINIMUM_OBL_DEPTH=2.5 m for tx2_3v2 * Enable bug fixes byt setting them to False * Update MLE front-length climatology * Update MLE_FL_FILE * Turn off backscatter and decrease biharmonic Leith constant * Revert back to LEITH_AH, MLE_FRONT_LENGTH=100, and Laplacian=False
1 parent 6443e9b commit 70a0af7

4 files changed

Lines changed: 171 additions & 31 deletions

File tree

‎param_templates/MOM_input.yaml‎

Lines changed: 87 additions & 16 deletions
Original file line numberDiff line numberDiff line change
@@ -220,7 +220,7 @@ Global:
220220
for backward compatibility but ideally should be 0."
221221
datatype: real
222222
units: PPT
223-
value:
223+
value:
224224
$COMP_ATM == "cam": 1.0e-6
225225
else: 0.0
226226
BOUND_SALINITY:
@@ -560,6 +560,14 @@ Global:
560560
value:
561561
$OCN_GRID == "tx2_3v2": "ocean_topo_tx2_3v2_240501.nc"
562562
$OCN_GRID == "tx0.25v1": "ocean_topog.nc"
563+
TOPO_EDITS_FILE:
564+
description: |
565+
"default = ''
566+
The file from which to read a list of i,j,z topography overrides."
567+
datatype: string
568+
value:
569+
$OCN_GRID == "tx2_3v2": "topo_edits_tx2_3v2_250107.nc"
570+
$OCN_GRID == "tx0.25v1": "All_edits.nc"
563571
MAXIMUM_DEPTH:
564572
description: |
565573
"[m]
@@ -990,7 +998,9 @@ Global:
990998
If true, use a Laplacian horizontal viscosity."
991999
datatype: logical
9921000
units: Boolean
993-
value: True
1001+
value:
1002+
$OCN_GRID == "tx2_3v2": False
1003+
else: True
9941004
KH:
9951005
description: |
9961006
"[m2 s-1] default = 0.0
@@ -1052,6 +1062,16 @@ Global:
10521062
value:
10531063
$OCN_GRID == "tx0.25v1": 0.01
10541064
$OCN_GRID == "MISOMIP": 0.001
1065+
USE_LEITHY:
1066+
description: |
1067+
"[Boolean] default = False
1068+
If true, use a modified version of the biharmonic Leith nonlinear
1069+
eddy viscosity. This implementation can include harmonic backscatter
1070+
when LEITHY_CK > 0."
1071+
datatype: logical
1072+
units: Boolean
1073+
value:
1074+
$OCN_GRID == "tx2_3v2": False
10551075
LEITH_AH:
10561076
description: |
10571077
"[Boolean] default = False
@@ -1069,6 +1089,14 @@ Global:
10691089
units: nondim
10701090
value:
10711091
$OCN_GRID == "tx2_3v2": 128.0
1092+
LEITHY_CK:
1093+
description: |
1094+
"[nondim] default = 1.0
1095+
Fraction of biharmonic dissipation that gets backscattered, in Leith+E."
1096+
datatype: real
1097+
units: nondim
1098+
value:
1099+
$OCN_GRID == "tx2_3v2": 0.0
10721100
USE_LAND_MASK_FOR_HVISC:
10731101
description: |
10741102
"[Boolean] default = False
@@ -1080,6 +1108,15 @@ Global:
10801108
datatype: logical
10811109
units: Boolean
10821110
value: True
1111+
FRICTWORK_BUG:
1112+
description: |
1113+
"[Boolean] default = True
1114+
If true, retain an answer-changing bug in calculating the FrictWork, which
1115+
cancels the h in thickness flux and the h at velocity point. This is not
1116+
recommended."
1117+
datatype: logical
1118+
units: Boolean
1119+
value: False
10831120
HMIX_FIXED:
10841121
description: |
10851122
"[m]
@@ -1504,6 +1541,21 @@ Global:
15041541
value:
15051542
$OCN_GRID == "tx2_3v2": 1000.0
15061543
$OCN_GRID == "tx0.25v1": 500.0
1544+
MLE_FRONT_LENGTH_FROM_FILE:
1545+
description: |
1546+
"[Boolean] default = False
1547+
If true, the MLE front-length scale is read from a file."
1548+
datatype: logical
1549+
units: Boolean
1550+
value:
1551+
$OCN_GRID == "tx2_3v2": False
1552+
MLE_FL_FILE:
1553+
description: |
1554+
"The path to the file containing the MLE
1555+
front-length scale.."
1556+
datatype: string
1557+
value:
1558+
$OCN_GRID == "tx2_3v2": "mle-lf-clim-tx2_3v2_20250115.nc"
15071559
MLE_MLD_DECAY_TIME:
15081560
description: |
15091561
"[s] default = 0.0
@@ -1587,7 +1639,7 @@ Global:
15871639
datatype: real
15881640
units: m2 s-1
15891641
value:
1590-
$OCN_GRID == "tx2_3v2": 0.0
1642+
$OCN_GRID == "tx2_3v2": 1.0E-07
15911643
$OCN_GRID == "tx0.25v1": 1.5E-05
15921644
$OCN_GRID == "MISOMIP": 5.0E-05
15931645
else: 2.0E-05
@@ -1609,7 +1661,7 @@ Global:
16091661
value: # KV/KD:
16101662
$OCN_GRID == "tx0.25v1": = 1.0E-04 / 1.5E-05
16111663
$OCN_GRID == "MISOMIP": = 1.0E-04 / 5.0E-05
1612-
else: = 1.0E-04 / 2.0E-05
1664+
else: = 1.0E-06 / 1.0E-07
16131665
KD_MIN:
16141666
description: |
16151667
"[m2 s-1] default = 2.0E-07
@@ -1620,7 +1672,7 @@ Global:
16201672
datatype: real
16211673
units: m2 s-1
16221674
value:
1623-
$OCN_GRID == "tx2_3v2": 2.0E-07
1675+
$OCN_GRID == "tx2_3v2": 1.0E-07
16241676
else: 2.0E-06
16251677
INT_TIDE_DECAY_SCALE:
16261678
description: |
@@ -2564,7 +2616,7 @@ Global:
25642616
datatype: string
25652617
value:
25662618
$OCN_GRID == "tx0.25v1": "MOM_channels_global_025"
2567-
$OCN_GRID == "tx2_3v2": "MOM_channels_global_tx2_3v2_240501"
2619+
$OCN_GRID == "tx2_3v2": "channels_tx2_3v2_250107.txt"
25682620
SMAG_BI_CONST:
25692621
description: |
25702622
"[nondim] default = 0.0
@@ -2697,13 +2749,6 @@ Global:
26972749
value:
26982750
$OCN_GRID == "tx2_3v2": True
26992751
$OCN_GRID == "tx0.25v1": True
2700-
TOPO_EDITS_FILE:
2701-
description: |
2702-
"default =
2703-
The file from which to read a list of i,j,z topography overrides."
2704-
datatype: string
2705-
value:
2706-
$OCN_GRID == "tx0.25v1": "All_edits.nc"
27072752
PARALLEL_RESTARTFILES:
27082753
description: |
27092754
"[Boolean] default = False
@@ -2730,6 +2775,16 @@ Global:
27302775
units: nondim
27312776
value:
27322777
$COMP_WAV == "ww3" and $MOM6_WW3_CPL_METHOD == "most": 2.0
2778+
VISC_REM_BUG:
2779+
description: |
2780+
"[Boolean] default = True
2781+
If true, visc_rem_[uv] in split mode is incorrectly calculated or accounted
2782+
for in two places. This parameter controls the defaults of two individual
2783+
flags, VISC_REM_TIMESTEP_BUG in MOM_dynamics_split_RK2(b) and
2784+
VISC_REM_BT_WEIGHT_BUG in MOM_barotropic."
2785+
datatype: logical
2786+
units: Boolean
2787+
value: False
27332788

27342789
# MISOMIP-only variables:
27352790

@@ -3150,7 +3205,7 @@ Global:
31503205
datatype: real
31513206
units: m2 s-1
31523207
value:
3153-
$OCN_GRID == "tx2_3v2": 0.0
3208+
$OCN_GRID == "tx2_3v2": 1.0E-06
31543209
else: 1.0E-04
31553210
KV_BBL_MIN:
31563211
description: |
@@ -3698,6 +3753,22 @@ KPP:
36983753
If True, use Stokes Similarity package.
36993754
datatype: logical
37003755
value:
3701-
$COMP_WAV == "ww3" and $MOM6_WW3_CPL_METHOD == "most": True
3702-
3756+
$COMP_WAV == "ww3" and $MOM6_WW3_CPL_METHOD == "most": True
3757+
COMPUTE_MONIN_OBUKHOV:
3758+
description: |
3759+
default = False
3760+
If True, limit the OBL depth to be no deeper than Monin-Obukhov depth.
3761+
datatype: logical
3762+
units: Boolean
3763+
value:
3764+
$OCN_GRID == "tx2_3v2": True
3765+
MINIMUM_OBL_DEPTH:
3766+
description: |
3767+
"[m] default = 0.0
3768+
If non-zero, a minimum depth to use for KPP OBL depth. Independent of this
3769+
parameter, the OBL depth is always at least as deep as the first layer."
3770+
datatype: real
3771+
value:
3772+
$OCN_GRID == "tx2_3v2": 2.5
37033773
...
3774+

‎param_templates/input_data_list.yaml‎

Lines changed: 4 additions & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -18,6 +18,7 @@ mom.input_data_list:
1818
$OCN_GRID == "tx2_3v2": "${INPUTDIR}/ocean_topo_tx2_3v2_240501.nc"
1919
$OCN_GRID == "tx0.25v1": "${INPUTDIR}/ocean_topog.nc"
2020
TOPO_EDITS_FILE:
21+
$OCN_GRID == "tx2_3v2": "${INPUTDIR}/topo_edits_tx2_3v2_250107.nc"
2122
$OCN_GRID == "tx0.25v1": "${INPUTDIR}/All_edits.nc"
2223
TEMP_SALT_Z_INIT_FILE:
2324
$OCN_GRID in ["tx2_3v2", "tx0.25v1"]:
@@ -34,14 +35,16 @@ mom.input_data_list:
3435
$OCN_GRID == "tx2_3v2": "${INPUTDIR}/energy_new_tx2_3_conserve_230415_cdf5.nc"
3536
$OCN_GRID == "tx0.25v1": "${INPUTDIR}/tidal_amplitude.v20140616.nc"
3637
CHANNEL_LIST_FILE:
37-
$OCN_GRID == "tx2_3v2": "${INPUTDIR}/MOM_channels_global_tx2_3v2_240501"
38+
$OCN_GRID == "tx2_3v2": "${INPUTDIR}/channels_tx2_3v2_250107.txt"
3839
$OCN_GRID == "tx0.25v1": "${INPUTDIR}/MOM_channels_global_025"
3940
GEOTHERMAL_FILE:
4041
$OCN_GRID == "tx2_3v2": "${INPUTDIR}/geothermal_davies2013_tx2_3_20240318_cdf5.nc"
4142
$OCN_GRID == "tx0.25v1": "${INPUTDIR}/geothermal_davies2013_v1.nc"
4243
CHL_FILE:
4344
$OCN_GRID == "tx0.25v1": "${INPUTDIR}/seawifs-clim-1997-2010.1440x1080.v20180328.nc"
4445
$OCN_GRID == "tx2_3v2": "${INPUTDIR}/seawifs-clim-1997-2010-tx2_3v2.230416.nc"
46+
MLE_FL_FILE:
47+
$OCN_GRID == "tx2_3v2": "${INPUTDIR}/mle-lf-clim-tx2_3v2_20250115.nc"
4548
CFC_BC_FILE: "${DIN_LOC_ROOT}/ocn/mom/grid_indpt/cfc_atm_20230310.nc"
4649
DIAG_COORD_DEF_RHO2:
4750
$OCN_GRID == "tx2_3v2": "${INPUTDIR}/ocean_rho2_190917.nc"

0 commit comments

Comments
 (0)