urbanopt-reporting 0.4.3 → 0.6.1

data/lib/measures/.rubocop.yml

@@ -1,5 +1,7 @@
 AllCops:
   Exclude:
-    - 'spec/test_measures/**/*'
+    - gems/**/*
+require: rubocop-performance
+
 inherit_from:
-  - http://s3.amazonaws.com/openstudio-resources/styles/rubocop_v3.yml
+  - http://s3.amazonaws.com/openstudio-resources/styles/rubocop_v4.yml

data/lib/measures/default_feature_reports/measure.rb

@@ -1,5 +1,5 @@
 # *********************************************************************************
-# URBANopt™, Copyright (c) 2019-2021, Alliance for Sustainable Energy, LLC, and other
+# URBANopt™, Copyright (c) 2019-2022, Alliance for Sustainable Energy, LLC, and other
 # contributors. All rights reserved.
 
 # Redistribution and use in source and binary forms, with or without modification,
@@ -43,7 +43,7 @@ require 'csv'
 require 'benchmark'
 require 'logger'
 
-@@logger = Logger.new(STDOUT)
+@@logger = Logger.new($stdout)
 
 # start the measure
 class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
@@ -215,6 +215,8 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
     result << OpenStudio::IdfObject.load("Output:Meter:MeterFileOnly,NaturalGas:Facility,#{reporting_frequency};").get
     result << OpenStudio::IdfObject.load("Output:Meter:MeterFileOnly,DistrictCooling:Facility,#{reporting_frequency};").get
     result << OpenStudio::IdfObject.load("Output:Meter:MeterFileOnly,DistrictHeating:Facility,#{reporting_frequency};").get
+    result << OpenStudio::IdfObject.load("Output:Meter:MeterFileOnly,Propane:Facility,#{reporting_frequency};").get
+    result << OpenStudio::IdfObject.load("Output:Meter:MeterFileOnly,FuelOilNo2:Facility,#{reporting_frequency};").get
 
     # result << OpenStudio::IdfObject.load("Output:Meter:MeterFileOnly,Cooling:Electricity,#{reporting_frequency};").get
     # result << OpenStudio::IdfObject.load("Output:Meter:MeterFileOnly,Heating:Electricity,#{reporting_frequency};").get
@@ -230,6 +232,15 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
     result << OpenStudio::IdfObject.load('Output:Variable,*,Heating Coil Heating Rate,hourly; !- HVAC Average [W];').get
     # result << OpenStudio::IdfObject.load("Output:Variable,*,Exterior Equipment:Electric Vehicles,#{reporting_frequency};").get
 
