Fixed issue where failed fragmentation didn't get the copy of the "af…

…ter" system for recording
MintonGroup · Dec 18, 2022 · 348b8de · 348b8de
1 parent 4336100
commit 348b8de
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 2 deletions.
diff --git a/src/collision/collision_io.f90 b/src/collision/collision_io.f90
@@ -227,9 +227,9 @@ module subroutine collision_io_write_frame_snapshot(self, nc, param)
                   if (allocated(pl)) deallocate(pl)
                   select case(stage)
                   case(1)
-                     allocate(pl, source=system%before%pl)
+                     allocate(pl, source=self%collision_system%before%pl)
                   case(2)
-                     allocate(pl, source=system%after%pl)
+                     allocate(pl, source=self%collision_system%after%pl)
                   end select
                   npl = pl%nbody
                   do i = 1, npl

diff --git a/src/symba/symba_collision.f90 b/src/symba/symba_collision.f90
@@ -60,6 +60,7 @@ module function symba_collision_casedisruption(system, param, t)  result(status)
                pl%ldiscard(impactors%idx(:)) = .false.
                pl%lcollision(impactors%idx(:)) = .false.
             end select
+            allocate(collision_system%after%pl, source=collision_system%before%pl) ! Be sure to save the pl so that snapshots still work 
          else
             ! Populate the list of new bodies
             nfrag = fragments%nbody