openstudio-standards 0.8.3 → 0.8.4

data/lib/openstudio-standards/standards/necb/README.md

@@ -0,0 +1,343 @@
+# NECB (National Energy Code for Buildings) Standards
+
+## Overview
+
+The NECB (National Energy Code for Buildings) standards module provides comprehensive implementation of Canadian building energy codes. This module represents one of the most complete implementations of a national building energy standard, covering multiple versions of the Canadian National Energy Code for Buildings, vintage building standards, and energy conservation measures.
+
+## Architecture
+
+The NECB module is organized into four main categories:
+
+1. **NECB Standard Versions** (2011, 2015, 2017, 2020)
+2. **BTAP Vintage Standards** (Pre-1980, 1980-2010)  
+3. **Energy Conservation Measures (ECMS)**
+4. **Shared Resources** (Common data, documentation, QAQC)
+
+## NECB Standard Versions
+
+### NECB2011 - Foundation Implementation
+
+#### Core Module: `necb_2011.rb`
+**Most comprehensive NECB implementation with full building system coverage**
+
+**Key Features:**
+- Complete climate zone implementation for Canada (8 climate zones)
+- All 8 HVAC system types with both single-speed and multi-speed variants
+- Comprehensive building envelope requirements
+- LED lighting compliance and occupancy sensor controls
+- Service water heating system requirements
+- BEPS (Building Energy Performance Standards) compliance pathways
+- Thermal bridging calculations (TBD integration)
+- Quality assurance and quality control (QAQC) framework
+
+**Major Subsystems:**
+
+##### HVAC Systems (`hvac_system_*.rb`)
+- **System 1**: Single Zone systems (`hvac_system_1_single_speed.rb`, `hvac_system_1_multi_speed.rb`)
+- **Systems 2 & 5**: Multi-zone VAV systems (`hvac_system_2_and_5.rb`)
+- **Systems 3 & 8**: Single zone heat pumps (`hvac_system_3_and_8_single_speed.rb`, `hvac_system_3_and_8_multi_speed.rb`)  
+- **System 4**: Multi-zone VAV with reheat (`hvac_system_4.rb`)
+- **System 6**: Multi-zone VAV with PFP boxes (`hvac_system_6.rb`)
+- **General**: `hvac_systems.rb`, `hvac_namer.rb`
+
+##### Building Systems
+- **Envelope**: `building_envelope.rb` - Comprehensive building envelope requirements
+- **Lighting**: `lighting.rb` - Interior lighting power density and controls
+- **Service Water Heating**: `service_water_heating.rb` - Hot water system efficiency requirements
+- **Electrical**: `electrical_power_systems_and_motors.rb` - Motor efficiency requirements
+
+##### Specialized Features  
+- **Auto-zoning**: `autozone.rb` - Automated thermal zone creation
+- **System Fuels**: `system_fuels.rb` - Fuel type selection and mapping
+- **BEPS Compliance**: `beps_compliance_path.rb` - Building Energy Performance Standards
+
+##### Data Resources (`data/`)
+**Extensive JSON database covering:**
+- **Climate & Geography**: `climate_zone_sets.json`, `province_map.json`
+- **Building Envelope**: `materials.json`, `surface_thermal_transmittance.json`
+- **HVAC Equipment**: `boilers.json`, `chillers.json`, `furnaces.json`, `heat_pumps.json`, `unitary_acs.json`
+- **System Selection**: `systems.json`, `necb_hvac_system_selection_type.json`
+- **Energy Efficiency**: `motors.json`, `heat_rejection.json`, `erv.json`
+- **Building Types**: `space_types.json`, `prototype_inputs.json`
+- **Lighting**: `led_lighting_data.json`
+- **Economics**: `regional_fuel_use.json`
+- **Reference Data**: `necb_2015_table_c1.json`, `constants.json`, `formulas.json`
+
+##### Quality Assurance (`qaqc/`)
+- **Main Module**: `necb_qaqc.rb` - Comprehensive model validation
+- **Resources**: `qaqc_resources/` - Reference data for validation
+- **Data**: `qaqc_data/` - Expected results for testing
+
+### NECB2015 - Enhanced Version  
+
+#### Core Module: `necb_2015.rb`
+**Updated NECB implementation with enhanced lighting and HVAC systems**
+
+**Key Improvements over 2011:**
+- Enhanced lighting power density requirements
+- Updated HVAC system efficiency requirements
+- Refined control strategies
+
+**Modules:**
+- **Lighting**: `lighting.rb` - Updated lighting requirements
+- **HVAC**: `hvac_systems.rb` - Enhanced system implementations
+- **QAQC**: `qaqc/necb_2015_qaqc.rb` - Version-specific validation
+
+### NECB2017 - System Updates
+
+#### Core Module: `necb_2017.rb`  
+**Newer version focusing on HVAC system improvements**
+
+**Key Features:**
+- Updated HVAC system requirements
+- Enhanced system control strategies
+- Refined equipment efficiency requirements
+
+**Modules:**
+- **HVAC**: `hvac_systems.rb` - Updated system implementations
+
+### NECB2020 - Latest Version
+
+#### Core Module: `necb_2020.rb`
+**Latest NECB implementation with building envelope and service water heating improvements**
+
+**Key Features:**
+- Enhanced building envelope performance requirements
+- Updated service water heating efficiency standards
+- Latest building code compliance requirements
+
+**Modules:**
+- **Envelope**: `building_envelope.rb` - Updated envelope requirements
+- **Service Water Heating**: `service_water_heating.rb` - Enhanced efficiency requirements
+
+## BTAP Vintage Standards
+
+### BTAPPRE1980 - Pre-1980 Building Standards
+
+#### Core Module: `btap_pre1980.rb`
+**Building performance standards for buildings constructed before 1980**
+
+**Purpose:** Historical building analysis and retrofit potential assessment
+
+**Key Features:**
+- Legacy building envelope performance levels
+- Historical HVAC system types and efficiencies  
+- Vintage building infiltration rates
+- Baseline performance for retrofit analysis
+
+**Modules:**
+- **Envelope**: `building_envelope.rb` - Pre-1980 construction standards
+- **HVAC Systems**: Multiple system files for historical equipment types
+- **Data**: `data/` - Historical performance data
+
+### BTAP1980TO2010 - Transitional Period Standards
+
+#### Core Module: `btap_1980to2010.rb`
+**Building standards for the 1980-2010 construction period**  
+
+**Purpose:** Bridge between historical and modern building standards
+
+**Key Features:**
+- Transitional building envelope requirements
+- Equipment efficiency improvements over pre-1980 era
+- Building code evolution tracking
+
+## Energy Conservation Measures (ECMS)
+
+### ECMS Framework (`ECMS/`)
+
+#### Core Module: `ecms.rb`
+**Energy Conservation Measures implementation framework**
+
+**Purpose:** Systematic evaluation and implementation of building energy efficiency measures
+
+**Key ECM Categories:**
+
+##### HVAC ECMs (`hvac_systems.rb`)
+- Advanced HVAC system configurations
+- High-efficiency equipment options
+- System optimization strategies
+
+##### Energy Recovery Ventilation (`erv.rb`)
+- Energy recovery ventilator implementation
+- Heat recovery system integration
+- Ventilation system optimization
+
+##### Natural Ventilation (`nv.rb`)  
+- Natural ventilation strategies
+- Mixed-mode system implementation
+- Climate-responsive ventilation control
+
+##### Renewable Energy (`pv_ground.rb`)
+- Ground-mounted photovoltaic systems
+- Solar energy integration
+- On-site renewable energy generation
+
+##### Load Optimization (`loads.rb`)
+- Building load reduction strategies
+- Equipment efficiency improvements
+- Operational optimization measures
+
+##### Data Resources (`data/`)
+- ECM performance data
+- Cost-benefit analysis parameters
+- Implementation guidelines
+
+## Shared Resources
+
+### Common Data and Utilities (`common/`)
+
+#### BTAP Integration
+- **`btap_datapoint.rb`**: Large-scale building analysis utilities
+- **`btap_data.rb`**: Common BTAP data handling functions
+
+#### Economic and Environmental Data
+- **`national_carbon_price_sched_feb_2025.json`**: Canadian carbon pricing schedules
+- **`utility_pricing_2025-02-20.csv`**: Canadian utility rate data
+- **`neb_end_use_prices.csv`**: End-use energy pricing
+- **`eccc_electric_grid_intensity_20250311.csv`**: Electric grid emissions factors
+- **`nir_gas_grid_intensity_20250311.csv`**: Natural gas grid emissions factors
+
+#### Reference Data
+- **`necb_reference_runs.csv`**: Reference building performance data
+- **`space_type_upgrade_map.json`**: Space type transformation mappings
+- **`system_types.yaml`**: HVAC system type definitions
+- **`construction_defaults.osm`**: Default construction assemblies
+
+#### Documentation (`docs/`)
+- **`air_system_names_method.md`**: Air system naming conventions
+
+#### Specialized Resources
+- **`phius.md`**: Passive House Institute US integration notes
+- **`bc_step_code_indicators.md`**: British Columbia Step Code compliance
+
+## Key Features and Capabilities
+
+### Climate Zone Implementation
+- **8 Canadian Climate Zones**: Complete implementation from Zone 4 (warmest) to Zone 8 (coldest)
+- **Provincial Mapping**: Automatic climate zone determination by province and city
+- **Weather Data Integration**: Seamless integration with Canadian weather files
+
+### HVAC System Types
+**Complete implementation of all 8 NECB system types:**
+
+1. **System 1**: Residential/Small Commercial (≤ 25,000 ft²)
+2. **System 2**: Multi-zone VAV with reheat
+3. **System 3**: Single-zone heat pump  
+4. **System 4**: VAV with reheat (large buildings)
+5. **System 5**: Packaged VAV with reheat
+6. **System 6**: VAV with PFP boxes and reheat
+7. **System 7**: VAV with hot water reheat
+8. **System 8**: VAV with electric reheat
+
+**Each system type includes:**
+- Single-speed and multi-speed equipment options
+- Appropriate control strategies
+- Equipment sizing and selection
+- Energy efficiency requirements
+
+### Building Performance Standards (BEPS)  
+- **Performance Path Compliance**: Alternative compliance using energy modeling
+- **Prescriptive Path Integration**: Seamless integration with prescriptive requirements
+- **Trade-off Analysis**: Building component performance trade-offs
+- **Compliance Reporting**: Automated compliance documentation
+
+### Advanced Features
+
+#### Thermal Bridging Integration
+- **TBD Integration**: Heat transfer through thermal bridges
+- **Detailed Heat Transfer**: Beyond simple assembly U-values
+- **Whole Building Performance**: Complete thermal performance assessment
+
+#### LED Lighting Compliance
+- **Power Density Requirements**: Automatic lighting power density calculations
+- **Control Integration**: Occupancy sensor and daylight dimming controls
+- **Compliance Verification**: Automated lighting compliance checking
+
+#### Quality Assurance Framework
+- **Comprehensive QAQC**: Model validation and quality checking
+- **Expected Results Testing**: Comparison against reference buildings
+- **Performance Verification**: Energy use intensity and component performance validation
+
+### Canadian-Specific Features
+
+#### Provincial Integration  
+- **Provincial Fuel Mapping**: Province-specific fuel type availability
+- **Utility Rate Integration**: Provincial electricity and gas rate structures
+- **Climate Data**: Canadian weather station integration
+
+#### Economic Analysis
+- **Regional Cost Factors**: Provincial and regional cost adjustments
+- **Carbon Pricing**: Integration with Canadian carbon pricing policies
+- **Lifecycle Cost Analysis**: Complete economic analysis framework
+
+## Usage Patterns
+
+### Typical Implementation Workflow
+
+1. **Standard Selection**: Choose appropriate NECB version (2011, 2015, 2017, 2020)
+2. **Climate Zone Determination**: Automatic based on building location
+3. **System Type Selection**: Based on building type and size
+4. **Component Specification**: Apply standard-specific requirements
+5. **Performance Verification**: QAQC validation and compliance checking
+
+### ECM Evaluation Process
+
+1. **Baseline Model**: Create NECB-compliant baseline
+2. **ECM Selection**: Choose from available energy conservation measures
+3. **Performance Modeling**: Model ECM implementation
+4. **Economic Analysis**: Cost-benefit analysis using Canadian data
+5. **Optimization**: Identify optimal ECM combinations
+
+### Integration Points
+
+- **BTAP Framework**: Complete integration with BTAP costing and analysis
+- **OpenStudio Standards**: Seamless integration with base standards framework
+- **Canadian Data**: Weather, utility, and economic data integration
+
+## Dependencies
+
+### Internal Dependencies
+- **Base Standards Framework**: Inherits from OpenStudio-Standards base classes
+- **BTAP Integration**: Extensive integration with BTAP utilities and costing
+- **Common Utilities**: Uses shared utility functions and data structures
+
+### External Dependencies  
+- **OpenStudio SDK**: Building energy modeling capabilities
+- **EnergyPlus**: Energy simulation engine
+- **TBD**: Thermal bridging analysis (optional)
+
+### Data Dependencies
+- **Canadian Weather Data**: Environment and Climate Change Canada weather files
+- **Economic Data**: Statistics Canada and provincial utility data
+- **Building Code Data**: National Research Council of Canada standards
+
+## Development and Maintenance
+
+### Version Evolution
+- **NECB2011**: Foundation implementation (most complete)
+- **NECB2015**: Lighting and HVAC updates
+- **NECB2017**: System control improvements  
+- **NECB2020**: Envelope and service water heating updates
+
+### Testing Framework
+- **Reference Buildings**: Comprehensive test suite using Canadian reference buildings
+- **Regression Testing**: Automated testing against expected results
+- **Compliance Verification**: Validation against published code requirements
+
+### Maintenance Considerations
+- **Code Updates**: Regular updates to reflect changes in Canadian building codes
+- **Data Updates**: Annual updates to utility rates, carbon pricing, and cost data
+- **Weather Data**: Updates to Canadian weather files as available
+
+## Canadian Building Code Context
+
+The NECB standards module provides the most comprehensive implementation of Canadian building energy requirements available in any building energy modeling tool. It reflects the unique aspects of Canadian building codes:
+
+- **Climate-Based Requirements**: Robust climate zone implementation
+- **System Type Flexibility**: Complete range of HVAC system options
+- **Performance-Based Compliance**: Both prescriptive and performance path options  
+- **Provincial Variations**: Accommodation of provincial building code differences
+- **Economic Integration**: Canadian-specific cost and economic data
+
+This makes the NECB module essential for any building energy analysis work in Canada.

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

@@ -20,6 +20,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 +161,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 +230,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 +250,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 +271,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 +741,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 +1592,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 +2227,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 },