urbanopt-reopt 1.0.0 → 1.2.0

data/lib/urbanopt/reopt/{reopt_ghp_adapter.rb → reopt_ghp_adapter_ghp.rb}

@@ -1,5 +1,5 @@
 # *********************************************************************************
-# URBANopt (tm), Copyright (c) Alliance for Sustainable Energy, LLC.
+# URBANopt (tm), Copyright (c) Alliance for Energy Innovation, LLC.
 # See also https://github.com/urbanopt/urbanopt-reopt-gem/blob/develop/LICENSE.md
 # *********************************************************************************
 
@@ -10,11 +10,16 @@ module URBANopt # :nodoc:
       def initialize
         # initialize @@logger
         @@logger ||= URBANopt::REopt.reopt_logger
-        # Define class variable
+        # Define class variables 
         @@hours_in_year = 8760
+        @@nat_gas_dollars_per_mmbtu = 13.5
+        @@year_of_simulation = 2023
+        @@year_of_simulation = 2023
+        @@small_multiplier = [0.00001]
+
       end
 
-      def create_reopt_input_building(run_dir, system_parameter_hash, reopt_ghp_assumptions_hash, building_id, modelica_result)
+      def create_reopt_input_building_ghp(run_dir, system_parameter_hash, reopt_ghp_assumptions_hash, building_id, modelica_result)
 
         # Define variables
         reopt_inputs_building = {}
@@ -28,9 +33,6 @@ module URBANopt # :nodoc:
             SpaceHeatingLoad: {},
             DomesticHotWaterLoad: {},
             ElectricLoad: {},
-            ElectricTariff: {
-              urdb_label: ""
-            },
             GHP: {},
             ExistingBoiler: {}
           }
@@ -38,7 +40,7 @@ module URBANopt # :nodoc:
 
         # The URDB label is required to be specified in the input assumption file
         if reopt_inputs_building[:ElectricTariff][:urdb_label].nil? || reopt_inputs_building[:ElectricTariff][:urdb_label].empty?
-          raise "Missing value for urdb_label - this is a required input"
+         raise "Missing value for urdb_label - this is a required input"
         end
 
         scenario_json_path = File.join(run_dir, "default_scenario_report.json")
@@ -50,7 +52,6 @@ module URBANopt # :nodoc:
             @longitude = scenario_json_data[:scenario_report][:location][:longitude_deg]
             reopt_inputs_building[:Site][:latitude] = @latitude
             reopt_inputs_building[:Site][:longitude] = @longitude
-
           end
         end
 
@@ -72,8 +73,8 @@ module URBANopt # :nodoc:
           end
             # Check if the total kBtu is zero
           if total_kbtu.zero?
-            # If zero, populate with hourly values meet reopts formatting requirements
-            reopt_inputs_building[:SpaceHeatingLoad][:fuel_loads_mmbtu_per_hour] = [0.000001] * @@hours_in_year
+            # If zero, populate with near zero hourly values to meet reopts formatting requirements
+            reopt_inputs_building[:SpaceHeatingLoad][:fuel_loads_mmbtu_per_hour] = @@small_multiplier * @@hours_in_year
           else
             # If not zero, convert and append to the array
             timeseries_data.each do |row|
@@ -86,8 +87,8 @@ module URBANopt # :nodoc:
           end
 
         else
-          # Calculate space heating load
-          reopt_inputs_building[:SpaceHeatingLoad][:fuel_loads_mmbtu_per_hour] = [0.000001] * @@hours_in_year
+          # populate with near zero hourly values to meet reopts formatting requirements
+          reopt_inputs_building[:SpaceHeatingLoad][:fuel_loads_mmbtu_per_hour] = @@small_multiplier * @@hours_in_year
           puts "Existing heating fuel cost was not taken into consideration in result calculations."
         end
 
@@ -103,12 +104,21 @@ module URBANopt # :nodoc:
 
         if File.exist?(@modelica_csv)
           modelica_data = CSV.read(@modelica_csv, headers: true)
