Fixed bug in angular momentum conservation in mergers

MintonGroup · Jan 17, 2023 · 9683841 · 9683841
1 parent 2defff3
commit 9683841
Showing 1 changed file with 1 addition and 1 deletion.
diff --git 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)