From a7ee2927c4a172f702c2429091692e1afa98773b Mon Sep 17 00:00:00 2001 From: David A Minton Date: Tue, 1 Mar 2022 17:25:51 -0500 Subject: [PATCH] Switched vh2vb coordinate conversion back to the original Swifter algorithm --- src/util/util_coord.f90 | 5 +++-- 1 file changed, 3 insertions(+), 2 deletions(-) diff --git a/src/util/util_coord.f90 b/src/util/util_coord.f90 index 41d230878..e261f3be7 100644 --- a/src/util/util_coord.f90 +++ b/src/util/util_coord.f90 @@ -198,9 +198,10 @@ module subroutine util_coord_vh2vb_pl(self, cb) associate(pl => self, npl => self%nbody) Gmtot = cb%Gmass + sum(pl%Gmass(1:npl)) cb%vb(:) = 0.0_DP - do i = npl, 1, -1 - cb%vb(:) = cb%vb(:) - pl%Gmass(i) * pl%vh(:, i) / Gmtot + do i = 1, npl + cb%vb(:) = cb%vb(:) - pl%Gmass(i) * pl%vh(:, i) end do + cb%vb(:) = cb%vb(:) / Gmtot do concurrent(i = 1:npl) pl%vb(:, i) = pl%vh(:, i) + cb%vb(:) end do