From fb9c96698a9e5ae2786e3a018ba108c8a7482a6f Mon Sep 17 00:00:00 2001 From: David A Minton Date: Sat, 7 Jan 2023 18:28:54 -0500 Subject: [PATCH] Reverted back to the old velocity vector. The "fix" made cascades --- src/fraggle/fraggle_generate.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/fraggle/fraggle_generate.f90 b/src/fraggle/fraggle_generate.f90 index 2c6cc9461..1883476e8 100644 --- a/src/fraggle/fraggle_generate.f90 +++ b/src/fraggle/fraggle_generate.f90 @@ -291,7 +291,7 @@ module subroutine fraggle_generate_pos_vec(collider) ! Shift to the cloud center coordinates fragments%rc(:,i) = fragments%rc(:,i) + fragment_cloud_center(:,j) - ! Make sure that the fragments are positioned away from the impact point + ! Make sure that the fragments are positioned away from the impact point direction = dot_product(fragments%rc(:,i) - impactors%rbimp(:), fragment_cloud_center(:,j) - impactors%rbimp(:)) if (direction < 0.0_DP) then fragments%rc(:,i) = fragments%rc(:,i) - fragment_cloud_center(:,j) @@ -448,7 +448,7 @@ module subroutine fraggle_generate_vel_vec(collider, nbody_system, param, lfailu vmag = vesc * vscale(i) rimp(:) = fragments%rc(:,i) - impactors%rbimp(:) vimp_unit(:) = .unit. rimp(:) - fragments%vc(:,i) = vimp_unit(:) * vsign(i) + fragments%vc(:,i) = (impactors%bounce_unit(:) + vimp_unit(:)) * vsign(i) fragments%vc(:,i) = vmag * .unit.fragments%vc(:,i) + vrot(:) end if end do