Skip to content

Allow tau OH H2 coefficient to vary with CO coefficient #787

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Merged
kdorheim merged 10 commits into from
Aug 18, 2025
Merged

Allow tau OH H2 coefficient to vary with CO coefficient #787

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
10 commits
Select commit Hold shift + click to select a range
57f1377
remove un needed interpolation call
kdorheim Jul 29, 2025
88b8ce6
simplify the zero H2 emissions
kdorheim Jul 29, 2025
3ed32d5
rho_h2o_h2 coefficent
kdorheim Jul 29, 2025
c6e0314
add H2 coeff for O3 RF to all input files
kdorheim Jul 29, 2025
b9abcd2
passing the ini tests
kdorheim Jul 29, 2025
f81bc3d
saving a wip, passing ini and h2 tests
kdorheim Jul 29, 2025
777954c
add coeff dervation notation
kdorheim Jul 30, 2025
9edc73f
spelling fix
kdorheim Aug 6, 2025
3f5a608
add doi
kdorheim Aug 6, 2025
2269895
change ini commnet
kdorheim Aug 6, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
495 changes: 249 additions & 246 deletions data-raw/input_params.csv
View file
Open in desktop

Large diffs are not rendered by default.

Binary file modified data/inputstable.rda
View file
Open in desktop
Binary file not shown.
2 changes: 1 addition & 1 deletion inst/include/component_data.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -285,7 +285,7 @@
#define D_COEFFICIENT_CH4 "CCH4"
#define D_COEFFICIENT_NMVOC "CNMVOC"
#define D_COEFFICIENT_CO "CCO"
#define D_COEFFICIENT_H2 "CH2"
#define D_COEFFICIENT_H2 "H2_CCO"
#define D_EMISSIONS_H2 "H2_emissions"


Expand Down
2 changes: 1 addition & 1 deletion inst/include/oh_component.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -71,7 +71,7 @@ class OHComponent : public IModelComponent {
double CNMVOC; // coefficient for NMVOC
double CNOX; // coefficient for NOX
double CCH4; // coefficient for CH4
unitval CH2; // coefficient for CH4
unitval H2_CCO; // adjust CO coefficient for H2


// logger
Expand Down
3 changes: 2 additions & 1 deletion inst/input/hector_picontrol.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,6 +130,7 @@ CCH4=-0.32 ; coefficent for CH4
;------------------------------------------------------------------------
[ozone]
PO3=30.0 ; preindustrial O3 concentration
CO3_H2=4.09e-3 ; coefficient for H2 contribution to trop. O3 RF (DU O3 Tg H2-1) SI Table 3 of Sand et al. 2023
NOX_emissions[1745]=3.877282 ; emissions time series
CO_emissions[1745]=348.5273588 ; emissions time series
NMVOC_emissions[1745]= 60.02182622 ; emissions time series
Expand Down Expand Up @@ -162,7 +163,7 @@ rho_bc=0.0508 ; (W yr m–2 C Tg–1) IPCC AR6 radiative efficiency BC
rho_oc=-.00621 ; (W yr m–2 C Tg–1) IPCC AR6 radiative efficiency OC (7.SM.1.3.1 of IPCC AR6)
rho_so2=-.00000724 ; (W yr m–2 S Gg-1) IPCC AR6 radiative efficiency SO2 (7.SM.1.3.1 of IPCC AR6)
rho_nh3=-.00208 ; (W yr m–2 NH3Tg–1) IPCC AR6 radiative efficiency NH3 (7.SM.1.3.1 of IPCC AR6)
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg−1) from Sand et al. 2023
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg-1) Sand et al. 2023, 10.1038/s43247-023-00857-8

; Miscellaneous radiative forcings default set to 0, or read in from a input table, may be used to read in
; additional forcings not modeled by Hector (i.e. solar, bc on snow , contrails from)
Expand Down
9 changes: 6 additions & 3 deletions inst/input/hector_ssp119.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,22 +130,25 @@ CH4_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time se
NOX_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
H2_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs


