docs: def qsat(T,p):
      """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).
  
      """
      es = clausius_clapeyron(T)
      q = eps * es / (p - (1 - eps) * es )
      return q