-          heating_power = "heating_electric_power_#{building_id}"
-          cooling_power = "cooling_electric_power_#{building_id}"
+          heating_power_header = modelica_data.headers.select{ |h| h.to_s.start_with?("heating_electric_power") }
+          prefix = ""
+          
+          if heating_power_header[0].split("_")[-1].include? ("B")
+            prefix = "B"
+          else
+            prefix = ""
+          end
+
+          heating_power = "heating_electric_power_#{prefix}#{building_id}"
+          cooling_power = "cooling_electric_power_#{prefix}#{building_id}"
+          ets_pump_power = "ets_pump_power_#{prefix}#{building_id}"
+          heating_system_capacity = "heating_system_capacity_#{prefix}#{building_id}"
+          cooling_system_capacity = "cooling_system_capacity_#{prefix}#{building_id}"
           pump_power = "pump_power_#{building_id}"
-          ets_pump_power = "ets_pump_power_#{building_id}"
-          heating_system_capacity = "heating_system_capacity_#{building_id}"
-          cooling_system_capacity = "cooling_system_capacity_#{building_id}"
 
           heating_power_values = cooling_power_values = pump_power_values = ets_pump_power_values = []
           total_electric_load_building = []
@@ -130,7 +140,6 @@ module URBANopt # :nodoc:
             # Convert watts to kilowatts
             elements.map { |e| e.to_f / 1000 }.sum
           end
-
           peak_combined_heatpump_thermal_ton = 0
 
           if modelica_data.headers.include?(heating_system_capacity)
@@ -139,13 +148,13 @@ module URBANopt # :nodoc:
           if modelica_data.headers.include?(cooling_system_capacity)
             cooling_system_capacity_value = modelica_data[cooling_system_capacity][0]
           end
-
+          
           watts_per_ton_cooling_capacity = 3517
           peak_combined_heatpump_thermal_ton = ([heating_system_capacity_value.to_f.abs, cooling_system_capacity_value.to_f.abs].max) / watts_per_ton_cooling_capacity
 
           # Store the result in reopt_inputs_building ElectricLoad
           reopt_inputs_building[:ElectricLoad][:loads_kw] = total_electric_load_building
-
+          reopt_inputs_building[:ElectricLoad][:year] = @@year_of_simulation
 
           domestic_hot_water = total_electric_load_building.map do |load|
             load * 0
@@ -156,16 +165,29 @@ module URBANopt # :nodoc:
 
           # Add GHP Fields
           reopt_inputs_building[:GHP] = {}
+
+          # Add avoided capital cost of all buildings
+          if reopt_inputs_building[:features] && !reopt_inputs_building[:features].empty?
+            # Find the feature matching this building_id
+            matching_feature = reopt_inputs_building[:features].find { |feature| feature[:feature_id] == building_id.to_i}
+
+            if matching_feature
+              # Set avoided capex value into GHP block
+              reopt_inputs_building[:GHP][:avoided_capex_by_ghp_present_value] = matching_feature[:avoided_capex_by_ghp_present_value]
+            end
+          end
           # REopt default
           reopt_inputs_building[:GHP][:require_ghp_purchase] = 1
           reopt_inputs_building[:GHP][:om_cost_per_sqft_year] = 0
           reopt_inputs_building[:GHP][:heatpump_capacity_sizing_factor_on_peak_load] = 1.0
           # Add the floor area
           building_json_path = File.join(run_dir, building_id.to_s, "feature_reports", "default_feature_report.json")
+
           if File.exist?(building_json_path)
+            puts building_json_path
             File.open(building_json_path, 'r') do |file|
               building_json_data = JSON.parse(file.read, symbolize_names: true)
-              reopt_inputs_building[:GHP][:building_sqft] = building_json_data[:program][:floor_area_sqft]
+              reopt_inputs_building[:GHP][:building_sqft] = building_json_data[:program][:footprint_area_sqft].to_f
             end
           else
             puts "File not found: #{building_json_path}"
@@ -173,8 +195,8 @@ module URBANopt # :nodoc:
 
           # Add existing boiler fuel cost
           # TODO : Add this as optional user input
-          nat_gas_dollars_per_mmbtu = 13.5
-          reopt_inputs_building[:ExistingBoiler][:fuel_cost_per_mmbtu] = nat_gas_dollars_per_mmbtu
+          # Cost of Natural Gas in $/mmbtu as per REopt Defaults
+          reopt_inputs_building[:ExistingBoiler][:fuel_cost_per_mmbtu] = @@nat_gas_dollars_per_mmbtu
 
           # Add ghpghx_responses
           ghpghx_output = {}
