From 1dc102414deb1c1f49676be5f5863b9fe3729382 Mon Sep 17 00:00:00 2001
From: James E McClure <jemcclur@gmail.com>
Date: Wed, 9 Jan 2019 11:08:24 -0500
Subject: [PATCH] prepare to skip 2nd distance map in morph

---
 models/ColorModel.cpp | 7 +++++--
 1 file changed, 5 insertions(+), 2 deletions(-)

diff --git a/models/ColorModel.cpp b/models/ColorModel.cpp
index f17e423b..fb38f81e 100644
--- a/models/ColorModel.cpp
+++ b/models/ColorModel.cpp
@@ -817,11 +817,13 @@ double ScaLBL_ColorModel::MorphInit(const double beta, const double target_delta
 	for (int k=0; k<Nz; k++){
 		for (int j=0; j<Ny; j++){
 			for (int i=0; i<Nx; i++){
-				if (phase_distance(i,j,k) < 0.0 ) phase_id(i,j,k) = 0;
-				else 		     				  phase_id(i,j,k) = 1;
+				//if (phase_distance(i,j,k) < 0.0 ) phase_id(i,j,k) = 0;
+				//else 		     				  phase_id(i,j,k) = 1;
+				if (phase_distance(i,j,k) < 0.0 ) phase(i,j,k) = 1.0;
 			}
 		}
 	}	
+	/*
 	CalcDist(phase_distance,phase_id,*Dm); // re-calculate distance
 
 	// 5. Update phase indicator field based on new distnace
@@ -839,6 +841,7 @@ double ScaLBL_ColorModel::MorphInit(const double beta, const double target_delta
 			} 
 		}
 	}
+	*/
 	fillDouble.fill(phase);
 
 	count = 0.f;