Fix radiated energy saved to statistics for EAGLE cooling.
(cherry picked from commit 340658c0)
I think the radiated energy stored in xp.cooling_data was incorrect for the EAGLE cooling scheme.
The variable cooling_du_dt actually stores the total du_dt (hydro + cooling), so the hydro_du_dt needs to be subtracted when computing the radiated energy.
Some examples (from CoolingSedovBlast3D):
| Without fix | With fix | |
|---|---|---|
| Sphenix: | ![]() |
![]() |
| Gizmo-MFV: | ![]() |
![]() |



