def estimated_inversion_strength(T0,T700):
<EOL>T850 = (T0+T700)/<NUM_LIT>;<EOL>LCL = lifting_condensation_level(T0, <NUM_LIT>)<EOL>LTS = potential_temperature(T700, <NUM_LIT>) - T0<EOL>Gammam = (g/cp*(<NUM_LIT:1.0> - (<NUM_LIT:1.0> + Lhvap*qsat(T850,<NUM_LIT>) / Rd / T850) /<EOL>(<NUM_LIT:1.0> + Lhvap**<NUM_LIT:2> * qsat(T850,<NUM_LIT>)/ cp/Rv/T850**<NUM_LIT:2>)))<EOL>z700 = (Rd*T0/g)*log(<NUM_LIT>/<NUM_LIT>)<EOL>return LTS - Gammam*(z700 - LCL)<EOL>
Compute the Estimated Inversion Strength or EIS,
    following Wood and Bretherton (2006, J. Climate)

    Inputs: T0 is surface temp in Kelvin
           T700 is air temperature at 700 hPa in Kelvin

    Output: EIS in Kelvin

    EIS is a normalized measure of lower tropospheric stability acccounting for
    temperature-dependence of the moist adiabat.