diff --git a/src/swiftest/swiftest_io.f90 b/src/swiftest/swiftest_io.f90
index 1924fc2e8..3347511f4 100644
--- a/src/swiftest/swiftest_io.f90
+++ b/src/swiftest/swiftest_io.f90
@@ -1157,7 +1157,6 @@ module subroutine swiftest_io_netcdf_read_hdr_system(self, nc, param)
          status = nf90_inq_varid(nc%id, nc%status_varname, nc%status_varid) 
          if (status == NF90_NOERR) then
             call netcdf_io_check( nf90_get_var(nc%id, nc%status_varid,  body_status, start=[1,tslot], count=[idmax,1]), "netcdf_io_read_hdr_system nf90_getvar status_varid"  )
-            where(body_status(:) /= INACTIVE) body_status(:) = ACTIVE
          else
             body_status(:) = ACTIVE
          end if
@@ -1167,8 +1166,8 @@ module subroutine swiftest_io_netcdf_read_hdr_system(self, nc, param)
          plmask(1) = .false. ! This is the central body
 
          ! Select only active bodies
-         plmask(:) = plmask(:) .and. (body_status(:) == ACTIVE)
-         tpmask(:) = tpmask(:) .and. (body_status(:) == ACTIVE)
+         plmask(:) = plmask(:) .and. (body_status(:) /= INACTIVE)
+         tpmask(:) = tpmask(:) .and. (body_status(:) /= INACTIVE)
 
          if (param%lmtiny_pl) then
             where(plmask(:))