diff --git a/sources/integrals.F90 b/sources/integrals.F90
index e7de725..c31eafd 100644
--- a/sources/integrals.F90
+++ b/sources/integrals.F90
@@ -179,9 +179,10 @@ module integrals
 
 ! write the integral file header
 !
-      write(funit,"('#',a8,10(1x,a18))") 'step', 'time', 'dt'                  &
+      write(funit,"('#',a8,13(1x,a18))") 'step', 'time', 'dt'                  &
                                        , 'mass', 'momx', 'momy', 'momz'        &
-                                       , 'ener', 'ekin', 'emag', 'eint'
+                                       , 'ener', 'ekin', 'emag', 'eint'        &
+                                       , 'einj', 'erat', 'arms'
       write(funit,"('#')")
 
 ! depending on the append parameter, choose the right file initialization for
@@ -321,6 +322,7 @@ module integrals
     use equations      , only : ien, imx, imy, imz
     use equations      , only : magnetized, gamma, csnd
     use evolution      , only : step, time, dtn
+    use forcing        , only : einj, rinj, arms
     use mpitools       , only : master
 #ifdef MPI
     use mpitools       , only : reduce_sum_real_array
@@ -458,6 +460,11 @@ module integrals
 #endif /* NDIMS == 3 */
       end if
 
+! sum up the injected energy and injection rate
+!
+      inarr( 9) = einj
+      inarr(10) = rinj
+
 ! get average, minimum and maximum values of density
 !
 #if NDIMS == 3
@@ -578,7 +585,7 @@ module integrals
 ! write down the integrals and statistics to appropriate files
 !
     if (master) then
-      write(funit,"(i9,10(1x,1es18.8e3))") step, time, dtn, inarr(1:8)
+      write(funit,"(i9,13(1x,1es18.8e3))") step, time, dtn, inarr(1:10), arms
       write(sunit,"(i9,23(1x,1es18.8e3))") step, time                          &
                                             , avarr(1), mnarr(1), mxarr(1)     &
                                             , avarr(2), mnarr(2), mxarr(2)     &