@@ -192,7 +214,8 @@ module URBANopt # :nodoc:
           ghpghx_output[:outputs][:yearly_heating_heatpump_electric_consumption_series_kw] = total_electric_load_building
           ghpghx_output[:outputs][:yearly_cooling_heatpump_electric_consumption_series_kw] = [0] * @@hours_in_year
           # This is not used in REopt calculation but required for formatting.
-          ghpghx_output[:inputs][:heating_thermal_load_mmbtu_per_hr] = [0.0001] * @@hours_in_year
+          # populate with near zero hourly values to meet reopts formatting requirements
+          ghpghx_output[:inputs][:heating_thermal_load_mmbtu_per_hr] = @@small_multiplier * @@hours_in_year
           # This is not used in REopt calculation but required for formatting.
           ghpghx_output[:inputs][:cooling_thermal_load_ton] = [0] * @@hours_in_year
 
@@ -210,7 +233,7 @@ module URBANopt # :nodoc:
 
       end
 
-      def create_reopt_input_district(run_dir, system_parameter_hash, reopt_ghp_assumptions_hash, ghp_id, modelica_result)
+      def create_reopt_input_district_ghp(run_dir, system_parameter_hash, reopt_ghp_assumptions_hash, ghp_id, modelica_result)
 
         reopt_inputs_district = {}
 
@@ -238,27 +261,29 @@ module URBANopt # :nodoc:
         # The URDB label is required to be specified in the input assumption file
         if reopt_inputs_district[:ElectricTariff][:urdb_label].nil? || reopt_inputs_district[:ElectricTariff][:urdb_label].empty?
 
-          raise "Missing value for urdb_label - this is a required input"
+         raise "Missing value for urdb_label - this is a required input"
 
         end
-
+        # populate with near zero hourly values to meet reopts formatting requirements
         # This is not used in REopt calculation but required for formatting.
-        reopt_inputs_district[:SpaceHeatingLoad][:fuel_loads_mmbtu_per_hour] = [0.000001]*@@hours_in_year
+        reopt_inputs_district[:SpaceHeatingLoad][:fuel_loads_mmbtu_per_hour] = @@small_multiplier*@@hours_in_year
+        # populate with near zero hourly values to meet reopts formatting requirements
         # This is not used in REopt calculation but required for formatting.
-        reopt_inputs_district[:DomesticHotWaterLoad][:fuel_loads_mmbtu_per_hour] = [0.0000001]*@@hours_in_year
+        reopt_inputs_district[:DomesticHotWaterLoad][:fuel_loads_mmbtu_per_hour] = @@small_multiplier*@@hours_in_year
 
-        # Adding year for ElectricLoad so district electric load can be calculated with REopt API v3.11
-        reopt_inputs_district[:ElectricLoad] = {:year => 2017}
-        #required for reopt formatting
-        reopt_inputs_district[:ElectricLoad][:loads_kw] = [0.00001]*@@hours_in_year
+        reopt_inputs_district[:ElectricLoad] = {}
+        # populate with near zero hourly values to meet reopts formatting requirements
+        # This is not used in REopt calculation but required for formatting.
+        reopt_inputs_district[:ElectricLoad][:loads_kw] = @@small_multiplier*@@hours_in_year
+        reopt_inputs_district[:ElectricLoad][:year] = @@year_of_simulation
 
         reopt_inputs_district[:ExistingBoiler] = {}
-        reopt_inputs_district[:ExistingBoiler][:fuel_cost_per_mmbtu] = 13.5
+        reopt_inputs_district[:ExistingBoiler][:fuel_cost_per_mmbtu] = @@nat_gas_dollars_per_mmbtu
 
         # GHP inputs
         reopt_inputs_district[:GHP] = {}
         reopt_inputs_district[:GHP][:require_ghp_purchase] = 1
