Merge branch 'debug'

src/symba/symba_fragmentation.f90

@@ -81,6 +81,12 @@ module function symba_fragmentation_casedisruption(system, param, family, x, v,
          write(*,*) 'No fragment solution found, so treat as a pure hit-and-run'
          status = ACTIVE 
          nfrag = 0
+         select type(pl => system%pl)
+         class is (symba_pl)
+            pl%status(family(:)) = status
+            pl%ldiscard(family(:)) = .false.
+            pl%lcollision(family(:)) = .false.
+         end select
       else
          ! Populate the list of new bodies
          write(*,'("Generating ",I2.0," fragments")') nfrag
@@ -360,6 +366,12 @@ module function symba_fragmentation_casesupercatastrophic(system, param, family,
          write(*,*) 'No fragment solution found, so treat as a pure hit-and-run'
          status = ACTIVE 
          nfrag = 0
+         select type(pl => system%pl)
+         class is (symba_pl)
+            pl%status(family(:)) = status
+            pl%ldiscard(family(:)) = .false.
+            pl%lcollision(family(:)) = .false.
+         end select
       else
          ! Populate the list of new bodies
          write(*,'("Generating ",I2.0," fragments")') nfrag

src/util/util_append.f90

@@ -177,10 +177,9 @@ module subroutine util_append_body(self, source, lsource_mask)
          call util_append(self%capom, source%capom, nold, nsrc, lsource_mask)
          call util_append(self%omega, source%omega, nold, nsrc, lsource_mask)
          call util_append(self%capm, source%capm, nold, nsrc, lsource_mask)
+         self%nbody = nold + count(lsource_mask(:))
       end associate
 
-      self%nbody = count(self%status(:) /= INACTIVE)
-
       return
    end subroutine util_append_body