docs: def clausius_clapeyron(T):
      """Compute saturation vapor pressure as function of temperature T.
  
      Input: T is temperature in Kelvin
      Output: saturation vapor pressure in mb or hPa
  
      Formula from Rogers and Yau "A Short Course in Cloud Physics" (Pergammon Press), p. 16
      claimed to be accurate to within 0.1% between -30degC and 35 degC
      Based on the paper by Bolton (1980, Monthly Weather Review).
  
      """
      Tcel = T - tempCtoK
      es = 6.112 * exp(17.67*Tcel/(Tcel+243.5))
      return es