We read every piece of feedback, and take your input very seriously.
To see all available qualifiers, see our documentation.
1 parent 71f2285 commit 5f4cad6Copy full SHA for 5f4cad6
1 file changed
t2thermo.py
@@ -50,6 +50,7 @@
50
from collections.abc import Iterable # Python >= 3.3
51
except ImportError:
52
from collections import Iterable
53
+import numpy as np
54
55
def cowat(t, p, bounds = False):
56
"""Density d and internal energy u of liquid water as a function of
@@ -273,7 +274,7 @@ def sat(t, bounds = False):
273
274
SC = SC * X1 + a[2]
275
SC = SC * X1 + a[1]
276
SC = SC * X1
- PC = exp(SC / (TC * (1.0 + a[6] * X1 + a[7] * X2)) - \
277
+ PC = np.exp(SC / (TC * (1.0 + a[6] * X1 + a[7] * X2)) - \
278
X1 / (a[8] * X2 + a[9]))
279
return PC * 2.212e7
280
else: return None
0 commit comments