def qsat(T,p):
es = clausius_clapeyron(T)<EOL>q = eps * es / (p - (<NUM_LIT:1> - eps) * es )<EOL>return q<EOL>
Compute saturation specific humidity as function of temperature and pressure.

    Input:  T is temperature in Kelvin
            p is pressure in hPa or mb
    Output: saturation specific humidity (dimensionless).