Put in the oblateness accelerations for central body for test particl…

…es in outer encounter of RMVS

src/rmvs/rmvs_step.f90

@@ -202,6 +202,10 @@ subroutine rmvs_step_out(cb, pl, tp, nbody_system, param, t, dt)
                call tp%step(nbody_system, param, outer_time, dto)
                tp%lfirst = lfirsttp
             else
+               if (param%loblatecb) then
+                  call swiftest_obl_acc(npl, cb%Gmass, cb%j2rp2, cb%j4rp4, pl%rbeg, pl%lmask, pl%outer(outer_index-1)%aobl, pl%Gmass, cb%aoblbeg)
+                  call swiftest_obl_acc(npl, cb%Gmass, cb%j2rp2, cb%j4rp4, pl%rend, pl%lmask, pl%outer(outer_index)%aobl, pl%Gmass, cb%aoblend)
+               end if
                call tp%step(nbody_system, param, outer_time, dto)
             end if
             do j = 1, npl

src/rmvs/rmvs_util.f90

@@ -298,13 +298,15 @@ module subroutine rmvs_util_setup_pl(self, n, param)
                allocate(pl%outer(i)%v(NDIM, n))
                pl%outer(i)%x(:,:) = 0.0_DP
                pl%outer(i)%v(:,:) = 0.0_DP
+               allocate(pl%outer(i)%aobl(NDIM, n))
+               pl%outer(i)%aobl(:,:) = 0.0_DP
             end do
             do i = 0, NTPHENC
                allocate(pl%inner(i)%x(NDIM, n))
                allocate(pl%inner(i)%v(NDIM, n))
-               allocate(pl%inner(i)%aobl(NDIM, n))
                pl%inner(i)%x(:,:) = 0.0_DP
                pl%inner(i)%v(:,:) = 0.0_DP
+               allocate(pl%inner(i)%aobl(NDIM, n))
                pl%inner(i)%aobl(:,:) = 0.0_DP
             end do
             ! if (param%ltides) then