diff --git a/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-imp.h b/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-imp.h
index f4f548601c37fccdba77eb12f0b0f0b0b1f801df..0b393f1e00518e861153af514cb2b70ea4fb5eac 100644
--- a/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-imp.h
+++ b/ProcessLib/LiquidFlow/LiquidFlowLocalAssembler-imp.h
@@ -47,6 +47,10 @@ void LiquidFlowLocalAssembler<ShapeFunction, IntegrationMethod, GlobalDim>::
     const Eigen::MatrixXd& perm =
         _material_properties.intrinsic_permeabiliy[mat_id];
 
+    // Note: For Inclined 1D in 2D/3D or 2D element in 3D, the first item in
+    //  the assert must be changed to perm.rows() == _element->getDimension()
+    assert(perm.rows() == GlobalDim || perm.rows() == 1);
+
     // TODO: The following two variables should be calculated inside the
     //       the integration loop for non-constant porosity and storage models.
     double porosity_variable = 0.;
diff --git a/Tests/Data b/Tests/Data
index a4b3ec0e41dbc4365025f98c1810a9986f3c1f28..80398320e2bc58005283f31b32dba67cd748e903 160000
--- a/Tests/Data
+++ b/Tests/Data
@@ -1 +1 @@
-Subproject commit a4b3ec0e41dbc4365025f98c1810a9986f3c1f28
+Subproject commit 80398320e2bc58005283f31b32dba67cd748e903