-        reopt_inputs_district[:GHP][:building_sqft] = 0.00001
+        reopt_inputs_district[:GHP][:building_sqft] = format('%.8f', @@small_multiplier[0]).to_f
         reopt_inputs_district[:GHP][:om_cost_per_sqft_year] = 0
         reopt_inputs_district[:GHP][:heatpump_capacity_sizing_factor_on_peak_load] = 1.0
 
@@ -272,15 +297,39 @@ module URBANopt # :nodoc:
 
 
         # Read GHX sizes from system parameter hash
-        ghe_specific_params = system_parameter_hash[:district_system][:fifth_generation][:ghe_parameters][:ghe_specific_params]
-        ghe_specific_params.each do |ghe_specific_param|
-          if ghe_specific_param[:ghe_id] = ghp_id
-            number_of_boreholes = ghe_specific_param[:borehole][:number_of_boreholes]
-            length_of_boreholes = ghe_specific_param[:borehole][:length_of_boreholes]
-            ghpghx_output[:outputs][:number_of_boreholes] = number_of_boreholes
-            # convert meters to feet
-            ghpghx_output[:outputs][:length_boreholes_ft] = (length_of_boreholes)*3.28084
-          end
+        ghe_specific_params = system_parameter_hash[:district_system][:fifth_generation][:ghe_parameters][:borefields]
+        
+        ghe_specific_params.each do |ghe|
+          if ghe[:ghe_id] == ghp_id
+            unless ghe[:pre_designed_borefield]
+              if ghe[:autosized_rectangle_borefield]
+                borefield = ghe[:autosized_rectangle_borefield]
+
+              elsif ghe[:autosized_rectangle_constrained_borefield]
+                borefield = ghe[:autosized_rectangle_constrained_borefield]
+
+              elsif ghe[:autosized_birectangle_borefield]
+                borefield = ghe[:autosized_birectangle_borefield]
+
+              elsif ghe[:autosized_birectangle_constrained_borefield]
+                borefield = ghe[:autosized_birectangle_constrained_borefield]
+
+              elsif ghe[:autosized_bizoned_rectangle_borefield]
+                borefield = ghe[:autosized_bizoned_rectangle_borefield]
+
+              elsif ghe[:autosized_near_square_borefield]
+                borefield = ghe[:autosized_near_square_borefield]
+
+              elsif ghe[:autosized_rowwise_borefield]
+                borefield = ghe[:autosized_rowwise_borefield]
+              end
+            end
+
+            ghpghx_output[:outputs][:number_of_boreholes] = borefield[:number_of_boreholes]
+            # convert meters to feet by multiplying with 3.28084
+            ghpghx_output[:outputs][:length_boreholes_ft] = (borefield[:borehole_length])*3.28084
+
+          end    
         end
 
         if File.exist?(@modelica_csv)
@@ -309,16 +358,19 @@ module URBANopt # :nodoc:
           # column_values = modelica_data.by_col[ghp_column]
 
           ghpghx_output[:outputs][:yearly_ghx_pump_electric_consumption_series_kw] = electrical_power_consumed_kw
-        else
-          ghpghx_output[:outputs][:yearly_ghx_pump_electric_consumption_series_kw] = [0.000000001]*@@hours_in_year
+        else        
+          # populate with near zero hourly values to meet reopts formatting requirements
+          # This is not used in REopt calculation but required for formatting.
+          ghpghx_output[:outputs][:yearly_ghx_pump_electric_consumption_series_kw] =  @@small_multiplier*@@hours_in_year
         end
 
         # This is not used in REopt calculation but required for formatting.
         ghpghx_output[:outputs][:peak_combined_heatpump_thermal_ton] = 0.000000001
 
         ghpghx_output[:outputs][:heat_pump_configuration] = "WSHP"
-        # Required for REpot formatting
-        ghpghx_output[:outputs][:yearly_total_electric_consumption_series_kw] = [0.00001] * @@hours_in_year
+        # populate with near zero hourly values to meet reopts formatting requirements
+        # This is not used in REopt calculation but required for formatting.
+        ghpghx_output[:outputs][:yearly_total_electric_consumption_series_kw] = @@small_multiplier * @@hours_in_year
         ghpghx_output[:outputs][:yearly_heating_heatpump_electric_consumption_series_kw] = [0] * @@hours_in_year
         ghpghx_output[:outputs][:yearly_cooling_heatpump_electric_consumption_series_kw] = [0] * @@hours_in_year
 
