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

added permeability log file

Browse Source
This commit is contained in:
James McClure 2018-11-24 14:02:31 -05:00
parent 3b6a7632ff
commit 61a580f44d
1 changed files with 13 additions and 6 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

19
models/MRTModel.cpp
View File

@ -146,8 +146,12 @@ void ScaLBL_MRTModel::Run(){
Minkowski Morphology(Mask); Minkowski Morphology(Mask);
int SIZE=Np*sizeof(double); int SIZE=Np*sizeof(double);
if (rank==0) printf("time Fx Fy Fz mu Vs As Js Xs vx vy vz\n"); 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\n");
fclose(log_file);
}
//.......create and start timer............ //.......create and start timer............
double starttime,stoptime,cputime; double starttime,stoptime,cputime;
ScaLBL_DeviceBarrier(); MPI_Barrier(comm); ScaLBL_DeviceBarrier(); MPI_Barrier(comm);
@ -206,11 +210,14 @@ void ScaLBL_MRTModel::Run(){
//if (rank==0) printf("Computing Minkowski functionals \n"); //if (rank==0) printf("Computing Minkowski functionals \n");
Morphology.ComputeScalar(Distance,0.f); Morphology.ComputeScalar(Distance,0.f);
Morphology.PrintAll(); //Morphology.PrintAll();
double mu = (tau-0.5)/3.f; double mu = (tau-0.5)/3.f;
if (rank==0) printf("%i %.8g %.8g %.8g %.8g %.8g %.8g %.8g %.8g %.8g %.8g %.8g\n",timestep, Fx, Fy, Fz, mu, if (rank==0) {
Morphology.V(),Morphology.A(),Morphology.J(),Morphology.X(),vax,vay,vaz); FILE * log_file = fopen("Permeability.csv","a");
fprintf(log_file,"%i %.8g %.8g %.8g %.8g %.8g %.8g %.8g %.8g %.8g %.8g %.8g\n",timestep, Fx, Fy, Fz, mu,
Morphology.V(),Morphology.A(),Morphology.J(),Morphology.X(),vax,vay,vaz);
fclose(log_file);
}
} }
} }
//************************************************************************/ //************************************************************************/