From 6fd0d8f0e4a5afe03f20aec50457d37d89208e2f Mon Sep 17 00:00:00 2001 From: David Minton Date: Tue, 18 May 2021 08:53:02 -0400 Subject: [PATCH] Fixed incorrect coordinate transform. --- src/symba/symba_frag_pos.f90 | 4 +--- 1 file changed, 1 insertion(+), 3 deletions(-) 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