@@ -332,6 +384,144 @@ module URBANopt # :nodoc:
         File.write(json_file_path, pretty_json)
       end
 
+      def create_reopt_input_building_bau(run_dir, system_parameter_hash, reopt_ghp_assumptions_hash, building_id, modelica_result)
+        reopt_inputs_building_bau = {}
+        if !reopt_ghp_assumptions_hash.nil?
+          reopt_inputs_building_bau = reopt_ghp_assumptions_hash
+        else
+          @@logger.info('Using default REopt assumptions')
+          # create a dictionary for REopt Inputs
+          reopt_inputs_building_bau = {
+            Site: {},
+            SpaceHeatingLoad: {},
+            DomesticHotWaterLoad: {},
+            ElectricLoad: {},
+            ElectricTariff: {
+              urdb_label: ""
+            },
+            ExistingBoiler: {}
+          }
+        end
+
+        # The URDB label is required to be specified in the input assumption file
+        if reopt_inputs_building_bau.nil? || reopt_inputs_building_bau[:ElectricTariff].nil? || reopt_inputs_building_bau[:ElectricTariff][:urdb_label].nil? || reopt_inputs_building_bau[:ElectricTariff][:urdb_label].empty?
+          raise "Missing value for urdb_label - this is a required input"
+        end
+
+        scenario_json_path = File.join(run_dir, "default_scenario_report.json")
+        if File.exist?(scenario_json_path)
+          File.open(scenario_json_path, 'r') do |file|
+            scenario_json_data = JSON.parse(file.read, symbolize_names: true)
+            # update site location
+            @latitude = scenario_json_data[:scenario_report][:location][:latitude_deg]
+            @longitude = scenario_json_data[:scenario_report][:location][:longitude_deg]
+            reopt_inputs_building_bau[:Site][:latitude] = @latitude
+            reopt_inputs_building_bau[:Site][:longitude] = @longitude
+
+          end
+        end
+
+        reopt_inputs_building_bau[:SpaceHeatingLoad][:fuel_loads_mmbtu_per_hour] = []
+        # Read the default csv report
+        default_feature_report_path = File.join(run_dir, building_id.to_s, "feature_reports", "default_feature_report.csv")
+        if File.exist?(default_feature_report_path)
+          timeseries_data = CSV.read(default_feature_report_path, headers: true)
+
+          # Initialize the total kBtu sum
+          heating_kbtu = 0.0
+          cooling_kbtu = 0.0
+          total_kwh_heating = 0.0
+          total_kwh_cooling = 0.0
+
+          # TODO : add other Heating Fuels if present (Heating:Propane(kBtu) etc.)
+          # Convert each value in "Heating:NaturalGas(kBtu)" to MMBtu and store in the array
+          timeseries_data.each do |row|
+            if row['Heating:NaturalGas(kBtu)'] # Ensure the value exists
+              heating_kBtu_value = row['Heating:NaturalGas(kBtu)'].to_f # Convert to float
+              heating_kbtu += heating_kBtu_value # Sum kBtu values
+            end
+            if row.headers.include?('Cooling:NaturalGas(kBtu)') && row['Cooling:NaturalGas(kBtu)']
+              cooling_kBtu_value = row['Cooling:NaturalGas(kBtu)'].to_f # Convert to float
+              cooling_kbtu += cooling_kBtu_value
+            end
+            if row['Heating:Electricity(kWh)']
+              heating_value = row['Heating:Electricity(kWh)'].to_f # Convert to float
+              total_kwh_heating += heating_value # Sum heating values
+            end
+            if row['Cooling:Electricity(kWh)']
+              cooling_value = row['Cooling:Electricity(kWh)'].to_f # Convert to float
+              total_kwh_cooling += cooling_value # Sum cooling values
+            end
+          end
+          # Check if the heating  kBtu is zero
+          if heating_kbtu.zero?
+            # If zero, populate with near zero hourly values to meet reopts formatting requirements
+            reopt_inputs_building_bau[:SpaceHeatingLoad][:fuel_loads_mmbtu_per_hour] = @@small_multiplier * @@hours_in_year
+          else
+            # If not zero, convert and append to the array
+            timeseries_data.each do |row|
+              if row['Heating:NaturalGas(kBtu)'] # Ensure the value exists
+                kBtu_value = row['Heating:NaturalGas(kBtu)'].to_f # Convert to float
+                mMBtu_value = kBtu_value / 1000 # Convert kBtu to MMBtu
+                reopt_inputs_building_bau[:SpaceHeatingLoad][:fuel_loads_mmbtu_per_hour] << mMBtu_value # Append to the array
+              end
+            end
+            # Add fuel cost for existing boiler
+            reopt_inputs_building_bau[:ExistingBoiler][:fuel_cost_per_mmbtu] = @@nat_gas_dollars_per_mmbtu
+          end
+
+
+          # if cooling_kbtu.zero?
+          #   # If zero, populate with near zero hourly values to meet reopts formatting requirements
+          #   reopt_inputs_building_bau[:CoolingLoad][:fuel_loads_mmbtu_per_hour] = @@small_multiplier * @@hours_in_year
+          # else
+
+          # if there is cooling through natural gas it needs to be added to coolingload. If there is electricity based cooling, it is included in electric load
+          if !cooling_kbtu.zero?
+            # Add fuel load values for cooling
+            reopt_inputs_building_bau[:CoolingLoad] = {}
+            # If not zero, convert and append to the array
+            timeseries_data.each do |row|
+              if row['Cooling:NaturalGas(kBtu)'] # Ensure the value exists
+                kBtu_value = row['Cooling:NaturalGas(kBtu)'].to_f # Convert to float
+                ton_value = kBtu_value / 12 # Convert kBtu to ton
+                reopt_inputs_building_bau[:CoolingLoad][:thermal_loads_ton] << ton_value # Append to the array
+              end
+            end
+            # Add fuel cost for existing chiller
+            reopt_inputs_building_bau[:ExistingChiller][:fuel_cost_per_mmbtu] = @@nat_gas_dollars_per_mmbtu
+          end
+
+          total_kwh_load = total_kwh_heating + total_kwh_cooling
+          reopt_inputs_building_bau[:ElectricLoad][:year] = @@year_of_simulation
+
+          # Check if the total kwh is zero
+          if total_kwh_load.zero?
+            # This is the total of heating and cooling electricity loads
+            reopt_inputs_building_bau[:ElectricLoad][:loads_kw] = @@small_multiplier * @@hours_in_year
+          else
+            total_value_kwh = [] # Initialize the array
+
+            timeseries_data.each do |row|
+              heating_value_kwh = row['Heating:Electricity(kWh)'] ? row['Heating:Electricity(kWh)'].to_f : 0.0
+              cooling_value_kwh = row['Cooling:Electricity(kWh)'] ? row['Cooling:Electricity(kWh)'].to_f : 0.0
+
+              total = heating_value_kwh + cooling_value_kwh # Sum the values
+              total_value_kwh << total # Append to the array
+              reopt_inputs_building_bau[:ElectricLoad][:loads_kw] = total_value_kwh
+            end
+          end
+
+        end
+
+        #save output report in reopt_ghp directory
+        reopt_ghp_dir = File.join(run_dir, "reopt_ghp", "reopt_ghp_inputs")
+        json_file_path = File.join(reopt_ghp_dir, "BAU_building_#{building_id}.json")
+        pretty_json = JSON.pretty_generate(reopt_inputs_building_bau)
+        File.write(json_file_path, pretty_json)
+      
+      end
+
     end
   end
 end

