diff --git a/src/collision/collision_generate.f90 b/src/collision/collision_generate.f90 index b3ef6303e..7bcf7bdbe 100644 --- a/src/collision/collision_generate.f90 +++ b/src/collision/collision_generate.f90 @@ -218,7 +218,7 @@ module subroutine collision_generate_merge(self, nbody_system, param, t) if (param%lrotation) then do concurrent(i = 1:NDIM) fragments%Ip(i,1) = sum(impactors%mass(:) * impactors%Ip(i,:)) - L_spin_new(i) = sum(impactors%L_orbit(i,:) + impactors%L_orbit(i,:)) + L_spin_new(i) = sum(impactors%L_orbit(i,:) + impactors%L_spin(i,:)) end do fragments%Ip(:,1) = fragments%Ip(:,1) / fragments%mass(1) fragments%rot(:,1) = L_spin_new(:) / (fragments%Ip(3,1) * fragments%mass(1) * fragments%radius(1)**2)