#Metview Macro

function specific_humidity_from_vapour_pressure(e: number, p: number)

    EPSILON = 0.621981
    v = p + (EPSILON - 1) * e
    if p - e < 1E-4 then
        return nil
    end if
    return EPSILON * e / v

end specific_humidity_from_vapour_pressure