data/lib/urbanopt/reopt/reopt_ghp_api.rb

@@ -1,5 +1,5 @@
 # *********************************************************************************
-# URBANopt (tm), Copyright (c) Alliance for Sustainable Energy, LLC.
+# URBANopt (tm), Copyright (c) Alliance for Energy Innovation, LLC.
 # See also https://github.com/urbanopt/urbanopt-reopt-gem/blob/develop/LICENSE.md
 # *********************************************************************************
 
@@ -22,7 +22,7 @@ module URBANopt # :nodoc:
                 if [nil, '', '<insert your key here>'].include? DEVELOPER_NREL_KEY
                     raise 'A developer.nrel.gov API key is required. Please see https://developer.nrel.gov/signup/ then update the file urbanopt-reopt-gem/developer_nrel_key.rb'
                 else
-                    #Store the NREL developer key
+                    #Store the NLR developer key
                     nrel_developer_key = DEVELOPER_NREL_KEY
                 end
             end
@@ -111,7 +111,7 @@ module URBANopt # :nodoc:
             run_id
         end
 
-        def reopt_request(results_url, poll_interval = 5, max_timeout = 300)
+        def reopt_request(results_url, poll_interval = 5, max_timeout = 600)
 
             key_error_count = 0
             key_error_threshold = 3

