diff --git a/ProcessLib/ThermoRichardsMechanics/Tests.cmake b/ProcessLib/ThermoRichardsMechanics/Tests.cmake
index 54c39047d50f2d26fa4c2d41498a546a6a968fc0..d08ec317ed2e86b61a69061afab7727022904add 100644
--- a/ProcessLib/ThermoRichardsMechanics/Tests.cmake
+++ b/ProcessLib/ThermoRichardsMechanics/Tests.cmake
@@ -448,6 +448,22 @@ AddTest(
     expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_ts_10_t_864000.000000.vtu sigma sigma 5e-2 1e-8
 )
 
+AddTest(
+    NAME ThermoRichardsMechanics_total_initial_stress_dont_apply_body_force_for_deformation
+    PATH ThermoRichardsMechanics/BodyForce
+    RUNTIME 1
+    EXECUTABLE ogs
+    EXECUTABLE_ARGS square_total_stress_test.xml
+    WRAPPER time
+    TESTER vtkdiff
+    DIFF_DATA
+    expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_total_stess0_test_ts_10_t_864000.000000.vtu displacement displacement 1e-10 1e-10
+    expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_total_stess0_test_ts_10_t_864000.000000.vtu pressure pressure 5e-2 1e-8
+    expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_total_stess0_test_ts_10_t_864000.000000.vtu temperature temperature 1e-10 1e-10
+    expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_total_stess0_test_ts_10_t_864000.000000.vtu epsilon epsilon 1e-10 1e-10
+    expected_dont_apply_body_force_for_deformation_ts_0_t_0.000000.vtu dont_apply_body_force_for_deformation_total_stess0_test_ts_10_t_864000.000000.vtu sigma sigma 5e-2 1e-8
+)
+
 if(OGS_USE_MFRONT)
     OgsTest(PROJECTFILE ThermoRichardsMechanics/MultiMaterialEhlers/square_1e1_2_matIDs.prj RUNTIME 1)
     OgsTest(PROJECTFILE ThermoRichardsMechanics/MultiMaterialEhlers/square_1e1_2_matIDs_restart.prj RUNTIME 1)
diff --git a/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj b/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
index 1b7bcb4e91200f3a227b0ff0245f67d1e304eafc..10d9d7bb72a4058f92fda1248fd29b7da53bc9ad 100644
--- a/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
+++ b/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
@@ -344,6 +344,14 @@
             <eigen>
                 <solver_type>SparseLU</solver_type>
             </eigen>
+            <petsc>
+                <parameters>
+                    -ksp_type bcgs
+                    -pc_type jacobi
+                    -ksp_rtol 1.e-20 -ksp_atol 1.e-10
+                    -ksp_max_it 4000
+                </parameters>
+            </petsc>
         </linear_solver>
     </linear_solvers>
 </OpenGeoSysProject>
diff --git a/Tests/Data/ThermoRichardsMechanics/BodyForce/square_total_stress_test.xml b/Tests/Data/ThermoRichardsMechanics/BodyForce/square_total_stress_test.xml
new file mode 100644
index 0000000000000000000000000000000000000000..d9110b3d59a8741f76eb1a414ea84495b6208723
--- /dev/null
+++ b/Tests/Data/ThermoRichardsMechanics/BodyForce/square_total_stress_test.xml
@@ -0,0 +1,17 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProjectDiff base_file="square.prj">
+    <remove sel="/*/processes/process[1]/initial_stress"/>
+    <add sel="/*/processes/process[1]">
+        <initial_stress type = "total">total_stress0</initial_stress>
+    </add>
+    <add sel="/*/parameters">
+        <parameter>
+            <name>total_stress0</name>
+            <type>Constant</type>
+            <values>-4.5e6 -2.5e6 -6.0e6 0.0</values>
+        </parameter>
+    </add>
+    <replace msel="/*/time_loop/output/prefix/text()">
+        dont_apply_body_force_for_deformation_total_stess0_test
+    </replace>
+</OpenGeoSysProjectDiff>