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

Fix to lbpm_nondarcy_simulator.cpp

Browse Source
This commit is contained in:
James E McClure
2016-11-22 13:40:25 -05:00
parent cbec621d7a
commit d486289a30
1 changed files with 5 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

11
tests/lbpm_nondarcy_simulator.cpp
View File

@@ -393,8 +393,7 @@ int main(int argc, char **argv)
starttime = MPI_Wtime();
//.........................................
double D32,Fo,Re,velocity,err1D,mag_force,vel_prev;
double D32,vawx,vawy,vawz,Fo,Re,velocity,err1D,mag_force,vel_prev;
FILE * NONDARCY;
if (rank == 0){
NONDARCY = fopen("nondarcy.csv","a");
@@ -470,9 +469,9 @@ int main(int argc, char **argv)
Averages.ComputeLocal();
Averages.Reduce();
double vawx = -Averages.vaw_global(0);
double vawy = -Averages.vaw_global(1);
double vawz = -Averages.vaw_global(2);
vawx = -Averages.vaw_global(0);
vawy = -Averages.vaw_global(1);
vawz = -Averages.vaw_global(2);
// Compute local measures
err = Re; // previous Reynolds number
@@ -509,7 +508,7 @@ int main(int argc, char **argv)
// Write steady state variables to csv file
if (rank==0){
fprintf(NONDARCY,"%.5g %.5g %.5g %.5g %.5g %.5g %.5g %.5g %.5g\n",D32,Fx,Fy,Fz,vx,vy,vz,Re,Fo);
fprintf(NONDARCY,"%.5g %.5g %.5g %.5g %.5g %.5g %.5g %.5g %.5g\n",D32,Fx,Fy,Fz,vawx,vawy,vawz,Re,Fo);
fflush(NONDARCY);
}
}