From 6ed5f5aea6e4fd714b09a9dd207a3aebfeb76df1 Mon Sep 17 00:00:00 2001
From: zherexli <zhe.rex.li@gmail.com>
Date: Wed, 2 Nov 2016 17:39:49 -0400
Subject: [PATCH] fix how the W phase saturation is calculated in morphological
 drainage

---
 tests/lbpm_morphdrain_pp.cpp | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/tests/lbpm_morphdrain_pp.cpp b/tests/lbpm_morphdrain_pp.cpp
index 11d484d8..2c2b5e25 100644
--- a/tests/lbpm_morphdrain_pp.cpp
+++ b/tests/lbpm_morphdrain_pp.cpp
@@ -409,7 +409,7 @@ int main(int argc, char **argv)
 		for (int j=1; j<Ny-1; j++){
 			for (int i=1; i<Nx-1; i++){
 				n=k*Nx*Ny+j*Nx+i;
-				if (id[n] == 1){
+				if (id[n] == 2){
 					count++;
 				}
 			}