Same thing as before but with less code

MintonGroup · Jan 5, 2023 · 96ab9f1 · 96ab9f1
1 parent 210e1e9
commit 96ab9f1
Showing 1 changed file with 2 additions and 3 deletions.
diff --git 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(:))