refactor morphopen

This commit is contained in:
James E McClure
2018-06-18 13:39:05 -04:00
parent e85f43ae5c
commit 89da203c2c

View File

@@ -74,7 +74,7 @@ int main(int argc, char **argv)
filename=argv[1];
Rcrit_new=0.f;
SW=strtod(argv[2],NULL);
printf("Target saturation %f \n",SW);
if (rank==0) printf("Target saturation %f \n",SW);
}
else ERROR("No input database provided\n");
// read the input database
@@ -133,6 +133,7 @@ int main(int argc, char **argv)
}
}
}
if (maxdist>50.f) maxdist=50.f;
// total Global is the number of nodes in the pore-space
MPI_Allreduce(&count,&totalGlobal,1,MPI_DOUBLE,MPI_SUM,comm);
MPI_Allreduce(&maxdist,&maxdistGlobal,1,MPI_DOUBLE,MPI_MAX,comm);
@@ -141,7 +142,6 @@ int main(int argc, char **argv)
if (rank==0) printf("Media Porosity: %f \n",porosity);
if (rank==0) printf("Maximum pore size: %f \n",maxdistGlobal);\
Dm->CommInit();
int iproc = Dm->iproc();
int jproc = Dm->jproc();
@@ -383,8 +383,6 @@ int main(int argc, char **argv)
}
sprintf(LocalRankFilename,"ID.%05i",rank);
FILE *ID = fopen(LocalRankFilename,"wb");
fwrite(id,1,N,ID);