Our calculation

Estimate pseudo adiabatic lapse rate 
By equilibrium cloud condensed model 
But with a method different from that of Atreya and Romani (1985)
      --> A new calculation method: 
          Minimize Gibbs free energy for a given 
          arbitrary elemental abundance

Estimate static stability 
By assuming that it is given as the difference between moist adiabat and dry adiabat