openstudio-standards 0.8.3 → 0.8.5.rc1

data/lib/openstudio-standards/standards/necb/common/btap_data.rb

@@ -5,10 +5,22 @@ class BTAPData
   attr_accessor :sqlite_file
   attr_accessor :btap_data
 
-  def initialize(model:, runner: nil, cost_result:, baseline_cost_equipment_total_cost_per_m_sq: -1.0,
-                 baseline_cost_utility_neb_total_cost_per_m_sq: -1.0, baseline_energy_eui_total_gj_per_m_sq: -1.0, qaqc:,
-                 npv_start_year: 2022, npv_end_year: 2041, npv_discount_rate: 0.03, npv_discount_rate_carbon: 0.03, oerd_utility_pricing: nil,
-                 utility_pricing_year: 2020)
+  def initialize(
+    model:, 
+    runner: nil, 
+    cost_result:, 
+    carbon_result:,
+    baseline_cost_equipment_total_cost_per_m_sq: -1.0,
+    baseline_cost_utility_neb_total_cost_per_m_sq: -1.0, 
+    baseline_energy_eui_total_gj_per_m_sq: -1.0, 
+    qaqc:,
+    npv_start_year: 2022, 
+    npv_end_year: 2041, 
+    npv_discount_rate: 0.03, 
+    npv_discount_rate_carbon: 0.03, 
+    oerd_utility_pricing: nil,
+    utility_pricing_year: 2020)
+
     @model = model
     @error_warning = []
     # sets sql file.
@@ -20,6 +32,8 @@ class BTAPData
     @neb_prices_csv_file_name = File.join(__dir__, 'neb_end_use_prices.csv')
     @neb_prices_csv_file_name = File.join(__dir__, 'utility_pricing_2025-02-20.csv') if oerd_utility_pricing.to_bool
     @necb_reference_runs_csv_file_name = File.join(__dir__, 'necb_reference_runs.csv')
+    @eccc_ghg_electricity = File.join(__dir__, 'eccc_electric_grid_intensity_20250311.csv') #REF (Retrieved March 11, 2025) for ECCC's emissions factors for the current and future years for GHG "without Biomass and RNG CO2 emissions": https://data-donnees.az.ec.gc.ca/data/substances/monitor/canada-s-greenhouse-gas-emissions-projections/Current-Projections-Actuelles/Energy-Energie/AM%20Scenario%20AMS/Grid-O%26G-Intensities-Intensites-Reseau-Delectricite-P%26G?lang=en
+    @nir_ghg_gas = File.join(__dir__, 'nir_gas_grid_intensity_20250311.csv') #REF (Retrieved March 11, 2025) for Year 2022 of "Table A6.1–1 CO2 Emission Factors for Marketable Natural Gas" of National Inventory Report: https://data-donnees.az.ec.gc.ca/api/file?path=%2Fsubstances%2Fmonitor%2Fcanada-s-official-greenhouse-gas-inventory%2FD-Emission-Factors%2FEN_Annex6_Emission_Factors.pdf
 
     # Conditioned floor area is used so much. May as well make it a object variable.
     # setup the queries
@@ -159,8 +173,8 @@ class BTAPData
 
   def building_costing_data(cost_result)
     building_data = {}
-    building_data['cost_rs_means_prov'] = cost_result['rs_means_prov']
-    building_data['cost_rs_means_city'] = cost_result['rs_means_city']
+    building_data['cost_province_state'] = cost_result['province_state']
+    building_data['cost_city'] = cost_result['city']
     building_data['cost_equipment_envelope_total_cost_per_m_sq'] = (cost_result['totals']['envelope']) / @conditioned_floor_area_m_sq
     building_data['cost_equipment_thermal_bridging_total_cost_per_m_sq'] = (cost_result['totals']['thermal_bridging']) / @conditioned_floor_area_m_sq
     building_data['cost_equipment_lighting_total_cost_per_m_sq'] = (cost_result['totals']['lighting']) / @conditioned_floor_area_m_sq
@@ -228,6 +242,9 @@ class BTAPData
     # Create a hash of the neb data.
     neb_data = CSV.parse(File.read(@neb_prices_csv_file_name), headers: true, converters: :numeric).map(&:to_h)
 
+    # Create a hash of the ECCC's electricity emissions factors for the current and future years for GHG "without Biomass and RNG CO2 emissions".
+    eccc_ghg_elec_emission_factors = CSV.parse(File.read(@eccc_ghg_electricity), headers: true, converters: :numeric).map(&:to_h)
+
     # Find which province the proposed building is located in
     building_type = 'Commercial'
     geography_data = climate_data
@@ -245,17 +262,14 @@ class BTAPData
     onsite_elec_generation = @btap_data['total_site_eui_gj_per_m_sq'] - @btap_data['net_site_eui_gj_per_m_sq']
     if onsite_elec_generation > 0.0
       eui_elec = @btap_data['energy_eui_electricity_gj_per_m_sq'] - onsite_elec_generation
