Skip to content

Commit

Permalink
Tune MARBL for CMIP6
Browse files Browse the repository at this point in the history 
Changes the default value of several MARBL settings parameters, and also
converts some Fortran parameters to MARBL settings parameters. Adds two new
options:

* lp_remin_scalef -- control whether particulate remin is scaled
* lo2_consumption_scalef -- control whether oxygen consumption is scaled

(For the CESM CMIP6 runs, we will scale oxygen consumption in the deep
Pacific.)

Several new diagnostics were added as well.
  • Loading branch information
@mnlevy1981
mnlevy1981 committed Jul 19, 2018
2 parents cb7576a + 56bc933 commit 3c7fa6a
Show file tree
Hide file tree
Showing 10 changed files with 770 additions and 113 deletions.
41 changes: 41 additions & 0 deletions autogenerated_src/default_diagnostics.json
View file
Original file line number Diff line number Diff line change
Expand Up @@ -1475,6 +1475,16 @@
"units": "mmol/m^3/s",
"vertical_grid": "layer_avg"
},
"O2_CONSUMPTION_SCALEF": {
"dependencies": {
"lo2_consumption_scalef": ".true."
},
"frequency": "never",
"longname": "O2 Consumption Scale Factor",
"operator": "average",
"units": 1,
"vertical_grid": "layer_avg"
},
"O2_PRODUCTION": {
"frequency": "medium",
"longname": "O2 Production",
Expand Down Expand Up @@ -1609,6 +1619,20 @@
"units": "mmol/m^3 cm/s",
"vertical_grid": "none"
},
"POC_hFLUX_IN": {
"frequency": "never",
"longname": "POC hFlux into Cell",
"operator": "average",
"units": "mmol/m^3 cm/s",
"vertical_grid": "layer_avg"
},
"POC_sFLUX_IN": {
"frequency": "never",
"longname": "POC sFlux into Cell",
"operator": "average",
"units": "mmol/m^3 cm/s",
"vertical_grid": "layer_avg"
},
"PON_REMIN_DONr": {
"frequency": "low",
"longname": "PON Remineralization routed to DONr",
Expand Down Expand Up @@ -1672,6 +1696,16 @@
"units": "cm/s",
"vertical_grid": "none"
},
"P_REMIN_SCALEF": {
"dependencies": {
"lp_remin_scalef": ".true."
},
"frequency": "never",
"longname": "Particulate Remin Scale Factor",
"operator": "average",
"units": 1,
"vertical_grid": "layer_avg"
},
"P_iron_FLUX_((particulate_flux_ref_depth_str))": {
"frequency": "medium",
"longname": "P_iron Flux at ((particulate_flux_ref_depth_str))",
Expand Down Expand Up @@ -2173,6 +2207,13 @@
"units": "mmol/m^3/s",
"vertical_grid": "layer_avg"
},
"pocToFloor": {
"frequency": "never",
"longname": "POC Flux Hitting Sea Floor",
"operator": "average",
"units": "nmol/cm^2/s",
"vertical_grid": "none"
},
"pocToSed": {
"frequency": "medium",
"longname": "POC Flux to Sediments",
Expand Down
176 changes: 149 additions & 27 deletions autogenerated_src/default_settings.json
View file
Original file line number Diff line number Diff line change
Expand Up @@ -100,7 +100,7 @@
"alphaPI_per_day": {
"datatype": "real",
"default_value": {
"((autotroph_sname)) == \"diat\"": 0.29,
"((autotroph_sname)) == \"diat\"": 0.28,
"((autotroph_sname)) == \"diaz\"": 0.39,
"((autotroph_sname)) == \"sp\"": 0.39,
"default": "1e34"
Expand All @@ -121,9 +121,9 @@
"gQfe_0": {
"datatype": "real",
"default_value": {
"((autotroph_sname)) == \"diat\"": "35e-6",
"((autotroph_sname)) == \"diaz\"": "70e-6",
"((autotroph_sname)) == \"sp\"": "35e-6",
"((autotroph_sname)) == \"diat\"": "30e-6",
"((autotroph_sname)) == \"diaz\"": "60e-6",
"((autotroph_sname)) == \"sp\"": "30e-6",
"default": "1e34"
},
"longname": "Initial Fe/C ratio for growth",
Expand All @@ -133,9 +133,9 @@
"gQfe_min": {
"datatype": "real",
"default_value": {
"((autotroph_sname)) == \"diat\"": 2.7e-06,
"((autotroph_sname)) == \"diaz\"": 5.4e-06,
"((autotroph_sname)) == \"sp\"": 2.7e-06,
"((autotroph_sname)) == \"diat\"": 2.5e-06,
"((autotroph_sname)) == \"diaz\"": 2.5e-06,
"((autotroph_sname)) == \"sp\"": 2.5e-06,
"default": "1e34"
},
"longname": "Minimum Fe/C ratio for growth",
Expand Down Expand Up @@ -167,7 +167,7 @@
"kFe": {
"datatype": "real",
"default_value": {
"((autotroph_sname)) == \"diat\"": 6e-05,
"((autotroph_sname)) == \"diat\"": 7e-05,
"((autotroph_sname)) == \"diaz\"": 4.5e-05,
"((autotroph_sname)) == \"sp\"": 3e-05,
"default": "1e34"
Expand Down Expand Up @@ -403,7 +403,7 @@
"graze_poc": {
"datatype": "real",
"default_value": {
"((grazing_sname)) == \"diat_zoo\"": 0.38,
"((grazing_sname)) == \"diat_zoo\"": 0.39,
"((grazing_sname)) == \"diaz_zoo\"": 0.1,
"((grazing_sname)) == \"sp_zoo\"": 0,
"default": "1e34"
Expand Down Expand Up @@ -471,8 +471,8 @@
"z_umax_0_per_day": {
"datatype": "real",
"default_value": {
"((grazing_sname)) == \"diat_zoo\"": 3.1,
"((grazing_sname)) == \"diaz_zoo\"": 3.25,
"((grazing_sname)) == \"diat_zoo\"": 3.15,
"((grazing_sname)) == \"diaz_zoo\"": 3.3,
"((grazing_sname)) == \"sp_zoo\"": 3.3,
"default": "1e34"
},
Expand Down Expand Up @@ -780,13 +780,27 @@
"user-specified"
]
},
"POM_bury_frac_max": {
"datatype": "real",
"default_value": 0.8,
"longname": "maximum bury fraction for POM",
"subcategory": "4. general parameters (bury coeffs)",
"units": "unitless"
},
"PON_bury_coeff": {
"datatype": "real",
"default_value": 0.5,
"longname": "Scale factor for burial of PON",
"subcategory": "4. general parameters (bury coeffs)",
"units": "unitless"
},
"bSi_bury_frac_max": {
"datatype": "real",
"default_value": 1.0,
"longname": "maximum bury fraction for bSi",
"subcategory": "4. general parameters (bury coeffs)",
"units": "unitless"
},
"bury_coeff_rmean_timescale_years": {
"datatype": "real",
"default_value": 10,
Expand All @@ -803,7 +817,7 @@
},
"caco3_bury_thres_omega_calc": {
"datatype": "real",
"default_value": 1.0,
"default_value": 0.89,
"longname": "omega calcite threshold for CaCO3 burial when opt = 'omega_calc'",
"subcategory": "4. general parameters",
"units": 1
Expand Down Expand Up @@ -847,6 +861,41 @@
"subcategory": "2. config flags",
"units": "unitless"
},
"gQ_Fe_kFe_thres": {
"datatype": "real",
"default_value": 10.0,
"longname": "Fe:kFe ratio threshold in uptake ratio computations",
"subcategory": "4. general parameters",
"units": 1
},
"gQ_Si_kSi_thres": {
"datatype": "real",
"default_value": 6.0,
"longname": "Si:kSi ratio threshold in uptake ratio computations",
"subcategory": "4. general parameters",
"units": 1
},
"gQsi_0": {
"datatype": "real",
"default_value": 0.137,
"longname": "initial Si/C ratio for growth",
"subcategory": "4. general parameters",
"units": 1
},
"gQsi_max": {
"datatype": "real",
"default_value": 0.822,
"longname": "max Si/C ratio for growth",
"subcategory": "4. general parameters",
"units": 1
},
"gQsi_min": {
"datatype": "real",
"default_value": 0.0457,
"longname": "min Si/C ratio for growth",
"subcategory": "4. general parameters",
"units": 1
},
"ladjust_bury_coeff": {
"_append_to_config_keywords": true,
"datatype": "logical",
Expand Down Expand Up @@ -883,6 +932,23 @@
"subcategory": "2. config flags",
"units": "unitless"
},
"lo2_consumption_scalef": {
"datatype": "logical",
"default_value": {
"GRID == \"CESM_x1\"": ".true.",
"default": ".false."
},
"longname": "Apply o2_consumption_scalef to o2 consumption (and request it as a forcing)",
"subcategory": "2. config flags",
"units": "unitless"
},
"lp_remin_scalef": {
"datatype": "logical",
"default_value": ".false.",
"longname": "Apply p_remin_scalef to particulate remin (and request it as a forcing)",
"subcategory": "2. config flags",
"units": "unitless"
},
"lsource_sink": {
"datatype": "logical",
"default_value": ".true.",
Expand All @@ -897,16 +963,44 @@
"subcategory": "2. config flags",
"units": "unitless"
},
"o2_sf_o2_range_hi": {
"datatype": "real",
"default_value": 45.0,
"longname": "o2_scalefactor is applied to diss length scales for O2 less than this",
"subcategory": "4. general parameters",
"units": "mmol/m^3"
},
"o2_sf_o2_range_lo": {
"datatype": "real",
"default_value": 5.0,
"longname": "o2_scalefactor is constant for O2 less than this",
"subcategory": "4. general parameters",
"units": "mmol/m^3"
},
"o2_sf_val_lo_o2": {
"datatype": "real",
"default_value": 2.6,
"longname": "o2_scalefactor constant for O2 less than o2_sf_o2_range_lo",
"subcategory": "4. general parameters",
"units": 1
},
"parm_CaCO3_diss": {
"datatype": "real",
"default_value": "500e2",
"longname": "Base CaCO3 dissipation length scale",
"subcategory": "4. general parameters (dissipation)",
"longname": "Base CaCO3 dissolution length scale",
"subcategory": "4. general parameters (dissolution)",
"units": "cm"
},
"parm_CaCO3_gamma": {
"datatype": "real",
"default_value": 0.02,
"longname": "CaCO3 gamma (fraction of production -> hard subclass)",
"subcategory": "4. general parameters",
"units": 1
},
"parm_FeLig_scavenge_rate0": {
"datatype": "real",
"default_value": 1.4,
"default_value": 1.2,
"longname": "Scavenging base rate for bound iron",
"subcategory": "4. general parameters (scavenging)",
"units": "unitless"
Expand All @@ -927,7 +1021,7 @@
},
"parm_Fe_scavenge_rate0": {
"datatype": "real",
"default_value": 18.0,
"default_value": 22.0,
"longname": "Scavenging base rate for Fe",
"subcategory": "4. general parameters (scavenging)",
"units": "unitless"
Expand All @@ -949,17 +1043,24 @@
"parm_POC_diss": {
"datatype": "real",
"default_value": "100e2",
"longname": "Base POC dissipation length scale",
"subcategory": "4. general parameters (dissipation)",
"longname": "Base POC dissolution length scale",
"subcategory": "4. general parameters (dissolution)",
"units": "cm"
},
"parm_SiO2_diss": {
"datatype": "real",
"default_value": "770e2",
"longname": "Base SiO2 dissipation length scale",
"subcategory": "4. general parameters (dissipation)",
"default_value": "650e2",
"longname": "Base SiO2 dissolution length scale",
"subcategory": "4. general parameters (dissolution)",
"units": "cm"
},
"parm_SiO2_gamma": {
"datatype": "real",
"default_value": 0.0,
"longname": "SiO2 gamma (fraction of production -> hard subclass)",
"subcategory": "4. general parameters",
"units": 1
},
"parm_f_prod_sp_CaCO3": {
"datatype": "real",
"default_value": {
Expand All @@ -970,11 +1071,32 @@
"subcategory": "4. general parameters",
"units": "unitless"
},
"parm_hPOC_CaCO3_ratio": {
"datatype": "real",
"default_value": 0.01,
"longname": "hPOC to CaCO3 ratio",
"subcategory": "4. general parameters",
"units": 1
},
"parm_hPOC_SiO2_ratio": {
"datatype": "real",
"default_value": 0.01,
"longname": "hPOC to SiO2 ratio",
"subcategory": "4. general parameters",
"units": 1
},
"parm_hPOC_dust_ratio": {
"datatype": "real",
"default_value": 0.01,
"longname": "hPOC to dust ratio",
"subcategory": "4. general parameters",
"units": 1
},
"parm_init_POC_bury_coeff": {
"datatype": "real",
"default_value": {
"GRID == \"CESM_x3\"": 5.0,
"default": 1.1
"default": 2.54
},
"longname": "Initial scale factor for burial of POC, PON",
"subcategory": "4. general parameters (bury coeffs)",
Expand All @@ -984,7 +1106,7 @@
"datatype": "real",
"default_value": {
"GRID == \"CESM_x3\"": 5.0,
"default": 1.1
"default": 0.36
},
"longname": "Initial scale factor for burial of POP",
"subcategory": "4. general parameters (bury coeffs)",
Expand All @@ -994,7 +1116,7 @@
"datatype": "real",
"default_value": {
"GRID == \"CESM_x3\"": 1.1,
"default": 1.0
"default": 1.53
},
"longname": "Initial scale factor for burial of bSi",
"subcategory": "4. general parameters (bury coeffs)",
Expand Down Expand Up @@ -1047,9 +1169,9 @@
],
"default": [
1,
3.0,
4.5,
5.5
3.6,
4.7,
4.8
]
},
"longname": "Prescribed scalelen values",
Expand Down
Loading

0 comments on commit 3c7fa6a

Please sign in to comment.