diff --git a/src/collision/collision_util.f90 b/src/collision/collision_util.f90 index 9898e3710..c2d3a1924 100644 --- a/src/collision/collision_util.f90 +++ b/src/collision/collision_util.f90 @@ -503,6 +503,7 @@ module subroutine collision_util_snapshot(self, param, nbody_system, t, arg) case("before") ! Saves the states of the bodies involved in the collision before the collision is resolved associate (idx => nbody_system%collider%impactors%id, ncoll => nbody_system%collider%impactors%ncoll) + allocate(pl, mold=nbody_system%pl) call pl%setup(ncoll, param) pl%id(:) = nbody_system%pl%id(idx(:)) pl%Gmass(:) = nbody_system%pl%Gmass(idx(:)) diff --git a/src/swiftest/swiftest_setup.f90 b/src/swiftest/swiftest_setup.f90 index 85084bfe0..a08daf29d 100644 --- a/src/swiftest/swiftest_setup.f90 +++ b/src/swiftest/swiftest_setup.f90 @@ -83,8 +83,8 @@ module subroutine swiftest_setup_construct_system(nbody_system, param) allocate(symba_pl :: nbody_system%pl_adds) allocate(symba_pl :: nbody_system%pl_discards) - allocate(collision_list_pltp :: nbody_system%pltp_encounter) - allocate(collision_list_plpl :: nbody_system%plpl_encounter) + allocate(symba_list_pltp :: nbody_system%pltp_encounter) + allocate(symba_list_plpl :: nbody_system%plpl_encounter) allocate(collision_list_plpl :: nbody_system%plpl_collision) if (param%lenc_save_trajectory .or. param%lenc_save_closest) then