Merge branch 'debug'

MintonGroup · Aug 18, 2021 · cdb972a · cdb972a
2 parents 97e4933 + 5fb1e63
commit cdb972a
Show file tree

Hide file tree

Showing 3 changed files with 5 additions and 2 deletions.
diff --git a/examples/symba_chambers_2013/param.in b/examples/symba_chambers_2013/param.in
@@ -18,14 +18,15 @@ 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)
 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
@@ -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
@@ -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