Fixed issue that was preventing discard information from being writte…

…n to NetCDF

src/swiftest/swiftest_discard.f90

@@ -25,7 +25,7 @@ module subroutine swiftest_discard_system(self, param)
       lpl_check = allocated(self%pl_discards)
       ltp_check = allocated(self%tp_discards)
 
-      associate(nbody_system => self, tp => self%tp, pl => self%pl, tp_discards => self%tp_discards, pl_discards => self%pl_discards)
+      associate(nbody_system => self, tp => self%tp, pl => self%pl, tp_discards => self%tp_discards, pl_discards => self%pl_discards, nc => param%system_history%nc)
          lpl_discards = .false.
          ltp_discards = .false.
          if (lpl_check) then
@@ -37,9 +37,11 @@ module subroutine swiftest_discard_system(self, param)
             call tp%discard(nbody_system, param)
             ltp_discards = (tp_discards%nbody > 0)
          end if
+         call nc%open(param)
          if (ltp_discards) call tp_discards%write_info(param%system_history%nc, param)
          if (lpl_discards) call pl_discards%write_info(param%system_history%nc, param)
          if (lpl_discards .and. param%lenergy) call self%conservation_report(param, lterminal=.false.)
+         call nc%close()
          if (lpl_check) call pl_discards%setup(0,param) 
          if (ltp_check) call tp_discards%setup(0,param)