Differentiating the
semi-grey irradiated atmosphere profile gives
4T3dτdT=a[Tint4+fTirr4(1−γ2)e−γcτ].
Hydrostatic
balance gives
dτ/dP=κIR/g when the thermal
mean opacity is locally constant. Hence
dPdT=4gT3aκIR[Tint4+fTirr4(1−γ2)e−γcτ].
For variable opacity, replace
κIR by its local value and integrate
dτ=κIR(P,T)dP/g.