From 34ec1644a4e0f4c6a9a00302a50ae89ae74a49f6 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Tue, 17 Jan 2023 19:53:37 -0500 Subject: [PATCH] Fixed issue where the rotation of impactors was not being saved properly --- src/collision/collision_resolve.f90 | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) 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.