From ec4d44b28a3a8fae3315c82cb654b7f0d34da6c9 Mon Sep 17 00:00:00 2001 From: David A Minton Date: Tue, 3 Jan 2023 10:42:56 -0500 Subject: [PATCH] Fixed issue that was preventing discard information from being written to NetCDF --- src/swiftest/swiftest_discard.f90 | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/src/swiftest/swiftest_discard.f90 b/src/swiftest/swiftest_discard.f90 index c9c6df340..bc78f8454 100644 --- a/src/swiftest/swiftest_discard.f90 +++ b/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)