Made relative velocity in symba_encounter_check to be heliocentric in…

…stead of barycentric to stay consistent with Swifter
MintonGroup · Feb 17, 2022 · 858003f · 858003f
1 parent af6c59b
commit 858003f
Showing 1 changed file with 2 additions and 2 deletions.
diff --git a/src/symba/symba_encounter_check.f90 b/src/symba/symba_encounter_check.f90
@@ -131,7 +131,7 @@ module function symba_encounter_check(self, param, system, dt, irec) result(lany
                   i = self%index1(k)
                   j = self%index2(k)
                   xr(:) = pl%xh(:,j) - pl%xh(:,i)
-                  vr(:) = pl%vb(:,j) - pl%vb(:,i)
+                  vr(:) = pl%vh(:,j) - pl%vh(:,i)
                   rcrit12 = pl%renc(i) + pl%renc(j)
                   call encounter_check_one(xr(1), xr(2), xr(3), vr(1), vr(2), vr(3), rcrit12, dt, lencounter(lidx), self%lvdotr(k))
                   if (lencounter(lidx)) then
@@ -146,7 +146,7 @@ module function symba_encounter_check(self, param, system, dt, irec) result(lany
                   i = self%index1(k)
                   j = self%index2(k)
                   xr(:) = tp%xh(:,j) - pl%xh(:,i)
-                  vr(:) = tp%vb(:,j) - pl%vb(:,i)
+                  vr(:) = tp%vh(:,j) - pl%vh(:,i)
                   call encounter_check_one(xr(1), xr(2), xr(3), vr(1), vr(2), vr(3), pl%renc(i), dt, &
                                            lencounter(lidx), self%lvdotr(k))
                   if (lencounter(lidx)) then