From 858003fd9c31f0430a57683c94a8296997151c7c Mon Sep 17 00:00:00 2001 From: David A Minton Date: Thu, 17 Feb 2022 16:34:49 -0500 Subject: [PATCH] Made relative velocity in symba_encounter_check to be heliocentric instead of barycentric to stay consistent with Swifter --- src/symba/symba_encounter_check.f90 | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/src/symba/symba_encounter_check.f90 b/src/symba/symba_encounter_check.f90 index 7281916b2..702bd09bd 100644 --- 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