From 4324e19dd5b7e9ebf5f1bf907ed5d03295ee8b5c Mon Sep 17 00:00:00 2001
From: James E McClure <jemcclur@gmail.com>
Date: Fri, 22 Mar 2019 16:15:29 -0400
Subject: [PATCH] radius exit criteria for morphdrain

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

diff --git a/analysis/morphology.cpp b/analysis/morphology.cpp
index e7bac352..af701770 100644
--- a/analysis/morphology.cpp
+++ b/analysis/morphology.cpp
@@ -450,8 +450,9 @@ double MorphDrain(DoubleArray &SignDist, char *id, std::shared_ptr<Domain> Dm, d
 	//	Rcrit_new = strtod(argv[2],NULL);
 	//	if (rank==0) printf("Max. distance =%f, Initial critical radius = %f \n",maxdistGlobal,Rcrit_new);
 	//}
-	
-	while (void_fraction_new > VoidFraction)
+	MPI_Barrier(Dm->Comm);
+
+	while (void_fraction_new > VoidFraction && Rcrit_new > 0.5)
 	{
 		void_fraction_diff_old = void_fraction_diff_new;
 		void_fraction_old = void_fraction_new;