From 3807039a602d02b3dbc3999d3a6d316839bb6ba1 Mon Sep 17 00:00:00 2001
From: James E McClure <mcclurej@vt.edu>
Date: Sun, 9 Jul 2017 18:56:52 -0400
Subject: [PATCH] Working on morphological initialization when wide scale
 variability is present

---
 tests/lbpm_morphdrain_pp.cpp |  2 +-
 tests/lbpm_morphopen_pp.cpp  | 11 +++++------
 2 files changed, 6 insertions(+), 7 deletions(-)

diff --git a/tests/lbpm_morphdrain_pp.cpp b/tests/lbpm_morphdrain_pp.cpp
index acc1e325..45b5f53e 100644
--- a/tests/lbpm_morphdrain_pp.cpp
+++ b/tests/lbpm_morphdrain_pp.cpp
@@ -288,7 +288,7 @@ int main(int argc, char **argv)
 	while (sw_new > SW && Rcrit_new > 0.99){
         
         Rcrit_old = Rcrit_new;
-		Rcrit_new -= deltaR;// decrease critical radius
+		Rcrit_new -= deltaR*Rcrit_old;// decrease critical radius
         sw_old = sw_new;
         sw_diff_old = sw_diff_new;
 
diff --git a/tests/lbpm_morphopen_pp.cpp b/tests/lbpm_morphopen_pp.cpp
index fb77a550..0c88176b 100644
--- a/tests/lbpm_morphopen_pp.cpp
+++ b/tests/lbpm_morphopen_pp.cpp
@@ -352,7 +352,7 @@ int main(int argc, char **argv)
         sw_diff_old = sw_diff_new;
         sw_old = sw_new;
         Rcrit_old = Rcrit_new;
-		Rcrit_new -= deltaR;
+		Rcrit_new -= deltaR*Rcrit_old;
 	    int Window=round(Rcrit_new);
         if (Window == 0) Window = 1; // If Window = 0 at the begining, after the following process will have sw=1.0
                                      // and sw<Sw will be immediately broken
@@ -483,11 +483,10 @@ int main(int argc, char **argv)
         sw_new = countGlobal/totalGlobal;
         sw_diff_new = abs(sw_new-SW);
         // for test only
-//        if (rank==0){
-//            printf("Final saturation=%f\n",sw_new);
-//            printf("Final critical radius=%f\n",Rcrit_new);
-//
-//        }
+        if (rank==0){
+            printf("Final saturation=%f\n",sw_new);
+            printf("Final critical radius=%f\n",Rcrit_new);
+        }
     }
 
     if (sw_diff_new<sw_diff_old){