diff --git a/ProcessLib/LiquidFlow/Tests.cmake b/ProcessLib/LiquidFlow/Tests.cmake
index 9feb5c6c64991252b9d0d999d7da1ec1a0a6508c..3938c2921918268e14a201350d400f0c50293bd5 100644
--- a/ProcessLib/LiquidFlow/Tests.cmake
+++ b/ProcessLib/LiquidFlow/Tests.cmake
@@ -560,7 +560,7 @@ AddTest(
REQUIREMENTS NOT OGS_USE_MPI
DIFF_DATA
transient_flow_in_inclined_2D_plane_ts_t_864000.000000.vtu transient_flow_in_inclined_2D_plane_ts_t_864000.000000.vtu pressure pressure 1e-14 1e-11
- transient_flow_in_inclined_2D_plane_ts_t_864000.000000.vtu transient_flow_in_inclined_2D_plane_ts_t_864000.000000.vtu v v 1e-14 1e-14
+ transient_flow_in_inclined_2D_plane_ts_t_864000.000000.vtu transient_flow_in_inclined_2D_plane_ts_t_864000.000000.vtu v v 1e-11 1e-14
)
AddTest(
diff --git a/Tests/Data/Parabolic/LiquidFlow/InclinedMeshElements/Inclined2DMesh/hydrostatic_flow_in_inclined_2D_plane.prj b/Tests/Data/Parabolic/LiquidFlow/InclinedMeshElements/Inclined2DMesh/hydrostatic_flow_in_inclined_2D_plane.prj
index 77184ce0bdd113eaa9df51a87ad8d31cc6208df4..805ba2fd59bc6e67fe15eef911ef718f4e072a38 100644
--- a/Tests/Data/Parabolic/LiquidFlow/InclinedMeshElements/Inclined2DMesh/hydrostatic_flow_in_inclined_2D_plane.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/InclinedMeshElements/Inclined2DMesh/hydrostatic_flow_in_inclined_2D_plane.prj
@@ -7,7 +7,7 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <process_variables>
+ <process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<specific_body_force>0 0 -9.8</specific_body_force>
@@ -104,7 +104,18 @@
<parameter>
<name>k</name>
<type>Constant</type>
- <values>1.e-12 0.0 0.0 5.e-12 </values>
+ <!--Permeability tensor in the local coordinate system on the
+ inclined plane. Its compononts associated with the off plane
+ direction have to be set to zero. In other words, the local
+ permeability tensor has to input as a diagonal block of the full
+ size tensor in the global coodinate system. The bases of the
+ local coordinate system has to be defined in order to rotate
+ this local tensor to the global coodinate system.-->
+ <values>1.e-12 0.0 0.0
+ 0.0 5.e-12 0.0
+ 0.0 0.0 0.0</values>
+ <use_local_coordinate_system>true
+ </use_local_coordinate_system>
</parameter>
<parameter>
<name>p0</name>
@@ -116,7 +127,28 @@
<type>Constant</type>
<value>1.0e+5</value>
</parameter>
+ <parameter>
+ <name>b0</name>
+ <type>Constant</type>
+ <values> -0.5000001748438417 0.0 0.8660253028382761</values>
+ </parameter>
+ <parameter>
+ <name>b1</name>
+ <type>Constant</type>
+ <values>0.0 -1.0 0.0</values>
+ </parameter>
+ <parameter>
+ <name>b2</name>
+ <!-- off plane direction -->
+ <type>Constant</type>
+ <values>0.8660253028382761 0.0 0.5000001748438417</values>
+ </parameter>
</parameters>
+ <local_coordinate_system>
+ <basis_vector_0>b0</basis_vector_0>
+ <basis_vector_1>b1</basis_vector_1>
+ <basis_vector_2>b2</basis_vector_2>
+ </local_coordinate_system>
<process_variables>
<process_variable>
<name>pressure</name>
@@ -126,7 +158,7 @@
<boundary_conditions>
<boundary_condition>
<type>Dirichlet</type>
- <mesh>geometry_top</mesh>
+ <mesh>geometry_top</mesh>
<parameter>p_top</parameter>
</boundary_condition>
</boundary_conditions>
diff --git a/Tests/Data/Parabolic/LiquidFlow/InclinedMeshElements/Inclined2DMesh/transient_flow_in_inclined_2D_plane.prj b/Tests/Data/Parabolic/LiquidFlow/InclinedMeshElements/Inclined2DMesh/transient_flow_in_inclined_2D_plane.prj
index 94061785920f8e9eff45aa08b012be5bc429426f..1ad729208705ef5b0b6c3ff00225240d0c34f55b 100644
--- a/Tests/Data/Parabolic/LiquidFlow/InclinedMeshElements/Inclined2DMesh/transient_flow_in_inclined_2D_plane.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/InclinedMeshElements/Inclined2DMesh/transient_flow_in_inclined_2D_plane.prj
@@ -7,7 +7,7 @@
<name>LiquidFlow</name>
<type>LIQUID_FLOW</type>
<integration_order>2</integration_order>
- <process_variables>
+ <process_variables>
<process_variable>pressure</process_variable>
</process_variables>
<specific_body_force>0 0 -9.8</specific_body_force>
@@ -33,11 +33,11 @@
<reference_value>1000.0</reference_value>
<independent_variable>
<variable_name>phase_pressure
- </variable_name>
+ </variable_name>
<reference_condition>1e5
- </reference_condition>
+ </reference_condition>
<slope>4.6511627906976743356e-10
- </slope>
+ </slope>
</independent_variable>
</property>
</properties>
@@ -112,7 +112,18 @@
<parameter>
<name>k</name>
<type>Constant</type>
- <values>1.e-12 0.0 0.0 5.e-12 </values>
+ <!--Permeability tensor in the local coordinate system on the
+ inclined plane. Its compononts associated with the off plane
+ direction have to be set to zero. In other words, the local
+ permeability tensor has to input as a diagonal block of the full
+ size tensor in the global coodinate system. The bases of the
+ local coordinate system has to be defined in order to rotate
+ this local tensor to the global coodinate system.-->
+ <values>1.e-12 0.0 0.0
+ 0.0 5.e-12 0.0
+ 0.0 0.0 0.0</values>
+ <use_local_coordinate_system>true
+ </use_local_coordinate_system>
</parameter>
<parameter>
<name>p0</name>
@@ -129,7 +140,28 @@
<type>Constant</type>
<value>1.0e+6</value>
</parameter>
+ <parameter>
+ <name>b0</name>
+ <type>Constant</type>
+ <values> -0.5000001748438417 0.0 0.8660253028382761</values>
+ </parameter>
+ <parameter>
+ <name>b1</name>
+ <type>Constant</type>
+ <values>0.0 -1.0 0.0</values>
+ </parameter>
+ <parameter>
+ <name>b2</name>
+ <!-- off plane direction -->
+ <type>Constant</type>
+ <values>0.8660253028382761 0.0 0.5000001748438417</values>
+ </parameter>
</parameters>
+ <local_coordinate_system>
+ <basis_vector_0>b0</basis_vector_0>
+ <basis_vector_1>b1</basis_vector_1>
+ <basis_vector_2>b2</basis_vector_2>
+ </local_coordinate_system>
<process_variables>
<process_variable>
<name>pressure</name>
@@ -139,12 +171,12 @@
<boundary_conditions>
<boundary_condition>
<type>Dirichlet</type>
- <mesh>geometry_top</mesh>
+ <mesh>geometry_top</mesh>
<parameter>p_top</parameter>
</boundary_condition>
<boundary_condition>
<type>Dirichlet</type>
- <mesh>geometry_bottom</mesh>
+ <mesh>geometry_bottom</mesh>
<parameter>p_bottom</parameter>
</boundary_condition>
</boundary_conditions>