diff --git a/ProcessLib/RichardsMechanics/Tests.cmake b/ProcessLib/RichardsMechanics/Tests.cmake
index f72b6458778898d1ecdb5b396ab23fcde82efeaa..4b55cd167fd482c676aa5054f7d2a1d8b0f4528e 100644
--- a/ProcessLib/RichardsMechanics/Tests.cmake
+++ b/ProcessLib/RichardsMechanics/Tests.cmake
@@ -244,3 +244,39 @@ AddTest(
GLOB orthotropic_swelling_xy_pcs_0_ts_*.vtu porosity porosity 1e-15 1e-15
GLOB orthotropic_swelling_xy_pcs_0_ts_*.vtu porosity_avg porosity_avg 1e-15 1e-15
)
+AddTest(
+ NAME RichardsMechanics_bishops_effective_stress_power_law
+ PATH RichardsMechanics
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS bishops_effective_stress_power_law.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT OGS_USE_MPI
+ DIFF_DATA
+ GLOB bishops_effective_stress_power_law_pcs_0_ts_*.vtu pressure pressure 1e-15 0
+ GLOB bishops_effective_stress_power_law_pcs_0_ts_*.vtu saturation saturation 1e-14 1e-15
+ GLOB bishops_effective_stress_power_law_pcs_0_ts_*.vtu displacement displacement 1e-15 0
+ GLOB bishops_effective_stress_power_law_pcs_0_ts_*.vtu sigma sigma 5e-15 0
+ GLOB bishops_effective_stress_power_law_pcs_0_ts_*.vtu epsilon epsilon 2e-15 0
+ GLOB bishops_effective_stress_power_law_pcs_0_ts_*.vtu velocity velocity 1e-15 0
+ GLOB bishops_effective_stress_power_law_pcs_0_ts_*.vtu porosity porosity 5e-15 1e-15
+ GLOB bishops_effective_stress_power_law_pcs_0_ts_*.vtu porosity_avg porosity_avg 5e-15 1e-15
+)
+AddTest(
+ NAME RichardsMechanics_bishops_effective_stress_saturation_cutoff
+ PATH RichardsMechanics
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS bishops_effective_stress_saturation_cutoff.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT OGS_USE_MPI
+ DIFF_DATA
+ GLOB bishops_effective_stress_saturation_cutoff_pcs_0_ts_*.vtu pressure pressure 1e-15 0
+ GLOB bishops_effective_stress_saturation_cutoff_pcs_0_ts_*.vtu saturation saturation 1e-14 1e-15
+ GLOB bishops_effective_stress_saturation_cutoff_pcs_0_ts_*.vtu displacement displacement 1e-15 0
+ GLOB bishops_effective_stress_saturation_cutoff_pcs_0_ts_*.vtu sigma sigma 1e-15 0
+ GLOB bishops_effective_stress_saturation_cutoff_pcs_0_ts_*.vtu epsilon epsilon 1e-15 0
+ GLOB bishops_effective_stress_saturation_cutoff_pcs_0_ts_*.vtu velocity velocity 1e-15 0
+ GLOB bishops_effective_stress_saturation_cutoff_pcs_0_ts_*.vtu porosity porosity 5e-15 1e-15
+ GLOB bishops_effective_stress_saturation_cutoff_pcs_0_ts_*.vtu porosity_avg porosity_avg 5e-15 1e-15
+)
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law.prj b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law.prj
new file mode 100644
index 0000000000000000000000000000000000000000..6e1ec61b27ab8ec280ff96dfc227bb413ec69816
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law.prj
@@ -0,0 +1,335 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>cube_1x1x1_hex20_1e0.vtu</mesh>
+ <geometry>cube_1x1x1.gml</geometry>
+ <processes>
+ <process>
+ <name>RM</name>
+ <type>RICHARDS_MECHANICS</type>
+ <integration_order>3</integration_order>
+ <dimension>3</dimension>
+ <!--
+ <jacobian_assembler>
+ <type>CentralDifferences</type>
+ <component_magnitudes>1e5 1e-5 1e-5 1e-5</component_magnitudes>
+ <relative_epsilons>1e-8 1e-8 1e-8 1e-8</relative_epsilons>
+ </jacobian_assembler>
+ -->
+ <constitutive_relation>
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <pressure>pressure</pressure>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+ <secondary_variable internal_name="velocity" output_name="velocity"/>
+ <secondary_variable internal_name="saturation" output_name="saturation"/>
+ <secondary_variable internal_name="porosity" output_name="porosity"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <media>
+ <medium>
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>bulk_modulus</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>bulk_modulus</name>
+ <type>Constant</type>
+ <value>10</value>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.5</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>reference_temperature</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </property>
+ <property>
+ <name>saturation</name>
+ <type>SaturationVanGenuchten</type>
+ <residual_liquid_saturation>0.1</residual_liquid_saturation>
+ <residual_gas_saturation>0.05</residual_gas_saturation>
+ <exponent>0.8</exponent>
+ <entry_pressure>0.5</entry_pressure>
+ </property>
+ <property>
+ <name>relative_permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>bishops_effective_stress</name>
+ <type>BishopsPowerLaw</type>
+ <exponent>3</exponent>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <processes>
+ <process ref="RM">
+ <nonlinear_solver>nonlinear_solver</nonlinear_solver>
+ <compensate_non_equilibrium_initial_residuum>true</compensate_non_equilibrium_initial_residuum>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstols>1e-13 1e-13 1e-13 1e-13</abstols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>3</t_end>
+ <timesteps>
+ <pair>
+ <repeat>50</repeat>
+ <delta_t>.02</delta_t>
+ </pair>
+
+ <pair>
+ <repeat>4</repeat>
+ <delta_t>.25</delta_t>
+ </pair>
+
+ <pair>
+ <repeat>50</repeat>
+ <delta_t>.02</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>bishops_effective_stress_power_law</prefix>
+ <timesteps>
+ <pair>
+ <repeat>1000</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>displacement</variable>
+ <variable>pressure</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ <variable>velocity</variable>
+ <variable>saturation</variable>
+ <variable>porosity</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <local_coordinate_system>
+ <basis_vector_0>e0</basis_vector_0>
+ <basis_vector_1>e1</basis_vector_1>
+ <basis_vector_2>e2</basis_vector_2>
+ </local_coordinate_system>
+ <parameters>
+ <parameter>
+ <name>e0</name>
+ <type>Constant</type>
+ <values>1 0 0</values>
+ </parameter>
+ <parameter>
+ <name>e1</name>
+ <type>Constant</type>
+ <values>0 1 0</values>
+ </parameter>
+ <parameter>
+ <name>e2</name>
+ <type>Constant</type>
+ <values>0 0 1</values>
+ </parameter>
+ <!-- Mechanics -->
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.3</value>
+ </parameter>
+ <!-- Model parameters -->
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>p0</name>
+ <type>Constant</type>
+ <value>0.1</value>
+ </parameter>
+ <parameter>
+ <name>dirichlet0</name>
+ <type>Constant</type>
+ <values>0</values>
+ </parameter>
+ <parameter>
+ <name>pressure_bc</name>
+ <type>CurveScaled</type>
+ <curve>pressure_ramp</curve>
+ <parameter>p0</parameter>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>pressure_ramp</name>
+ <coords>0 1 2 3</coords>
+ <values>1 -10 -10 1</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>2</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>p0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>back</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>nonlinear_solver</name>
+ <type>Newton</type>
+ <max_iter>100</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ <scaling>true</scaling>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_104_t_3.000000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_104_t_3.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..6968e497ff126cbd1398dadf11e4548c7f878f0b
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_104_t_3.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f02b79a6560a16d684770632812fe55379674caedfee6a3f5092211c5696f937
+size 8187
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_25_t_0.500000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_25_t_0.500000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..5e149d5aa840af32d5fe887ff4464318cfd3b59e
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_25_t_0.500000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:974584d21ffcdabc7d84ad77b22a4006613ee22b719a552a403b1786f74408b1
+size 7423
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_50_t_1.000000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_50_t_1.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..39d85628b4a62dfa38e5bd2a93103991736fe4c8
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_50_t_1.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:de0db17f4546c25722993b64087a6b41168be7bbd8c1ab996e68a631ef44470d
+size 7299
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_52_t_1.500000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_52_t_1.500000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..10bd69b5350bf54f6db098a62df0169658042deb
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_52_t_1.500000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:338240d63d18184e0297419a33ecfa2308f866eb705978307334c41dd90fd8e7
+size 7283
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_54_t_2.000000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_54_t_2.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..ca8ef43f40e1267b0635d5002f3e852d8ebf6a42
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_54_t_2.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:abbb64b4f5da0d1468f73a898b4023150e39e88f748f38b71b512b81c2f29537
+size 7447
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_79_t_2.500000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_79_t_2.500000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..0b5677ebfeb15ce41a3cb485663c130908d4cb23
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law_pcs_0_ts_79_t_2.500000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:4d17a6b7751c56d1fda352d13f8211618422b791a6a993025594d4e3f703c5cc
+size 7331
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff.prj b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff.prj
new file mode 100644
index 0000000000000000000000000000000000000000..fe65fa784251bf5070761b66459f363cb14d0688
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff.prj
@@ -0,0 +1,335 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>cube_1x1x1_hex20_1e0.vtu</mesh>
+ <geometry>cube_1x1x1.gml</geometry>
+ <processes>
+ <process>
+ <name>RM</name>
+ <type>RICHARDS_MECHANICS</type>
+ <integration_order>3</integration_order>
+ <dimension>3</dimension>
+ <!--
+ <jacobian_assembler>
+ <type>CentralDifferences</type>
+ <component_magnitudes>1e5 1e-5 1e-5 1e-5</component_magnitudes>
+ <relative_epsilons>1e-8 1e-8 1e-8 1e-8</relative_epsilons>
+ </jacobian_assembler>
+ -->
+ <constitutive_relation>
+ <type>LinearElasticIsotropic</type>
+ <youngs_modulus>E</youngs_modulus>
+ <poissons_ratio>nu</poissons_ratio>
+ </constitutive_relation>
+ <process_variables>
+ <displacement>displacement</displacement>
+ <pressure>pressure</pressure>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="sigma" output_name="sigma"/>
+ <secondary_variable internal_name="epsilon" output_name="epsilon"/>
+ <secondary_variable internal_name="velocity" output_name="velocity"/>
+ <secondary_variable internal_name="saturation" output_name="saturation"/>
+ <secondary_variable internal_name="porosity" output_name="porosity"/>
+ </secondary_variables>
+ <specific_body_force>0 0 0</specific_body_force>
+ </process>
+ </processes>
+ <media>
+ <medium>
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>bulk_modulus</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>bulk_modulus</name>
+ <type>Constant</type>
+ <value>10</value>
+ </property>
+ <property>
+ <name>biot_coefficient</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.5</value>
+ </property>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>reference_temperature</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </property>
+ <property>
+ <name>saturation</name>
+ <type>SaturationVanGenuchten</type>
+ <residual_liquid_saturation>0.1</residual_liquid_saturation>
+ <residual_gas_saturation>0.05</residual_gas_saturation>
+ <exponent>0.8</exponent>
+ <entry_pressure>0.5</entry_pressure>
+ </property>
+ <property>
+ <name>relative_permeability</name>
+ <type>Constant</type>
+ <value>1</value>
+ </property>
+ <property>
+ <name>bishops_effective_stress</name>
+ <type>BishopsSaturationCutoff</type>
+ <cutoff_value>0.95</cutoff_value>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <time_loop>
+ <processes>
+ <process ref="RM">
+ <nonlinear_solver>nonlinear_solver</nonlinear_solver>
+ <compensate_non_equilibrium_initial_residuum>true</compensate_non_equilibrium_initial_residuum>
+ <convergence_criterion>
+ <type>PerComponentDeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstols>1e-13 1e-13 1e-13 1e-13</abstols>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0</t_initial>
+ <t_end>3</t_end>
+ <timesteps>
+ <pair>
+ <repeat>50</repeat>
+ <delta_t>.02</delta_t>
+ </pair>
+
+ <pair>
+ <repeat>4</repeat>
+ <delta_t>.25</delta_t>
+ </pair>
+
+ <pair>
+ <repeat>50</repeat>
+ <delta_t>.02</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>bishops_effective_stress_saturation_cutoff</prefix>
+ <timesteps>
+ <pair>
+ <repeat>1000</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable>displacement</variable>
+ <variable>pressure</variable>
+ <variable>sigma</variable>
+ <variable>epsilon</variable>
+ <variable>velocity</variable>
+ <variable>saturation</variable>
+ <variable>porosity</variable>
+ </variables>
+ </output>
+ </time_loop>
+ <local_coordinate_system>
+ <basis_vector_0>e0</basis_vector_0>
+ <basis_vector_1>e1</basis_vector_1>
+ <basis_vector_2>e2</basis_vector_2>
+ </local_coordinate_system>
+ <parameters>
+ <parameter>
+ <name>e0</name>
+ <type>Constant</type>
+ <values>1 0 0</values>
+ </parameter>
+ <parameter>
+ <name>e1</name>
+ <type>Constant</type>
+ <values>0 1 0</values>
+ </parameter>
+ <parameter>
+ <name>e2</name>
+ <type>Constant</type>
+ <values>0 0 1</values>
+ </parameter>
+ <!-- Mechanics -->
+ <parameter>
+ <name>E</name>
+ <type>Constant</type>
+ <value>1</value>
+ </parameter>
+ <parameter>
+ <name>nu</name>
+ <type>Constant</type>
+ <value>.3</value>
+ </parameter>
+ <!-- Model parameters -->
+ <parameter>
+ <name>displacement0</name>
+ <type>Constant</type>
+ <values>0 0 0</values>
+ </parameter>
+ <parameter>
+ <name>p0</name>
+ <type>Constant</type>
+ <value>0.1</value>
+ </parameter>
+ <parameter>
+ <name>dirichlet0</name>
+ <type>Constant</type>
+ <values>0</values>
+ </parameter>
+ <parameter>
+ <name>pressure_bc</name>
+ <type>CurveScaled</type>
+ <curve>pressure_ramp</curve>
+ <parameter>p0</parameter>
+ </parameter>
+ </parameters>
+ <curves>
+ <curve>
+ <name>pressure_ramp</name>
+ <coords>0 1 2 3</coords>
+ <values>1 -10 -10 1</values>
+ </curve>
+ </curves>
+ <process_variables>
+ <process_variable>
+ <name>displacement</name>
+ <components>3</components>
+ <order>2</order>
+ <initial_condition>displacement0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <component>0</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <component>1</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <component>2</component>
+ <parameter>dirichlet0</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>pressure</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>p0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>bottom</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>front</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>top</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>back</geometry>
+ <type>Dirichlet</type>
+ <parameter>pressure_bc</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>nonlinear_solver</name>
+ <type>Newton</type>
+ <max_iter>100</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ <scaling>true</scaling>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_104_t_3.000000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_104_t_3.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..42658f65522ffdd625f12045287b3b3e1c3a1352
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_104_t_3.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:4ea27efd0ddc782cf3679849c7bd64f229573166b0a977d00956db5fd1e93a0b
+size 8171
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_25_t_0.500000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_25_t_0.500000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..06fc2256d7a78c760ded73dcfa799d8866acf107
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_25_t_0.500000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:e324c4251c2d820142b1cd9a35b2cf0186123e4c60eeed810f03c5b515519a3e
+size 7379
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_50_t_1.000000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_50_t_1.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..785f5b776b55079e525bd786ad7d297f4a81a2b7
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_50_t_1.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:e0edbb427dd2793f32ebcebaf8995a20a4e6283da05e4ff02e81f9bedd4499c3
+size 7251
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_52_t_1.500000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_52_t_1.500000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..bc2577ac6537d64072031cb11ef931a60a5030f8
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_52_t_1.500000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:719db305516fec867f2ad7220b2ee75b7d4ac390e173bbea0e7773b294b092e4
+size 7235
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_54_t_2.000000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_54_t_2.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..cea5428a833297926ff15ae58618804ef7bd0da7
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_54_t_2.000000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f1136f1ff8cf7c03b9763863b03c7907e753b4add265600285c1f68e2610af43
+size 7395
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_79_t_2.500000.vtu b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_79_t_2.500000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..8f3ab459a4f9e13934d2afae910b66175348824a
--- /dev/null
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff_pcs_0_ts_79_t_2.500000.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:ac0aab0bc0b91b2028da9ff6d239cfe0aeaa6a1b48e5125da1004804541f67b4
+size 7291