From 9683841a137092e2534650e8c136ff74d4f5df07 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Tue, 17 Jan 2023 10:47:58 -0500 Subject: [PATCH] Fixed bug in angular momentum conservation in mergers --- src/collision/collision_generate.f90 | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) 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)