diff --git a/src/collision/collision_resolve.f90 b/src/collision/collision_resolve.f90 index ae928b39b..b50b21334 100644 --- a/src/collision/collision_resolve.f90 +++ b/src/collision/collision_resolve.f90 @@ -134,7 +134,10 @@ module subroutine collision_resolve_consolidate_impactors(self, nbody_system, pa end if density(j) = impactors%mass(j) / volume(j) impactors%radius(j) = (3 * volume(j) / (4 * PI))**(1.0_DP / 3.0_DP) - if (param%lrotation) impactors%Ip(:, j) = impactors%Ip(:, j) / impactors%mass(j) + if (param%lrotation) then + impactors%Ip(:, j) = impactors%Ip(:, j) / impactors%mass(j) + impactors%rot(:,j) = impactors%L_spin(:, j) / (impactors%Ip(3,j) * impactors%mass(j) * impactors%radius(j)**2) + end if end do lflag = .true.