diff --git a/src/symba/symba_collision.f90 b/src/symba/symba_collision.f90 index f813e068a..7c49009f2 100644 --- a/src/symba/symba_collision.f90 +++ b/src/symba/symba_collision.f90 @@ -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)