-      # Calculate the saved GHG emissions from onsite electricity generatation (assume it does not emmit GHGs)
-      ghg_saved = get_utility_ghg_kg_per_gj(province: province_abbreviation, fuel_type: "Electricity") * onsite_elec_generation
-      # Convert to tonnes of ghg/m^2
-      ghg_elec = (@btap_data['cost_utility_ghg_electricity_kg_per_m_sq'].to_f - ghg_saved.to_f)/1000.0
     else
       eui_elec = @btap_data['energy_eui_electricity_gj_per_m_sq']
-      # Convert to tonnes of ghg/m^2
-      ghg_elec = @btap_data['cost_utility_ghg_electricity_kg_per_m_sq'].to_f/1000.0
     end
-    # puts "onsite_elec_generation is #{onsite_elec_generation}"
-    # puts "eui_elec is #{eui_elec}"
+    # Find the appropriate row for electricity GHG emmission factors
+    row_emm_factor = eccc_ghg_elec_emission_factors.detect do |data|
+      (data['province'] == province)
+    end
+    # Find the appropriate row for electricity pricing
     row = neb_data.detect do |data|
       (data['building_type'] == building_type) && (data['province'] == province) && (data['fuel_type'] == 'Electricity')
     end
@@ -269,16 +283,15 @@ class BTAPData
     year_index = 1.0
     npv_elec_ghg = 0.0
     npv_elec_peak = 0.0
+    emm_factor_conv = 3600.0 #GWH to GJ
     if eui_elec > 0.0
       for year in npv_start_year.to_int..npv_end_year.to_int
-        # puts "year, #{year}, #{row[year.to_s]}, year_index, #{year_index}"
         npv_elec += (eui_elec * row[year.to_s]) / (1+npv_discount_rate)**year_index
-        npv_elec_ghg += (ghg_elec * get_national_ghg_cost(year: year.to_i)) / (1+npv_discount_rate_carbon)**year_index
-        npv_elec_peak += (annual_peak_kW_per_m_sq * row_peak[year.to_s]) / (1+npv_discount_rate)**year_index unless row_peak.nil? || annual_peak_kW_per_m_sq.nil?
+        npv_elec_ghg += (((eui_elec * row_emm_factor[year.to_i.to_s].to_f) / emm_factor_conv) * get_national_ghg_cost(year: year.to_i)) / (1+npv_discount_rate_carbon)**year_index
+        npv_elec_peak += (annual_peak_kW_per_m_sq * row_peak[year.to_i.to_s]).to_f / (1+npv_discount_rate)**year_index unless row_peak.nil? || annual_peak_kW_per_m_sq.nil?
         year_index += 1.0
       end
     end
-    # puts "npv_elec is #{npv_elec}"
 
     # Calculate npv of natural gas
     eui_ngas= @btap_data['energy_eui_natural_gas_gj_per_m_sq']
@@ -740,7 +753,16 @@ class BTAPData
       economics_data["cost_utility_neb_#{neb_fuel.downcase}_cost_per_m_sq"] = fuel_consumption_gj * neb_fuel_cost.to_f / @conditioned_floor_area_m_sq
       economics_data['cost_utility_neb_total_cost_per_m_sq'] += economics_data["cost_utility_neb_#{neb_fuel.downcase}_cost_per_m_sq"]
       # Determine cost in GHG kg of CO2
-      economics_data["cost_utility_ghg_#{neb_fuel.downcase}_kg_per_m_sq"] = fuel_consumption_gj * get_utility_ghg_kg_per_gj(province: model.getWeatherFile.stateProvinceRegion, fuel_type: ep_fuel[:eplus_fuel_name]) / @conditioned_floor_area_m_sq
+      if neb_fuel.downcase == 'electricity' # Use ECCC data for utility emission factors rather than NEB data.
+        # Create a hash of the ECCC's electricity emissions factors for the current and future years for GHG "without Biomass and RNG CO2 emissions".
+        eccc_ghg_elec_emission_factors = CSV.parse(File.read(@eccc_ghg_electricity), headers: true, converters: :numeric).map(&:to_h)
+        # Find the appropriate row for electricity GHG emmission factors
+        row_emm_factor = eccc_ghg_elec_emission_factors.detect { |data| (data['province'] == province) }
+        emm_factor_conv = 3600.0
+        economics_data["cost_utility_ghg_#{neb_fuel.downcase}_kg_per_m_sq"] = ((fuel_consumption_gj * row_emm_factor[utility_pricing_year.to_i.to_s].to_f) / emm_factor_conv) * 1000 / @conditioned_floor_area_m_sq
+      else
+        economics_data["cost_utility_ghg_#{neb_fuel.downcase}_kg_per_m_sq"] = fuel_consumption_gj * get_utility_ghg_kg_per_gj(province: model.getWeatherFile.stateProvinceRegion, fuel_type: ep_fuel[:eplus_fuel_name]) / @conditioned_floor_area_m_sq
+      end
       economics_data['cost_utility_ghg_total_kg_per_m_sq'] += economics_data["cost_utility_ghg_#{neb_fuel.downcase}_kg_per_m_sq"]
     end
     # Commenting out block charge rates for now....
