regions retain mutual exclusivity

This commit is contained in:
James E McClure 2019-04-24 14:18:45 -04:00
parent 0169f68999
commit 4e4d5e3a8c

View File

@ -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;