Fixing velocity vector

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