TOH0=6.6 ; inital OH lifetime (years)
CNOX=0.0042 ; coefficient for NOX
CCO=-0.000105 ; coefficient for CO
CNMVOC=-0.000315 ; coefficient for NMVOC (non methane VOC)
CCH4=-0.32 ; coefficient for CH4
CH2=-0.00044625 ; coefficient for H2 TODO-H2
H2_CCO=4.25 ; scaling factor for the CO coefficient that affects OH lifetime; see PR787

;------------------------------------------------------------------------
[ozone]
PO3=30.0 ; preindustrial O3 concentration
CO3_H2=4.09e-3 ; coefficient for H2 contribution to trop. O3 RF (DU O3 Tg H2-1) SI Table 3 of Sand et al. 2023
NOX_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs


;------------------------------------------------------------------------
[N2O]
Expand Down Expand Up @@ -174,7 +177,7 @@ rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg−1) from Sand et al. 2023
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg-1) Sand et al. 2023, 10.1038/s43247-023-00857-8


; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
Expand Down
8 changes: 5 additions & 3 deletions inst/input/hector_ssp126.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,22 +130,24 @@ CH4_emissions=csv:tables/ssp126_emiss-constraints_rf.csv ; emissions time se
NOX_emissions=csv:tables/ssp126_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp126_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp126_emiss-constraints_rf.csv ; emissions time series
H2_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs

TOH0=6.6 ; inital OH lifetime (years)
CNOX=0.0042 ; coefficient for NOX
CCO=-0.000105 ; coefficient for CO
CNMVOC=-0.000315 ; coefficient for NMVOC (non methane VOC)
CCH4=-0.32 ; coefficient for CH4
CH2=-0.00044625 ; coefficient for H2 TODO-H2
H2_CCO=4.25 ; scaling factor for the CO coefficient that affects OH lifetime; see PR787


;------------------------------------------------------------------------
[ozone]
PO3=30.0 ; preindustrial O3 concentration
CO3_H2=4.09e-3 ; coefficient for H2 contribution to trop. O3 RF (DU O3 Tg H2-1) SI Table 3 of Sand et al. 2023
NOX_emissions=csv:tables/ssp126_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp126_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp126_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs

;------------------------------------------------------------------------
[N2O]
Expand Down Expand Up @@ -174,7 +176,7 @@ rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg−1) from Sand et al. 2023
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg-1) Sand et al. 2023, 10.1038/s43247-023-00857-8


; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
Expand Down
12 changes: 6 additions & 6 deletions inst/input/hector_ssp245.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,25 +130,25 @@ CH4_emissions=csv:tables/ssp245_emiss-constraints_rf.csv ; emissions time se
NOX_emissions=csv:tables/ssp245_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp245_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp245_emiss-constraints_rf.csv ; emissions time series
H2_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs


TOH0=6.6 ; inital OH lifetime (years)
CNOX=0.0042 ; coefficient for NOX
CCO=-0.000105 ; coefficient for CO
CNMVOC=-0.000315 ; coefficient for NMVOC (non methane VOC)
CCH4=-0.32 ; coefficient for CH4
CH2=-0.00044625 ; coefficient for H2 TODO-H2
H2_CCO=4.25 ; scaling factor for the CO coefficient that affects OH lifetime; see PR787


;------------------------------------------------------------------------
[ozone]
PO3=30.0 ; preindustrial O3 concentration
CO3_H2=4.09e-3 ; coefficient for H2 TODO-H2
PO3=30.0 ; preindustrial O3 concentration
CO3_H2=4.09e-3 ; coefficient for H2 contribution to trop. O3 RF (DU O3 Tg H2-1) SI Table 3 of Sand et al. 2023
NOX_emissions=csv:tables/ssp245_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp245_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp245_emiss-constraints_rf.csv ; emissions time series
H2_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs


