Fix incorrect method of setting royalties output param value

Author: softwareengineerprogrammer

src/geophires_x/Economics.py

@@ -2507,7 +2507,9 @@ def Calculate(self, model: Model) -> None:
             self.CCap.value = (self.sam_economics_calculations.capex.quantity()
                                .to(self.CCap.CurrentUnits.value).magnitude)
 
-            average_annual_royalties = np.average(self.sam_economics_calculations.royalties_opex[1:])  # ignore Year 0
+            average_annual_royalties = np.average(
+                self.sam_economics_calculations.royalties_opex.value[1:]  # ignore pre-revenue year(s) (Year 0)
+            )
             if average_annual_royalties > 0:
                 self.royalties_opex.value = average_annual_royalties
                 self.Coam.value += self.royalties_opex.quantity().to(self.Coam.CurrentUnits.value).magnitude

src/geophires_x/EconomicsSam.py

@@ -176,7 +176,7 @@ def sf(_v: float, num_sig_figs: int = 5) -> float:
     royalty_rate = model.economics.royalty_rate.quantity().to('dimensionless').magnitude
     ppa_revenue_row = _cash_flow_profile_row(cash_flow, 'PPA revenue ($)')
     royalties_unit = sam_economics.royalties_opex.CurrentUnits.value.replace('/yr', '')
-    sam_economics.royalties_opex = [
+    sam_economics.royalties_opex.value = [
         quantity(x * royalty_rate, 'USD').to(royalties_unit).magnitude for x in ppa_revenue_row
     ]
 

tests/geophires_x_tests/test_economics_sam.py

@@ -655,7 +655,7 @@ def get_row(name: str):
         expected_royalties_USD = [x * royalty_rate for x in ppa_revenue_row_USD]
         expected_royalties_MUSD = [x * 1e-6 for x in expected_royalties_USD]
 
-        self.assertListEqual(expected_royalties_MUSD, sam_econ.royalties_opex)
+        self.assertListEqual(expected_royalties_MUSD, sam_econ.royalties_opex.value)
 
         om_prod_based_expense_row = get_row('O&M production-based expense ($)')
         self.assertListAlmostEqual(expected_royalties_USD, om_prod_based_expense_row, places=0)