diff --git a/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/Decovalex-0-TRF.prj b/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/Decovalex-0-TRF.prj
index 8865cb39463ca9c37838e6d2d5d1ebe42c6933f9..140edaf9951cdc5eed8b896a05682174923b9943 100644
--- a/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/Decovalex-0-TRF.prj
+++ b/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/Decovalex-0-TRF.prj
@@ -559,12 +559,14 @@
         <parameter>
             <name>b0</name>
             <type>Constant</type>
-            <values> 0.829037572555042 0.559192903470747</values>
+            <!-- 34 degrees bedding angle -->
+            <values> 0.829037572555041692 0.55919290347074683016</values>
         </parameter>
         <parameter>
             <name>b1</name>
             <type>Constant</type>
-            <values>-0.559192903470747 0.829037572555042</values>
+            <!-- 34 degrees bedding angle -->
+            <values>-0.55919290347074683016 0.829037572555041692</values>
         </parameter>
         <parameter>
             <name>ThermalConductivityClay</name>
diff --git a/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/WithPicardNonLinearSolverAndPETSc/Decovalex-0-TRF.prj b/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/WithPicardNonLinearSolverAndPETSc/Decovalex-0-TRF.prj
index 08853826e3d4939fbf22692d07b33337f7b1e4d7..92f99f0549dba8f18d8652e3b52de44c6ddb2488 100644
--- a/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/WithPicardNonLinearSolverAndPETSc/Decovalex-0-TRF.prj
+++ b/Tests/Data/ThermoRichardsFlow/TaskCDECOVALEX2023/WithPicardNonLinearSolverAndPETSc/Decovalex-0-TRF.prj
@@ -559,12 +559,14 @@
         <parameter>
             <name>b0</name>
             <type>Constant</type>
-            <values> 0.829037572555042 0.559192903470747</values>
+            <!-- 34 degrees bedding angle -->
+            <values> 0.829037572555041692 0.55919290347074683016</values>
         </parameter>
         <parameter>
             <name>b1</name>
             <type>Constant</type>
-            <values>-0.559192903470747 0.829037572555042</values>
+            <!-- 34 degrees bedding angle -->
+            <values>-0.55919290347074683016 0.829037572555041692</values>
         </parameter>
         <parameter>
             <name>ThermalConductivityClay</name>
diff --git a/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj b/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
index 10d9d7bb72a4058f92fda1248fd29b7da53bc9ad..b429403bc6f081465a9ba23c0298558a88fae326 100644
--- a/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
+++ b/Tests/Data/ThermoRichardsMechanics/BodyForce/square.prj
@@ -194,12 +194,14 @@
         <parameter>
             <name>b0</name>
             <type>Constant</type>
-            <values> 0.829037572555042 0.559192903470747</values>
+            <!-- 34 degrees bedding angle -->
+            <values> 0.829037572555041692 0.55919290347074683016</values>
         </parameter>
         <parameter>
             <name>b1</name>
             <type>Constant</type>
-            <values>-0.559192903470747 0.829037572555042</values>
+            <!-- 34 degrees bedding angle -->
+            <values>-0.55919290347074683016 0.829037572555041692</values>
         </parameter>
         <parameter>
             <name>ThermalConductivityClay</name>
diff --git a/Tests/Data/ThermoRichardsMechanics/TaskCDECOVALEX2023/Decovalex-0.prj b/Tests/Data/ThermoRichardsMechanics/TaskCDECOVALEX2023/Decovalex-0.prj
index daff056884bc3c2b4950de708a906498b934bc0f..5420e1018ebd950cc97ceb2e2375829ebe1db9ee 100644
--- a/Tests/Data/ThermoRichardsMechanics/TaskCDECOVALEX2023/Decovalex-0.prj
+++ b/Tests/Data/ThermoRichardsMechanics/TaskCDECOVALEX2023/Decovalex-0.prj
@@ -37,16 +37,11 @@
                 <displacement>displacement</displacement>
             </process_variables>
             <secondary_variables>
-                <secondary_variable internal_name="saturation"
-                                    output_name="saturation"/>
-                <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="liquid_density"
-                                    output_name="liquid_density"/>
+                <secondary_variable internal_name="saturation" output_name="saturation"/>
+                <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="liquid_density" output_name="liquid_density"/>
             </secondary_variables>
             <specific_body_force>0 0</specific_body_force>
             <initial_stress>Initial_stress</initial_stress>
@@ -552,14 +547,15 @@
         </output>
     </time_loop>
     <local_coordinate_system>
-        <basis_vector_0 implicit = "true" />
+        <basis_vector_0 implicit="true"/>
         <basis_vector_1>b1</basis_vector_1>
     </local_coordinate_system>
     <parameters>
         <parameter>
             <name>b1</name>
             <type>Constant</type>
-            <values>-0.559192903470747 0.829037572555042</values>
+            <!-- 34 degrees bedding angle -->
+            <values>-0.55919290347074683016 0.829037572555041692</values>
         </parameter>
         <parameter>
             <name>ThermalConductivityClay</name>