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