urbanopt-reopt 0.12.0 → 1.0.0

checksums.yaml

@@ -1,7 +1,7 @@
 ---
 SHA256:
-  metadata.gz: ebc2f3c6f40a363e9f234917e55bd165eab74734c2ae477dd393893fb6210a3d
-  data.tar.gz: a2502bf5a657e7d256e39be74426f8eb91d7774f64851274a98c67347a33de12
+  metadata.gz: d7585e601e08ef8ae4b677124c3f6790129b44dbfe7da4c701168698437df94c
+  data.tar.gz: fb375bcd95645f6821ee8f06fce5d27eb682f1253db82d71ada98b9f043c921f
 SHA512:
-  metadata.gz: dfad70c6283e27679164c2b6c14dcbb7a009643a838614a4d31f04aca864f548324cea4dfa0978e02eba29bd7705ea32a31867763aca16120a90205cf10ad241
-  data.tar.gz: 243aca41f2c5f45970f346cb7d430a7443fd8f0a5db08b59e43db23ef8a6063d9aefd56e4a5b45a54cb89d708c147531e34424f8172c23276943faa987169b74
+  metadata.gz: a8d65764ca0eb6ec58a63f188882c2e70ba5457299f7ddbac75c7a98cc5d2762f62a82f1ecb7cb5d61c330d9ad29d6e7a38d1b24ee429c5d30ef910199f1919c
+  data.tar.gz: c91ec78859161c11e381771b0a5f231e7ee168dc758f49f9e80035e15483ca8c1b95c42a32334dca830f60c67d167c5016cc6ea776be2cff39180bde42091852

data/.github/workflows/nightly_ci_build.yml

@@ -2,13 +2,11 @@ name: REopt-gem CI
 
 on:
   workflow_dispatch:
-  # push:
+  push:
   schedule:
     # https://docs.github.com/en/actions/using-workflows/events-that-trigger-workflows#schedule
     # 5:23 am UTC (11:23pm MDT the day before) every weekday night in MDT
     - cron: '23 5 * * 2-6'
-  pull_request:
-    types: [review_requested]
 
 env:
   # This env var should enforce develop branch of all dependencies
@@ -21,17 +19,16 @@ jobs:
     # https://github.com/rbenv/ruby-build/discussions/1940
     runs-on: ubuntu-latest
     container:
-      image: docker://nrel/openstudio:3.7.0
+      image: docker://nrel/openstudio:3.9.0
     steps:
       - uses: actions/checkout@v4
       - name: set git config options
         shell: bash
-        run: |
-          git config --global --add safe.directory '*'
+        run: git config --global --add safe.directory '*'
       - name: Update gems
         run: |
+          bundle install
           bundle update
-          bundle exec certified-update
       - name: Run Rspec
         run: bundle exec rspec
       - name: Coveralls

data/.gitignore

@@ -16,7 +16,7 @@
 /lib/measures/test_results
 /lib/measures/.rubocop*
 # test files
-spec/run/
+spec/files/run/baseline_scenario_ghe/reopt_ghp/
 
 developer_nrel_key.rb
 # rspec failure tracking

data/CHANGELOG.md

@@ -1,12 +1,29 @@
 # URBANopt REopt Gem
 
+## Version 1.0.0
+
+## What's Changed
+
+* Upgrade to OpenStudio 3.9 and Ruby 3.2.2 by @vtnate in https://github.com/urbanopt/urbanopt-reopt-gem/pull/155
+* Updates for REopt v3.11 by @vtnate in https://github.com/urbanopt/urbanopt-reopt-gem/pull/157
+
+### Other Changes
+
+* expose error messages for certificate verify failed by @kflemin in https://github.com/urbanopt/urbanopt-reopt-gem/pull/151
+* REopt GHP LCCA Analysis by @tanushree04 in https://github.com/urbanopt/urbanopt-reopt-gem/pull/153
+
+**Full Changelog**: https://github.com/urbanopt/urbanopt-reopt-gem/compare/v0.12.0...v1.0.0
+
 ## Version 0.12.0
+
 Date Range: 01/11/24 - 05/06/24
+
 * Use Reopt v3 by @vtnate in https://github.com/urbanopt/urbanopt-reopt-gem/pull/149
 
 **Full Changelog**: https://github.com/urbanopt/urbanopt-reopt-gem/compare/v0.11.0...v0.12.0
 
 ## Version 0.11.0
