Fixed additional memory errors when there are 0 of one kind of body i…

…n the system

Makefile.Defines

@@ -65,13 +65,13 @@ GPAR = -fopenmp -ftree-parallelize-loops=4
 GMEM = -fsanitize=undefined -fsanitize=address -fsanitize=leak 
 GWARNINGS = -Wall -Warray-bounds  -Wimplicit-interface -Wextra -Warray-temporaries
 
-FFLAGS  = $(IDEBUG) $(HEAPARR)
+#FFLAGS  = $(IDEBUG) $(HEAPARR)
 #FFLAGS  = -init=snan,arrays -no-wrap-margin -O3 $(STRICTREAL) $(SIMDVEC) $(PAR) 
-FORTRAN = ifort
+#FORTRAN = ifort
 #AR     = xiar
 
-#FORTRAN = gfortran
-#FFLAGS =  -ffree-line-length-none $(GDEBUG) $(GMEM)
+FORTRAN = gfortran
+FFLAGS =  -ffree-line-length-none $(GDEBUG) #$(GMEM)
 AR     = ar
 
 # DO NOT include in CFLAGS the "-c" option to compile object only

src/discard/discard.f90

@@ -8,13 +8,19 @@ module subroutine discard_system(self, param)
       !! Calls the discard methods for each body class and then the write method if any discards were detected
       !!
       implicit none
+      ! Arguments
       class(swiftest_nbody_system), intent(inout) :: self   !! Swiftest system object
       class(swiftest_parameters),   intent(in)    :: param  !! Current run configuration parameters
+      ! Internals
+      logical :: lany_discards
 
       associate(system => self, tp => self%tp, pl => self%pl)
-         call tp%discard(system, param)
+         lany_discards = .false.
          call pl%discard(system, param)
-         if (any(tp%ldiscard(:) .or. any(pl%ldiscard(:)))) call system%write_discard(param)
+         call tp%discard(system, param)
+         if (tp%nbody > 0) lany_discards = lany_discards .or.  any(tp%ldiscard(:))
+         if (pl%nbody > 0) lany_discards = lany_discards .or.  any(pl%ldiscard(:))
+         if (lany_discards) call system%write_discard(param)
       end associate
 
       return
@@ -31,6 +37,7 @@ module subroutine discard_pl(self, system, param)
       class(swiftest_pl),           intent(inout) :: self   !! Swiftest massive body object
       class(swiftest_nbody_system), intent(inout) :: system !! Swiftest nbody system object
       class(swiftest_parameters),   intent(in)    :: param  !! Current run configuration parameter
+
       self%ldiscard(:) = .false.
 
       return

src/modules/symba_classes.f90

@@ -387,9 +387,10 @@ module subroutine symba_step_interp_system(self, param, t, dt)
          real(DP),                   intent(in)    :: dt    !! Current stepsize
       end subroutine symba_step_interp_system
 
-      module subroutine symba_step_set_recur_levels_system(self)
+      module subroutine symba_step_set_recur_levels_system(self, ireci)
          implicit none
-         class(symba_nbody_system),  intent(inout) :: self !! SyMBA nbody system objec
+         class(symba_nbody_system),  intent(inout) :: self  !! SyMBA nbody system objec
+         integer(I4B),               intent(in)    :: ireci !! Input recursion level
       end subroutine symba_step_set_recur_levels_system
 
       module recursive subroutine symba_step_recur_system(self, param, t, ireci)

src/symba/symba_step.f90

@@ -94,14 +94,15 @@ module subroutine symba_step_interp_system(self, param, t, dt)
    end subroutine symba_step_interp_system
 
 
-   module subroutine symba_step_set_recur_levels_system(self)
+   module subroutine symba_step_set_recur_levels_system(self, ireci)
       !! author: David A. Minton
       !!
       !! Resets pl, tp,and encounter structures at the start of a new step
       !! 
       implicit none
       ! Arguments
-      class(symba_nbody_system),  intent(inout) :: self !! SyMBA nbody system object
+      class(symba_nbody_system),  intent(inout) :: self  !! SyMBA nbody system object
+      integer(I4B),               intent(in)    :: ireci !! Input recursion level 
       ! Internals
       integer(I4B) :: i, irecp
 
@@ -115,20 +116,21 @@ module subroutine symba_step_set_recur_levels_system(self)
                           plind3 => pltpenc_list%index1(1:pltpenc_list%nenc), &
                           tpind  => pltpenc_list%index2(1:pltpenc_list%nenc))
 
-                  irecp = system%irec + 1
+                  irecp = ireci + 1
 
                   do i = 1, plplenc_list%nenc
-                     if (pl%levelg(plind1(i)) == irecp) pl%levelg(plind1(i)) = system%irec
-                     if (pl%levelg(plind2(i)) == irecp) pl%levelg(plind2(i)) = system%irec
+                     if (pl%levelg(plind1(i)) == irecp) pl%levelg(plind1(i)) = ireci
+                     if (pl%levelg(plind2(i)) == irecp) pl%levelg(plind2(i)) = ireci
                   end do
                   do i = 1, pltpenc_list%nenc
-                     if (pl%levelg(plind3(i)) == irecp) pl%levelg(plind3(i)) = system%irec
-                     if (tp%levelg(tpind(i))  == irecp) tp%levelg(tpind(i))  = system%irec
+                     if (pl%levelg(plind3(i)) == irecp) pl%levelg(plind3(i)) = ireci
+                     if (tp%levelg(tpind(i))  == irecp) tp%levelg(tpind(i))  = ireci
                   end do
                end associate
 
-               where(plplenc_list%level(1:plplenc_list%nenc) == irecp) plplenc_list%level(:) = system%irec
-               where(pltpenc_list%level(1:pltpenc_list%nenc) == irecp) pltpenc_list%level(:) = system%irec
+               if (plplenc_list%nenc > 0) where(plplenc_list%level(1:plplenc_list%nenc) == irecp) plplenc_list%level(:) = ireci
+               if (pltpenc_list%nenc > 0) where(pltpenc_list%level(1:pltpenc_list%nenc) == irecp) pltpenc_list%level(:) = ireci
+               system%irec = ireci
             end select
          end select
       end associate
@@ -205,7 +207,7 @@ module recursive subroutine symba_step_recur_system(self, param, t, ireci)
                      call pltpenc_list%collision_check(system, param, t+dtl, dtl, ireci) 
                   end if
 
-                  call self%set_recur_levels()
+                  call self%set_recur_levels(ireci)
 
                end do
             end select