Fixed status variable output

MintonGroup · Jan 5, 2023 · 64fbf6d · 64fbf6d
1 parent 825c2d9
commit 64fbf6d
Show file tree

Hide file tree

Showing 3 changed files with 2 additions and 6 deletions.
diff --git a/python/swiftest/swiftest/io.py b/python/swiftest/swiftest/io.py
@@ -63,7 +63,7 @@
 # handles strings differently than Python's Xarray.
 string_varnames = ["name", "particle_type", "status", "origin_type", "stage", "regime"]
 char_varnames = ["space"]
-int_varnames = ["id", "ntp", "npl", "nplm", "discard_body_id", "collision_id", "loopnum"]
+int_varnames = ["id", "ntp", "npl", "nplm", "discard_body_id", "collision_id", "status"]
 
 def bool2yesno(boolval):
     """

diff --git a/src/swiftest/swiftest_io.f90 b/src/swiftest/swiftest_io.f90
@@ -1557,7 +1557,7 @@ module subroutine swiftest_io_netcdf_write_frame_cb(self, nc, param)
 
          idslot = self%id + 1
          call netcdf_io_check( nf90_put_var(nc%id, nc%id_varid, self%id, start=[idslot]), "netcdf_io_write_frame_cb nf90_put_var cb id_varid"  )
-         call netcdf_io_check( nf90_put_var(nc%id, nc%status_varid, ACTIVE, start=[idslot]), "netcdf_io_write_frame_cb nf90_put_var cb id_varid"  )
+         call netcdf_io_check( nf90_put_var(nc%id, nc%status_varid, ACTIVE, start=[idslot, tslot]), "netcdf_io_write_frame_cb nf90_put_var cb id_varid"  )
 
          call netcdf_io_check( nf90_put_var(nc%id, nc%Gmass_varid, self%Gmass, start=[idslot, tslot]), "netcdf_io_write_frame_cb nf90_put_var cb Gmass_varid"  )
          if (param%lclose) call netcdf_io_check( nf90_put_var(nc%id, nc%radius_varid, self%radius, start=[idslot, tslot]), "netcdf_io_write_frame_cb nf90_put_var cb radius_varid"  )

diff --git a/src/swiftest/swiftest_util.f90 b/src/swiftest/swiftest_util.f90
@@ -1717,10 +1717,6 @@ module subroutine swiftest_util_rearray_pl(self, nbody_system, param)
             end where
          end if
 
-         call nc%open(param)
-         call pl%write_info(nc, param)
-         call nc%close()
-
          ! Reindex the new list of bodies 
          call pl%sort("mass", ascending=.false.)
          call pl%flatten(param)