;------------------------------------------------------------------------
Expand Down Expand Up @@ -178,7 +178,7 @@ rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg−1) from Sand et al. 2023
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg-1) Sand et al. 2023, 10.1038/s43247-023-00857-8



Expand Down
8 changes: 5 additions & 3 deletions inst/input/hector_ssp370.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,23 +130,25 @@ CH4_emissions=csv:tables/ssp370_emiss-constraints_rf.csv ; emissions time se
NOX_emissions=csv:tables/ssp370_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp370_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp370_emiss-constraints_rf.csv ; emissions time series
H2_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs


TOH0=6.6 ; inital OH lifetime (years)
CNOX=0.0042 ; coefficient for NOX
CCO=-0.000105 ; coefficient for CO
CNMVOC=-0.000315 ; coefficient for NMVOC (non methane VOC)
CCH4=-0.32 ; coefficient for CH4
CH2=-0.00044625 ; coefficient for H2 TODO-H2
H2_CCO=4.25 ; scaling factor for the CO coefficient that affects OH lifetime; see PR787


;------------------------------------------------------------------------
[ozone]
PO3=30.0 ; preindustrial O3 concentration
CO3_H2=4.09e-3 ; coefficient for H2 contribution to trop. O3 RF (DU O3 Tg H2-1) SI Table 3 of Sand et al. 2023, 10.1038/s43247-023-00857-8
NOX_emissions=csv:tables/ssp370_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp370_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp370_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs

;------------------------------------------------------------------------
[N2O]
Expand Down Expand Up @@ -175,7 +177,7 @@ rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg−1) from Sand et al. 2023
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg-1) Sand et al. 2023, 10.1038/s43247-023-00857-8


; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
Expand Down
8 changes: 5 additions & 3 deletions inst/input/hector_ssp434.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,23 +130,25 @@ CH4_emissions=csv:tables/ssp434_emiss-constraints_rf.csv ; emissions time se
NOX_emissions=csv:tables/ssp434_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp434_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp434_emiss-constraints_rf.csv ; emissions time series
H2_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs


TOH0=6.6 ; inital OH lifetime (years)
CNOX=0.0042 ; coefficient for NOX
CCO=-0.000105 ; coefficient for CO
CNMVOC=-0.000315 ; coefficient for NMVOC (non methane VOC)
CCH4=-0.32 ; coefficient for CH4
CH2=-0.00044625 ; coefficient for H2 TODO-H2
H2_CCO=4.25 ; scaling factor for the CO coefficient that affects OH lifetime; see PR787


;------------------------------------------------------------------------
[ozone]
PO3=30.0 ; preindustrial O3 concentration
CO3_H2=4.09e-3 ; coefficient for H2 contribution to trop. O3 RF (DU O3 Tg H2-1) SI Table 3 of Sand et al. 2023
NOX_emissions=csv:tables/ssp434_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp434_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp434_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs

;------------------------------------------------------------------------
[N2O]
Expand Down Expand Up @@ -175,7 +177,7 @@ rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg−1) from Sand et al. 2023
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg-1) Sand et al. 2023, 10.1038/s43247-023-00857-8


; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
Expand Down
8 changes: 5 additions & 3 deletions inst/input/hector_ssp460.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,23 +130,25 @@ CH4_emissions=csv:tables/ssp460_emiss-constraints_rf.csv ; emissions time se
NOX_emissions=csv:tables/ssp460_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp460_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp460_emiss-constraints_rf.csv ; emissions time series
H2_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs


TOH0=6.6 ; inital OH lifetime (years)
CNOX=0.0042 ; coefficient for NOX
CCO=-0.000105 ; coefficient for CO
CNMVOC=-0.000315 ; coefficient for NMVOC (non methane VOC)
CCH4=-0.32 ; coefficient for CH4
CH2=-0.00044625 ; coefficient for H2 TODO-H2
H2_CCO=4.25 ; scaling factor for the CO coefficient that affects OH lifetime; see PR787


