Skip to content
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

Jump to bottom

Update Aerosols #724

Merged
merged 3 commits into from
Jan 10, 2024
Merged

Update Aerosols #724

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
3 commits
Select commit Hold shift + click to select a range
aab43a8
New Aerosol factors
kdorheim Jan 9, 2024
d09a494
update old new comparison data
kdorheim Jan 10, 2024
5e6426d
Fix typo
kdorheim Jan 10, 2024
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
2 changes: 1 addition & 1 deletion inst/include/forcing_component.hpp
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -120,7 +120,7 @@ class ForcingComponent : public IModelComponent {

// Aerosol parameters for aerosol-cloud interactions (RFaci) see equation
// Equation 7.SM.1.2 of IPCC AR6
double const aci_beta = 2.09841432; // (W m-2) IPCC AR6 7.SM.1.3.1
double const aci_beta = 2.279759; // see Dorheim et al. in prep
kdorheim marked this conversation as resolved.
Show resolved Hide resolved
double const s_BCOC = 111.05064063; // (Tg C yr-1) IPCC AR6 7.SM.1.3.1
double const s_SO2 =
(260.34644166 * 1000) *
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp119.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year
delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6)
delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6)
delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6)
rho_bc=0.0508 ; (W yr m–2 C Tg–1) IPCC AR6 radiative efficiency BC (7.SM.1.3 of IPCC AR6)
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_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep

; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails)
Expand Down
14 changes: 7 additions & 7 deletions inst/input/hector_ssp126.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -37,7 +37,7 @@ preind_interdeep_c=37100 ; intermediate and deep ocean C, IPCC AR6 Figure 5.12
atmos_co2=590.33 ; initial atmospheric carbon pool Pg C of CO2 size used in spin up, consistent with C0
C0=277.15 ; preindustrial CO2 from table 7.SM.1 IPCC AR6; CO2 radiative forcing is calibrated to this value
npp_flux0=56.2 ; preindustrial net primary production, Pg C/yr from Ito et al. 2011
;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978
;pf.npp_flux0=5.6 ; cf. Figure 3 in Chen et al. 2017 10.1088/1748-9326/aa8978
;nonpf.npp_flux0=50.6

veg_c=550 ; initial vegetation pool, Pg C
Expand Down Expand Up @@ -73,7 +73,7 @@ f_litterd=0.98 ; Fraction of litter flux to detritus (balance to soil)
;nonpf.f_litterd=0.98

; Anthropogenic contributions: emissions, direct air carbon capture and storage, and land use change, Pg C/yr
ffi_emissions=csv:tables/ssp126_emiss-constraints_rf.csv
ffi_emissions=csv:tables/ssp126_emiss-constraints_rf.csv
daccs_uptake=csv:tables/ssp126_emiss-constraints_rf.csv
luc_emissions=csv:tables/ssp126_emiss-constraints_rf.csv
luc_uptake=csv:tables/ssp126_emiss-constraints_rf.csv
Expand Down Expand Up @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year
delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6)
delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6)
delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6)
rho_bc=0.0508 ; (W yr m–2 C Tg–1) IPCC AR6 radiative efficiency BC (7.SM.1.3 of IPCC AR6)
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_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep

; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails)
Expand Down Expand Up @@ -232,7 +232,7 @@ delta_HFC32=0.0 ; tropospheric adjustments scalar unitless (7.3.2.4 IPCC AR6
molarMass=52.0 ; grams
HFC32_emissions=csv:tables/ssp126_emiss-constraints_rf.csv

[HFC4310_halocarbon]
[HFC4310_halocarbon]
; This is equivalent to HFC-43-10mee
tau=17.0 ; lifetime in years (Table 7.SM.7 IPPC AR6)
rho_HFC4310=0.000357 ; radiative efficiencies W/m2/ppt (Table 7.SM.7 IPPC AR6)
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp245.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year
delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6)
delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6)
delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6)
rho_bc=0.0508 ; (W yr m–2 C Tg–1) IPCC AR6 radiative efficiency BC (7.SM.1.3 of IPCC AR6)
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_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep

; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails)
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp370.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year
delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6)
delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6)
delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6)
rho_bc=0.0508 ; (W yr m–2 C Tg–1) IPCC AR6 radiative efficiency BC (7.SM.1.3 of IPCC AR6)
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_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep

; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails)
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp434.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year
delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6)
delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6)
delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6)
rho_bc=0.0508 ; (W yr m–2 C Tg–1) IPCC AR6 radiative efficiency BC (7.SM.1.3 of IPCC AR6)
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_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep

; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails)
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp460.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year
delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6)
delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6)
delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6)
rho_bc=0.0508 ; (W yr m–2 C Tg–1) IPCC AR6 radiative efficiency BC (7.SM.1.3 of IPCC AR6)
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_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep

; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails)
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp534-over.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year
delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6)
delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6)
delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6)
rho_bc=0.0508 ; (W yr m–2 C Tg–1) IPCC AR6 radiative efficiency BC (7.SM.1.3 of IPCC AR6)
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_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep

; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails)
Expand Down
8 changes: 4 additions & 4 deletions inst/input/hector_ssp585.ini
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -165,10 +165,10 @@ baseyear=1750 ; when to start reporting; by definition, all F=0 in this year
delta_co2=0.05 ; (unitless) forcing tropospheric adjustments for CO2 (7.3.2.1 of IPCC AR6)
delta_ch4=-.14 ; (unitless) forcing tropospheric adjustments for CH4 (7.3.2.2 of IPCC AR6)
delta_n2o=0.07 ; (unitless) forcing tropospheric adjustments for N2O (7.3.2.3 of IPCC AR6)
rho_bc=0.0508 ; (W yr m–2 C Tg–1) IPCC AR6 radiative efficiency BC (7.SM.1.3 of IPCC AR6)
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_bc=0.06386286 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_oc=-0.006407143 ; (W yr m–2 C Tg–1) see Dorheim et al. in prep
rho_so2=-7.469841e-06 ; (W yr m–2 S Gg-1) see Dorheim et al. in prep
rho_nh3=-0.002146032 ; (W yr m–2 NH3Tg–1) see Dorheim et al. in prep

; Miscellaneous radiative forcings are by default zero, but can read in from a input table to
; represent additional forcings not modeled by Hector (e.g. solar, bc on snow, jet contrails)
Expand Down
Loading
Loading