Horizontally and temporally averaged heating and cooling contributions to the internal energy equation as a function of height. The plotted terms are the radiative terms Q_loss(orange), Q_rad(blue), divergence of heat flux (∇ ∙ q) (green), joule heating from numerical resistivity (η_num J², (red)), dissipation from numerical viscosity (∇ ∙ (Π_num ∙ v)), (violet), adiabatic heating (p∇ ∙ v), (navy blue), and flux of internal energy (∇ ∙ (e_intv)), (light green). Left panels show heating terms, and right panels show cooling terms. The top row corresponds to the simulation with a vertical magnetic field (Simulation 2), while the bottom row corresponds to the simulation with an arcade-like magnetic field (Simulation 3). The black dotted line in each panel represents the total contribution (sum of all terms). For clarity, values smaller than 1 × 10⁻¹⁰ are omitted from the plots.