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

check header for Peremeability.csv

Browse Source
This commit is contained in:
JamesEMcclure
2020-01-13 11:36:29 -05:00
parent 62674e4e0c
commit b69b505b2f
1 changed files with 13 additions and 3 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
models/MRTModel.cpp
View File

@@ -194,11 +194,21 @@ void ScaLBL_MRTModel::Run(){
int SIZE=Np*sizeof(double);
if (rank==0){
FILE * log_file = fopen("Permeability.csv","a");
fprintf(log_file,"time Fx Fy Fz mu Vs As Js Xs vx vy vz k\n");
fclose(log_file);
bool WriteHeader=false;
log_file = fopen("Permeability.csv","r");
if (log_file != NULL)
fclose(log_file);
else
WriteHeader=true;
if (WriteHeader){
FILE * log_file = fopen("Permeability.csv","a");
fprintf(log_file,"time Fx Fy Fz mu Vs As Js Xs vx vy vz k\n");
fclose(log_file);
}
}
//.......create and start timer............
double starttime,stoptime,cputime;
ScaLBL_DeviceBarrier(); MPI_Barrier(comm);