Decided it's best to take binding energy out of the orbit energy calc

MintonGroup · Jan 11, 2023 · 5fa50c0 · 5fa50c0
1 parent 310598f
commit 5fa50c0
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/swiftest/swiftest_io.f90 b/src/swiftest/swiftest_io.f90
@@ -141,7 +141,7 @@ module subroutine swiftest_io_conservation_report(self, param, lterminal)
          te_now = nbody_system%te
          L_orbit_now(:) = nbody_system%L_orbit(:)
          L_spin_now(:) = nbody_system%L_spin(:)
-         E_orbit_now = ke_orbit_now + pe_now + be_now
+         E_orbit_now = ke_orbit_now + pe_now
          L_total_now(:) = nbody_system%L_total(:) + nbody_system%L_escape(:)
          GMtot_now = nbody_system%GMtot + nbody_system%GMescape