Skip to content
This repository was archived by the owner on Aug 28, 2024. It is now read-only.

Commit

Permalink
Fixed issue where failed fragmentation event lead to lost particles i…
Browse files Browse the repository at this point in the history
…nstead of pure hit and run
  • Loading branch information
daminton committed Aug 12, 2021
1 parent f13b6c7 commit 9d7a741
Showing 1 changed file with 12 additions and 0 deletions.
12 changes: 12 additions & 0 deletions src/symba/symba_fragmentation.f90
Original file line number Diff line number Diff line change
Expand Up @@ -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
Expand Down Expand Up @@ -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
Expand Down

0 comments on commit 9d7a741

Please sign in to comment.