output interest rate (tangentially related to WIP to pass it to SAM)

Author: softwareengineerprogrammer

src/geophires_x/EconomicsSam.py

@@ -138,7 +138,8 @@ def _get_single_owner_parameters(model: Model) -> dict[str, Any]:
     # Debt/equity ratio ('Fraction of Investment in Bonds' parameter)
     ret['debt_percent'] = econ.FIB.quantity().to(convertible_unit('%')).magnitude
 
-    # TODO interest rate
+    # TODO/WIP interest rate
+    # ret['real_discount_rate']
 
     # 'Property Tax Rate'
     ret['property_tax_rate'] = econ.PTR.quantity().to(convertible_unit('%')).magnitude

src/geophires_x/Outputs.py

@@ -250,6 +250,13 @@ def PrintOutputs(self, model: Model):
 
                 elif model.economics.econmodel.value in (EconomicModel.BICYCLE, EconomicModel.SAM_SINGLE_OWNER_PPA):
                     f.write(f'      Economic Model = {model.economics.econmodel.value.value}\n')
+
+                if model.economics.econmodel.value in EconomicModel.SAM_SINGLE_OWNER_PPA:
+                    # TODO disambiguate interest rate for all economic models - see
+                    #  535c02d4adbeeeca553b61e9b996fccf00016529
+                    irfl = Outputs._field_label(econ.interest_rate.Name, 49)
+                    f.write(f'      {irfl}{econ.interest_rate.value:10.2f} {econ.interest_rate.CurrentUnits.value}\n')
+
                 f.write(f'      Accrued financing during construction:            {model.economics.inflrateconstruction.value*100:10.2f} {model.economics.inflrateconstruction.CurrentUnits.value}\n')
                 f.write(f'      Project lifetime:                              {model.surfaceplant.plant_lifetime.value:10.0f} {model.surfaceplant.plant_lifetime.CurrentUnits.value}\n')
                 f.write(f'      Capacity factor:                                 {model.surfaceplant.utilization_factor.value * 100:10.1f} %\n')

tests/geophires_x_tests/generic-egs-case-2.txt renamed to tests/geophires_x_tests/generic-egs-case-2_sam-single-owner-ppa.txt

@@ -25,7 +25,7 @@
 class EconomicsSamTestCase(BaseTestCase):
 
     def _egs_test_file_path(self) -> str:
-        return self._get_test_file_path('generic-egs-case-2.txt')
+        return self._get_test_file_path('generic-egs-case-2_sam-single-owner-ppa.txt')
 
     def _get_result(self, _params) -> GeophiresXResult:
         return GeophiresXClient().get_geophires_result(
@@ -46,6 +46,10 @@ def _npv(r: GeophiresXResult) -> float:
         base_lcoe = _lcoe(base_result)
         self.assertGreater(base_lcoe, 7)
 
+        ir = base_result.result['ECONOMIC PARAMETERS']['Interest Rate']
+        self.assertEqual(ir['value'], 7.0)
+        self.assertEqual(ir['unit'], '%')
+
         npvs = [_npv(self._get_result({'Starting Electricity Sale Price': x / 100.0})) for x in range(1, 20, 4)]
         for i in range(len(npvs) - 1):
             self.assertLess(npvs[i], npvs[i + 1])