thermo.py 文件源码

def svpice( t) :
    '''
    Returns saturation vapor pressure over ice, in hPa, given temperature in K.
    The Goff-Gratch equation (Smithsonian Met. Tables,  5th ed., pp. 350, 1984)
    '''
    a = 273.16 / t
    exponent = -9.09718 * (a - 1.) - 3.56654 * log10(a) + 0.876793 * (1. - 1./a) + log10(6.1071)

    return 10.0**exponent