diff --git a/src/swiftest/swiftest_driver.f90 b/src/swiftest/swiftest_driver.f90 index eed4fe446..8cc7f9a2a 100644 --- a/src/swiftest/swiftest_driver.f90 +++ b/src/swiftest/swiftest_driver.f90 @@ -106,7 +106,6 @@ program swiftest_driver if (idump == dump_cadence) then idump = 0 call nbody_system%dump(param) - end if call integration_timer%report(message="Integration steps:", unit=display_unit, nsubsteps=istep_out) diff --git a/src/swiftest/swiftest_io.f90 b/src/swiftest/swiftest_io.f90 index a2b604357..af434b2c5 100644 --- a/src/swiftest/swiftest_io.f90 +++ b/src/swiftest/swiftest_io.f90 @@ -129,7 +129,11 @@ module subroutine swiftest_io_conservation_report(self, param, lterminal) associate(nbody_system => self, pl => self%pl, cb => self%cb, npl => self%pl%nbody, display_unit => param%display_unit, nc => self%system_history%nc) - call pl%vb2vh(cb) + select type(self) + class is (helio_nbody_system) ! Don't convert vh to vb for Helio-based integrators, because they are already have that calculated + class is (whm_nbody_system) + call pl%vh2vb(cb) + end select call pl%rh2rb(cb) call nbody_system%get_energy_and_momentum(param)