def pmelt_T_iceV(T):
    """
    EQ 3 / Melting pressure of ice V
    """
    T_star = 256.164
    p_star = 350.1
    theta = T / T_star
    pi_melt = 1 - 1.18721 * (1.0 - theta ** 8)
    return pi_melt * p_star