diff --git a/analysis/SubPhase.cpp b/analysis/SubPhase.cpp
index 4cbc5337..9c2f4529 100644
--- a/analysis/SubPhase.cpp
+++ b/analysis/SubPhase.cpp
@@ -514,7 +514,7 @@ void SubPhase::Full(){
 						double nz = 0.5*(Phi(i,j,k+1)-Phi(i,j,k-1));
 						InterfaceTransportMeasures(  beta,  rho_w,  rho_n,  nA, nB, nx, ny, nz, ux, uy, uz, iwn);
 					}
-					if ( phi > 0.0){
+					else if ( phi > 0.0){
 						if (morph_n->label(i,j,k) > 0 ){
 							vol_nd_bulk += 1.0;
 							nd.M += nA*rho_n;