diff --git a/tests/lbpm_uCT_pp.cpp b/tests/lbpm_uCT_pp.cpp
new file mode 100644
index 00000000..cf8e55aa
--- /dev/null
+++ b/tests/lbpm_uCT_pp.cpp
@@ -0,0 +1,57 @@
+// Sequential blob analysis 
+// Reads parallel simulation data and performs connectivity analysis
+// and averaging on a blob-by-blob basis
+// James E. McClure 2014
+
+#include "IO/netcdf.h"
+
+#include "ProfilerApp.h"
+
+
+void load( const std::string filename )
+{
+    int fid = netcdf::open( filename );
+
+    std::vector<std::string> vars = netcdf::getVarNames( fid );
+    for (size_t i=0; i<vars.size(); i++) {
+        printf("Reading variable %s\n",vars[i].c_str());
+        netcdf::VariableType type = netcdf::getVarType( fid, vars[i] );
+        if ( type == netcdf::STRING )
+            Array<std::string> tmp = netcdf::getVar<std::string>( fid, vars[i] );
+        else if ( type == netcdf::SHORT )
+            Array<short> tmp = netcdf::getVar<short>( fid, vars[i] );
+        else
+            Array<double> tmp = netcdf::getVar<double>( fid, vars[i] );
+    }
+
+    std::vector<std::string> attr = netcdf::getAttNames( fid );
+    for (size_t i=0; i<attr.size(); i++) {
+        printf("Reading attribute %s\n",attr[i].c_str());
+        netcdf::VariableType type = netcdf::getAttType( fid, attr[i] );
+        if ( type == netcdf::STRING )
+            Array<std::string> tmp = netcdf::getAtt<std::string>( fid, attr[i] );
+        else
+            Array<double> tmp = netcdf::getAtt<double>( fid, attr[i] );
+    }
+    netcdf::close( fid );
+}
+
+
+int main(int argc, char **argv)
+{
+    PROFILE_START("Main");
+
+    std::vector<std::string> filenames;
+    
+    if ( argc==0 ) {
+        printf("At least one filename must be specified\n");
+        return 1;
+    }
+
+    for (int i=1; i<argc; i++)
+        load( argv[i] );
+
+    PROFILE_SAVE("TestNetcdf");
+    return 0;  
+}
+