def pvap_water(TC):
    """Return vapor pressure of water (Pa) for temperature in C"""

    # Antoine equation
    a, b, c = 8.07131, 1730.63, 233.426
    p_mmHg = 10**(a - (b/(c+TC)))
    p_Pa = p_mmHg * 133.3

    return p_Pa