def p_sat(T=293.15):
    """
    Water vapor saturation pressure [Pa]
    
    input : temperature T [K]
    """
    T_celsius = T - 273.15
    log_p_sat = 2.7858 + 7.5 * T_celsius / (237.3 + T_celsius)
    return 10 ** log_p_sat
