Updates to TwoPhase:AssignComponentLabels

2015-08-30 08:22:05 -04:00
parent de51a11dfd
commit e3961007eb
1 changed files with 3 additions and 0 deletions
--- a/common/TwoPhase.cpp
+++ b/common/TwoPhase.cpp
@@ -464,6 +464,9 @@ void TwoPhase::AssignComponentLabels()
 		for (int j=0; j<Ny; j++){
 			for (int i=0; i<Nx; i++){
 				Label_WP(i,j,k) = 0;
+				if (SDs(i,j,k) > 0.0) PhaseID(i,j,k) = 1;
+				else if (Phase(i,j,k) > 0.0) PhaseID(i,j,k) = 1;
+				else PhaseID(i,j,l) = 2;
 			}
 		}
 	}