Skip to content

Commit

Permalink
Fix units for unaddressable heating emissions
Browse files Browse the repository at this point in the history
  • Loading branch information
Bill-Becker committed Sep 23, 2024
1 parent 8ef48e4 commit 78b0429
Show file tree
Hide file tree
Showing 2 changed files with 2 additions and 2 deletions.
2 changes: 1 addition & 1 deletion src/results/heating_cooling_load.jl
Original file line number Diff line number Diff line change
Expand Up @@ -95,7 +95,7 @@ function add_heating_load_results(m::JuMP.AbstractModel, p::REoptInputs, d::Dict
p.s.space_heating_load.unaddressable_annual_fuel_mmbtu +
p.s.process_heat_load.unaddressable_annual_fuel_mmbtu)

r["annual_emissions_from_unaddressable_heating_load_mmbtu"] = r["annual_total_unaddressable_heating_load_mmbtu"] * p.s.existing_boiler.emissions_factor_lb_CO2_per_mmbtu * TONNE_PER_LB
r["annual_emissions_from_unaddressable_heating_load_tonnes_CO2"] = r["annual_total_unaddressable_heating_load_mmbtu"] * p.s.existing_boiler.emissions_factor_lb_CO2_per_mmbtu * TONNE_PER_LB

d["HeatingLoad"] = r
nothing
Expand Down
2 changes: 1 addition & 1 deletion test/runtests.jl
Original file line number Diff line number Diff line change
Expand Up @@ -677,7 +677,7 @@ else # run HiGHS tests
weighted_avg_addressable_fraction = (0.5 * 0.6 + 0.5 * 0.8 + 0.3 * 0.7) / (0.5 + 0.5 + 0.3)
@test round(abs(addressable / total - weighted_avg_addressable_fraction), digits=3) == 0

unaddressable_emissions = results["HeatingLoad"]["annual_emissions_from_unaddressable_heating_load_mmbtu"]
unaddressable_emissions = results["HeatingLoad"]["annual_emissions_from_unaddressable_heating_load_tonnes_CO2"]
addressable_site_fuel_emissions = results["Site"]["annual_emissions_from_fuelburn_tonnes_CO2"]
total_site_emissions = unaddressable_emissions + addressable_site_fuel_emissions
@test round(abs(addressable_site_fuel_emissions / total_site_emissions - weighted_avg_addressable_fraction), digits=3) == 0
Expand Down

0 comments on commit 78b0429

Please sign in to comment.