Switched vh2vb coordinate conversion back to the original Swifter alg…

…orithm
MintonGroup · Mar 1, 2022 · a7ee292 · a7ee292
1 parent b1b7bc8
commit a7ee292
Showing 1 changed file with 3 additions and 2 deletions.
diff --git 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