From 23d0e09956e2032d54be92a327c231f17814b055 Mon Sep 17 00:00:00 2001 From: James McClure Date: Sat, 24 Nov 2018 15:40:44 -0500 Subject: [PATCH] added absperm --- models/MRTModel.cpp | 3 +++ 1 file changed, 3 insertions(+) diff --git a/models/MRTModel.cpp b/models/MRTModel.cpp index d0d61957..6bab022b 100644 --- a/models/MRTModel.cpp +++ b/models/MRTModel.cpp @@ -213,6 +213,9 @@ void ScaLBL_MRTModel::Run(){ //Morphology.PrintAll(); double mu = (tau-0.5)/3.f; if (rank==0) { + double h = Lz/double(Nz); + double absperm = h*h*mu*sqrt(vax*vax+vay*vay+vaz*vaz)/sqrt(Fx*Fx+Fy*Fy+Fz*Fz); + printf(" %f\n",absperm); 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);