#Metview Macro

function saturation_specific_humidity(t, p)
    return saturation_specific_humidity(t, p, "water")
end saturation_specific_humidity

function saturation_specific_humidity(t, p, phase)
    e = saturation_vapour_pressure(t, phase)
    return specific_humidity_from_vapour_pressure(e, p)
end saturation_specific_humidity