diff --git a/examples/symba_chambers_2013/param.in b/examples/symba_chambers_2013/param.in
index 52b09a882..28500d32b 100644
--- a/examples/symba_chambers_2013/param.in
+++ b/examples/symba_chambers_2013/param.in
@@ -18,6 +18,7 @@ OUT_STAT       REPLACE
 CHK_CLOSE      yes                ! check for planetary close encounters
 CHK_RMAX       100000.0            ! discard outside of
 EXTRA_FORCE    no                 ! no extra user-defined forces
+DISCARD_OUT    discard.out
 BIG_DISCARD    no                 ! output all planets if anything discarded
 RHILL_PRESENT  yes                ! Hill's sphere radii in input file
 MU2KG          1.98847e30         ! (M_sun-> kg)
@@ -25,7 +26,7 @@ DU2M           1.495979e11        ! distance unit to meters (AU --> m)
 TU2S           3.15569259747e7    ! time unit to seconds (years --> seconds)
 GMTINY         1e-10		  
 ENERGY         yes
+ENERGY_OUT     energy.dat
 ROTATION       yes
 FRAGMENTATION  yes
-DISCARD_OUT    discard.out
 SEED           8 12261555 871132 92734722 21132443 36344777 4334443 219291656 3848566
diff --git a/examples/symba_clement_2018/param.in b/examples/symba_clement_2018/param.in
index f90f0f89e..79890595e 100644
--- a/examples/symba_clement_2018/param.in
+++ b/examples/symba_clement_2018/param.in
@@ -19,13 +19,14 @@ CHK_CLOSE      yes                ! check for planetary close encounters
 CHK_RMAX       100000.0            ! discard outside of
 EXTRA_FORCE    no                 ! no extra user-defined forces
 BIG_DISCARD    no                 ! output all planets if anything discarded
+DISCARD_OUT    discard.out
 RHILL_PRESENT  yes                ! Hill's sphere radii in input file
 MU2KG          1.98847e30         ! (M_sun-> kg)
 DU2M           1.495979e11        ! distance unit to meters (AU --> m)
 TU2S           3.15569259747e7    ! time unit to seconds (years --> seconds)
 GMTINY         1e-10		  
 ENERGY         yes
+ENERGY_OUT     energy.dat
 ROTATION       yes
 FRAGMENTATION  yes
-DISCARD_OUT    discard.out
 SEED           8 12261555 871132 92734722 21132443 36344777 4334443 219291656 3848566
diff --git a/src/io/io.f90 b/src/io/io.f90
index 3f7a98cfc..4aa3b88a3 100644
--- a/src/io/io.f90
+++ b/src/io/io.f90
@@ -624,6 +624,7 @@ module subroutine io_param_reader(self, unit, iotype, v_list, iostat, iomsg)
          write(*,*) "ROTATION      = ", param%lrotation
          write(*,*) "TIDES         = ", param%ltides
          write(*,*) "ENERGY         = ",param%lenergy
+         if (param%lenergy) write(*,*) "ENERGY_OUT     = ",trim(adjustl(param%energy_out))
          write(*,*) "MU2KG          = ",param%MU2KG       
          write(*,*) "TU2S           = ",param%TU2S        
          write(*,*) "DU2M           = ",param%DU2M