support Inflation Rate During Construction

Author: softwareengineerprogrammer

docs/SAM-Economic-Models.md

@@ -28,7 +28,7 @@ The following table describes how GEOPHIRES parameters are transformed into SAM
 | `Max Net Electricity Production`                                                                                                                | Generation Profile                                     | `Nameplate capacity`                                                             | `Singleowner`                     | `system_capacity`                                            | .. N/A                                                                                                                                                                                                                                                                                                                    |
 | `Utilizaton Factor`                                                                                                                             | Generation Profile                                     | `Nominal capacity factor`                                                        | `Singleowner`                     | `user_capacity_factor`                                       | .. N/A                                                                                                                                                                                                                                                                                                                    |
 | `Net Electricity Production`                                                                                                                    | AC Degradation                                         | `Annual AC degradation rate` schedule                                            | `Utilityrate5`                    | `degradation`                                                | Percentage difference of each year's `Net Electricity Production` from `Max Net Electricity Production` is input as SAM as the degradation rate schedule in order to match SAM's generation profile to GEOPHIRES                                                                                                          |
-| `Total Capital Cost` + `Investment Tax Credit Value`                                                                                            | Installation Costs                                     | `Total Installed Cost`                                                           | `Singleowner`                     | `total_installed_cost`                                       | .. N/A                                                                                                                                                                                                                                                                                                                    |
+| `Total Capital Cost` + `Investment Tax Credit Value` + (`Inflation Rate During Construction` × `Total Capital Cost`)                            | Installation Costs                                     | `Total Installed Cost`                                                           | `Singleowner`                     | `total_installed_cost`                                       | ITC, if present, is added to GEOPHIRES total capital cost since SAM handles ITC credits in cash flow analysis. Inflation during construction is treated as an indirect EPC capital cost percentage.                                                                                                                       |
 | `Total O&M Cost`, `Inflation Rate`                                                                                                              | Operating Costs                                        | `Fixed operating cost`, `Escalation rate` set to `Inflation Rate` × -1           | `Singleowner`                     | `om_fixed`, `om_fixed_escal`                                 | .. N/A                                                                                                                                                                                                                                                                                                                    |
 | `Plant Lifetime`                                                                                                                                | Financial Parameters → Analysis Parameters             | `Analysis period`                                                                | `CustomGeneration`, `Singleowner` | `CustomGeneration.analysis_period`, `Singleowner.term_tenor` | .. N/A                                                                                                                                                                                                                                                                                                                    |
 | `Inflation Rate`                                                                                                                                | Financial Parameters → Analysis Parameters             | `Inflation rate`                                                                 | `Utilityrate5`                    | `inflation_rate`                                             | .. N/A                                                                                                                                                                                                                                                                                                                    |

src/geophires_x/EconomicsSam.py

@@ -181,9 +181,15 @@ def _get_single_owner_parameters(model: Model) -> dict[str, Any]:
 
     ret['analysis_period'] = model.surfaceplant.plant_lifetime.value
 
-    itc = econ.RITCValue.value
-    total_capex_musd = econ.CCap.value + itc
-    ret['total_installed_cost'] = total_capex_musd * 1e6
+    itc = econ.RITCValue.quantity()
+    total_capex = econ.CCap.quantity() + itc
+    # ret['total_installed_cost'] = total_capex_musd * 1e6
+
+    # 'Inflation Rate During Construction'
+    construction_additional_cost = econ.inflrateconstruction.value * total_capex
+
+    ret['total_installed_cost'] = (total_capex + construction_additional_cost).to('USD').magnitude
+
     # TODO break out indirect costs (instead of lumping all into direct cost)
 
     opex_musd = econ.Coam.value
@@ -193,8 +199,6 @@ def _get_single_owner_parameters(model: Model) -> dict[str, Any]:
 
     # TODO construction years
 
-    # TODO 'Inflation Rate During Construction'
-
     # Note generation profile is generated relative to the max in _get_utility_rate_parameters
     ret['system_capacity'] = _get_max_net_generation_MW(model) * 1e3
 

tests/geophires_x_tests/test_economics_sam.py

@@ -175,10 +175,16 @@ def get_single_value(name: str) -> list[float]:
     @staticmethod
     def _get_cash_flow_row(cash_flow, name):
 
+        def r_0(r):
+            if r is not None and len(r) > 0:
+                return r[0]
+
+            return None
+
         return next(
             [GeophiresXResult._get_sam_cash_flow_profile_entry_display_to_entry_val(ed) for ed in r]
             for r in cash_flow
-            if r[0] == name
+            if r_0(r) == name
         )[1:]
 
     def test_only_electricity_end_use_supported(self):
@@ -258,6 +264,20 @@ def get_row(name: str):
 
         assert_incentives({'One-time Grants Etc': 100, 'Other Incentives': 0.1}, 100_100_000)
 
+    def test_inflation_rate_during_construction(self):
+        infl_rate = 0.05
+        r_no_infl = self._get_result({'Inflation Rate During Construction': 0})
+        r_infl = self._get_result({'Inflation Rate During Construction': infl_rate})
+
+        cash_flow_no_infl = r_no_infl.result['SAM CASH FLOW PROFILE']
+        cash_flow_infl = r_infl.result['SAM CASH FLOW PROFILE']
+
+        tic = 'Total installed cost ($)'
+        tic_no_infl = EconomicsSamTestCase._get_cash_flow_row(cash_flow_no_infl, tic)[0]
+        tic_infl = EconomicsSamTestCase._get_cash_flow_row(cash_flow_infl, tic)[0]
+
+        self.assertAlmostEqual(tic_no_infl * (1 + infl_rate), tic_infl, places=0)
+
     def test_ptc(self):
         def assert_ptc(params, expected_ptc_usd_per_kWh):
             m: Model = EconomicsSamTestCase._new_model(self._egs_test_file_path(), additional_params=params)