porosity factor in effperm

This commit is contained in:
James McClure 2022-08-12 18:34:50 -04:00
parent d50c98c71a
commit 4a016eee6c

View File

@ -411,7 +411,7 @@ void SubPhase::Basic() {
dir_z = 1.0;
force_mag = 1.0;
}
double Porosity = (gwb.V + gnb.V)/Dm->Volume();
double Porosity = (gwb.V + gnb.V)/Dm->Volume;
double saturation = gwb.V / (gwb.V + gnb.V);
double water_flow_rate =
gwb.V * (gwb.Px * dir_x + gwb.Py * dir_y + gwb.Pz * dir_z) / gwb.M /