data/lib/urbanopt/reopt/reopt_ghp_files/reopt_ghp_assumption.json

@@ -14,9 +14,9 @@
         "urdb_label": "594976725457a37b1175d089"
     },
     "GHP":{
-        "installed_cost_heatpump_per_ton": 1075,
-        "installed_cost_ghx_per_ft": 14,
-        "installed_cost_building_hydronic_loop_per_sqft": 1.7,
+        "installed_cost_heatpump_per_ton": 1904,
+        "installed_cost_ghx_per_ft": 17,
+        "installed_cost_building_hydronic_loop_per_sqft": 0,
         "om_cost_per_sqft_year": 0,
         "macrs_bonus_fraction": 0.6,
         "macrs_itc_reduction": 0.5,

data/lib/urbanopt/reopt/reopt_ghp_post_processor.rb

@@ -1,11 +1,12 @@
 # *********************************************************************************
-# URBANopt (tm), Copyright (c) Alliance for Sustainable Energy, LLC.
+# URBANopt (tm), Copyright (c) Alliance for Energy Innovation, LLC.
 # See also https://github.com/urbanopt/urbanopt-reopt-gem/blob/develop/LICENSE.md
 # *********************************************************************************
 
 require 'bundler/setup'
 require 'urbanopt/reopt/reopt_logger'
 require 'urbanopt/reopt/reopt_ghp_api'
+require 'urbanopt/reopt/reopt_ghp_result'
 require 'csv'
 require 'json'
 require 'fileutils'
@@ -115,12 +116,14 @@ module URBANopt # :nodoc:
         end
 
         building_ids.each do |building_id|
-          # create REopt building input file for all buildings in loop order list
-          reopt_input_building = adapter.create_reopt_input_building(@run_dir, @system_parameter_input_hash, @reopt_ghp_assumptions_input_hash, building_id, @modelica_result_input)
+          # create REopt building input file for all buildings in loop order list in GHP scenario
+          reopt_input_building = adapter.create_reopt_input_building_ghp(@run_dir, @system_parameter_input_hash, @reopt_ghp_assumptions_input_hash, building_id, @modelica_result_input)
+          #create REopt building input file for all buildings in loop order list in BAU scenario
+          reopt_input_building_bau = adapter.create_reopt_input_building_bau(@run_dir, @system_parameter_input_hash, @reopt_ghp_assumptions_input_hash, building_id, @modelica_result_input)
         end
         ghp_ids.each do |ghp_id|
           # create REopt district input file
-          reopt_input_district = adapter.create_reopt_input_district(@run_dir, @system_parameter_input_hash, @reopt_ghp_assumptions_input_hash, ghp_id, @modelica_result_input)
+          reopt_input_district = adapter.create_reopt_input_district_ghp(@run_dir, @system_parameter_input_hash, @reopt_ghp_assumptions_input_hash, ghp_id, @modelica_result_input)
         end
 
         Dir.foreach(reopt_ghp_input) do |input_file|
@@ -142,8 +145,14 @@ module URBANopt # :nodoc:
           # call the REopt API
           api = URBANopt::REopt::REoptLiteGHPAPI.new(reopt_input_data, DEVELOPER_NREL_KEY, reopt_output_file, @localhost)
           api.get_api_results
+
         end
+
+        ## POST PROCESS RESULTS
+        ghp_results = URBANopt::REopt::REoptGHPResult.new
+        results = ghp_results.result_calculate(reopt_ghp_dir)
       end
+
     end # REoptGHPPostProcessor
   end # REopt
 end # URBANopt