Added status check to vb2vh conversion

MintonGroup · Mar 20, 2023 · 2739286 · 2739286
1 parent cbece6c
commit 2739286
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/swiftest/swiftest_util.f90 b/src/swiftest/swiftest_util.f90
@@ -494,7 +494,7 @@ module subroutine swiftest_util_coord_vb2vh_pl(self, cb)
       associate(pl => self, npl => self%nbody)
          cb%vb(:) = 0.0_DP
          do i = npl, 1, -1
-            cb%vb(:) = cb%vb(:) - pl%Gmass(i) * pl%vb(:, i) / cb%Gmass
+            if (pl%status(i) /= INACTIVE) cb%vb(:) = cb%vb(:) - pl%Gmass(i) * pl%vb(:, i) / cb%Gmass
          end do
          do concurrent(i = 1:npl)
             pl%vh(:, i) = pl%vb(:, i) - cb%vb(:)