changed to barycentric velocity in encounter output

MintonGroup · Dec 5, 2022 · ae92326 · ae92326
1 parent 73b72c8
commit ae92326
Showing 1 changed file with 1 addition and 1 deletion.
diff --git a/src/symba/symba_util.f90 b/src/symba/symba_util.f90
@@ -1356,7 +1356,7 @@ module subroutine symba_util_take_encounter_snapshot(self, param, t)
                   allocate(snapshot%pl%vh(NDIM,npl_snap))
                   do i = 1, NDIM
                      snapshot%pl%rh(i,:) = pack(pl%rh(i,1:npl), pl%lmask(1:npl))
-                     snapshot%pl%vh(i,:) = pack(pl%vh(i,1:npl), pl%lmask(1:npl))
+                     snapshot%pl%vh(i,:) = pack(pl%vb(i,1:npl), pl%lmask(1:npl))
                   end do
                   if (param%lclose) then
                      allocate(snapshot%pl%radius(npl_snap))