diff --git a/src/symba/symba_fragmentation.f90 b/src/symba/symba_fragmentation.f90 index b36c54e9a..430b8f81d 100644 --- a/src/symba/symba_fragmentation.f90 +++ b/src/symba/symba_fragmentation.f90 @@ -30,7 +30,9 @@ module function symba_fragmentation_casedisruption(system, param, family, x, v, real(DP), dimension(:), allocatable :: m_frag, rad_frag integer(I4B), dimension(:), allocatable :: id_frag logical :: lfailure - + + write(*, '("Disruption between bodies ",I8,99(:,",",I8))') system%pl%id(family(:)) + ! Collisional fragments will be uniformly distributed around the pre-impact barycenter nfrag = NFRAG_DISRUPT allocate(m_frag(nfrag)) @@ -117,6 +119,8 @@ module function symba_fragmentation_casehitandrun(system, param, family, x, v, m logical :: lpure logical, dimension(system%pl%nbody) :: lmask + write(*, '("Hit and run between bodies ",I8,99(:,",",I8))') system%pl%id(family(:)) + mtot = sum(mass(:)) xcom(:) = (mass(1) * x(:,1) + mass(2) * x(:,2)) / mtot vcom(:) = (mass(1) * v(:,1) + mass(2) * v(:,2)) / mtot @@ -217,7 +221,7 @@ module function symba_fragmentation_casemerge(system, param, family, x, v, mass, select type(pl => system%pl) class is (symba_pl) - write(*, '("Merging bodies ",99(I8,",",:))') pl%id(family(:)) + write(*, '("Merging bodies ",I8,99(:,",",I8))') pl%id(family(:)) ibiggest = maxloc(pl%Gmass(family(:)), dim=1) id_frag(1) = pl%id(family(ibiggest)) @@ -304,6 +308,8 @@ module function symba_fragmentation_casesupercatastrophic(system, param, family, logical :: lfailure logical, dimension(system%pl%nbody) :: lmask + write(*, '("Supercatastrophic disruption between bodies ",I8,99(:,",",I8))') system%pl%id(family(:)) + ! Collisional fragments will be uniformly distributed around the pre-impact barycenter nfrag = NFRAG_SUPERCAT allocate(m_frag(nfrag))