@@ -247,6 +247,9 @@ def write_MARBL_diagnostics_file(
247247 full_diag_dict ["O2" ]["properties" ]["has surface flux" ] = True
248248 # DIC
249249 if "DIC" in full_diag_dict .keys ():
250+ # STF_SALT_DIC is special case
251+ if valid_diag_modes .index (diag_mode ) >= valid_diag_modes .index ("minimal" ):
252+ full_diag_dict ["DIC" ]["diags" ]["STF_SALT_DIC" ] = "medium_average"
250253 full_diag_dict ["DIC" ]["diags" ]["DIC_RIV_FLUX" ] = freq_op
251254 full_diag_dict ["DIC" ]["diags" ]["J_DIC" ] = freq_op
252255 full_diag_dict ["DIC" ]["diags" ]["Jint_100m_DIC" ] = Jint_100m_freq_op ["DIC" ]
@@ -265,6 +268,9 @@ def write_MARBL_diagnostics_file(
265268 full_diag_dict ["DIC_ALT_CO2" ]["properties" ]["has surface flux" ] = True
266269 # ALK
267270 if "ALK" in full_diag_dict .keys ():
271+ # STF_SALT_ALK is special case
272+ if valid_diag_modes .index (diag_mode ) >= valid_diag_modes .index ("minimal" ):
273+ full_diag_dict ["ALK" ]["diags" ]["STF_SALT_ALK" ] = "medium_average"
268274 full_diag_dict ["ALK" ]["diags" ]["ALK_RIV_FLUX" ] = freq_op
269275 full_diag_dict ["ALK" ]["diags" ]["STF_ALK" ] = freq_op
270276 full_diag_dict ["ALK" ]["diags" ]["J_ALK" ] = low_freq_op
@@ -353,6 +359,9 @@ def write_MARBL_diagnostics_file(
353359 full_diag_dict ["DO14Ctot" ]["properties" ]["has surface flux" ] = True
354360 # ABIO_DIC
355361 if "ABIO_DIC" in full_diag_dict .keys ():
362+ # STF_SALT_ABIO_DIC is special case
363+ if valid_diag_modes .index (diag_mode ) >= valid_diag_modes .index ("minimal" ):
364+ full_diag_dict ["ABIO_DIC" ]["diags" ]["STF_SALT_ABIO_DIC" ] = "medium_average"
356365 # full_diag_dict['ABIO_DIC']['diags']['J_ABIO_DIC'] = freq_op
357366 full_diag_dict ["ABIO_DIC" ]["diags" ]["STF_ABIO_DIC" ] = freq_op
358367 # full_diag_dict['ABIO_DIC']['diags']['FvPER_ABIO_DIC'] = freq_op
0 commit comments