Fixed a problem with id values getting shifted by one. We count from …

…0 because the central body is always 0

src/netcdf/netcdf.f90

@@ -808,7 +808,7 @@ module subroutine netcdf_read_particle_info_system(self, nc, param, plmask, tpma
                call pl%info(i)%set_value(status="ACTIVE")
             end do
             allocate(plind(npl))
-            plind(:) = pack([(i, i = 1, idmax)], plmask(:))
+            plind(:) = pack([(i, i = 0, idmax-1)], plmask(:))
          end if
          if (ntp > 0) then
             tp%status(:) = ACTIVE
@@ -817,16 +817,16 @@ module subroutine netcdf_read_particle_info_system(self, nc, param, plmask, tpma
                call tp%info(i)%set_value(status="ACTIVE")
             end do
             allocate(tpind(ntp))
-            tpind(:) = pack([(i, i = 1, idmax)], tpmask(:))
+            tpind(:) = pack([(i, i = 0, idmax-1)], tpmask(:))
          end if
 
          call check( nf90_get_var(nc%id, nc%id_varid, itemp), "netcdf_read_particle_info_system nf90_getvar id_varid"  )
          cb%id = itemp(1)
          pl%id(:) = pack(itemp, plmask)
          tp%id(:) = pack(itemp, tpmask)
          cb%id = 0
-         pl%id(:) = pack([(i,i=1,idmax)],plmask)
-         tp%id(:) = pack([(i,i=1,idmax)],tpmask)
+         pl%id(:) = pack([(i,i=0,idmax-1)],plmask)
+         tp%id(:) = pack([(i,i=0,idmax-1)],tpmask)
 
          call check( nf90_get_var(nc%id, nc%name_varid, ctemp, count=[NAMELEN, idmax]), "netcdf_read_particle_info_system nf90_getvar name_varid"  )
          call cb%info%set_value(name=ctemp(1))