Output segment thickness values as given in order to avoid dropping precision

Author: softwareengineerprogrammer

src/geophires_x/Outputs.py

@@ -1621,7 +1621,7 @@ def PrintOutputs(self, model: Model):
                 else:
                     for i in range(1, model.reserv.numseg.value):
                         f.write(f'      Segment {str(i):s}   Geothermal gradient:                    {model.reserv.gradient.value[i-1]:10.4g} ' + model.reserv.gradient.CurrentUnits.value +NL)
-                        f.write(f'      Segment {str(i):s}   Thickness:                         {model.reserv.layerthickness.value[i-1]:10.0f} ' + model.reserv.layerthickness.CurrentUnits.value + NL)
+                        f.write(f'      Segment {str(i):s}   Thickness:                         {model.reserv.layerthickness.value[i-1]} {model.reserv.layerthickness.CurrentUnits.value}\n')
                     f.write(f'      Segment {str(i+1):s}   Geothermal gradient:                    {model.reserv.gradient.value[i]:10.4g} ' + model.reserv.gradient.CurrentUnits.value + NL)
                 if model.economics.DoCarbonCalculations.value:
                     f.write(f'      Total Avoided Carbon Emissions:                       {model.economics.CarbonThatWouldHaveBeenProducedTotal.value:10.2f} '
@@ -1702,7 +1702,7 @@ def PrintOutputs(self, model: Model):
                 else:
                     for i in range(1, model.reserv.numseg.value):
                         f.write(f'      Segment {str(i):s}   Geothermal gradient:                    {model.reserv.gradient.value[i-1]:10.4g} ' + model.reserv.gradient.CurrentUnits.value +NL)
-                        f.write(f'      Segment {str(i):s}   Thickness:                         {model.reserv.layerthickness.value[i-1]:10.0f} ' + model.reserv.layerthickness.CurrentUnits.value + NL)
+                        f.write(f'      Segment {str(i):s}   Thickness:                         {model.reserv.layerthickness.value[i-1]} {model.reserv.layerthickness.CurrentUnits.value}\n')
                     f.write(f'      Segment {str(i+1):s}   Geothermal gradient:                    {model.reserv.gradient.value[i]:10.4g} ' + model.reserv.gradient.CurrentUnits.value + NL)
 
                 f.write(NL)

tests/test_geophires_x.py

@@ -819,3 +819,25 @@ def test_drilling_cost_curves(self):
                     'value'
                 ]
                 self.assertAlmostEqual(indirect_cost_factor * expected_cost_musd, cost_per_well_val, delta=0.1)
+
+    def test_segment_thickness_output(self):
+        thickness_1 = 0.793
+        thickness_2 = 1.646
+        result = GeophiresXClient().get_geophires_result(
+            GeophiresInputParameters(
+                from_file_path=self._get_test_file_path('geophires_x_tests/generic-egs-case.txt'),
+                params={
+                    'Number of Segments': 3,
+                    'Gradient 1': 42.69972,
+                    'Gradient 2': 51.66667,
+                    'Thickness 1': thickness_1,
+                    'Gradient 3': 46.9697,
+                    'Thickness 2': thickness_2,
+                },
+            )
+        )
+
+        self.assertEqual(thickness_1, result.result['SUMMARY OF RESULTS']['Segment 1   Thickness']['value'])
+        self.assertAlmostEqual(
+            thickness_2, result.result['SUMMARY OF RESULTS']['Segment 2   Thickness']['value'], places=3
+        )