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), )