Adding to energy eq.Chtmultiregion
 

I want to model a structure in the weather temperature inversion,kind of atmospheric boundary layer under thermal stratification, I need to change the energy equation and add -gi/cp which is adiabatic lapse rate, i want to modify energy equation instead of solving potential temperature equation
because i want to consider radiation model (FVDOM- part of the structure is made of glass and i want to see the green house effect) and solar radiation

How can i modify the energy eq?
i want to follow the compressible equations of the bellow paper: