set sam ppa price from geophires electricity price

Author: softwareengineerprogrammer

src/geophires_x/EconomicsSam.py

@@ -23,40 +23,8 @@
 import geophires_x.Model as Model
 
 
-def get_single_owner_parameters(model: Model) -> dict[str, Any]:
-    econ = model.economics
-
-    ret: dict[str, Any] = {}
-
-    itc = econ.CCap.value * econ.RITC.value
-    total_capex_musd = (
-        econ.CCap.value - itc
-    )
-    ret['total_installed_cost'] = total_capex_musd * 1e6
-
-    opex_musd = econ.Coam.value
-    ret['om_fixed'] = [opex_musd * 1e6]
-
-    average_net_generation_MW = get_average_net_generation_MW(model)
-    ret['system_capacity'] = average_net_generation_MW * 1e3
-
-    geophires_ctr_tenths = Decimal(econ.CTR.value)
-    fed_rate_tenths = geophires_ctr_tenths * (Decimal(0.7))
-    state_rate_tenths = geophires_ctr_tenths - fed_rate_tenths
-    ret['federal_tax_rate'] = [float(fed_rate_tenths * Decimal(100))]
-    ret['state_tax_rate'] = [float(state_rate_tenths * Decimal(100))]
-
-    geophires_itc_tenths = Decimal(econ.RITC.value)
-    ret['itc_fed_percent'] = [float(geophires_itc_tenths * Decimal(100))]
-
-    geophires_ptr_tenths = Decimal(econ.PTR.value)
-    ret['property_tax_rate'] = float(geophires_ptr_tenths * Decimal(100))
-
-    return ret
-
-
 def calculate_sam_economics(
-    model: Model
+        model: Model
 ) -> dict[str, dict[str, Any]]:
     custom_gen = CustomGeneration.new()
     grid = Grid.from_existing(custom_gen)
@@ -72,12 +40,11 @@ def calculate_sam_economics(
     for module_file, module in zip(file_names, modules):
         with open(Path(project_dir, f'{module_file}.json'), encoding='utf-8') as file:
             data = json.load(file)
-            # loop through each key-value pair
             for k, v in data.items():
                 if k != 'number_inputs':
                     module.value(k, v)
 
-    for k, v in get_single_owner_parameters(model).items():
+    for k, v in _get_single_owner_parameters(model).items():
         single_owner.value(k, v)
 
     for module in modules:
@@ -97,17 +64,51 @@ def calculate_sam_economics(
     ret = {}
     for e in display_data:
         field_display = e[0] + ':' + ' ' * (max_field_name_len - len(e[0]) - 1)
-        print(f'{field_display}\t{sig_figs(e[1], 5)} {e[2]}')
-        ret[e[0]] = {'value': sig_figs(e[1], 5), 'unit': e[2]}
+        # print(f'{field_display}\t{sig_figs(e[1], 5)} {e[2]}')
+        ret[e[0]] = {'value': _sig_figs(e[1], 5), 'unit': e[2]}
+
+    return ret
+
+
+def _get_single_owner_parameters(model: Model) -> dict[str, Any]:
+    econ = model.economics
+
+    ret: dict[str, Any] = {}
+
+    itc = econ.CCap.value * econ.RITC.value
+    total_capex_musd = (
+            econ.CCap.value - itc
+    )
+    ret['total_installed_cost'] = total_capex_musd * 1e6
+
+    opex_musd = econ.Coam.value
+    ret['om_fixed'] = [opex_musd * 1e6]
+
+    average_net_generation_MW = _get_average_net_generation_MW(model)
+    ret['system_capacity'] = average_net_generation_MW * 1e3
+
+    geophires_ctr_tenths = Decimal(econ.CTR.value)
+    fed_rate_tenths = geophires_ctr_tenths * (Decimal(0.7))
+    state_rate_tenths = geophires_ctr_tenths - fed_rate_tenths
+    ret['federal_tax_rate'] = [float(fed_rate_tenths * Decimal(100))]
+    ret['state_tax_rate'] = [float(state_rate_tenths * Decimal(100))]
+
+    geophires_itc_tenths = Decimal(econ.RITC.value)
+    ret['itc_fed_percent'] = [float(geophires_itc_tenths * Decimal(100))]
+
+    geophires_ptr_tenths = Decimal(econ.PTR.value)
+    ret['property_tax_rate'] = float(geophires_ptr_tenths * Decimal(100))
+
+    ret['ppa_price_input'] = [econ.ElecStartPrice.value]
 
     return ret
 
 
-def get_average_net_generation_MW(model: Model) -> float:
+def _get_average_net_generation_MW(model: Model) -> float:
     return np.average(model.surfaceplant.NetElectricityProduced.value)
 
 
-def sig_figs(val: float, num_sig_figs: int) -> float:
+def _sig_figs(val: float, num_sig_figs: int) -> float:
     if val is None:
         return None
 

tests/test_geophires_x.py

@@ -915,13 +915,27 @@ def test_negative_electricity_production_raises_error(self):
         self.assertIn('Electricity production calculated as negative', str(e.exception))
 
     def test_economic_model_single_owner_ppa_sam(self):
-        result = GeophiresXClient().get_geophires_result(
-            GeophiresInputParameters(
-                from_file_path=self._get_test_file_path('geophires_x_tests/generic-egs-case.txt'),
-                params={
-                    'Economic Model': 5,
-                },
+        def _get_result(_params) -> GeophiresXResult:
+            return GeophiresXClient().get_geophires_result(
+                GeophiresInputParameters(
+                    from_file_path=self._get_test_file_path('geophires_x_tests/generic-egs-case.txt'),
+                    params={
+                        'Economic Model': 5,
+                    }
+                    | _params,
+                )
             )
-        )
 
-        self.assertGreater(result.result['SUMMARY OF RESULTS']['Electricity breakeven price']['value'], 6)
+        def _lcoe(r: GeophiresXResult) -> float:
+            return r.result['SUMMARY OF RESULTS']['Electricity breakeven price']['value']
+
+        def _npv(r: GeophiresXResult) -> float:
+            return r.result['ECONOMIC PARAMETERS']['Project NPV']['value']
+
+        base_result = _get_result({})
+        base_lcoe = _lcoe(base_result)
+        self.assertGreater(base_lcoe, 6)
+
+        npvs = [_npv(_get_result({'Starting Electricity Sale Price': x / 100.0})) for x in range(1, 20, 2)]
+        for i in range(len(npvs) - 1):
+            self.assertLess(npvs[i], npvs[i + 1])