From 8695744b644e707f5348311cada2402b261ae742 Mon Sep 17 00:00:00 2001
From: James E McClure <jemcclur@gmail.com>
Date: Wed, 21 Feb 2018 20:39:27 -0500
Subject: [PATCH] debugging perm sim

---
 tests/lbpm_permeability_simulator.cpp | 10 +++++-----
 1 file changed, 5 insertions(+), 5 deletions(-)

diff --git a/tests/lbpm_permeability_simulator.cpp b/tests/lbpm_permeability_simulator.cpp
index 0c22d9de..5591e3dc 100644
--- a/tests/lbpm_permeability_simulator.cpp
+++ b/tests/lbpm_permeability_simulator.cpp
@@ -186,7 +186,7 @@ int main(int argc, char **argv)
 		id = new char[N];
 		double sum, sum_local;
 		double iVol_global = 1.0/(1.0*(Nx-2)*(Ny-2)*(Nz-2)*nprocs);
-		if (BoundaryCondition > 0) iVol_global = 1.0/(1.0*(Nx-2)*nprocx*(Ny-2)*nprocy*((Nz-2)*nprocz-6));
+		//if (BoundaryCondition > 0) iVol_global = 1.0/(1.0*(Nx-2)*nprocx*(Ny-2)*nprocy*((Nz-2)*nprocz-6));
 		double porosity, pore_vol;
 		//...........................................................................
 		if (rank == 0) cout << "Reading in domain from signed distance function..." << endl;
@@ -195,7 +195,7 @@ int main(int argc, char **argv)
 		// Read the signed distance
 		sprintf(LocalRankString,"%05d",rank);
 		sprintf(LocalRankFilename,"%s%s","SignDist.",LocalRankString);
-		ReadBinaryFile(LocalRankFilename, Averages->SDs.data(), N);
+		ReadBinaryFile(LocalRankFilename, Averages.SDs.data(), N);
 		MPI_Barrier(comm);
 		if (rank == 0) cout << "Domain set." << endl;
 
@@ -217,11 +217,11 @@ int main(int argc, char **argv)
 			for ( j=0;j<Ny;j++){
 				for ( i=0;i<Nx;i++){
 					int n = k*Nx*Ny+j*Nx+i;
-					if (Averages->SDs(n) > 0.0){
+					if (Averages.SDs(n) > 0.0){
 						id[n] = 2;	
 					}
 					// compute the porosity (actual interface location used)
-					if (Averages->SDs(n) > 0.0){
+					if (Averages.SDs(n) > 0.0){
 						sum++;	
 					}
 				}
@@ -241,7 +241,7 @@ int main(int argc, char **argv)
 		// Compute the media porosity, assign phase labels and solid composition
 		//.......................................................................
 		sum_local=0.0;
-		Np=0;  // number of local pore nodes
+		int Np=0;  // number of local pore nodes
 		//.......................................................................
 		for (k=1;k<Nz-1;k++){
 			for (j=1;j<Ny-1;j++){