diff --git a/Tests/Data/TH2M/H2M/StrainDependentPermeability/Strain_Dependent_Permeability_Test.prj b/Tests/Data/TH2M/H2M/StrainDependentPermeability/Strain_Dependent_Permeability_Test.prj
index 034a310523d52d8d99257e18242e6ba30795bcc5..8a84de4804071439db11ec5955d6248feb8d4fad 100644
--- a/Tests/Data/TH2M/H2M/StrainDependentPermeability/Strain_Dependent_Permeability_Test.prj
+++ b/Tests/Data/TH2M/H2M/StrainDependentPermeability/Strain_Dependent_Permeability_Test.prj
@@ -304,6 +304,7 @@
                 <variable>velocity_gas</variable>
                 <variable>velocity_liquid</variable>
                 <variable>liquid_pressure</variable>
+                <variable>liquid_pressure_interpolated</variable>
                 <variable>liquid_density</variable>
                 <variable>gas_density</variable>
                 <variable>porosity</variable>
diff --git a/Tests/Data/ThermoRichardsMechanics/A2/A2.prj b/Tests/Data/ThermoRichardsMechanics/A2/A2.prj
index f97d31b350b409caf9431ff9a8121701d474a33b..af78397f558a88fd29e3ce7d02bad4a6145420ed 100644
--- a/Tests/Data/ThermoRichardsMechanics/A2/A2.prj
+++ b/Tests/Data/ThermoRichardsMechanics/A2/A2.prj
@@ -245,6 +245,11 @@
                 <variable>velocity</variable>
                 <variable>liquid_density</variable>
                 <variable>viscosity</variable>
+                <variable>HeatFlowRate</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>liquid_density_avg</variable>
+                <variable>viscosity_avg</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/ThermoRichardsMechanics/FullySaturatedFlowMechanics/flow_fully_saturated.prj b/Tests/Data/ThermoRichardsMechanics/FullySaturatedFlowMechanics/flow_fully_saturated.prj
index 06fafa874af2a1b8c76bb996579851ed2c1dc425..b5585d0921a1e3af01cf0e3d5e77dd944b305fea 100644
--- a/Tests/Data/ThermoRichardsMechanics/FullySaturatedFlowMechanics/flow_fully_saturated.prj
+++ b/Tests/Data/ThermoRichardsMechanics/FullySaturatedFlowMechanics/flow_fully_saturated.prj
@@ -164,6 +164,8 @@
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
                 <variable>saturation</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/ThermoRichardsMechanics/LinearMechanics/mechanics_linear.prj b/Tests/Data/ThermoRichardsMechanics/LinearMechanics/mechanics_linear.prj
index 16da0150f0f22f9f468de587551fe399877d6707..55bd36253a0d3fbd34fed69f88cdd16edafbd72c 100644
--- a/Tests/Data/ThermoRichardsMechanics/LinearMechanics/mechanics_linear.prj
+++ b/Tests/Data/ThermoRichardsMechanics/LinearMechanics/mechanics_linear.prj
@@ -184,6 +184,7 @@
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
                 <variable>saturation</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/ThermoRichardsMechanics/MFront/ThermoPoroElasticity/uniaxial_isothermal_drainage_imbibition_extended_mfront_model.prj b/Tests/Data/ThermoRichardsMechanics/MFront/ThermoPoroElasticity/uniaxial_isothermal_drainage_imbibition_extended_mfront_model.prj
index 09344c2fc4ac91018dce25489b84e7697adb4741..4687b388bc523a2daae139d36649395e7d91aa23 100644
--- a/Tests/Data/ThermoRichardsMechanics/MFront/ThermoPoroElasticity/uniaxial_isothermal_drainage_imbibition_extended_mfront_model.prj
+++ b/Tests/Data/ThermoRichardsMechanics/MFront/ThermoPoroElasticity/uniaxial_isothermal_drainage_imbibition_extended_mfront_model.prj
@@ -195,6 +195,9 @@
                 <variable>saturation</variable>
                 <variable>sigma_total</variable>
                 <variable>epsilon</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>temperature_interpolated</variable>
+                <variable>saturation_avg</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/ThermoRichardsMechanics/PointHeatSource/point_heat_source_2D.prj b/Tests/Data/ThermoRichardsMechanics/PointHeatSource/point_heat_source_2D.prj
index 91d13d8bf7ab47f7fcba58c28ff28b443b9e4bda..9cf5da9d693aefa7525c91c69686a9482738f5a3 100644
--- a/Tests/Data/ThermoRichardsMechanics/PointHeatSource/point_heat_source_2D.prj
+++ b/Tests/Data/ThermoRichardsMechanics/PointHeatSource/point_heat_source_2D.prj
@@ -172,9 +172,14 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
                 <variable>temperature</variable>
+                <variable>temperature_interpolated</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>HeatFlowRate</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
@@ -193,6 +198,9 @@
                 <variable>temperature</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>HeatFlowRate</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
             <meshes>
diff --git a/Tests/Data/ThermoRichardsMechanics/RichardsFlow2D/RichardsFlow_2d_small.prj b/Tests/Data/ThermoRichardsMechanics/RichardsFlow2D/RichardsFlow_2d_small.prj
index da99593e28c9dd4b8581e4152a727295236a8fc6..b550bdce11f036ae70b783952afa0db14b57a620 100644
--- a/Tests/Data/ThermoRichardsMechanics/RichardsFlow2D/RichardsFlow_2d_small.prj
+++ b/Tests/Data/ThermoRichardsMechanics/RichardsFlow2D/RichardsFlow_2d_small.prj
@@ -201,6 +201,10 @@
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
                 <variable>saturation</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>saturation_avg</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/ThermoRichardsMechanics/Simple3DThermoMechanicsFromTM/cube_1e3.prj b/Tests/Data/ThermoRichardsMechanics/Simple3DThermoMechanicsFromTM/cube_1e3.prj
index 4d7774421bae1e7179001d68cc73b96760845ab1..a1dd762634560927974919372baa06566bf6a3cc 100644
--- a/Tests/Data/ThermoRichardsMechanics/Simple3DThermoMechanicsFromTM/cube_1e3.prj
+++ b/Tests/Data/ThermoRichardsMechanics/Simple3DThermoMechanicsFromTM/cube_1e3.prj
@@ -161,6 +161,7 @@
                 <variable>temperature</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>