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: |