OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

adjustment to erodelabel in morph open

Browse Source
This commit is contained in:
James E McClure 2019-04-27 10:50:24 -04:00
parent d93fe1f176
commit 5dba2fe702
1 changed files with 3 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
analysis/morphology.cpp
View File

@ -449,6 +449,8 @@ double MorphDrain(DoubleArray &SignDist, signed char *id, std::shared_ptr<Domain
// if (rank==0) printf("Max. distance =%f, Initial critical radius = %f \n",maxdistGlobal,Rcrit_new);
//}
MPI_Barrier(Dm->Comm);
if (ErodeLabel == 1) VoidFraction = 1-VoidFraction;
while (void_fraction_new > VoidFraction && Rcrit_new > 0.5)
{
@ -641,7 +643,7 @@ double MorphDrain(DoubleArray &SignDist, signed char *id, std::shared_ptr<Domain
for (int j=1; j<ny-1; j++){
for (int i=1; i<nx-1; i++){
n=k*nx*ny+j*nx+i;
if (id[n] > 1){
if (id[n] == ErodeLabel){
count+=1.0;
}
}