Left: cooling powers of the included cooling molecules as functions of radius. The total line cooling power (Λ_line) is also shown in black. Right: breakdown of the cooling and heating terms (indicated in the figure; see also Sect. 2.2) as a function of radius. Γ_comp denotes the first term in Eq. (5), that is, heating/cooling by gas compression/expansion, converted to the same units as the other heating and cooling terms (positive values indicate heating). The plots correspond to t = 5 × 10⁵ yr.