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

Fixed bug in pressure boundary conditions lbpm_color_simulator

Browse Source
This commit is contained in:
James E McClure
2015-08-13 02:23:27 -04:00
parent f580de19c2
commit f978ce4eb7
1 changed files with 0 additions and 15 deletions
Download Patch File Download Diff File Expand all files Collapse all files

15
tests/lbpm_color_simulator.cpp
View File

@@ -375,21 +375,6 @@ int main(int argc, char **argv)
fread(id,1,N,IDFILE);
fclose(IDFILE);
for ( k=0;k<Nz;k++){
for ( j=0;j<Ny;j++){
for ( i=0;i<Nx;i++){
n = k*Nx*Ny+j*Nx+i;
// The following turns off communication if external BC are being set
if (BoundaryCondition > 0){
if (Dm.kproc==0 && k==0) id[n]=0;
if (Dm.kproc==0 && k==1) id[n]=0;
if (Dm.kproc==nprocz-1 && k==Nz-2) id[n]=0;
if (Dm.kproc==nprocz-1 && k==Nz-1) id[n]=0;
}
}
}
}
// Set up kstart, kfinish so that the reservoirs are excluded from averaging
int kstart,kfinish;
kstart = 1;