From 999f9e737c21f68b2957ee3be00b2ac93fa1baca Mon Sep 17 00:00:00 2001
From: James McClure <mcclurej@vt.edu>
Date: Thu, 5 Jun 2014 15:03:25 -0400
Subject: [PATCH] fixed porosity bug

---
 tests/lb2_Color_wia_mpi.cpp | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/tests/lb2_Color_wia_mpi.cpp b/tests/lb2_Color_wia_mpi.cpp
index 12b50ec1..b6a423fd 100644
--- a/tests/lb2_Color_wia_mpi.cpp
+++ b/tests/lb2_Color_wia_mpi.cpp
@@ -452,7 +452,7 @@ int main(int argc, char **argv)
 					id[n] = 2;	
 				}
 				// compute the porosity (actual interface location used)
-				if (SignDist.data[n] > -1.0){ 
+				if (SignDist.data[n] > 1.0){ 
 					sum++;	
 				}
 			}
@@ -499,6 +499,7 @@ int main(int argc, char **argv)
 		phi_s = -phi_s;
 	 	das = (phi_s+1.0)*0.5;
 		dbs = 1.0 - das;
+		printf("Resetting phi_s = %f, das = %f, dbs = %f \n", phi_s, das, dbs);
 		FlipID(id,Nx*Ny*Nz);
 	}
 #endif