Fixed incorrect coordinate transform.

MintonGroup · May 18, 2021 · 6fd0d8f · 6fd0d8f
1 parent e1ec26d
commit 6fd0d8f
Showing 1 changed file with 1 addition and 3 deletions.
diff --git 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