...
  View open merge request
Commits (1)
......@@ -80,11 +80,23 @@ public:
//! Function for defining the (intrinsic) permeability \f$[m^2]\f$.
PermeabilityType permeabilityAtPos(const GlobalPosition& globalPos) const
{ return globalPos[1] > 35.0 ? permeabilityOverburden_ : permeability_; }
{
if (globalPos[1] > 35.0)
return permeabilityOverburden_;
else if ((globalPos[1] > 32)||(globalPos[1] < 30))
return permeability_;
else return permeability_/10000.;
}
//! Return the porosity
Scalar porosityAtPos(const GlobalPosition& globalPos) const
{ return globalPos[1] > 35.0 ? porosityOverburden_ : porosity_; }
{
if (globalPos[1] > 35.0)
return porosityOverburden_;
else if ((globalPos[1] > 32)||(globalPos[1] < 30))
return porosity_;
else return 0.05;
}
//! Return the material law parameters
const MaterialLawParams& materialLawParamsAtPos(const GlobalPosition& globalPos) const
......