From e0c6f7d33edf94763011f0ea41842385c87143ea Mon Sep 17 00:00:00 2001
From: James E McClure <jemcclur@gmail.com>
Date: Wed, 25 Jun 2014 12:58:48 -0400
Subject: [PATCH] updated lb2_Blob_wia_mpi

---
 tests/lb2_Blob_wia_mpi.cpp | 13 +++++++++++--
 1 file changed, 11 insertions(+), 2 deletions(-)

diff --git a/tests/lb2_Blob_wia_mpi.cpp b/tests/lb2_Blob_wia_mpi.cpp
index 6895abf5..abfd08f3 100644
--- a/tests/lb2_Blob_wia_mpi.cpp
+++ b/tests/lb2_Blob_wia_mpi.cpp
@@ -1322,6 +1322,7 @@ int main(int argc, char **argv)
 			}
 		}
 	}
+	int nblobs = 0;
 	IntArray BlobList(3,ncubes);	// store indices for blobs (cubes)
 	IntArray temp(3,ncubes);		// temporary storage array
 	IntArray BlobSizes(MAX_LOCAL_BLOB_COUNT);			// number of nodes in each blob
@@ -2308,13 +2309,21 @@ int main(int argc, char **argv)
 			Jwn = Kwn = efawns = 0.0;
 			trJwn = trawn = trRwn = 0.0;
 			
+			for (k=0;k<Nz;k++){
+				for (j=0;j<Ny;j++){
+					for (i=0;i<Nx;i++){
+						LocalBlobID(i,j,k) = -1;	// Initialize each time
+					}
+				}
+			}
+			
 			nblobs=0;
 			double vF,vS;
 			vF = vS = 0.0;
 			for (k=0;k<Nz;k++){
 				for (j=0;j<Ny;j++){
-					for (i=1;i<Nx;i++){
-						if ( indicator(i,j,k) == -1 ){
+					for (i=0;i<Nx;i++){
+						if ( LocalBlobID(i,j,k) == -1 ){
 							if ( Phase(i,j,k) > 0.0 ){
 								if ( SignDist(i,j,k) > 0.0 ){
 									// node i,j,k is in the porespace