diff --git a/src/symba/symba_util.f90 b/src/symba/symba_util.f90 index a7c4c6ae7..50521f7a8 100644 --- 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))