def clausius_clapeyron(T):
Tcel = T - tempCtoK<EOL>es = <NUM_LIT> * exp(<NUM_LIT>*Tcel/(Tcel+<NUM_LIT>))<EOL>return es<EOL>
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).