def printheatcap():
S = array([0., 30., 35., 40.])
entry = " %6.1f"
print
print " SPECIFIC HEAT SEAWATER Cp [ J/(Kg K) ]"
print "PRESSURE TEMPERATURE",
print deg+"C IPTS-68 SALINTY:", "%2.0f" % S[0]
print "DECIBARS %8d %11d %11d %11d %11d" % tuple(T)
for p in P:
print "%6.0f" % p, 5*entry % tuple(heatcap(S[0],T,p))
for s in S[1:]:
print 68*"-" + "SALINTIY:", "%2.0f" % s
for p in P:
print "%6.0f" % p, 5*entry % tuple(heatcap(s,T,p))
def values(S, T, P):
"""Compute a row in the table"""
return (
P / 10,
S,
T,
dens(S, T, P) - 1000,
drhods(S, T, P),
10 ** 7 * alpha(S, T, P),
heatcap(S, T, P),
1000 * temppot0(S, T, P),
soundvel(S, T, P),
)
