Incompressible entropy

[jowr]

I must have broken the entropy function for incompressible fluids.

from __future__ import print_function, division
import numpy as np
from CoolProp.CoolProp import PropsSI
fluids = ['INCOMP::TX22','INCOMP::Water']
p = np.array([13961853.79075649, 15362882.97044676, 2438758.77856947])
T = np.array([     413.20095544,      408.06711477,     457.06607800])
for fluid in fluids:
    print(fluid)
    for _ in range(3):
        #print(PropsSI('S','T',T,'P',p,fluid))
        print("[",end="")
        for i in range(len(p)):
            print(PropsSI('S','T',T[i],'P',p[i],fluid),end=" ")
        print("]")
    print("")

INCOMP::TX22
[-4797.62947591 -15856.3798244 -37644.8817001 ]
[-82006.4531606 -170274.027194 -346480.176439 ]
[-699677.042638 -1405615.20615 -2817162.53435 ]

INCOMP::Water
[-22111.2063579 -69264.7382502 -162975.733605 ]
[-351813.118495 -728668.562524 -1481783.38215 ]
[-2989428.41559 -6003899.15672 -12032244.5705 ]

[jowr]

Looks like it works:

from __future__ import print_function, division
import numpy as np
from CoolProp.CoolProp import PropsSI
fluids = ['HEOS::Water','INCOMP::Water']
p = np.array([13961853.79075649, 15362882.97044676, 2438758.77856947])
T = np.array([     413.20095544,      408.06711477,     457.06607800])
for fluid in fluids:
    print(fluid)
    for out in ['H','S']:
        for _ in range(3):
            #print(PropsSI('S','T',T,'P',p,fluid))
            print("[",end="")
            for i in range(len(p)):
                print(PropsSI(out,'T',T[i],'P',p[i],fluid),end=" ")
            print("]")
    print("")

HEOS::Water
[598278.521655 577459.606388 781035.062803 ]
[598278.521655 577459.606388 781035.062803 ]
[598278.521655 577459.606388 781035.062803 ]
[1725.88039873 1671.53476145 2175.19280846 ]
[1725.88039873 1671.53476145 2175.19280846 ]
[1725.88039873 1671.53476145 2175.19280846 ]

INCOMP::Water
[517406.179903 496937.605375 695615.549984 ]
[517406.179903 496937.605375 695615.549984 ]
[517406.179903 496937.605375 695615.549984 ]
[24989.0668045 24935.8080746 25425.3590451 ]
[24989.0668045 24935.8080746 25425.3590451 ]
[24989.0668045 24935.8080746 25425.3590451 ]