diff --git a/src/symba/symba_frag_pos.f90 b/src/symba/symba_frag_pos.f90 index 3e3bca5d4..165995f54 100644 --- a/src/symba/symba_frag_pos.f90 +++ b/src/symba/symba_frag_pos.f90 @@ -581,10 +581,8 @@ subroutine set_fragment_radial_velocities(lmerge) ! Recombine the tangential and radial components into the final velocity vector do i = 1, nfrag - v_r_mag(i) = dot_product(v_r(:,i), v_r_unit(:, i)) - v_frag(:, i) = v_r_mag(i) * v_r_unit(:, i) + v_t_mag(i) * v_t_unit(:, i) + v_frag(:, i) = v_r(:, i) + v_t_mag(i) * v_t_unit(:, i) end do - call shift_vector_to_origin(m_frag, v_frag) end if do i = 1, nfrag