;------------------------------------------------------------------------
[ozone]
PO3=30.0 ; preindustrial O3 concentration
CO3_H2=4.09e-3 ; coefficient for H2 contribution to trop. O3 RF (DU O3 Tg H2-1) SI Table 3 of Sand et al. 2023
NOX_emissions=csv:tables/ssp460_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp460_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp460_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs

;------------------------------------------------------------------------
[N2O]
Expand Down Expand Up @@ -175,7 +177,7 @@ rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg−1) from Sand et al. 2023
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg-1) Sand et al. 2023, 10.1038/s43247-023-00857-8


; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
Expand Down
8 changes: 5 additions & 3 deletions inst/input/hector_ssp534-over.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,23 +130,25 @@ CH4_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv ; emissions ti
NOX_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv ; emissions time series
H2_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs


TOH0=6.6 ; inital OH lifetime (years)
CNOX=0.0042 ; coefficient for NOX
CCO=-0.000105 ; coefficient for CO
CNMVOC=-0.000315 ; coefficient for NMVOC (non methane VOC)
CCH4=-0.32 ; coefficient for CH4
CH2=-0.00044625 ; coefficient for H2 TODO-H2
H2_CCO=4.25 ; scaling factor for the CO coefficient that affects OH lifetime; see PR787


;------------------------------------------------------------------------
[ozone]
PO3=30.0 ; preindustrial O3 concentration
CO3_H2=4.09e-3 ; coefficient for H2 contribution to trop. O3 RF (DU O3 Tg H2-1) SI Table 3 of Sand et al. 2023
NOX_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp534-over_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs

;------------------------------------------------------------------------
[N2O]
Expand Down Expand Up @@ -175,7 +177,7 @@ rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg−1) from Sand et al. 2023
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg-1) Sand et al. 2023, 10.1038/s43247-023-00857-8


; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
Expand Down
8 changes: 5 additions & 3 deletions inst/input/hector_ssp585.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -130,22 +130,24 @@ CH4_emissions=csv:tables/ssp585_emiss-constraints_rf.csv ; emissions time se
NOX_emissions=csv:tables/ssp585_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp585_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp585_emiss-constraints_rf.csv ; emissions time series
H2_emissions=csv:tables/ssp119_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs

TOH0=6.6 ; inital OH lifetime (years)
CNOX=0.0042 ; coefficient for NOX
CCO=-0.000105 ; coefficient for CO
CNMVOC=-0.000315 ; coefficient for NMVOC (non methane VOC)
CCH4=-0.32 ; coefficient for CH4
CH2=-0.00044625 ; coefficient for H2 TODO-H2
H2_CCO=4.25 ; scaling factor for the CO coefficient that affects OH lifetime; see PR787


;------------------------------------------------------------------------
[ozone]
PO3=30.0 ; preindustrial O3 concentration
CO3_H2=4.09e-3 ; coefficient for H2 contribution to trop. O3 RF (DU O3 Tg H2-1) SI Table 3 of Sand et al. 2023
NOX_emissions=csv:tables/ssp585_emiss-constraints_rf.csv ; emissions time series
CO_emissions=csv:tables/ssp585_emiss-constraints_rf.csv ; emissions time series
NMVOC_emissions=csv:tables/ssp585_emiss-constraints_rf.csv ; emissions time series
H2_emissions[1745]=0 ; assumed 0 emissions for the RCMIP era SSPs

;------------------------------------------------------------------------
[N2O]
Expand Down Expand Up @@ -174,7 +176,7 @@ rho_bc=0.06386286 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) Dorheim et al. 2024
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) Dorheim et al. 2024
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) Dorheim et al. 2024
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg−1) from Sand et al. 2023
rho_h2o_h2=0.00019 ; (W m −2 H2 Tg-1) Sand et al. 2023, 10.1038/s43247-023-00857-8


; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
Expand Down
Loading