+
 Date Range: 07/06/23 - 01/11/24
 
 * Use different error messages in different places by @vtnate in https://github.com/urbanopt/urbanopt-reopt-gem/pull/145
@@ -17,6 +34,7 @@ Date Range: 07/06/23 - 01/11/24
 **Full Changelog**: https://github.com/urbanopt/urbanopt-reopt-gem/compare/v0.10.0...v0.11.0
 
 ## Version 0.10.0
+
 Date Range: 12/13/22 - 7/6/23
 
 - Update dependencies for OpenStudio 3.6.1

data/Gemfile

@@ -17,7 +17,9 @@ allow_local = ENV['FAVOR_LOCAL_GEMS']
 # if allow_local && File.exist?('../OpenStudio-extension-gem')
 #   gem 'openstudio-extension', path: '../OpenStudio-extension-gem'
 # elsif allow_local
-#   gem 'openstudio-extension', github: 'NREL/OpenStudio-extension-gem', branch: 'develop'
+# gem 'openstudio-extension', github: 'NREL/OpenStudio-extension-gem', branch: 'develop'
+# else
+# gem 'openstudio-extension', '~> 0.8.1'
 # end
 #
 # if allow_local && File.exist?('../openstudio-common-measures-gem')
@@ -41,14 +43,17 @@ allow_local = ENV['FAVOR_LOCAL_GEMS']
 #   gem 'openstudio-model-articulation', '0.1.0'
 # end
 
-# if allow_local && File.exist?('../urbanopt-scenario-gem')
-#  gem 'urbanopt-scenario', path: '../urbanopt-scenario-gem'
-# elsif allow_local
-# gem 'urbanopt-scenario', github: 'URBANopt/urbanopt-scenario-gem', branch: 'reopt-v3'
-# end
+if allow_local && File.exist?('../urbanopt-scenario-gem')
+  gem 'urbanopt-scenario', path: '../urbanopt-scenario-gem'
+elsif allow_local
+  gem 'urbanopt-scenario', github: 'URBANopt/urbanopt-scenario-gem', branch: 'develop'
+end
 
-# if allow_local && File.exists?('../urbanopt-geojson-gem')
-#   gem 'urbanopt-geojson', path: '../urbanopt-geojson-gem'
-# elsif allow_local
-#   gem 'urbanopt-geojson', github: 'URBANopt/urbanopt-geojson-gem', branch: 'develop'
-# end
+# Temporary! Remove this once reporting-gem is merged/released
+# gem 'urbanopt-reporting', github: 'URBANopt/urbanopt-reporting-gem', branch: 'develop'
+
+if allow_local && File.exist?('../urbanopt-geojson-gem')
+  gem 'urbanopt-geojson', path: '../urbanopt-geojson-gem'
+elsif allow_local
+  gem 'urbanopt-geojson', github: 'URBANopt/urbanopt-geojson-gem', branch: 'develop'
+end

data/LICENSE.md

@@ -1,4 +1,4 @@
-URBANopt (tm), Copyright (c) 2019-2024, Alliance for Sustainable Energy, LLC, and other
+URBANopt (tm), Copyright (c) 2019-2025, Alliance for Sustainable Energy, LLC, and other
 contributors. All rights reserved.
 
 Redistribution and use in source and binary forms, with or without modification,

data/doc_templates/LICENSE.md

@@ -1,4 +1,4 @@
-URBANopt (tm), Copyright (c) 2019-2024, Alliance for Sustainable Energy, LLC, and other
+URBANopt (tm), Copyright (c) 2019-2025, Alliance for Sustainable Energy, LLC, and other
 contributors. All rights reserved.
 
 Redistribution and use in source and binary forms, with or without modification,

data/doc_templates/copyright_erb.txt

@@ -1,6 +1,6 @@
 <%
   # *********************************************************************************
-  # URBANopt (tm), Copyright (c) 2019-2024, Alliance for Sustainable Energy, LLC, and other
+  # URBANopt (tm), Copyright (c) 2019-2025, Alliance for Sustainable Energy, LLC, and other
   # contributors. All rights reserved.
 
   # Redistribution and use in source and binary forms, with or without modification,

data/docs/schemas/reopt-output-schema.md

@@ -1,6 +1,6 @@
 # REopt Lite Outputs Schema
 
