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

Commit

Permalink
Switched to adding new id to largest hit and run body to prevent id d…
Browse files Browse the repository at this point in the history
…uplication
  • Loading branch information
daminton committed Aug 18, 2021
1 parent 98a4d5d commit 10726e9
Showing 1 changed file with 1 addition and 1 deletion.
2 changes: 1 addition & 1 deletion src/symba/symba_collision.f90
Original file line number Diff line number Diff line change
Expand Up @@ -169,7 +169,7 @@ module function symba_collision_casehitandrun(system, param, family, x, v, mass,
m_frag(2:nfrag) = (mtot - m_frag(1)) / (nfrag - 1)
rad_frag(2:nfrag) = (3 * m_frag(2:nfrag) / (4 * PI * avg_dens))**(1.0_DP / 3.0_DP)
m_frag(nfrag) = m_frag(nfrag) + (mtot - sum(m_frag(:)))
id_frag(2:nfrag) = [(i, i = system%maxid + 1, system%maxid + nfrag - 1)]
id_frag(1:nfrag) = [(i, i = system%maxid + 1, system%maxid + nfrag)]

do i = 1, nfrag
Ip_frag(:, i) = Ip(:, jproj)
Expand Down

0 comments on commit 10726e9

Please sign in to comment.