  
  	  // set up the comment
  	  if (libMesh::processor_id() == 0)
-@@ -1632,6 +1692,43 @@
+@@ -1632,6 +1692,53 @@
  	    }
  	}
      }

+    io_blksize = std::min(max_io_blksize, n_objects),
+    num_blks   = std::ceil(static_cast<double>(n_objects)/static_cast<double>(io_blksize));
+
+  std::cout << "io_blksize = "    << io_blksize
+	    << ", num_objects = " << n_objects
+	    << ", num_blks = "    << num_blks
+	    << std::endl;
+  
+  unsigned int written_length=0;                       // The numer of values written.  This will be returned
+

 {
   parallel_only();
 
@@ -1900,8 +2307,58 @@
 
   if (libMesh::processor_id() == 0)
     libmesh_assert_equal_to (written_length, vec_length);
