diff --git a/ProcessLib/RichardsMechanics/Tests.cmake b/ProcessLib/RichardsMechanics/Tests.cmake
index 39ec1743fd21c5479833ebe5b7d89ff64e590e28..9b3a6fee10bd8bb89b5e4bbe8bbbd9c30e7cc885 100644
--- a/ProcessLib/RichardsMechanics/Tests.cmake
+++ b/ProcessLib/RichardsMechanics/Tests.cmake
@@ -182,3 +182,17 @@ AddTest(
     GLOB RichardsFlow_2d_richardsflow_pcs_0_ts_*.vtu pressure pressure 5e-11 1e-15
     GLOB RichardsFlow_2d_richardsflow_pcs_0_ts_*.vtu saturation saturation 1e-14 1e-15
 )
+
+#PETSc
+AddTest(
+    NAME Parallel_RichardsMechanics_RichardsFlow_2d_richardsflow
+    PATH RichardsMechanics
+    EXECUTABLE_ARGS RichardsFlow_2d_richardsflow.prj
+    WRAPPER mpirun
+    WRAPPER_ARGS -np 1
+    TESTER vtkdiff
+    REQUIREMENTS OGS_USE_MPI
+    DIFF_DATA
+    RichardsFlow_2d_richardsflow_pcs_0_ts_99_t_1900.000000.vtu RichardsFlow_2d_richardsflow_pcs_0_ts_99_t_1900_000000_0.vtu pressure pressure 5e-8 1e-10
+    RichardsFlow_2d_richardsflow_pcs_0_ts_99_t_1900.000000.vtu RichardsFlow_2d_richardsflow_pcs_0_ts_99_t_1900_000000_0.vtu saturation saturation 1e-10 1e-11
+)
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_richardsflow.prj b/Tests/Data/RichardsMechanics/RichardsFlow_2d_richardsflow.prj
index 64d61581d3855f1331c748a657e0e4af50b2ac75..6179be63dbee1e8671c5554ea6aa4539cbb98908 100644
--- a/Tests/Data/RichardsMechanics/RichardsFlow_2d_richardsflow.prj
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_richardsflow.prj
@@ -195,6 +195,10 @@
                 <solver_type>PardisoLU</solver_type>
                 -->
             </eigen>
+            <petsc>
+                <prefix>RM</prefix>
+                <parameters>-RM_ksp_type bcgs -RM_pc_type bjacobi -RM_ksp_rtol 1e-16 -RM_ksp_max_it 10000</parameters>
+            </petsc>
         </linear_solver>
     </linear_solvers>
 </OpenGeoSysProject>