-When the gem calls the REopt Lite APUI it recieves the following complete set of results described in the data dictionary below. Only those needed to update a Feature or Scenario Report's distributed_generation attibute set and timeseries CSV are pulled from the reponse and transferred to the Feature or Scenario Report. You may choose to modify the code to include more or less of the full REopt Lite response.
+When the gem calls the REopt Lite APUI it receives the following complete set of results described in the data dictionary below. Only those needed to update a Feature or Scenario Report's distributed_generation attribute set and timeseries CSV are pulled from the response and transferred to the Feature or Scenario Report. You may choose to modify the code to include more or less of the full REopt Lite response.
 
 ## Data Dictionary
 <ReoptOutputSchema />

data/lib/urbanopt/reopt/feature_report_adapter.rb

@@ -41,9 +41,9 @@ module URBANopt # :nodoc:
         else
           @@logger.info('Using default REopt assumptions')
           reopt_inputs = {
-            Settings:{},
+            Settings: {},
             Site: {},
-            Financial:{},
+            Financial: {},
             ElectricTariff: {
               monthly_demand_rates: [0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0],
               monthly_energy_rates: [0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13, 0.13]
@@ -151,14 +151,13 @@ module URBANopt # :nodoc:
       #
       # [*parameters:*]
       #
-      # * +feature_report+ - _URBANopt::Reporting::DefaultReports::FeatureReport_ - FeatureReport to update from a \REopt reponse hash.
-      # * +reopt_output+ - _Hash_ - A reponse hash from the \REopt API to use in overwriting FeatureReport technology sizes, costs and dispatch strategies.
+      # * +feature_report+ - _URBANopt::Reporting::DefaultReports::FeatureReport_ - FeatureReport to update from a \REopt response hash.
+      # * +reopt_output+ - _Hash_ - A response hash from the \REopt API to use in overwriting FeatureReport technology sizes, costs and dispatch strategies.
       # * +timeseries_csv_path+ - _String_ - Optional. The path to a file at which a new timeseries CSV will be written. If not provided a file is created based on the run_uuid of the \REopt optimization task.
       #
       # [*return:*] _URBANopt::Reporting::DefaultReports::FeatureReport_ - Returns an updated FeatureReport.
       ##
       def update_feature_report(feature_report, reopt_output, timeseries_csv_path = nil, resilience_stats = nil)
-
         # Check if the \REopt response is valid
         if reopt_output['status'] != 'optimal'
           @@logger.error("ERROR cannot update Feature Report #{feature_report.name} #{feature_report.id}  - REopt optimization was non-optimal")
@@ -439,7 +438,7 @@ module URBANopt # :nodoc:
           x[$storage_to_grid_col] = $storage_to_grid[i] || 0 if defined?(storage)
           x[$storage_soc_col] = $storage_soc[i] || 0 if defined?(storage)
           x[$generator_total_col] = $generator_total[i] || 0 if defined?(generator)
-          x[$generator_to_battery_col] = $generator_to_battery[i] || 0 if (defined?(generator) && defined?(storage))
+          x[$generator_to_battery_col] = $generator_to_battery[i] || 0 if defined?(generator) && defined?(storage)
           x[$generator_to_load_col] = $generator_to_load[i] || 0 if defined?(generator)
           x[$generator_to_grid_col] = $generator_to_grid[i] || 0 if defined?(generator)
           x[$pv_total_col] = $pv_total[i] || 0
@@ -447,7 +446,7 @@ module URBANopt # :nodoc:
           x[$pv_to_load_col] = $pv_to_load[i] || 0
           x[$pv_to_grid_col] = $pv_to_grid[i] || 0
           x[$wind_total_col] = $wind_total[i] || 0 if defined?(wind)
-          x[$wind_to_battery_col] = $wind_to_battery[i] || 0 if (defined?(wind) && defined?(storage))
+          x[$wind_to_battery_col] = $wind_to_battery[i] || 0 if defined?(wind) && defined?(storage)
           x[$wind_to_load_col] = $wind_to_load[i] || 0 if defined?(wind)
           x[$wind_to_grid_col] = $wind_to_grid[i] || 0 if defined?(wind)
           return x

data/lib/urbanopt/reopt/reopt_ghp_adapter.rb

@@ -0,0 +1,337 @@
+# *********************************************************************************
+# URBANopt (tm), Copyright (c) Alliance for Sustainable Energy, LLC.
+# See also https://github.com/urbanopt/urbanopt-reopt-gem/blob/develop/LICENSE.md
+# *********************************************************************************
+
+module URBANopt # :nodoc:
+  module REopt # :nodoc:
+    class REoptGHPAdapter
+
+      def initialize
+        # initialize @@logger
+        @@logger ||= URBANopt::REopt.reopt_logger
+        # Define class variable
+        @@hours_in_year = 8760
+      end
+
+      def create_reopt_input_building(run_dir, system_parameter_hash, reopt_ghp_assumptions_hash, building_id, modelica_result)
+
+        # Define variables
+        reopt_inputs_building = {}
+        if !reopt_ghp_assumptions_hash.nil?
+          reopt_inputs_building = reopt_ghp_assumptions_hash
+        else
+          @@logger.info('Using default REopt assumptions')
+          # create a dictionary for REopt Inputs
+          reopt_inputs_building = {
+            Site: {},
+            SpaceHeatingLoad: {},
+            DomesticHotWaterLoad: {},
+            ElectricLoad: {},
+            ElectricTariff: {
+              urdb_label: ""
+            },
+            GHP: {},
+            ExistingBoiler: {}
+          }
+        end
+
+        # 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"
+        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[:Site][:latitude] = @latitude
+            reopt_inputs_building[:Site][:longitude] = @longitude
+
+          end
+        end
+
+        reopt_inputs_building[: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
+          total_kbtu = 0.0
+
+          # 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
+              kBtu_value = row['Heating:NaturalGas(kBtu)'].to_f # Convert to float
+              total_kbtu += kBtu_value # Sum kBtu values
+            end
+          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
+          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[:SpaceHeatingLoad][:fuel_loads_mmbtu_per_hour] << mMBtu_value # Append to the array
+              end
+            end
+          end
+
+        else
+          # Calculate space heating load
+          reopt_inputs_building[:SpaceHeatingLoad][:fuel_loads_mmbtu_per_hour] = [0.000001] * @@hours_in_year
+          puts "Existing heating fuel cost was not taken into consideration in result calculations."
+        end
+
+
+        # read_modelica_result
+        modelica_project = File.expand_path(modelica_result)
+        project_name = File.basename(modelica_project)
+        @modelica_csv = File.join(
+          modelica_project,
+          "#{project_name}.Districts.DistrictEnergySystem_results",
+          "#{project_name}.Districts.DistrictEnergySystem_result.csv"
+        )
+
+        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}"
+          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 = []
+          # Ensure the column exists
+          if modelica_data.headers.include?(heating_power)
+            heating_power_values = modelica_data[heating_power]
+          end
+          if modelica_data.headers.include?(cooling_power)
+            cooling_power_values = modelica_data[cooling_power]
+          end
+          if modelica_data.headers.include?(pump_power)
+            pump_power_values = modelica_data[pump_power]
+          end
+          if modelica_data.headers.include?(ets_pump_power)
+            ets_pump_power_values = modelica_data[ets_pump_power]
+          end
+
+          total_electric_load_building = heating_power_values.zip(cooling_power_values, pump_power_values, ets_pump_power_values).map do |elements|
+            # 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)
+            heating_system_capacity_value = modelica_data[heating_system_capacity][0]
+          end
+          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
+
+
+          domestic_hot_water = total_electric_load_building.map do |load|
+            load * 0
+          end
+
+          # This is not used in REopt calculation but required for formatting.
+          reopt_inputs_building[:DomesticHotWaterLoad][:fuel_loads_mmbtu_per_hour] = domestic_hot_water
+
+          # Add GHP Fields
+          reopt_inputs_building[:GHP] = {}
+          # 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)
+            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]
+            end
+          else
+            puts "File not found: #{building_json_path}"
+          end
+
+          # 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
+
+          # Add ghpghx_responses
+          ghpghx_output = {}
+          ghpghx_output[:outputs] = {}
+          ghpghx_output[:inputs] = {}
+          ghpghx_output[:outputs][:heat_pump_configuration] = "WSHP"
+          # This is not used in REopt calculation but required for formatting.
+          ghpghx_output[:outputs][:yearly_ghx_pump_electric_consumption_series_kw] = [0] * @@hours_in_year
+          ghpghx_output[:outputs][:number_of_boreholes] = 0
+          ghpghx_output[:outputs][:length_boreholes_ft] = 0
+
+          ghpghx_output[:outputs][:peak_combined_heatpump_thermal_ton] = peak_combined_heatpump_thermal_ton
+          ghpghx_output[:outputs][:yearly_total_electric_consumption_kwh] = total_electric_load_building.sum
+          ghpghx_output[:outputs][:yearly_total_electric_consumption_series_kw] = total_electric_load_building
+          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
+          # This is not used in REopt calculation but required for formatting.
+          ghpghx_output[:inputs][:cooling_thermal_load_ton] = [0] * @@hours_in_year
+
+          ghpghx_output_all = [ghpghx_output]
+          reopt_inputs_building[:GHP][:ghpghx_responses] = {}
+          reopt_inputs_building[:GHP][:ghpghx_responses] = ghpghx_output_all
+
+        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, "GHP_building_#{building_id}.json")
+        pretty_json = JSON.pretty_generate(reopt_inputs_building)
+        File.write(json_file_path, pretty_json)
+
+      end
+
+      def create_reopt_input_district(run_dir, system_parameter_hash, reopt_ghp_assumptions_hash, ghp_id, modelica_result)
+
+        reopt_inputs_district = {}
+
+        if !reopt_ghp_assumptions_hash.nil?
+          reopt_inputs_district = reopt_ghp_assumptions_hash
+        else
+          @@logger.info('Using default REopt assumptions')
+          # create a dictionary for REopt Inputs
+          reopt_inputs_district = {
+            Site: {},
+            SpaceHeatingLoad: {},
+            DomesticHotWaterLoad: {},
+            ElectricLoad: {},
+            ElectricTariff: {
+              "urdb_label": ""
+            },
+            GHP: {},
+            ExistingBoiler: {}
+          }
+        end
+
+        reopt_inputs_district[:Site] = {}
+        reopt_inputs_district[:Site][:latitude] = @latitude
+        reopt_inputs_district[:Site][:longitude] = @longitude
+        # 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"
+
+        end
+
+        # 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
+        # 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
+
+        # 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[:ExistingBoiler] = {}
+        reopt_inputs_district[:ExistingBoiler][:fuel_cost_per_mmbtu] = 13.5
+
+        # 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][:om_cost_per_sqft_year] = 0
+        reopt_inputs_district[:GHP][:heatpump_capacity_sizing_factor_on_peak_load] = 1.0
+
+        # Add ghpghx outputs
+        ghpghx_output = {}
+        ghpghx_output[:outputs] = {}
+        ghpghx_output[:inputs] = {}
+
+        ghpghx_output[:inputs][:heating_thermal_load_mmbtu_per_hr] = [0]*@@hours_in_year
+        ghpghx_output[:inputs][:cooling_thermal_load_ton] = [0] * @@hours_in_year
+
+
+        # 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
+        end
+
+        if File.exist?(@modelica_csv)
+
+          modelica_data = CSV.read(@modelica_csv, headers: true)
+
+          electrical_power_consumed = modelica_data["electrical_power_consumed"]
+          # Convert watts to kilowatts
+          electrical_power_consumed_kw = electrical_power_consumed.map { |e| e.to_f / 1000 }
+          # if ghp_id.include?('-')
+          #   # Note: For some reason when reading columns, '-' from the column headers are removed, whereas ghp_id has -
+          #   ghp_id_formatted = ghp_id.delete('-')
+          #   ghp_column = "electrical_power_consumed_#{ghp_id_formatted}".to_sym
+
+          # else
+          #   # Note: For some reason when reading columns, '-' from the column headers are removed, whereas ghp_id has -
+          #   ghp_column = "electrical_power_consumed_#{ghp_id}".to_sym
+          # end
+
+          # # Ensure the column exists
+          # unless modelica_data.headers.include?(ghp_column)
+          #   puts "Column #{ghp_column} does not exist in the CSV file."
+          # end
+
+          # # Access values from the column
+          # 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
+        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
+        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
+
+        ghpghx_output_all = [ghpghx_output, ghpghx_output]
+        reopt_inputs_district[:GHP][:ghpghx_responses] = ghpghx_output_all
+
+        #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, "GHX_#{ghp_id}.json")
+        pretty_json = JSON.pretty_generate(reopt_inputs_district)
+        File.write(json_file_path, pretty_json)
+      end
+
+    end
+  end
+end