From e7987c9cc9f9cede608ddd49123e8d2cae60b827 Mon Sep 17 00:00:00 2001
From: rexzheli <rexzheli@thor.sv.vt.edu>
Date: Fri, 2 Sep 2016 14:38:34 -0400
Subject: [PATCH] Fixed the setup for the square tube test

---
 tests/lbpm_squaretube_pp.cpp | 5 +++--
 1 file changed, 3 insertions(+), 2 deletions(-)

diff --git a/tests/lbpm_squaretube_pp.cpp b/tests/lbpm_squaretube_pp.cpp
index 02c4ae76..672f894b 100644
--- a/tests/lbpm_squaretube_pp.cpp
+++ b/tests/lbpm_squaretube_pp.cpp
@@ -149,13 +149,14 @@ int main(int argc, char **argv)
 	double pore_vol;
 
 	sum=0;
+        //double wall_width = (Nx - TubeWidth)/2;
 	for (k=0;k<Nz;k++){
 		for (j=0;j<Ny;j++){
 			for (i=0;i<Nx;i++){
 				n = k*Nx*Ny + j*Nz + i;
 				// square capillary tube aligned with the z direction
-				Averages.SDs(i,j,k) = TubeWidth - fabs(i-0.5*Nx);
-				Averages.SDs(i,j,k) = min(Averages.SDs(i,j,k),TubeWidth-fabs(j-0.5*Ny));
+				Averages.SDs(i,j,k) = TubeWidth/2 - fabs(i-0.5*Nx);
+				Averages.SDs(i,j,k) = min(Averages.SDs(i,j,k),TubeWidth/2-fabs(j-0.5*Ny));
 
 				// Initialize phase positions
 				if (Averages.SDs(i,j,k) < 0.0){