@@ -1582,6 +1604,142 @@ class BTAPData
     cooling_peak_w = [cooling_peak_w_electricity.to_f, cooling_peak_w_gas.to_f].max
     data['cooling_peak_w_per_m_sq'] = cooling_peak_w / @conditioned_floor_area_m_sq
 
+    # Peak water system load  # @todo IMPORTANT NOTE: Peak water system load must be updated if a combination of fuel types is used in a building model.
+    command = "SELECT Value
+               FROM TabularDataWithStrings
+               WHERE ReportName='EnergyMeters'
+               AND ReportForString='Entire Facility'
+               AND TableName='Annual and Peak Values - Electricity'
+               AND RowName='WaterSystems:Electricity'
+               AND ColumnName='Electricity Maximum Value'
+               AND Units='W'"
+    water_heating_peak_w_electricity = @sqlite_file.get.execAndReturnFirstDouble(command)
+    command = "SELECT Value
+              FROM TabularDataWithStrings
+              WHERE ReportName='EnergyMeters'
+              AND ReportForString='Entire Facility'
+              AND TableName='Annual and Peak Values - Natural Gas'
+              AND RowName='WaterSystems:NaturalGas'
+              AND ColumnName='Natural Gas Maximum Value'
+              AND Units='W'"
+    water_heating_peak_w_gas = @sqlite_file.get.execAndReturnFirstDouble(command)
+    water_heating_peak_w = [water_heating_peak_w_electricity.to_f, water_heating_peak_w_gas.to_f].max
+    data['energy_peak_water_systems_w_per_m_sq'] = water_heating_peak_w / @conditioned_floor_area_m_sq
+
+    #===================================================================================================================
+    # Winter peak values
+    ### Get number of timesteps from the model
+    number_of_timesteps_per_hour = 1
+    timestep_second = 3600 / number_of_timesteps_per_hour
+    #===================================================================================================================
+    ### Calculate number of timesteps of the whole year, and create dates/times for the whole year
+    number_of_timesteps_of_year = 365 * 24 * number_of_timesteps_per_hour
+    timesteps_of_year = []
+    d = Time.new(2006, 1, 1, 0)
+    (0...number_of_timesteps_of_year).each do |increment|
+      timesteps_of_year << (d + (60 * 60 / number_of_timesteps_per_hour) * increment).strftime('%Y-%m-%d %H:%M')
+    end
+    # Find month of each timestep
+    keys_month = timesteps_of_year.map{ |timestep| DateTime.parse(timestep).to_date.month }
+    # Peak electricity in winter
+    query = "
+    SELECT ReportDataDictionaryIndex
+    FROM ReportDataDictionary
+    WHERE Name=='Electricity:Facility' AND ReportingFrequency == 'Hourly' AND Units == 'J'
+    "
+    index_timestep_electricity = @model.sqlFile.get.execAndReturnVectorOfInt(query).get[0]
+    number_of_timesteps_per_hour = 1
+    if index_timestep_electricity.nil?
+      puts "No hourly or timestep data available.  Peak data not available"
+      data['energy_peak_electric_w_per_m_sq_winter'] = 0
+    else
+      ### Get timestep output for Electricity:Facility output.
+      query = "
+      Select VariableValue
+      FROM ReportVariableData
+      WHERE
+      ReportVariableDataDictionaryIndex = #{index_timestep_electricity}
+      "
+      timestep_values = @model.sqlFile.get.execAndReturnVectorOfDouble(query).get
+          # Create a hash where its keys are months and values for eahc month are electricity energy use (J) of all timesteps for that month
+      hash_months_electricity = keys_month.zip(timestep_values).group_by(&:first).map{|key, value| [key, value.map(&:last)]}.to_h
+      # Find maximum of electricity energy use (J) among all timesteps for each month of the year
+      monthly_peak_J = hash_months_electricity.map { |key, value| [key, value.max()] }.to_h
+      ### Convert J to W of maximum electricity energy use (J)
+      monthly_peak_W = monthly_peak_J.map { |key, value| [key, value / timestep_second] }.to_h
+      # Get winter peak
+      winter_peak_W = [monthly_peak_W[12].to_f, monthly_peak_W[1].to_f, monthly_peak_W[2].to_f].max
+      ### Normalize by building floor area
+      annual_peak_W_per_m_sq = winter_peak_W / @btap_data['bldg_conditioned_floor_area_m_sq'].to_f
+      data['energy_peak_electric_w_per_m_sq_winter'] = annual_peak_W_per_m_sq
+    end
+
+    # Peak electricity heating in winter
+    query = "
+    SELECT ReportDataDictionaryIndex
+    FROM ReportDataDictionary
+    WHERE Name=='Heating:Electricity' AND ReportingFrequency == 'Hourly' AND Units == 'J'
+    "
+    index_timestep_electricity = @model.sqlFile.get.execAndReturnVectorOfInt(query).get[0]
+    number_of_timesteps_per_hour = 1
+    if index_timestep_electricity.nil?
+      puts "No hourly or timestep data available.  Peak data not available"
+      data['energy_peak_electric_heating_w_per_m_sq_winter'] = 0
+    else
+      ### Get timestep output for Electricity:Facility output.
+      query = "
+      Select VariableValue
+      FROM ReportVariableData
+      WHERE
+      ReportVariableDataDictionaryIndex = #{index_timestep_electricity}
+      "
+      timestep_values = @model.sqlFile.get.execAndReturnVectorOfDouble(query).get
+      # Create a hash where its keys are months and values for eahc month are electricity energy use (J) of all timesteps for that month
+      hash_months_electricity = keys_month.zip(timestep_values).group_by(&:first).map{|key, value| [key, value.map(&:last)]}.to_h
+      # Find maximum of electricity energy use (J) among all timesteps for each month of the year
+      monthly_peak_J = hash_months_electricity.map { |key, value| [key, value.max()] }.to_h
+      ### Convert J to W of maximum electricity energy use (J)
+      monthly_peak_W = monthly_peak_J.map { |key, value| [key, value / timestep_second] }.to_h
+      # Get winter peak
+      winter_peak_W = [monthly_peak_W[12].to_f, monthly_peak_W[1].to_f, monthly_peak_W[2].to_f].max
+      ### Normalize by building floor area
+      annual_peak_W_per_m_sq = winter_peak_W / @btap_data['bldg_conditioned_floor_area_m_sq'].to_f
+      data['energy_peak_electric_heating_w_per_m_sq_winter'] = annual_peak_W_per_m_sq
+    end
+
+    # Peak water system electricity heating in winter
+    query = "
+    SELECT ReportDataDictionaryIndex
+    FROM ReportDataDictionary
+    WHERE Name=='WaterSystems:Electricity' AND ReportingFrequency == 'Hourly' AND Units == 'J'
+    "
+    index_timestep_electricity = @model.sqlFile.get.execAndReturnVectorOfInt(query).get[0]
+    number_of_timesteps_per_hour = 1
+    if index_timestep_electricity.nil?
+      puts "No hourly or timestep data available.  Peak data not available"
+      data['energy_peak_electric_water_systems_w_per_m_sq_winter'] = 0
+    else
+      ### Get timestep output for Electricity:Facility output.
+      query = "
+      Select VariableValue
+      FROM ReportVariableData
+      WHERE
+      ReportVariableDataDictionaryIndex = #{index_timestep_electricity}
+      "
+      timestep_values = @model.sqlFile.get.execAndReturnVectorOfDouble(query).get
+          # Create a hash where its keys are months and values for eahc month are electricity energy use (J) of all timesteps for that month
+      hash_months_electricity = keys_month.zip(timestep_values).group_by(&:first).map{|key, value| [key, value.map(&:last)]}.to_h
+      # Find maximum of electricity energy use (J) among all timesteps for each month of the year
+      monthly_peak_J = hash_months_electricity.map { |key, value| [key, value.max()] }.to_h
+      ### Convert J to W of maximum electricity energy use (J)
+      monthly_peak_W = monthly_peak_J.map { |key, value| [key, value / timestep_second] }.to_h
+      # Get winter peak
+      winter_peak_W = [monthly_peak_W[12].to_f, monthly_peak_W[1].to_f, monthly_peak_W[2].to_f].max
+      ### Normalize by building floor area
+      annual_peak_W_per_m_sq = winter_peak_W / @btap_data['bldg_conditioned_floor_area_m_sq'].to_f
+      data['energy_peak_electric_water_systems_w_per_m_sq_winter'] = annual_peak_W_per_m_sq
+    end
+
     return data
   end
 
