diff --git a/src/cooling/const_du/cooling.h b/src/cooling/const_du/cooling.h
index 8f005ee0bf3d866ed196bd7b6a1e1094c39f22ef..40bf2e919d84c278b1860dcc3360e887aa3be71c 100644
--- a/src/cooling/const_du/cooling.h
+++ b/src/cooling/const_du/cooling.h
@@ -79,6 +79,9 @@ __attribute__((always_inline)) INLINE static void cooling_cool_part(
 
   /* Update the internal energy */
   hydro_set_internal_energy(p, u_new);
+
+  /* Store the radiated energy */
+  xp->cooling_data.radiated_energy += u_new - u_old;
 }
 
 /**
diff --git a/src/cooling/const_lambda/cooling.h b/src/cooling/const_lambda/cooling.h
index 458babf540a40db1f42da4750e0bb4c9a6ed8ce2..6d3f904d15ac8eb693175037294fccf73bbf1b9b 100644
--- a/src/cooling/const_lambda/cooling.h
+++ b/src/cooling/const_lambda/cooling.h
@@ -114,6 +114,9 @@ __attribute__((always_inline)) INLINE static void cooling_cool_part(
    * " */
   /*       "%g, du_dt*dt = %g, u_old + du_dt*dt = %g, u_new = %g\n", */
   /*       u_old, du_dt, dt, du_dt * dt, u_new, u_new_test); */
+
+  /* Store the radiated energy */
+  xp->cooling_data.radiated_energy += u_new - u_old;
 }
 
 /**