Merge branch 'master' into flux-tubes

sources/io.F90

@@ -4029,12 +4029,11 @@ module io
       allocate(fields(16,nm), children(nc,nm), bounds(3,2,nm), &
 #if NDIMS == 3
                faces(NDIMS,ns,ns,ns,nm), &
-               edges(NDIMS,ns,ns,ns,nm), corners(nm,ns,ns,ns), &
+               edges(NDIMS,ns,ns,ns,nm), corners(ns,ns,ns,nm), &
 #else /* NDIMS == 3 */
                edges(NDIMS,ns,ns,nm), corners(ns,ns,nm), &
 #endif /* NDIMS == 3 */
                stat = status)
-!                block_array(get_last_id()), stat = status)
 
       if (status == 0) then