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

fixing morphopen

Browse Source
This commit is contained in:
James E McClure 2018-11-02 14:33:48 -04:00
parent f9d5c2965a
commit d7ada6f10c
1 changed files with 3 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

13
tests/lbpm_morphopen_pp.cpp
View File

@ -165,7 +165,9 @@ int main(int argc, char **argv)
for (int j=0; j<ny; j++){ for (int j=0; j<ny; j++){
for (int i=0; i<nx; i++){ for (int i=0; i<nx; i++){
n = k*nx*ny+j*nx+i; n = k*nx*ny+j*nx+i;
if (SignDist(i,j,k) < 0.0) id[n] = 0; if (SignDist(i,j,k) < 0.f){
// don't do anything
}
else{ else{
// initially saturated with wetting phase // initially saturated with wetting phase
id[n] = 2; id[n] = 2;
@ -423,15 +425,6 @@ int main(int argc, char **argv)
if (rank==0) printf("Writing ID file \n"); if (rank==0) printf("Writing ID file \n");
sprintf(LocalRankFilename,"ID.%05i",rank); sprintf(LocalRankFilename,"ID.%05i",rank);
// Preserve mineral labels
for (int k=0; k<nz; k++){
for (int j=0; j<ny; j++){
for (int i=0; i<nx; i++){
n = k*nx*ny+j*nx+i;
if (SignDist(i,j,k) < 0.0) id[n] = Dm->id[n];
}
}
}
FILE *ID = fopen(LocalRankFilename,"wb"); FILE *ID = fopen(LocalRankFilename,"wb");
fwrite(id,1,N,ID); fwrite(id,1,N,ID);
fclose(ID); fclose(ID);