Fix testing disabling saturarion properties if convergence problem

jjgomera · Apr 20, 2021 · f252e01 · f252e01
1 parent 0a687f6
commit f252e01
Showing 1 changed file with 8 additions and 2 deletions.
diff --git a/iapws/humidAir.py b/iapws/humidAir.py
@@ -686,8 +686,14 @@ def f(T):
 
         # Saturation related properties
         A_sat = self._eq(self.T, self.P)
-        self.xa_sat = A_sat*MW/Ma/(1-A_sat*(1-MW/Ma))
-        self.RH = (1-self.xa)/(1-self.xa_sat)
+        if A_sat:
+            self.xa_sat = A_sat*MW/Ma/(1-A_sat*(1-MW/Ma))
+            self.RH = (1-self.xa)/(1-self.xa_sat)
+        else:
+            self.xa_sat = None
+            self.RH = None
+            self.msg = "Saturation state don't converge"
+            self.status = 3
 
     def derivative(self, z, x, y):
         """