Try using np.errstate = warn

softwareengineerprogrammer · Oct 4, 2023 · c10525e · c10525e
1 parent 3e4464b
commit c10525e
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/geophires_x/AGSOutputs.py b/src/geophires_x/AGSOutputs.py
@@ -52,7 +52,11 @@ def PrintOutputs(self, model: Model):
                     f = scipy.interpolate.interp1d(np.arange(0, len(model.surfaceplant.NetElectricityProduced.value)),
                                                    model.surfaceplant.NetElectricityProduced.value, fill_value="extrapolate")
                     aranged = np.arange(0, len(model.wellbores.ProducedTemperature.value), 1.0)
-                    model.surfaceplant.NetElectricityProduced.value = f(aranged)
+
+                    # FIXME WIP (warn errstate may or may not be appropriate)
+                    with np.errstate(divide='warn'):
+                        model.surfaceplant.NetElectricityProduced.value = f(aranged)
+
                 if len(model.surfaceplant.FirstLawEfficiency.value) != len(model.wellbores.ProducedTemperature.value):
                     f = scipy.interpolate.interp1d(np.arange(0, len(model.surfaceplant.FirstLawEfficiency.value)),
                                                    model.surfaceplant.FirstLawEfficiency.value, fill_value="extrapolate")