@@ -2081,20 +2239,36 @@ class BTAPData
   def get_utility_ghg_kg_per_gj(province:, fuel_type:)
     ghg_data = [
       # Obtained from Portfolio Manager https://portfoliomanager.energystar.gov/pdf/reference/Emissions.pdf 10/10/2020
-      { "province": 'AB', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 53.24, "CO2eq Emissions (g/m3)": 1939.0 },
-      { "province": 'BC', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 53.19, "CO2eq Emissions (g/m3)": 1937.0 },
-      { "province": 'MB', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.09, "CO2eq Emissions (g/m3)": 1897.0 },
-      { "province": 'NB', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
-      { "province": 'NL', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
-      { "province": 'NT', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
-      { "province": 'NS', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
-      { "province": 'NU', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
-      { "province": 'ON', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.14, "CO2eq Emissions (g/m3)": 1912.0 },
-      { "province": 'PE', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
-      { "province": 'QC', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.12, "CO2eq Emissions (g/m3)": 1898.0 },
-      { "province": 'SK', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 50.53, "CO2eq Emissions (g/m3)": 1840.0 },
-      { "province": 'YT', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
+      #{ "province": 'AB', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 53.24, "CO2eq Emissions (g/m3)": 1939.0 },
+      #{ "province": 'BC', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 53.19, "CO2eq Emissions (g/m3)": 1937.0 },
+      #{ "province": 'MB', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.09, "CO2eq Emissions (g/m3)": 1897.0 },
+      #{ "province": 'NB', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
+      #{ "province": 'NL', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
+      #{ "province": 'NT', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
+      #{ "province": 'NS', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
+      #{ "province": 'NU', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
+      #{ "province": 'ON', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.14, "CO2eq Emissions (g/m3)": 1912.0 },
+      #{ "province": 'PE', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
+      #{ "province": 'QC', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.12, "CO2eq Emissions (g/m3)": 1898.0 },
+      #{ "province": 'SK', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 50.53, "CO2eq Emissions (g/m3)": 1840.0 },
+      #{ "province": 'YT', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.50, "CO2eq Emissions (g/m3)": 1912.0 },
+
+      #REF (Retrieved March 11, 2025) for Year 2022 of "Table A6.1–1 CO2 Emission Factors for Marketable Natural Gas" of National Inventory Report: https://data-donnees.az.ec.gc.ca/api/file?path=%2Fsubstances%2Fmonitor%2Fcanada-s-official-greenhouse-gas-inventory%2FD-Emission-Factors%2FEN_Annex6_Emission_Factors.pdf
+      { "province": 'AB', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 53.87, "CO2eq Emissions (g/m3)": 1962.0 },
+      { "province": 'BC', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 53.98, "CO2eq Emissions (g/m3)": 1966.0 },
+      { "province": 'MB', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.58, "CO2eq Emissions (g/m3)": 1915.0 },
+      { "province": 'NB', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.69, "CO2eq Emissions (g/m3)": 1919.0 },
+      { "province": 'NL', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.69, "CO2eq Emissions (g/m3)": 1919.0 },
+      { "province": 'NT', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 53.98, "CO2eq Emissions (g/m3)": 1966.0 },
+      { "province": 'NS', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.69, "CO2eq Emissions (g/m3)": 1919.0 },
+      { "province": 'NU', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 53.98, "CO2eq Emissions (g/m3)": 1966.0 },
+      { "province": 'ON', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.75, "CO2eq Emissions (g/m3)": 1921.0 },
+      { "province": 'PE', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.69, "CO2eq Emissions (g/m3)": 1919.0 },
+      { "province": 'QC', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.88, "CO2eq Emissions (g/m3)": 1926.0 },
+      { "province": 'SK', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 52.72, "CO2eq Emissions (g/m3)": 1920.0 },
+      { "province": 'YT', "fuel_type": 'NaturalGas', "CO2eq Emissions (kg/MBtu)": 56.98, "CO2eq Emissions (g/m3)": 1966.0 },
 
+      # Obtained from Portfolio Manager https://portfoliomanager.energystar.gov/pdf/reference/Emissions.pdf 10/10/2020
       { "province": 'AB', "fuel_type": 'FuelOilNo2', "CO2eq Emissions (kg/MBtu)": 75.13, "CO2eq Emissions (g/m3)": 2763.0 },
       { "province": 'BC', "fuel_type": 'FuelOilNo2', "CO2eq Emissions (kg/MBtu)": 75.13, "CO2eq Emissions (g/m3)": 2763.0 },
       { "province": 'MB', "fuel_type": 'FuelOilNo2', "CO2eq Emissions (kg/MBtu)": 75.13, "CO2eq Emissions (g/m3)": 2763.0 },

data/lib/openstudio-standards/standards/necb/common/btap_datapoint.rb

@@ -2,7 +2,7 @@ require 'openstudio'
 require 'securerandom'
 require 'optparse'
 require 'yaml'
-#require 'git-revision'
+# require 'git-revision'
 # resource_folder = File.join(__dir__, '..', '..', 'measures/btap_results/resources')
 # OpenSSL::SSL::VERIFY_PEER = OpenSSL::SSL::VERIFY_NONE
 
@@ -77,7 +77,6 @@ class BTAPDatapoint
     @dp_output_folder = File.join(output_folder, @options[:datapoint_id])
 
     # Show versions in yml file.
-    @options[:btap_costing_git_revision] = Git::Revision.commit_short
     @options[:os_git_revision] = OpenstudioStandards.git_revision
 
     # Get users' inputs for the parameters needed for the calculation af net present value
@@ -125,7 +124,13 @@ class BTAPDatapoint
         @standard.model_apply_standard(model: model,
                                        epw_file: @options[:epw_file],
                                        custom_weather_folder: weather_folder,
+                                       btap_weather: @options[:btap_weather],
                                        sizing_run_dir: File.join(@dp_temp_folder, 'sizing_folder'),
+                                       hvac_system_primary: @options[:hvac_system_primary],
+                                       hvac_system_dwelling_units: @options[:hvac_system_dwelling_units],
+                                       hvac_system_washrooms: @options[:hvac_system_washrooms],
+                                       hvac_system_corridor: @options[:hvac_system_corridor],
+                                       hvac_system_storage: @options[:hvac_system_storage],
                                        primary_heating_fuel: @options[:primary_heating_fuel],
                                        necb_reference_hp: @options[:necb_reference_hp],
                                        necb_reference_hp_supp_fuel: @options[:necb_reference_hp_supp_fuel],
@@ -185,6 +190,7 @@ class BTAPDatapoint
                                        boiler_fuel: @options[:boiler_fuel],
                                        boiler_cap_ratio: @options[:boiler_cap_ratio],
                                        swh_fuel: @options[:swh_fuel],
+                                       airloop_fancoils_heating: @options[:airloop_fancoils_heating],
                                        oerd_utility_pricing: @oerd_utility_pricing
                                        )
       end
@@ -201,30 +207,7 @@ class BTAPDatapoint
         #Check if sql file was create.. if not we can't get much output.
         raise('no sql file found,simulation probably could not start due to error...see error logs.') if model.sqlFile.empty?
         # Create qaqc file and save it.
-        @qaqc = @standard.init_qaqc(model)
-        command = "SELECT Value
-                  FROM TabularDataWithStrings
-                  WHERE ReportName='LEEDsummary'
-                  AND ReportForString='Entire Facility'
-                  AND TableName='Sec1.1A-General Information'
-                  AND RowName = 'Principal Heating Source'
-                  AND ColumnName='Data'"
-        value = model.sqlFile.get.execAndReturnFirstString(command)
-        # make sure all the data are available
-        @qaqc[:building][:principal_heating_source] = 'unknown'
-        unless value.empty?
-          @qaqc[:building][:principal_heating_source] = value.get
-        end
-
-        if @qaqc[:building][:principal_heating_source] == 'Additional Fuel'
-          model.getPlantLoops.sort.each do |iplantloop|
-            boilers = iplantloop.components.select { |icomponent| icomponent.to_BoilerHotWater.is_initialized }
-            @qaqc[:building][:principal_heating_source] = 'FuelOilNo2' unless boilers.select { |boiler| boiler.to_BoilerHotWater.get.fuelType.to_s == 'FuelOilNo2' }.empty?
-          end
-        end
-
-        @qaqc[:aws_datapoint_id] = @options[:datapoint_id]
-        @qaqc[:aws_analysis_id] = @options[:analysis_id]
+        @qaqc = BTAPDatapoint.build_qaqc(model, @standard, @options[:datapoint_id], @options[:analysis_id])
 
         # Load the sql file into model
         sql_path = OpenStudio::Path.new(File.join(@run_dir, 'run/eplusout.sql'))
@@ -240,15 +223,26 @@ class BTAPDatapoint
         # Attach the sql file from the run to the sizing model
         model.setSqlFile(sql)
 
+        if @options[:enable_costing] or @options[:enable_carbon]
+          @cp = CommonPaths.instance
+          post_analysis = BTAPDatapointAnalysis.new(
+            model: model,
+            output_folder: @dp_temp_folder,
+            template: @options[:template],
+            standard: @standard,
+            qaqc: @qaqc)
+        end
+
         @cost_result = nil
-        if defined?(BTAPCosting)
-          # Perform costing
-          costing = BTAPCosting.new
-          costing.load_database
-          @cost_result, @btap_items = costing.cost_audit_all(model: model, prototype_creator: @standard, template_type: @options[:template])
-          @qaqc[:costing_information] = @cost_result
-          File.open(File.join(@dp_temp_folder, 'cost_results.json'), 'w') { |f| f.write(JSON.pretty_generate(@cost_result, allow_nan: true)) }
-          puts "Wrote File cost_results.json in #{Dir.pwd} "
+        if @options[:enable_costing]
+          @cost_result = post_analysis.run_costing(
+            costs_csv: @cp.costs_path,
+            factors_csv: @cp.costs_local_factors_path)
+        end
+
+        @carbon_result = nil
+        if @options[:enable_carbon]
+          @carbon_result = post_analysis.run_carbon
         end
 
         @qaqc[:options] = @options # This is options sent on the command line
@@ -256,6 +250,7 @@ class BTAPDatapoint
         @btap_data = BTAPData.new(model: model,
                                   runner: nil,
                                   cost_result: @cost_result,
+                                  carbon_result: @carbon_result,
                                   qaqc: @qaqc,
                                   npv_start_year: @npv_start_year,
                                   npv_end_year: @npv_end_year,
@@ -300,25 +295,62 @@ class BTAPDatapoint
       # clean temp/cache folder up.
       FileUtils.rm_rf(input_folder_cache)
       FileUtils.rm_rf(@dp_temp_folder)
-      if @failed == true
+      # Do not fail container if running on AWS, handle error on AWS instead to avoid isuses with large analyses.
+      unless @failed == false || (@dp_output_folder.start_with?('s3://'))
+      #if @failed == true
         raise(@bang)
       end
     end
   end
 
+  class << self
+
+    # Initializes the qaqc data structure.
+    # Scoped inside of the class so that it can be used in the intialization of
+    # this class as well as in BTAPAnalysis.
+    def build_qaqc(model, standard, datapoint_id, analysis_id)
+      qaqc = standard.init_qaqc(model)
+      command = "SELECT Value
+                FROM TabularDataWithStrings
+                WHERE ReportName='LEEDsummary'
+                AND ReportForString='Entire Facility'
+                AND TableName='Sec1.1A-General Information'
+                AND RowName = 'Principal Heating Source'
+                AND ColumnName='Data'"
+      value = model.sqlFile.get.execAndReturnFirstString(command)
+      # make sure all the data are available
+      qaqc[:building][:principal_heating_source] = 'unknown'
+      unless value.empty?
+        qaqc[:building][:principal_heating_source] = value.get
+      end
+
+      if qaqc[:building][:principal_heating_source] == 'Additional Fuel'
+        model.getPlantLoops.sort.each do |iplantloop|
+          boilers = iplantloop.components.select { |icomponent| icomponent.to_BoilerHotWater.is_initialized }
+          qaqc[:building][:principal_heating_source] = 'FuelOilNo2' unless boilers.select { |boiler| boiler.to_BoilerHotWater.get.fuelType.to_s == 'FuelOilNo2' }.empty?
+        end
+      end
+
+      qaqc[:aws_datapoint_id] = datapoint_id
+      qaqc[:aws_analysis_id]  = analysis_id
+
+      return qaqc
+    end
+  end
+
   def s3_copy_file_to_s3(bucket_name:, source_file:, target_file:, n: 0)
-    #require 'aws-sdk-core'
-    #require 'aws-sdk-s3'
+    # require 'aws-sdk-core'
+    # require 'aws-sdk-s3'
     Aws.use_bundled_cert!
-    s3_resource = Aws::S3::Resource.new(region: 'ca-central-1')
+    s3_client = Aws::S3::Client.new(region: 'ca-central-1')
+    # Using transfer manager class instead of depricated method
+    transfer_manager = Aws::S3::TransferManager.new(client: s3_client)
 
     puts("Copying File to S3. source_file:#{source_file} bucket:#{bucket_name} target_folder:#{target_file}")
     response = nil
     begin
-      obj = s3_resource.bucket(bucket_name).object(target_file)
-
-      # passing the TempFile object's path is massively faster than passing the TempFile object itself
-      result = obj.upload_file(source_file)
+      # Using upload_file from transfer manager class instead of depricated method from older class structure
+      result = transfer_manager.upload_file(source_file, bucket: bucket_name, key: target_file)
 
       if result == true
         puts "Object '#{source_file}' uploaded to bucket '#{bucket_name}'."
@@ -618,4 +650,4 @@ class BTAPDatapoint
     #===================================================================================================================
   end #output_timestep_data
   #=====================================================================================================================
-end
+end

data/lib/openstudio-standards/standards/necb/common/eccc_electric_grid_intensity_20250311.csv

@@ -0,0 +1,14 @@
+province,units,2005,2006,2007,2008,2009,2010,2011,2012,2013,2014,2015,2016,2017,2018,2019,2020,2021,2022,2023,2024,2025,2026,2027,2028,2029,2030,2031,2032,2033,2034,2035,2036,2037,2038,2039,2040,2041,2042,2043,2044,2045,2046,2047,2048,2049,2050
+Alberta,tCO2e/GWh,958.7,991.3,925.2,916.8,822.5,843.4,794.1,754.5,785.7,789.4,744.9,768.4,765.3,604.3,568.6,530.2,458.3,382.1,304.5,301.4,265.2,226.8,210.2,172.1,165,170.9,166.2,153.9,145,138.4,118.8,106.5,100.3,82.3,77.2,71.4,69.9,70.3,70.9,71.5,71.6,70.6,70.3,68.6,43.9,19.2
+British Columbia,tCO2e/GWh,20.7,20,17,25.3,20.9,21.6,9.4,10.1,10.6,6.9,4.6,2.9,2,5.4,11.6,5.3,8.2,5.7,5.1,5.4,5.3,4.2,3.9,3.7,1.5,1.4,1.4,1.3,1.2,1.1,1,0.6,0.6,0.5,0.5,0.5,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1
+Manitoba,tCO2e/GWh,19.4,24,26.7,21,9.6,4.5,5.8,5.4,5.4,5.6,5.6,3.2,3.1,1.8,1.8,1.9,2.6,2.1,1.8,1.1,0.5,0.7,0.7,0.7,0.7,0.7,0.7,0.7,0.7,0.7,0.7,0.7,0.7,0.7,0.6,0.6,0.6,0.6,0.6,0.6,0.6,0.6,0.6,0.6,0.6,0.6
+New Brunswick,tCO2e/GWh,613,595.4,539,546.9,549.6,379.4,333.2,281.4,236.6,264.7,259.4,333.3,261.1,281.6,263,182.5,264.5,291.2,298.4,295.7,319.1,306,278.1,266.1,252,37,31.8,24.8,29.1,29,46.9,52.2,65,64.4,27.8,27.7,2.9,2.5,2.5,2.4,1.9,14.3,12.5,10.6,1.2,0.9
+Newfoundland and Labrador,tCO2e/GWh,85.8,64.8,106.4,90.3,89.7,76.6,81.7,79.3,86,119.8,148.3,151.8,145.1,105.9,104.7,90.3,74.2,81.4,7.6,6.3,6.3,6.2,6,5.8,6.4,5.6,5.5,5.3,5.2,5.1,5,4.8,5,5.3,5.3,5.2,5.1,4.9,4.5,4.2,4.2,4.1,4.1,4,4,3.9
+Northwest Territories,tCO2e/GWh,303.6,275.4,304.7,263.2,205.1,236.3,226.1,115.4,190.2,241.7,352.2,196.2,189.8,202.7,186.1,191.8,192.6,215.5,234.9,175,1,0,0,5.2,4.6,4.1,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
+Nova Scotia,tCO2e/GWh,927.2,974.9,979,803.5,817.8,741.4,743,726.4,713.3,689.9,643.5,656.2,651.2,684.4,645.1,619.3,582.9,542.4,658.1,605.3,529.8,477.9,461.9,448.1,303.9,88.2,89.7,112,105,98.4,70.6,66.7,63.2,60.5,59,57.2,55.4,46.7,45.9,63.3,60.2,57.1,53.7,50.3,47.1,38.6
+Nunavut,tCO2e/GWh,606,507.6,510.1,539,566.8,566.8,482.4,484.2,453.9,567.5,481.8,607.5,574.8,675.3,673.3,616.2,639.8,628.2,788.4,724.3,494.4,456.6,442.9,444.2,456.9,454.3,466.7,410.7,425,449,401,411.4,423.9,378.1,389.3,403.4,361.2,373.9,339.2,304.3,318.7,288.9,257.7,267.9,237.4,208.7
+Ontario,tCO2e/GWh,224.3,206,228.7,189,112.2,143.4,91.2,90.9,67.3,40,42.1,36,16.4,25.5,24.8,24.2,26.3,30.2,23.9,55.5,65.3,76.1,56.9,50.4,50,65,55.9,65.8,68.4,76.1,30.5,23.3,19.1,13.6,10.7,13.1,12.9,11.9,11.7,11.5,14,13.8,13.6,13.4,12.5,3.4
+Prince Edward Island,tCO2e/GWh,5,8.4,4.8,4,5.3,1.6,1.1,9,3.2,3.4,11.3,3.5,6.5,2.1,0.8,0.2,1.4,0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.8,0.7,0.7,0.7,0.7,0.7,0.6,0.6,0.6,0.6,0.6,0.6,0.6,0.6,0.5,0.5,0.5,0.5,0.5
+Quebec,tCO2e/GWh,3.2,5.1,8.7,2.5,3.5,2.4,2.7,2.5,1.4,1.5,1.3,1.5,1.4,1.4,1.5,1.8,1.5,1.3,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0,0
+Saskatchewan,tCO2e/GWh,860.6,813.3,803.3,769.7,872,846.4,797.5,807.8,685.2,700.2,736.8,721.4,703.2,720.6,698.5,615.7,692.9,648.2,542.8,560.6,522.9,515.3,499.4,341.2,330.8,115.3,117.1,108.6,104.3,100.3,48.7,40.4,37,34.9,52.4,58.6,62.6,57.6,53.2,40.8,23.8,23.2,22.9,22.3,21.8,12.6
+Yukon,tCO2e/GWh,70,53.4,53.2,51.2,44.6,48.6,68.7,42.8,41.4,41.4,44.6,47.1,53.5,76.2,113.1,111.3,79.3,73.3,120,107.6,7,6.9,6.8,7,8.3,10.2,41.9,48,52.1,52.3,58.9,57.4,55.2,53.2,52.3,46.3,29.8,28.7,27.3,25.2,24.5,24.6,28.5,27.2,27.3,27.6

data/lib/openstudio-standards/standards/necb/common/nir_gas_grid_intensity_20250311.csv

@@ -0,0 +1,14 @@
+province,units,2022
+Alberta,g/m3,1962
+British Columbia,g/m3,1966
+Manitoba,g/m3,1915
+New Brunswick,g/m3,1919
+Newfoundland and Labrador,g/m3,1919
+Northwest Territories,g/m3,1966
+Nova Scotia,g/m3,1919
+Nunavut,g/m3,1966
+Ontario,g/m3,1921
+Prince Edward Island,g/m3,1919
+Quebec,g/m3,1926
+Saskatchewan,g/m3,1920
+Yukon,g/m3,1966