+    ## add environmental factor outputs
+    #result << OpenStudio::IdfObject.load("Output:Meter:MeterFileOnly,Output:EnvironmentalImpactFactors,#{reporting_frequency};").get
+    # result << OpenStudio::IdfObject.load("Output:Variable,*,Environmental Impact Total N2O Emissions Carbon Equivalent Mass,#{reporting_frequency}; !- HVAC Sum [kg];").get
+    # result << OpenStudio::IdfObject.load("Output:Variable,*,Environmental Impact Total CH4 Emissions Carbon Equivalent Mass,#{reporting_frequency}; !- HVAC Sum [kg];").get
+    # result << OpenStudio::IdfObject.load("Output:Variable,*,Environmental Impact Total CO2 Emissions Carbon Equivalent Mass,#{reporting_frequency}; !- HVAC Sum [kg];").get
+    # result << OpenStudio::IdfObject.load("Output:Variable,*,Environmental Impact NaturalGas CO2 Emissions Mass,#{reporting_frequency}; !- HVAC Sum [kg];").get
+    # result << OpenStudio::IdfObject.load("Output:Variable,*,Environmental Impact NaturalGas CH4 Emissions Mass,#{reporting_frequency}; !- HVAC Sum [kg];").get
+    # result << OpenStudio::IdfObject.load("Output:Variable,*,Environmental Impact NaturalGas N2O Emissions Mass,#{reporting_frequency}; !- HVAC Sum [kg];").get
+
     timeseries_data = ['District Cooling Chilled Water Rate', 'District Cooling Mass Flow Rate',
                        'District Cooling Inlet Temperature', 'District Cooling Outlet Temperature',
                        'District Heating Hot Water Rate', 'District Heating Mass Flow Rate',
@@ -240,7 +251,15 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
 
     ev_timeseries_data = ['Exterior Equipment:Electric Vehicles']
 
+    emissions_timeseries_data = ['Future_Annual_Electricity_Emissions', 'Future_Hourly_Electricity_Emissions',
+                                  'Historical_Annual_Electricity_Emissions', 'Historical_Hourly_Electricity_Emissions',
+                                  'Future_Annual_Electricity_Emissions_Intensity', 'Future_Hourly_Electricity_Emissions_Intensity',
+                                  'Historical_Annual_Electricity_Emissions_Intensity', 'Historical_Hourly_Electricity_Emissions_Intensity']
+
+
+
     timeseries_data += tes_timeseries_data
+    timeseries_data += emissions_timeseries_data
 
     timeseries_data.each do |ts|
       result << OpenStudio::IdfObject.load("Output:Variable,*,#{ts},#{reporting_frequency};").get
@@ -277,6 +296,7 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
     if value.nil?
       return nil
     end
+
     if from_units.nil? || to_units.nil?
       @runner.registerError("Cannot convert units...from_units: #{from_units} or to_units: #{to_units} left blank.")
       return nil
@@ -343,6 +363,20 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
     sql_file = sql_file.get
     model.setSqlFile(sql_file)
 
+    # Get the weather file run period (as opposed to design day run period)
+    ann_env_pd = nil
+    sql_file.availableEnvPeriods.each do |env_pd|
+      env_type = sql_file.environmentType(env_pd)
+      if env_type.is_initialized && (env_type.get == OpenStudio::EnvironmentType.new('WeatherRunPeriod'))
+        ann_env_pd = env_pd
+      end
+    end
+
+    if ann_env_pd == false
+      runner.registerError("Can't find a weather runperiod, make sure you ran an annual simulation, not just the design days.")
+      return false
+    end
+
     # get building from model
     building = model.getBuilding
 
@@ -408,8 +442,8 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
     # unconditioned_area
     unconditioned_area = sql_query(runner, sql_file, 'AnnualBuildingUtilityPerformanceSummary', "TableName='Building Area' AND RowName='Unconditioned Building Area' AND ColumnName='Area'")
     feature_report.program.unconditioned_area_sqft = convert_units(unconditioned_area, 'm^2', 'ft^2')
-    if building.standardsBuildingType.is_initialized
-      floor_area -= unconditioned_area if ['Residential'].include?(building.standardsBuildingType.get) # conditioned floor area only
+    if building.standardsBuildingType.is_initialized && ['Residential'].include?(building.standardsBuildingType.get)
+      floor_area -= unconditioned_area # conditioned floor area only
     end
 
     # maximum_number_of_stories
@@ -429,7 +463,7 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
 
     # footprint_area
     if building.standardsBuildingType.is_initialized
-      if not ['Residential'].include?(building.standardsBuildingType.get)
+      if !['Residential'].include?(building.standardsBuildingType.get)
         feature_report.program.footprint_area_sqft = feature_report.program.floor_area_sqft / number_of_stories
       else
         feature_report.program.footprint_area_sqft = convert_units(floor_area, 'm^2', 'ft^2') / building.additionalProperties.getFeatureAsInteger('NumberOfConditionedStories').get
@@ -459,6 +493,7 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
         building_type = building_type.get
       end
       next if ['Residential'].include?(building_type) # space types with empty building type fields will inherit from the building object
+
       space_type_areas[building_type] = 0 if space_type_areas[building_type].nil?
       space_type_areas[building_type] += convert_units(space_type.floorArea, 'm^2', 'ft^2')
     end
@@ -471,6 +506,7 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
       else
         building_type = space.spaceType.get.standardsBuildingType
       end
+
       if building_type.empty?
         building_type = 'unknown'
       else
@@ -689,15 +725,15 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
     # district_cooling
     district_cooling = sql_query(runner, sql_file, 'AnnualBuildingUtilityPerformanceSummary', "TableName='End Uses' AND RowName='Total End Uses' AND ColumnName='District Cooling'")
     feature_report.reporting_periods[0].district_cooling_kwh = convert_units(district_cooling, 'GJ', 'kWh')
-    if building.standardsBuildingType.is_initialized
-      feature_report.reporting_periods[0].district_cooling_kwh = 0.0 if ['Residential'].include?(building.standardsBuildingType.get)
+    if building.standardsBuildingType.is_initialized && ['Residential'].include?(building.standardsBuildingType.get)
+      feature_report.reporting_periods[0].district_cooling_kwh = 0.0
     end
 
     # district_heating
     district_heating = sql_query(runner, sql_file, 'AnnualBuildingUtilityPerformanceSummary', "TableName='End Uses' AND RowName='Total End Uses' AND ColumnName='District Heating'")
     feature_report.reporting_periods[0].district_heating_kwh = convert_units(district_heating, 'GJ', 'kWh')
-    if building.standardsBuildingType.is_initialized
-      feature_report.reporting_periods[0].district_heating_kwh = 0.0 if ['Residential'].include?(building.standardsBuildingType.get)
+    if building.standardsBuildingType.is_initialized && ['Residential'].include?(building.standardsBuildingType.get)
+      feature_report.reporting_periods[0].district_heating_kwh = 0.0
     end
 
     # water
@@ -725,10 +761,10 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
         # report each query in its corresponding feature report obeject
         x = ft.tr(' ', '_').downcase
         if x.include? 'water'
-          x_u = x + '_qbft'
-        else  
+          x_u = "#{x}_qbft"
+        else
           x = x.gsub('_#2', '')
-          x_u = x + '_kwh'
+          x_u = "#{x}_kwh"
         end
         m = feature_report.reporting_periods[0].end_uses.send(x_u)
 
@@ -739,8 +775,8 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
         end
         sql_r = convert_units(sql_r, 'GJ', 'kWh')
 
-        if building.standardsBuildingType.is_initialized
-          sql_r = 0.0 if ['Residential'].include?(building.standardsBuildingType.get) && x_u.include?('district')
+        if building.standardsBuildingType.is_initialized && (['Residential'].include?(building.standardsBuildingType.get) && x_u.include?('district'))
+          sql_r = 0.0
         end
         m.send("#{y}=", sql_r)
       end
@@ -755,7 +791,7 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
         sql = sql_query(runner, sql_file, 'AnnualBuildingUtilityPerformanceSummary', "TableName='End Uses' AND RowName='#{eu}' AND ColumnName='#{ft}'")
 
         # ensure not nil so the equations below don't error out
-        if not sql.nil?
+        if !sql.nil?
           sql_r += convert_units(sql, 'GJ', 'kWh')
         end
       end
@@ -798,24 +834,64 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
     time_setpoint_not_met_during_occupied_hours = time_setpoint_not_met_during_occupied_heating + time_setpoint_not_met_during_occupied_cooling
     feature_report.reporting_periods[0].comfort_result[:time_setpoint_not_met_during_occupied_hours] = time_setpoint_not_met_during_occupied_hours
 
-    ######################################## Reporting TImeseries Results FOR CSV File ######################################
+    # electricity emissions
+    begin
 
-    # Get the weather file run period (as opposed to design day run period)
-    ann_env_pd = nil
-    sql_file.availableEnvPeriods.each do |env_pd|
-      env_type = sql_file.environmentType(env_pd)
-      if env_type.is_initialized
-        if env_type.get == OpenStudio::EnvironmentType.new('WeatherRunPeriod')
-          ann_env_pd = env_pd
-        end
-      end
-    end
+      # future_annual_emissions
+      future_annual_emissions_ts = sql_file.timeSeries(ann_env_pd.to_s, reporting_frequency.to_s, 'Future_Annual_Electricity_Emissions', 'EMS')
+      feature_report.reporting_periods[0].emissions[:future_annual_electricity_emissions_mt] = future_annual_emissions_ts.get.values.sum
 
-    if ann_env_pd == false
-      runner.registerError("Can't find a weather runperiod, make sure you ran an annual simulation, not just the design days.")
-      return false
+      # future_hourly_emissions
+      future_hourly_emissions_ts = sql_file.timeSeries(ann_env_pd.to_s, reporting_frequency.to_s, 'Future_Hourly_Electricity_Emissions', 'EMS')
+      feature_report.reporting_periods[0].emissions[:future_hourly_electricity_emissions_mt] = future_hourly_emissions_ts.get.values.sum
+
+      # historical_annual_emissions
+      historical_annual_emissions_ts = sql_file.timeSeries(ann_env_pd.to_s, reporting_frequency.to_s, 'Historical_Annual_Electricity_Emissions', 'EMS')
+      feature_report.reporting_periods[0].emissions[:historical_annual_electricity_emissions_mt] = historical_annual_emissions_ts.get.values.sum
+
+      # historical_hourly_emissions
+      historical_hourly_emissions_ts = sql_file.timeSeries(ann_env_pd.to_s, reporting_frequency.to_s, 'Historical_Hourly_Electricity_Emissions', 'EMS')
+      feature_report.reporting_periods[0].emissions[:historical_hourly_electricity_emissions_mt] = historical_hourly_emissions_ts.get.values.sum
+
+      # future_annual_emissions_intensity
+      future_annual_emissions_intensity_ts = sql_file.timeSeries(ann_env_pd.to_s, reporting_frequency.to_s, 'Future_Annual_Electricity_Emissions_Intensity', 'EMS')
+      feature_report.reporting_periods[0].emissions[:future_annual_electricity_emissions_intensity_kg_per_ft2] = future_annual_emissions_intensity_ts.get.values.sum
+
+      # future_hourly_emissions_intensity
+      future_hourly_emissions_intensity_ts = sql_file.timeSeries(ann_env_pd.to_s, reporting_frequency.to_s, 'Future_Hourly_Electricity_Emissions_Intensity', 'EMS')
+      feature_report.reporting_periods[0].emissions[:future_hourly_electricity_emissions_intensity_kg_per_ft2] = future_hourly_emissions_intensity_ts.get.values.sum
+
+      # historical_annual_emissions_intensity
+      historical_annual_emissions_intensity_ts = sql_file.timeSeries(ann_env_pd.to_s, reporting_frequency.to_s, 'Historical_Annual_Electricity_Emissions_Intensity', 'EMS')
+      feature_report.reporting_periods[0].emissions[:historical_annual_electricity_emissions_intensity_kg_per_ft2] = historical_annual_emissions_intensity_ts.get.values.sum
+
+      # historical_hourly_emissions_intensity
+      historical_hourly_emissions_intensity_ts = sql_file.timeSeries(ann_env_pd.to_s, reporting_frequency.to_s, 'Historical_Hourly_Electricity_Emissions_Intensity', 'EMS')
+      feature_report.reporting_periods[0].emissions[:historical_hourly_electricity_emissions_intensity_kg_per_ft2] = historical_hourly_emissions_intensity_ts.get.values.sum
+
+    rescue StandardError
+      @@logger.info('Emissions are not reported for this feature')
     end
 
+    ##########################################################################################################################
+    # set conversion variables
+    conv_J_mwh = 1000000 * 60 * 60 # J to MWh (1000000J/MJ * 60hr/min * 60 min/sec)
+    conv_kg_mt = 0.001 # kg to metric ton
+    conv_kbtu_J = 1054852.32 # KBtu to J (1kBtu = 1054852.32 J)
+
+    ##### Emisison factors for natural gas, propane, and fuel oil based on EPA eGRID data and calculated using 20-year GWP horizon based on ASHRAE 189.1
+    ## natural gas :  277.358126 KG/MWH
+    ## propane : 323.896704 KG/MWH
+    ## Fuel oil : 294.962046 KG/MWH
+    nat_gas_val = 277.358126
+    lpg_val = 323.896704
+    fo1_val = 294.962046
+    fo2_val = 294.962046
+
+
+    ##########################################################################################################################
+    ######################################## Reporting TImeseries Results FOR CSV File #######################################
+
     # timeseries we want to report
     requested_timeseries_names = [
       'Electricity:Facility',
@@ -861,7 +937,31 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
       'District Heating Inlet Temperature',
       'District Heating Outlet Temperature',
       'Cooling Coil Total Cooling Rate',
-      'Heating Coil Heating Rate'
+      'Heating Coil Heating Rate',
+      'Future_Annual_Electricity_Emissions',
+      'Future_Hourly_Electricity_Emissions',
+      'Historical_Annual_Electricity_Emissions',
+      'Historical_Hourly_Electricity_Emissions',
+      'Future_Annual_Electricity_Emissions_Intensity',
+      'Future_Hourly_Electricity_Emissions_Intensity',
+      'Historical_Annual_Electricity_Emissions_Intensity',
+      'Historical_Hourly_Electricity_Emissions_Intensity',
+      'Natural_Gas_Emissions',
+      'Natural_Gas_Emissions_Intensity',
+      'Propane_Emissions', 
+      'Propane_Emissions_Intensity',
+      'FuelOilNo2_Emissions',
+      'FuelOilNo2_Emissions_Intensity',
+      'Curtailed EV Power',
+      'Daily EV Charge Energy Capacity',
+      'EV Charge Ratio',
+      'Total Charged EV Energy',
+      'Total Curtailed EV Energy',
+      'Total Scheduled EV Energy',
+      'Emission Intensity Schedule Output',
+      'EV Charging Effective Schedule',
+      'EV Charging Original Schedule',
+      'EV Charging Original Load'
     ]
 
     # add thermal comfort timeseries
@@ -880,7 +980,7 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
     runner.registerInfo("All timeseries: #{requested_timeseries_names}")
 
     # timeseries variables to keep to calculate power
-    tsToKeep = ['Electricity:Facility', 'ElectricityProduced:Facility']
+    tsToKeep = ['Electricity:Facility', 'ElectricityProduced:Facility','Propane:Facility', 'NaturalGas:Facility', 'FuelOilNo2:Facility', 'FuelOilNo1:Facility']
     tsToKeepIndexes = {}
 
     ### powerFactor ###
@@ -954,7 +1054,7 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
           # use key_value name
           # special case for Zone Thermal Comfort: use both timeseries_name and key_value
           if timeseries_name.include? 'Zone Thermal Comfort'
-            new_timeseries_name = timeseries_name + ' ' + key_value
+            new_timeseries_name = "#{timeseries_name} #{key_value}"
           else
             new_timeseries_name = key_value
           end
@@ -982,8 +1082,8 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
         end
 
         # residential considerations
-        if building.standardsBuildingType.is_initialized
-          values[key_cnt] = Array.new(n, 0) if ['DistrictCooling:Facility', 'DistrictHeating:Facility'].include?(timeseries_name) && ['Residential'].include?(building.standardsBuildingType.get)
+        if building.standardsBuildingType.is_initialized && (['DistrictCooling:Facility', 'DistrictHeating:Facility'].include?(timeseries_name) && ['Residential'].include?(building.standardsBuildingType.get))
+          values[key_cnt] = Array.new(n, 0)
         end
 
         # unit conversion
@@ -1001,7 +1101,13 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
                           'm3'
                         when 'W'
                           'W'
-                      end
+                        when 'kg'
+                          'kg'
+                        when 'MT'
+                          'MT'
+                        when 'KG/FT2'
+                          'KG/FT2'
+                     end
         end
 
         # loop through each value and apply unit conversion
@@ -1021,6 +1127,103 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
           tsToKeepIndexes[timeseries_name] = key_cnt
         end
 
+        ### add emissions for natural gas, propane and fuel oil
+        # # set conversion variables
+        # conv_J_mwh = 1000000 * 60 * 60 # J to MWh (1000000J/MJ * 60hr/min * 60 min/sec)
+        # conv_kg_mt = 0.001 # kg to metric ton
+        # conv_kbtu_J = 1054852.32 # KBtu to J (1kBtu = 1054852.32 J)
+
+        # ##### Emisison factors for natural gas, propane, and fuel oil based on EPA eGRID data and calculated using 20-year GWP horizon based on ASHRAE 189.1
+        # ## natural gas :  277.358126 KG/MWH
+        # ## propane : 323.896704 KG/MWH
+        # ## Fuel oil : 294.962046 KG/MWH
+        # nat_gas_val = 277.358126
+        # lpg_val = 323.896704
+        # fo1_val = 294.962046
+        # fo2_val = 294.962046
+
+        if timeseries_name == 'Natural_Gas_Emissions'
+          newVals = Array.new(n, 0)
+          (0..n - 1).each do |j|
+            newVals[j] = (nat_gas_val * (values[tsToKeepIndexes['NaturalGas:Facility']][j].to_f * conv_kbtu_J.to_f) / conv_J_mwh.to_f ) * conv_kg_mt.to_f 
+            j += 1
+          end
+          new_unit = 'MT'
+          values[key_cnt] = newVals
+
+          # add emissions sum to feature report
+          feature_report.reporting_periods[0].emissions[:natural_gas_emissions_mt] = newVals.sum
+        end
+
+        if timeseries_name == 'Propane_Emissions'
+          newVals = Array.new(n, 0)
+          (0..n - 1).each do |j|
+            newVals[j] = (lpg_val * (values[tsToKeepIndexes['Propane:Facility']][j].to_f * conv_kbtu_J.to_f) / conv_J_mwh.to_f ) * conv_kg_mt.to_f 
+            j += 1
+          end
+          new_unit = 'MT'
+          values[key_cnt] = newVals
+
+          # add emissions sum to feature report
+          feature_report.reporting_periods[0].emissions[:propane_emissions_mt] = newVals.sum
+        end
+
+        if timeseries_name == 'FuelOilNo2_Emissions'
+          newVals = Array.new(n, 0)
+          (0..n - 1).each do |j|
+            newVals[j] = (fo2_val * (values[tsToKeepIndexes['FuelOilNo2:Facility']][j].to_f * conv_kbtu_J.to_f) / conv_J_mwh.to_f ) * conv_kg_mt.to_f 
+            j += 1
+          end
+          new_unit = 'MT'
+          values[key_cnt] = newVals
+
+          # add emissions sum to feature report
+          feature_report.reporting_periods[0].emissions[:fueloil_no2_emissions_mt] = newVals.sum
+        end
+
+        ### calculate emissions intensity metric
+        # get flr_area
+        flr_area = building.floorArea * 10.764 #change from m2 to ft2
+
+        if timeseries_name == 'Natural_Gas_Emissions_Intensity'
+          newVals = Array.new(n, 0)
+          (0..n - 1).each do |j|
+            newVals[j] = (((nat_gas_val * (values[tsToKeepIndexes['NaturalGas:Facility']][j].to_f * conv_kbtu_J.to_f) / conv_J_mwh.to_f ) * conv_kg_mt.to_f ) * 1000 / flr_area) # unit: kg/ft2 - changed mt to kg
+            j += 1
+          end
+          new_unit = 'KG/FT2'
+          values[key_cnt] = newVals
+
+          # add emissions sum to feature report
+          feature_report.reporting_periods[0].emissions[:natural_gas_emissions_intensity_kg_per_ft2] = newVals.sum
+        end
+
+        if timeseries_name == 'Propane_Emissions_Intensity'
+          newVals = Array.new(n, 0)
+          (0..n - 1).each do |j|
+            newVals[j] = (((lpg_val * (values[tsToKeepIndexes['Propane:Facility']][j].to_f * conv_kbtu_J.to_f) / conv_J_mwh.to_f ) * conv_kg_mt.to_f ) * 1000 / flr_area) # unit: kg/ft2 - changed mt to kg
+            j += 1
+          end
+          new_unit = 'KG/FT2'
+          values[key_cnt] = newVals
+
+          # add emissions sum to feature report
+          feature_report.reporting_periods[0].emissions[:propane_emissions_intensity_kg_per_ft2] = newVals.sum
+        end
+
+        if timeseries_name == 'FuelOilNo2_Emissions_Intensity'
+          newVals = Array.new(n, 0)
+          (0..n - 1).each do |j|
+            newVals[j] = (((fo2_val * (values[tsToKeepIndexes['FuelOilNo2:Facility']][j].to_f * conv_kbtu_J.to_f) / conv_J_mwh.to_f ) * conv_kg_mt.to_f ) * 1000 / flr_area) # unit: kg/ft2 - changed mt to kg
+            j += 1
+          end
+          new_unit = 'KG/FT2'
+          values[key_cnt] = newVals
+
+          # add emissions sum to feature report
+          feature_report.reporting_periods[0].emissions[:fueloil_no2_emissions_intensity_kg_per_ft2] = newVals.sum
+        end
+
         # special processing: power
         if powerTimeseries.include? timeseries_name
           # special case: net series (subtract generation from load)
@@ -1077,7 +1280,8 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
               end
             end
           end
-        end
+        end          
+          
 
         # append units to headers
         new_timeseries_name += "(#{new_unit})"
@@ -1210,6 +1414,8 @@ class DefaultFeatureReports < OpenStudio::Measure::ReportingMeasure
       end
     end
 
+    puts "values = #{values}"
+
     # closing the sql file
     sql_file.close
 

data/lib/measures/export_modelica_loads/measure.rb

@@ -1,5 +1,5 @@
 # *********************************************************************************
-# URBANopt™, Copyright (c) 2019-2021, Alliance for Sustainable Energy, LLC, and other
+# URBANopt™, Copyright (c) 2019-2022, Alliance for Sustainable Energy, LLC, and other
 # contributors. All rights reserved.
 
 # Redistribution and use in source and binary forms, with or without modification,
@@ -337,6 +337,7 @@ class ExportModelicaLoads < OpenStudio::Measure::ReportingMeasure
       f << "double tab1(8760,4)\n"
       modelica_data.each_with_index do |row, index|
         next if index.zero?
+
         f << row.join(';') << "\n"
       end
     end