Fixed issue where the undefined nplm was being reported in non-SyMBA …

…runs
MintonGroup · May 18, 2023 · db6ad98 · db6ad98
1 parent 29afacf
commit db6ad98
Show file tree

Hide file tree

Showing 2 changed files with 2 additions and 2 deletions.
diff --git a/src/swiftest/swiftest_io.f90 b/src/swiftest/swiftest_io.f90
@@ -285,7 +285,7 @@ module subroutine swiftest_io_display_run_information(self, param, integration_t
          call self%compact_output(param,integration_timer)
       end if
 
-      if (self%pl%nplm > 0) then
+      if (param%lmtiny_pl) then
 #ifdef COARRAY
          if (param%lcoarray) then
             write(param%display_unit, co_symbastatfmt) this_image(),self%t, tfrac, self%pl%nbody, self%pl%nplm, self%tp%nbody

diff --git a/src/swiftest/swiftest_module.f90 b/src/swiftest/swiftest_module.f90
@@ -248,7 +248,7 @@ module swiftest
       integer(I8B)                                         :: nplpl   !! Number of body-body comparisons in the flattened upper triangular matrix
       type(swiftest_kinship),  dimension(:),   allocatable :: kin        !! Array of merger relationship structures that can account for multiple pairwise mergers in a single step
       logical,                 dimension(:),   allocatable :: lmtiny     !! flag indicating whether this body is below the GMTINY cutoff value
-      integer(I4B)                                         :: nplm       !! number of bodies above the GMTINY limit
+      integer(I4B)                                         :: nplm = 0   !! number of bodies above the GMTINY limit
       integer(I8B)                                         :: nplplm     !! Number of body (all massive)-body (only those above GMTINY) comparisons in the flattened upper triangular matrix 
       integer(I4B),            dimension(:),   allocatable :: nplenc     !! number of encounters with other planets this time step
       integer(I4B),            dimension(:),   allocatable :: ntpenc     !! number of encounters with test particles this time step