We use the standard mixing-length theory of Böhm-Vitense ([1]}) to deal with
the convection. We set the mixing-length \(l\)  to be equal to one local pressure scale
height, i.e.,
\(l=H_{\rm P}=-\cfrac{dr}{d\ln P}=\cfrac{P}{\rho g}~.\) 
