introduced a quasi-impermeable layer in the shale

I have discussed with Dennis and we think this is an improvement which is worth being part of the release

lecture/mm/fractures/matrixspatialparams.hh

@@ -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