Parse S_DAC_GT PROFILE

Author: softwareengineerprogrammer

src/geophires_x_client/geophires_x_result.py

@@ -371,6 +371,10 @@ def __init__(self, output_file_path, logger_name=None):
         if ccus_profile is not None:
             self.result['CCUS PROFILE'] = ccus_profile
 
+        sdacgt_profile = self._get_sdacgt_profile()
+        if sdacgt_profile is not None:
+            self.result['S_DAC_GT PROFILE'] = sdacgt_profile
+
         self.result['metadata'] = {'output_file_path': self.output_file_path}
         for metadata_field in GeophiresXResult._METADATA_FIELDS:
             self.result['metadata'][metadata_field] = self._get_equal_sign_delimited_field(metadata_field)
@@ -595,6 +599,27 @@ def extract_table_header(lines: list) -> list:
             self._logger.debug(f'Failed to get extended economic profile: {e}')
             return None
 
+    def _get_sdacgt_profile(self):
+        def extract_table_header(lines: list) -> list:
+            # Tried various regexy approaches to extract this programmatically but landed on hard-coding.
+            return [
+                'Year Since Start',
+                'Carbon Captured (tonne/yr)',
+                'Cumm. Carbon Captured (tonne)',
+                'S_DAC_GT Annual Cost (USD/yr)',
+                'S_DAC_GT Cumm. Cash Flow (USD)',
+                'Cumm. Cost Per Tonne (USD/tonne)',
+            ]
+
+        try:
+            lines = self._get_profile_lines('S_DAC_GT PROFILE')
+            profile = [extract_table_header(lines)]
+            profile.extend(self._extract_addons_style_table_data(lines))
+            return profile
+        except BaseException as e:
+            self._logger.debug(f'Failed to get S-DAC-GT profile: {e}')
+            return None
+
     def _get_ccus_profile(self):
         """
         FIXME TODO - transform from revenue & cashflow if present (CCUS profile replaced by revenue & cashflow

tests/test_geophires_x_client.py

@@ -532,3 +532,29 @@ def test_parse_annualized_capital_costs(self):
     def test_parse_number_with_commas(self):
         result = GeophiresXResult(self._get_test_file_path('examples/S-DAC-GT.out'))
         self.assertEqual(result.result['S_DAC_GT ECONOMICS']['Total Cost of Capture']['value'], 499_311_405.59)
+
+    def test_parse_sdacgt_profile(self):
+        result = GeophiresXResult(self._get_test_file_path('examples/S-DAC-GT.out'))
+        sdacgt_profile = result.result['S_DAC_GT PROFILE']
+        self.assertIsNotNone(sdacgt_profile)
+        self.assertEqual(
+            sdacgt_profile[0],
+            [
+                'Year Since Start',
+                'Carbon Captured (tonne/yr)',
+                'Cumm. Carbon Captured (tonne)',
+                'S_DAC_GT Annual Cost (USD/yr)',
+                'S_DAC_GT Cumm. Cash Flow (USD)',
+                'Cumm. Cost Per Tonne (USD/tonne)',
+            ],
+        )
+
+        # Values below need to be synchronized if S-DAC-GT example output values change.
+        self.assertEqual(sdacgt_profile[1], [1, 78330.8, 78330.8, 17411627.98, 17411627.98, 222.28])
+
+        self.assertEqual(
+            sdacgt_profile[15],
+            [15, 76263.89, 1167207.48, 16952186.81, 259450710.33, 222.28],
+        )
+
+        self.assertEqual(sdacgt_profile[30], [30, 61974.61, 2246284.1, 13775920.11, 499311405.59, 222.28])