Cleaned up code a bit.

MintonGroup · Jul 7, 2021 · 20dac50 · 20dac50
1 parent 986cfaf
commit 20dac50
Show file tree

Hide file tree

Showing 2 changed files with 5 additions and 10 deletions.
diff --git a/src/helio/helio_getacch.f90 b/src/helio/helio_getacch.f90
@@ -15,14 +15,8 @@ module subroutine helio_getacch_pl(self, system, param, t)
       class(swiftest_parameters),   intent(in)    :: param  !! Current run configuration parameters of 
       real(DP),                     intent(in)    :: t      !! Current simulation time
       ! Internals
-      logical, save                    :: lmalloc = .true.
-      integer(I4B)                     :: i
-      real(DP)                       :: r2
-      real(DP), dimension(:), allocatable, save    :: irh
-      real(DP), dimension(:, :), allocatable, save :: xh_loc, aobl
 
       associate(cb => system%cb, pl => self, npl => self%nbody)
-         pl%ahi(:,:) = 0.0_DP
          call helio_getacch_int_pl(pl, t)
          pl%ah(:,:) = pl%ahi(:,:)
          if (param%loblatecb) call pl%obl_acc(cb)
@@ -53,10 +47,9 @@ module subroutine helio_getacch_tp(self, system, param, t, xhp)
          real(DP)                                     :: r2, mu
          real(DP), dimension(:), allocatable, save    :: irh, irht
 
-         associate(tp => self, ntp => self%nbody, cb => system%cb, pl => system%pl, npl => system%pl%nbody)
+         associate(tp => self, ntp => self%nbody, cb => system%cb, npl => system%pl%nbody)
             select type(pl => system%pl)
             class is (helio_pl)
-               self%ahi(:,:) = 0.0_DP
                call helio_getacch_int_tp(tp, pl, t, xhp)
                tp%ah(:,:) = tp%ahi(:,:)
                if (param%loblatecb) call tp%obl_acc(cb) 
@@ -84,6 +77,7 @@ subroutine helio_getacch_int_pl(pl, t)
          real(DP), dimension(NDIM) :: dx
 
          associate(npl => pl%nbody)
+            pl%ahi(:,:) = 0.0_DP
             do i = 1, npl - 1
                do j = i + 1, npl
                   dx(:) = pl%xh(:,j) - pl%xh(:,i)
@@ -119,9 +113,10 @@ subroutine helio_getacch_int_tp(tp, pl, t, xhp)
          real(DP), dimension(NDIM) :: dx
 
          associate(ntp => tp%nbody, npl => pl%nbody)
+            tp%ahi(:,:) = 0.0_DP
             do i = 1, ntp
             if (tp%status(i) == ACTIVE) then
-               do j = 2, npl
+               do j = 1, npl
                   dx(:) = tp%xh(:,i) - xhp(:,j)
                   r2 = dot_product(dx(:), dx(:))
                   fac = pl%Gmass(j) / (r2 * sqrt(r2))

diff --git a/src/whm/whm_step.f90 b/src/whm/whm_step.f90
@@ -16,7 +16,7 @@ module subroutine whm_step_system(self, param, t, dt)
       real(DP),                   intent(in)    :: t     !! Current simulation time
       real(DP),                   intent(in)    :: dt    !! Current stepsize
 
-      associate(system => self, cb => self%cb,  pl => self%pl, tp => self%tp, ntp => self%tp%nbody, npl => self%pl%nbody)
+      associate(system => self, cb => self%cb,  pl => self%pl, tp => self%tp, ntp => self%tp%nbody)
          call pl%set_rhill(cb)
          call self%set_beg_end(xbeg = pl%xh)
          call pl%step(system, param, t, dt)