I think I finally got the velocity vectors right

MintonGroup · Jan 7, 2023 · 8b841fa · 8b841fa
1 parent fb9c966
commit 8b841fa
Showing 1 changed file with 1 addition and 2 deletions.
diff --git a/src/fraggle/fraggle_generate.f90 b/src/fraggle/fraggle_generate.f90
@@ -448,8 +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) = (impactors%bounce_unit(:) + vimp_unit(:)) * vsign(i) 
-                  fragments%vc(:,i) = vmag * .unit.fragments%vc(:,i) + vrot(:)
+                  fragments%vc(:,i) = vmag * vscale(i) * impactors%bounce_unit(:) * vsign(i) + vrot(:) + vmag * vimp_unit(:)
                end if
             end do