diff --git a/MeshLib/IO/VtkIO/VtuInterface-impl.h b/MeshLib/IO/VtkIO/VtuInterface-impl.h
index 5b0a3df70437cdb6bb63f243987ed91e29d3d784..55e13cdeda7c4188527b3811670b7dd19a5312ad 100644
--- a/MeshLib/IO/VtkIO/VtuInterface-impl.h
+++ b/MeshLib/IO/VtkIO/VtuInterface-impl.h
@@ -21,6 +21,7 @@
 #include <vtkUnstructuredGrid.h>
 
 #include "BaseLib/Logging.h"
+#include "MeshLib/Mesh.h"
 #include "MeshLib/Vtk/VtkMappedMeshSource.h"
 #include "VtuInterface.h"
 
@@ -48,6 +49,26 @@ bool VtuInterface::writeVTU(std::string const& file_name,
         return false;
     }
 
+#ifdef USE_PETSC
+    if (_mesh->getProperties().existsPropertyVector<unsigned char>(
+            "vtkGhostType", MeshLib::MeshItemType::Cell, 1))
+    {
+        auto* ghost_cell_property =
+            _mesh->getProperties().getPropertyVector<unsigned char>(
+                "vtkGhostType", MeshLib::MeshItemType::Cell, 1);
+        if (ghost_cell_property)
+        {
+            const_cast<MeshLib::PropertyVector<unsigned char>*>(
+                ghost_cell_property)
+                ->is_for_output = true;
+        }
+    }
+    else
+    {
+        DBUG("No vtkGhostType data in mesh '{}'.", _mesh->getName());
+    }
+#endif
+
     vtkNew<MeshLib::VtkMappedMeshSource> vtkSource;
     vtkSource->SetMesh(_mesh);
 
diff --git a/MeshLib/IO/XDMF/transformData.cpp b/MeshLib/IO/XDMF/transformData.cpp
index 77df3c0bd8917717db8062f892219ebfb89b0791..b3d0809f8f3438934336e5fcdc2196a46f6c534f 100644
--- a/MeshLib/IO/XDMF/transformData.cpp
+++ b/MeshLib/IO/XDMF/transformData.cpp
@@ -216,6 +216,11 @@ std::vector<XdmfHdfData> transformAttributes(MeshLib::Mesh const& mesh,
             continue;
         }
 
+        if (!property_base->is_for_output)
+        {
+            continue;
+        }
+
         if (auto const attribute = transformAttribute(
                 std::pair(std::string(name), property_base), n_files))
         {
diff --git a/MeshLib/PropertyVector.h b/MeshLib/PropertyVector.h
index df8d9991ce65f4589109c2baf6368756c57a1936..58b605e58d03d444d843b0747981bcb070e811f4 100644
--- a/MeshLib/PropertyVector.h
+++ b/MeshLib/PropertyVector.h
@@ -30,6 +30,7 @@ public:
     MeshItemType getMeshItemType() const { return _mesh_item_type; }
     std::string const& getPropertyName() const { return _property_name; }
     int getNumberOfGlobalComponents() const { return _n_components; }
+    bool is_for_output = true;
 
 protected:
     PropertyVectorBase(std::string property_name,
diff --git a/MeshLib/Vtk/VtkMappedMeshSource.cpp b/MeshLib/Vtk/VtkMappedMeshSource.cpp
index 739aee68d26aa58a577eb186ba65fbb1e4a127f3..435c24d90cbe8246a195265ffa2134f5b71cb255 100644
--- a/MeshLib/Vtk/VtkMappedMeshSource.cpp
+++ b/MeshLib/Vtk/VtkMappedMeshSource.cpp
@@ -119,6 +119,10 @@ int VtkMappedMeshSource::RequestData(vtkInformation* /*request*/,
     // Arrays
     for (auto [name, property] : _mesh->getProperties())
     {
+        if (!property->is_for_output)
+        {
+            continue;
+        }
         if (auto p = dynamic_cast<PropertyVector<double>*>(property))
         {
             addProperty(*p);
@@ -166,6 +170,10 @@ int VtkMappedMeshSource::RequestData(vtkInformation* /*request*/,
         {
             addProperty(*p);
         }
+        else if (auto p = dynamic_cast<PropertyVector<uint8_t>*>(property))
+        {
+            addProperty(*p);
+        }
         else
         {
             OGS_FATAL(
@@ -180,7 +188,8 @@ int VtkMappedMeshSource::RequestData(vtkInformation* /*request*/,
                 "\n\t long long,"
                 "\n\t unsigned long long,"
                 "\n\t char,",
-                "\n\t unsigned char.",
+                "\n\t unsigned char,",
+                "\n\t uint8_t.",
                 property->getPropertyName(),
                 typeid(*property).name());
         }
diff --git a/ProcessLib/LiquidFlow/Tests.cmake b/ProcessLib/LiquidFlow/Tests.cmake
index 4613b470d9199bc10b4f7a240a9576c4f95c3bb2..4b56987952842d7bb27ef7e744acd14ed84df80f 100644
--- a/ProcessLib/LiquidFlow/Tests.cmake
+++ b/ProcessLib/LiquidFlow/Tests.cmake
@@ -523,6 +523,7 @@ if(NOT OGS_USE_MPI)
         PROPERTIES WILL_FAIL TRUE
     )
 
+    OgsTest(PROJECTFILE Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet_multiple_outputs_different_variables.xml)
 endif()
 
 AddTest(
@@ -550,7 +551,7 @@ AddTest(
 )
 
 AddTest(
-    NAME SimpleSynthetics_XDMF
+    NAME LiquidFlow_SimpleSynthetics_XDMF
     PATH Parabolic/LiquidFlow/SimpleSynthetics/XDMF
     EXECUTABLE ogs
     EXECUTABLE_ARGS FunctionParameterTest_XDMF.prj
@@ -571,7 +572,7 @@ AddTest(
 )
 
 AddTest(
-    NAME SimpleSynthetics_XDMF_compression_off
+    NAME LiquidFlow_SimpleSynthetics_XDMF_compression_off
     PATH Parabolic/LiquidFlow/SimpleSynthetics/XDMF_compression_off
     EXECUTABLE ogs
     EXECUTABLE_ARGS FunctionParameterTest_XDMF.prj
diff --git a/ProcessLib/Output/Output.cpp b/ProcessLib/Output/Output.cpp
index 19ec56f654a5a625acf157e7779a852fd6c664cf..d6e057a40688bbd4fb595b7eff41690048547d40 100644
--- a/ProcessLib/Output/Output.cpp
+++ b/ProcessLib/Output/Output.cpp
@@ -174,6 +174,28 @@ void Output::outputMeshes(
     std::vector<std::reference_wrapper<const MeshLib::Mesh>> const& meshes)
     const
 {
+    if (_output_data_specification.output_variables.empty())
+    {
+        // special case: no output properties specified => output all properties
+        for (auto const& mesh : meshes)
+        {
+            for (auto [name, property] : mesh.get().getProperties())
+            {
+                property->is_for_output = true;
+            }
+        }
+    }
+    else
+    {
+        for (auto const& mesh : meshes)
+        {
+            for (auto [name, property] : mesh.get().getProperties())
+            {
+                property->is_for_output =
+                    _output_data_specification.output_variables.contains(name);
+            }
+        }
+    }
     _output_format->outputMeshes(timestep, t, iteration, meshes,
                                  _output_data_specification.output_variables);
 }
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_calculatesurfaceflux.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_calculatesurfaceflux.prj
index fc279d7290f39a627d4d7077f373eed70ffc05c7..e5d6f088ba8aa241060da9cf430df02800a779d3 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_calculatesurfaceflux.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_calculatesurfaceflux.prj
@@ -47,6 +47,7 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v      </variable>
+                <variable>specific_flux</variable>
             </variables>
             <meshes>
                 <mesh>cube_1x1x1_hex_1e3</mesh>
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_neumann_calculatesurfaceflux.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_neumann_calculatesurfaceflux.prj
index 064020339fd58dde72462360d8ec2711af957072..5dc7abcf614cd4a36fa693d041659e34073567dc 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_neumann_calculatesurfaceflux.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_1e3_neumann_calculatesurfaceflux.prj
@@ -48,6 +48,7 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v      </variable>
+                <variable>specific_flux</variable>
             </variables>
             <meshes>
                 <mesh>cube_1x1x1_hex_1e3</mesh>
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_front_back.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_front_back.prj
index 273dbc7dde33243517f880ad5c5e63291e7ff265..b3f748ddf48f2b14fb4430e2680268c04d3cce1e 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_front_back.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_front_back.prj
@@ -47,6 +47,7 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v      </variable>
+                <variable>specific_flux</variable>
             </variables>
             <meshes>
                 <mesh>cube_1x1x1_prism_2e3</mesh>
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_left_right.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_left_right.prj
index 404df23b082e0823b541119e2dffb5270409459c..6581369c1a11d3a9084af8cda9b6eb9e0af8e802 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_left_right.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_left_right.prj
@@ -47,6 +47,7 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v      </variable>
+                <variable>specific_flux</variable>
             </variables>
             <meshes>
                 <mesh>cube_1x1x1_prism_2e3</mesh>
diff --git a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_top_bottom.prj b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_top_bottom.prj
index f4ab25d7714f889d9ff50ba677403a76371364da..5756e71d8305986020c47d663a0454d6b41232cf 100644
--- a/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_top_bottom.prj
+++ b/Tests/Data/Elliptic/cube_1x1x1_SteadyStateDiffusion/cube_2e3_prism_surfaceflux_top_bottom.prj
@@ -47,6 +47,7 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v      </variable>
+                <variable>specific_flux</variable>
             </variables>
             <meshes>
                 <mesh>cube_1x1x1_prism_2e3</mesh>
diff --git a/Tests/Data/Elliptic/wedge_1x1x1_SteadyStateDiffusion/wedge_1e3_prism_surfaceflux_diagonal.prj b/Tests/Data/Elliptic/wedge_1x1x1_SteadyStateDiffusion/wedge_1e3_prism_surfaceflux_diagonal.prj
index ebcb68f548dff2599987d24531e2d5f58953c5ac..4581860c3fea332f4e4e73ae7f95749588c403af 100644
--- a/Tests/Data/Elliptic/wedge_1x1x1_SteadyStateDiffusion/wedge_1e3_prism_surfaceflux_diagonal.prj
+++ b/Tests/Data/Elliptic/wedge_1x1x1_SteadyStateDiffusion/wedge_1e3_prism_surfaceflux_diagonal.prj
@@ -47,6 +47,7 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v      </variable>
+                <variable>specific_flux</variable>
             </variables>
             <meshes>
                 <mesh>wedge_1x1x1_1e3_prism</mesh>
diff --git a/Tests/Data/HydroMechanics/A2/A2.prj b/Tests/Data/HydroMechanics/A2/A2.prj
index 5bc8aad234713f86716214d1f14401ee8bb4890b..535f14dbb7fd533f724c4cf16a4316d3e4dc099b 100644
--- a/Tests/Data/HydroMechanics/A2/A2.prj
+++ b/Tests/Data/HydroMechanics/A2/A2.prj
@@ -222,6 +222,8 @@
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/HydroMechanics/ExcavationNiches/excavation_niches.prj b/Tests/Data/HydroMechanics/ExcavationNiches/excavation_niches.prj
index 51ba724810c4f252ba7ccc8638a47bc08da4e5c2..96356fd85c8cb7fb1a9db5f0ca97d7bf66ab801c 100644
--- a/Tests/Data/HydroMechanics/ExcavationNiches/excavation_niches.prj
+++ b/Tests/Data/HydroMechanics/ExcavationNiches/excavation_niches.prj
@@ -125,6 +125,10 @@
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>displacement_active</variable>
+                <variable>pressure_active</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/HydroMechanics/ExcavationNiches/excavation_niches2.prj b/Tests/Data/HydroMechanics/ExcavationNiches/excavation_niches2.prj
index bb0c0ba2acdbde16d6721f31f39fc447ff222a16..2c902ba37ef4e50ad981470be1ab8485d531a25f 100644
--- a/Tests/Data/HydroMechanics/ExcavationNiches/excavation_niches2.prj
+++ b/Tests/Data/HydroMechanics/ExcavationNiches/excavation_niches2.prj
@@ -172,6 +172,10 @@
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>displacement_active</variable>
+                <variable>pressure_active</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/HydroMechanics/FailureIndexDependentPermeability/quad_with_half_hole.prj b/Tests/Data/HydroMechanics/FailureIndexDependentPermeability/quad_with_half_hole.prj
index bbff892cc9354e3801006812b558e7ebacd76895..67ebcf53f1bce53bcc1115c0b7e26983257548ef 100644
--- a/Tests/Data/HydroMechanics/FailureIndexDependentPermeability/quad_with_half_hole.prj
+++ b/Tests/Data/HydroMechanics/FailureIndexDependentPermeability/quad_with_half_hole.prj
@@ -69,6 +69,7 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>permeability</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
diff --git a/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground.prj b/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground.prj
index ca7b58305f95d03aeb82e010a3330c851f6dba9a..142d9608dc34bba86a5d2971aac9cdd371295dd9 100644
--- a/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground.prj
+++ b/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground.prj
@@ -265,6 +265,8 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>sigma</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>DarcyVelocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_python.prj b/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_python.prj
index d8dbc38c08640d3b3a75419596bfeaa69011aa01..3aae782a0ba4e3be7dec7e7027fe38b08bc183be 100644
--- a/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_python.prj
+++ b/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_python.prj
@@ -266,6 +266,8 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>sigma</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>DarcyVelocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_quadBCu.prj b/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_quadBCu.prj
index 1a6ad07013015c1ef4592652b0a937b3a9c44285..4e39065b5b3fc815a3bf879b539d81dba01e0ea5 100644
--- a/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_quadBCu.prj
+++ b/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_quadBCu.prj
@@ -275,6 +275,8 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>sigma</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>DarcyVelocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_quadBCu_python.prj b/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_quadBCu_python.prj
index eb0dfa2795ab76df875390971dc7f501a5ba8653..711559e8f2635255db5633e89bbd09d3bee73f87 100644
--- a/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_quadBCu_python.prj
+++ b/Tests/Data/HydroMechanics/GroundEquilibrium/simHM_ground_quadBCu_python.prj
@@ -279,6 +279,8 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>sigma</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>DarcyVelocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/HydroMechanics/Linear/Confined_Compression/square_1e2.prj b/Tests/Data/HydroMechanics/Linear/Confined_Compression/square_1e2.prj
index 43d12156401334c99d18d70f5ef6e15dbf3b49e7..66fa4cf72461fe60c9e176b809601e71a9c0f5f9 100644
--- a/Tests/Data/HydroMechanics/Linear/Confined_Compression/square_1e2.prj
+++ b/Tests/Data/HydroMechanics/Linear/Confined_Compression/square_1e2.prj
@@ -134,6 +134,9 @@
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>pressure_interpolated</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/HydroMechanics/Linear/Confined_Compression/square_1e2_linear.prj b/Tests/Data/HydroMechanics/Linear/Confined_Compression/square_1e2_linear.prj
index 9730e54326e76c155a350d0b99c9912930da519f..f12b75096deb76811d7e79ea3c1178277aa4c337 100644
--- a/Tests/Data/HydroMechanics/Linear/Confined_Compression/square_1e2_linear.prj
+++ b/Tests/Data/HydroMechanics/Linear/Confined_Compression/square_1e2_linear.prj
@@ -134,6 +134,9 @@
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>pressure_interpolated</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_anisotropic.prj b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_anisotropic.prj
index 14a39f22bb1c1bc1bd0b3748f844d9c15516467a..2f46a88d04d3e2c288ed8ff91cab77021140533b 100644
--- a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_anisotropic.prj
+++ b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_anisotropic.prj
@@ -38,6 +38,7 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>permeability</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
diff --git a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_anisotropic_rotated.prj b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_anisotropic_rotated.prj
index 10058b5a7d0fb291b6110fc6dee131bbdc5b582b..e50a5dcf5af83a580befd2af9af5c0d8b68f12bd 100644
--- a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_anisotropic_rotated.prj
+++ b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_anisotropic_rotated.prj
@@ -38,6 +38,7 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>permeability</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
diff --git a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_quasiisotropic.prj b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_quasiisotropic.prj
index a88e05af3e4c0c6922d7958526d2f1736cac53a2..8357d391be69d686bac708d0c4adb81427ee6f1c 100644
--- a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_quasiisotropic.prj
+++ b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/disc_with_hole_quasiisotropic.prj
@@ -38,6 +38,7 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>permeability</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
diff --git a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/unconfined_biot.prj b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/unconfined_biot.prj
index 918b8fa167b7c46260ef2d0d5f537f4f07fdfc61..c813f3bbbfa2f355909e6fcbd1d0c4035178dfcd 100644
--- a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/unconfined_biot.prj
+++ b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/unconfined_biot.prj
@@ -121,6 +121,7 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>permeability</variable>
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
                 <variable>sigma_zz</variable>
@@ -129,6 +130,7 @@
                 <variable>epsilon_yy</variable>
                 <variable>epsilon_zz</variable>
                 <variable>epsilon_xy</variable>
+                <variable>epsilon</variable>
                 <variable>velocity</variable>
             </variables>
         </output>
diff --git a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/x_strain_y_flow.prj b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/x_strain_y_flow.prj
index b5a64de135316163e019c6b59066e001d5369557..70270cf014b2a24e6f47f48386ec6c6e1ee03256 100644
--- a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/x_strain_y_flow.prj
+++ b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/x_strain_y_flow.prj
@@ -38,10 +38,12 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>permeability</variable>
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
                 <variable>sigma_zz</variable>
                 <variable>sigma_xy</variable>
+                <variable>epsilon</variable>
                 <variable>epsilon_xx</variable>
                 <variable>epsilon_yy</variable>
                 <variable>epsilon_zz</variable>
diff --git a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/y_strain_z_flow.prj b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/y_strain_z_flow.prj
index 9262762e8aa2d3351e18a906027a20144f46a7d2..73fc2f37b8f9ccc9b2cd3091207042749e05ce82 100644
--- a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/y_strain_z_flow.prj
+++ b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/y_strain_z_flow.prj
@@ -38,10 +38,12 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>permeability</variable>
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
                 <variable>sigma_zz</variable>
                 <variable>sigma_xy</variable>
+                <variable>epsilon</variable>
                 <variable>epsilon_xx</variable>
                 <variable>epsilon_yy</variable>
                 <variable>epsilon_zz</variable>
diff --git a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/z_strain_x_flow.prj b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/z_strain_x_flow.prj
index d28e656b5a2157204dfcb2056458d1c9317a6b7f..9ac6d0d13ecc2bd917b7862db18d272eca397b3f 100644
--- a/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/z_strain_x_flow.prj
+++ b/Tests/Data/HydroMechanics/OrthotropicEmbeddedFracturePermeability/z_strain_x_flow.prj
@@ -38,6 +38,7 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>permeability</variable>
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
                 <variable>sigma_zz</variable>
@@ -46,6 +47,7 @@
                 <variable>epsilon_yy</variable>
                 <variable>epsilon_zz</variable>
                 <variable>epsilon_xy</variable>
+                <variable>epsilon</variable>
                 <variable>velocity</variable>
             </variables>
         </output>
diff --git a/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/MultiMesh/drainage.prj b/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/MultiMesh/drainage.prj
index 2003e12615a73afddc9bb543e48471b9dff476c7..5b3cbe0f4ce9ca3cc198631eeab2e7bb049c7ede 100644
--- a/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/MultiMesh/drainage.prj
+++ b/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/MultiMesh/drainage.prj
@@ -125,6 +125,7 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>q</variable>
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
diff --git a/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/MultiMesh/drainage_staggered.prj b/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/MultiMesh/drainage_staggered.prj
index 43bb9349c228a0f4254d74b3131b3dc54accd068..37ad0609f2c3e40af999329ea30d89d9739ae858 100644
--- a/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/MultiMesh/drainage_staggered.prj
+++ b/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/MultiMesh/drainage_staggered.prj
@@ -168,6 +168,7 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>q</variable>
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
diff --git a/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/SingleMesh/drainage.prj b/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/SingleMesh/drainage.prj
index db1915e7a12b1b371e4d690a31d4d5e72b401798..db5f63d8005a409227f8870f95ac127db244ebaa 100644
--- a/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/SingleMesh/drainage.prj
+++ b/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/SingleMesh/drainage.prj
@@ -120,6 +120,7 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>q</variable>
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
diff --git a/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/SingleMesh/drainage_staggered.prj b/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/SingleMesh/drainage_staggered.prj
index e8c52d47823b255d61c03f68c31ed7f59e62397b..b4f1103dcad1fc5ef26f1423cafcfe4c5103a25a 100644
--- a/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/SingleMesh/drainage_staggered.prj
+++ b/Tests/Data/HydroMechanics/ParallelComputing/SimpleHM/SingleMesh/drainage_staggered.prj
@@ -163,6 +163,7 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>q</variable>
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
diff --git a/Tests/Data/HydroMechanics/Principal_Stress/Hollow_Sphere/sphere.prj b/Tests/Data/HydroMechanics/Principal_Stress/Hollow_Sphere/sphere.prj
index 8b6db25f6cc1038174d80c8440cf29d7eee3603d..334c9b11843558f2d79f2af863f1885310a9fbdd 100644
--- a/Tests/Data/HydroMechanics/Principal_Stress/Hollow_Sphere/sphere.prj
+++ b/Tests/Data/HydroMechanics/Principal_Stress/Hollow_Sphere/sphere.prj
@@ -70,6 +70,7 @@
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
+                <variable>principal_stress_values</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/HydroMechanics/Principal_Stress/Tube/tube.prj b/Tests/Data/HydroMechanics/Principal_Stress/Tube/tube.prj
index 83736ae6e245711ccaf4dafcd7f83a13d2202e7e..0a7fee40b2414f8818fd2d907010002d52b7e2dd 100644
--- a/Tests/Data/HydroMechanics/Principal_Stress/Tube/tube.prj
+++ b/Tests/Data/HydroMechanics/Principal_Stress/Tube/tube.prj
@@ -70,6 +70,7 @@
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
+                <variable>principal_stress_values</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/HydroMechanics/SeabedResponse/seabed_response_200x100.prj b/Tests/Data/HydroMechanics/SeabedResponse/seabed_response_200x100.prj
index 7a0c24620195861c4059a44fae0d7a121a49dc62..ec59913c79ddb5002027e61c246f8aeee85d7c14 100644
--- a/Tests/Data/HydroMechanics/SeabedResponse/seabed_response_200x100.prj
+++ b/Tests/Data/HydroMechanics/SeabedResponse/seabed_response_200x100.prj
@@ -118,6 +118,8 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>sigma</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/HydroMechanics/StrainDependentPermeability/gas_loading.prj b/Tests/Data/HydroMechanics/StrainDependentPermeability/gas_loading.prj
index f606389e80dd4cb8a91718b6fb8431188937e4e6..4a938ac6b653e8292cfb6dc2e85cd7f16ea85d71 100644
--- a/Tests/Data/HydroMechanics/StrainDependentPermeability/gas_loading.prj
+++ b/Tests/Data/HydroMechanics/StrainDependentPermeability/gas_loading.prj
@@ -102,6 +102,7 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>permeability</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
diff --git a/Tests/Data/LIE/HydroMechanics/TaskB.prj b/Tests/Data/LIE/HydroMechanics/TaskB.prj
index de6faf8211c296eb789f75726431d2d25dfca847..6a8b7691105e1666f3889087b78630bfd2312419 100644
--- a/Tests/Data/LIE/HydroMechanics/TaskB.prj
+++ b/Tests/Data/LIE/HydroMechanics/TaskB.prj
@@ -90,8 +90,20 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>nodal_w</variable>
+                <variable>nodal_aperture</variable>
+                <variable>strain_xx</variable>
+                <variable>strain_yy</variable>
+                <variable>strain_xy</variable>
+                <variable>stress_xx</variable>
+                <variable>stress_yy</variable>
+                <variable>velocity</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>NodalForcesJump</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture.prj b/Tests/Data/LIE/HydroMechanics/single_fracture.prj
index 87ac762f823b2dfc8e68b1f8d245ac2ad6e2e1d0..e0ae4e5fa2cc852cf4c2a355a586e4a886006ade 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture.prj
@@ -84,8 +84,20 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>nodal_w</variable>
+                <variable>nodal_aperture</variable>
+                <variable>strain_xx</variable>
+                <variable>strain_yy</variable>
+                <variable>strain_xy</variable>
+                <variable>stress_xx</variable>
+                <variable>stress_yy</variable>
+                <variable>velocity</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>NodalForcesJump</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_3D.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_3D.prj
index e1ca893fefc94319fa6289ff438199f63ae4e924..2e9b0d4884c2bfbe9855a40b187f406ca21e5112 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture_3D.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_3D.prj
@@ -84,8 +84,27 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>nodal_w</variable>
+                <variable>nodal_aperture</variable>
+                <variable>strain_xx</variable>
+                <variable>strain_yy</variable>
+                <variable>strain_zz</variable>
+                <variable>strain_xy</variable>
+                <variable>strain_xz</variable>
+                <variable>strain_yz</variable>
+                <variable>stress_xx</variable>
+                <variable>stress_yy</variable>
+                <variable>stress_zz</variable>
+                <variable>stress_xy</variable>
+                <variable>stress_xz</variable>
+                <variable>stress_yz</variable>
+                <variable>velocity</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>NodalForcesJump</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow.prj
index 65ba2022d03d526bd75cb0cafc41e890563fc322..29d301421ada472053265d3febda189137c8485d 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow.prj
@@ -93,8 +93,17 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>nodal_w</variable>
+                <variable>nodal_aperture</variable>
+                <variable>strain_xx</variable>
+                <variable>strain_yy</variable>
+                <variable>strain_xy</variable>
+                <variable>stress_xx</variable>
+                <variable>stress_yy</variable>
+                <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ.prj
index 5116bd181f988b99c0a5dc58c4990ea3c144b873..7927b24d43657cf47cfe54ac32632e9c71648387 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ.prj
@@ -102,8 +102,17 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>nodal_w</variable>
+                <variable>nodal_aperture</variable>
+                <variable>strain_xx</variable>
+                <variable>strain_yy</variable>
+                <variable>strain_xy</variable>
+                <variable>stress_xx</variable>
+                <variable>stress_yy</variable>
+                <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ_sigma0.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ_sigma0.prj
index 960c1c572b54721aab19fa25a63a32ad722fe4bb..3a229d011ea5bd2a5d7299f85b9216059b30c922 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ_sigma0.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_CHZ_sigma0.prj
@@ -102,8 +102,17 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>nodal_w</variable>
+                <variable>nodal_aperture</variable>
+                <variable>strain_xx</variable>
+                <variable>strain_yy</variable>
+                <variable>strain_xy</variable>
+                <variable>stress_xx</variable>
+                <variable>stress_yy</variable>
+                <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_linear_aperture0.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_linear_aperture0.prj
index da1dc33c9759da89a998bc3fbf42809606b0506e..ca154824d4f9de111c6279a329acf85d93cf7311 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_linear_aperture0.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_linear_aperture0.prj
@@ -93,8 +93,18 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>nodal_w</variable>
+                <variable>nodal_aperture</variable>
+                <variable>strain_xx</variable>
+                <variable>strain_yy</variable>
+                <variable>strain_xy</variable>
+                <variable>stress_xx</variable>
+                <variable>stress_yy</variable>
+                <variable>stress_xy</variable>
+                <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_linear_aperture0_e.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_linear_aperture0_e.prj
index efcf2bdc646235b55efdf27f00a3e0841de33f8c..72443a68c92465b7e4063b991cb599e4b490a609 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_linear_aperture0_e.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_3compartments_flow_linear_aperture0_e.prj
@@ -93,8 +93,18 @@
             </timesteps>
             <variables>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>nodal_w</variable>
+                <variable>nodal_aperture</variable>
+                <variable>strain_xx</variable>
+                <variable>strain_yy</variable>
+                <variable>strain_xy</variable>
+                <variable>stress_xx</variable>
+                <variable>stress_yy</variable>
+                <variable>stress_xy</variable>
+                <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_CZ_kf_const.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_CZ_kf_const.prj
index f2e77b6e2d93d23bc7606b78e55c2ab040bf917a..9174e8add5dbda804f74f33d35a45b9329d97f4b 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture_CZ_kf_const.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_CZ_kf_const.prj
@@ -102,6 +102,9 @@
                 <variable>pressure</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>k_f</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_CZ_kf_cubic.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_CZ_kf_cubic.prj
index 44fb73e44084f6301830cd24949655b52970229c..4c16d42e2bc1ec6b93e65409cf8864d54b57da15 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture_CZ_kf_cubic.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_CZ_kf_cubic.prj
@@ -101,6 +101,9 @@
                 <variable>pressure</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>k_f</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_LE_kf_const.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_LE_kf_const.prj
index 3bf845c8e4c189c41f8aa4f7dfed9e22833454b9..44427cd4df1ba6dbbc2d4817bb068b963f6aa0d1 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture_LE_kf_const.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_LE_kf_const.prj
@@ -100,6 +100,9 @@
                 <variable>pressure</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>k_f</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_LE_kf_cubic.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_LE_kf_cubic.prj
index cb98d6231f5f81d99ba378b96602260891c8a1b0..2e526106483f8dba7f297ca6ccbdd60b7049dbba 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture_LE_kf_cubic.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_LE_kf_cubic.prj
@@ -99,6 +99,9 @@
                 <variable>pressure</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>k_f</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/HydroMechanics/single_fracture_constK.prj b/Tests/Data/LIE/HydroMechanics/single_fracture_constK.prj
index 0fca2f17057b749f5b71c3a341498c077098480d..4b1a9ec5693b8075bc9274b6cc1a3f144605dde5 100644
--- a/Tests/Data/LIE/HydroMechanics/single_fracture_constK.prj
+++ b/Tests/Data/LIE/HydroMechanics/single_fracture_constK.prj
@@ -87,6 +87,18 @@
                 <variable>pressure</variable>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>nodal_aperture</variable>
+                <variable>nodal_w</variable>
+                <variable>pressure_interpolated</variable>
+                <variable>strain_xx</variable>
+                <variable>strain_xy</variable>
+                <variable>strain_yy</variable>
+                <variable>strain_zz</variable>
+                <variable>stress_xx</variable>
+                <variable>stress_xy</variable>
+                <variable>stress_yy</variable>
+                <variable>stress_zz</variable>
+                <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.prj b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.prj
index 3d43c4e6e59b53e2affd8d170e408576b03584f5..76b03992ff957f98b07710ccd388eb1b63ff3530 100644
--- a/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.prj
+++ b/Tests/Data/LIE/Mechanics/cohesive_zone_load_path.prj
@@ -79,6 +79,7 @@
                 <variable>sigma_yy</variable>
                 <variable>sigma_zz</variable>
                 <variable>sigma_xy</variable>
+                <variable>f_stress_n</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/Mechanics/coulomb_load_path.prj b/Tests/Data/LIE/Mechanics/coulomb_load_path.prj
index cfaeda08190b18aaf6f97629960c1702209f9fba..1ec9bdf84ca1bb8b7a6717562e27afa3810c018b 100644
--- a/Tests/Data/LIE/Mechanics/coulomb_load_path.prj
+++ b/Tests/Data/LIE/Mechanics/coulomb_load_path.prj
@@ -83,6 +83,7 @@
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
                 <variable>sigma_xy</variable>
+                <variable>aperture</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures.prj b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures.prj
index 554a4b1303b6ffed1e513c9e9771314f54abfb1b..53600927925a45a4629ee2defffbe94d7339aec1 100644
--- a/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures.prj
+++ b/Tests/Data/LIE/Mechanics/elastic_push_pull_two_fractures.prj
@@ -82,6 +82,8 @@
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
                 <variable>sigma_xy</variable>
+                <variable>f_stress_n</variable>
+                <variable>f_stress_s</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/LIE/Mechanics/mohr_coulomb_load_path_nu0p3.prj b/Tests/Data/LIE/Mechanics/mohr_coulomb_load_path_nu0p3.prj
index 7a6e2dca43ce39ad5dc47da245b213c263031257..73b0731b64ebb74b058f3e2854f664ce30f3432c 100644
--- a/Tests/Data/LIE/Mechanics/mohr_coulomb_load_path_nu0p3.prj
+++ b/Tests/Data/LIE/Mechanics/mohr_coulomb_load_path_nu0p3.prj
@@ -83,6 +83,8 @@
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
                 <variable>sigma_xy</variable>
+                <variable>f_stress_s</variable>
+                <variable>stress_xy</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/Mechanics/sfrac.q.prj b/Tests/Data/LIE/Mechanics/sfrac.q.prj
index 0ef5c242b52c4e095cde5f23b4eeba0253ba5655..e18bd8d7fcbae6e1064cc595825c4b47c26dfdfd 100644
--- a/Tests/Data/LIE/Mechanics/sfrac.q.prj
+++ b/Tests/Data/LIE/Mechanics/sfrac.q.prj
@@ -88,6 +88,7 @@
                 <variable>displacement_jump2</variable>
                 <variable>displacement_jump3</variable>
                 <variable>displacement_jump4</variable>
+                <variable>f_stress_n</variable>
                 <!--
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
diff --git a/Tests/Data/LIE/Mechanics/single_joint_displacement_controlled.prj b/Tests/Data/LIE/Mechanics/single_joint_displacement_controlled.prj
index 39b42df974319eb9e3fc9c39e99d381db0bee321..1f4c2b7b41f9c56467910ccd1add57fa3514d61b 100644
--- a/Tests/Data/LIE/Mechanics/single_joint_displacement_controlled.prj
+++ b/Tests/Data/LIE/Mechanics/single_joint_displacement_controlled.prj
@@ -80,6 +80,7 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>aperture</variable>
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
                 <variable>sigma_xy</variable>
diff --git a/Tests/Data/LIE/Mechanics/single_joint_negative_aperture.prj b/Tests/Data/LIE/Mechanics/single_joint_negative_aperture.prj
index ff000ed23bc025c4fbcd6d1ccd78c9a2d3ee1fc0..c5767a8a54c1a9e23898464af57ef229e9365bfb 100644
--- a/Tests/Data/LIE/Mechanics/single_joint_negative_aperture.prj
+++ b/Tests/Data/LIE/Mechanics/single_joint_negative_aperture.prj
@@ -72,6 +72,7 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>displacement_jump1</variable>
+                <variable>aperture</variable>
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
                 <variable>sigma_xy</variable>
diff --git a/Tests/Data/LIE/Mechanics/two_cracks_branch_pull.prj b/Tests/Data/LIE/Mechanics/two_cracks_branch_pull.prj
index 8a0be9b1a1870e841625c95a846a28ce0bfe87d9..fe6293d00e1ce0de606f2ff796054666e96939a0 100644
--- a/Tests/Data/LIE/Mechanics/two_cracks_branch_pull.prj
+++ b/Tests/Data/LIE/Mechanics/two_cracks_branch_pull.prj
@@ -81,6 +81,7 @@
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
                 <variable>sigma_xy</variable>
+                <variable>f_stress_n</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/LIE/Mechanics/two_cracks_junction_pull.prj b/Tests/Data/LIE/Mechanics/two_cracks_junction_pull.prj
index 1959920ff514cc347bd2cd1370470af6e78ef145..4c5d7b57b33933a9bc2f21c1921044ba95cd6a00 100644
--- a/Tests/Data/LIE/Mechanics/two_cracks_junction_pull.prj
+++ b/Tests/Data/LIE/Mechanics/two_cracks_junction_pull.prj
@@ -83,6 +83,7 @@
                 <variable>sigma_xx</variable>
                 <variable>sigma_yy</variable>
                 <variable>sigma_xy</variable>
+                <variable>f_stress_n</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Ehlers/MFront/square_1e1_2_matIDs.prj b/Tests/Data/Mechanics/Ehlers/MFront/square_1e1_2_matIDs.prj
index 918ae3836ed61202ea9464d68901b4e2362214cd..f71bebbd7fa423d93e4ca20e817b1c408049a0c5 100644
--- a/Tests/Data/Mechanics/Ehlers/MFront/square_1e1_2_matIDs.prj
+++ b/Tests/Data/Mechanics/Ehlers/MFront/square_1e1_2_matIDs.prj
@@ -106,6 +106,8 @@
                 <variable>eps_p.D</variable>
                 <variable>eps_p.V</variable>
                 <variable>eps_p.eff</variable>
+                <variable>MaterialForces</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Mechanics/Ehlers/MFront/square_1e1_2_matIDs_restart.prj b/Tests/Data/Mechanics/Ehlers/MFront/square_1e1_2_matIDs_restart.prj
index a17eba3837c56aa37b1fbf7f8c8315e301178ac7..93538854808a26034ccffb69c07011797d5118a2 100644
--- a/Tests/Data/Mechanics/Ehlers/MFront/square_1e1_2_matIDs_restart.prj
+++ b/Tests/Data/Mechanics/Ehlers/MFront/square_1e1_2_matIDs_restart.prj
@@ -106,6 +106,8 @@
                 <variable>eps_p.D</variable>
                 <variable>eps_p.V</variable>
                 <variable>eps_p.eff</variable>
+                <variable>MaterialForces</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Mechanics/Ehlers/MFront/two_material_ids_single_solid.prj b/Tests/Data/Mechanics/Ehlers/MFront/two_material_ids_single_solid.prj
index 65eda141c427c042af67244c2225f3b2a753e36b..6ece6d08ebd51d7cdd1f5abbd851d482397fbba6 100644
--- a/Tests/Data/Mechanics/Ehlers/MFront/two_material_ids_single_solid.prj
+++ b/Tests/Data/Mechanics/Ehlers/MFront/two_material_ids_single_solid.prj
@@ -75,6 +75,8 @@
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>EquivalentPlasticStrain</variable>
+                <variable>MaterialForces</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_xyz.prj b/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_xyz.prj
index eaeb84496fe0f3707ab860b11ae18e636d9e17e6..0ec91df0f8fc3a15f5f122bd8bca9938b46dcf70 100644
--- a/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_xyz.prj
+++ b/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_xyz.prj
@@ -72,6 +72,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_yzx.prj b/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_yzx.prj
index df77b546b65f9ebf071d41e4349f663bcbaeb1ba..314a66fd942a7ce26cbc22ca46b5d0f0549b1792 100644
--- a/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_yzx.prj
+++ b/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_yzx.prj
@@ -72,6 +72,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_zxy.prj b/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_zxy.prj
index fbd6a4d38ea96a7eb59168ea6044d905e9837162..15c770a2291f0ba64362b212b98bcb3664e66f9b 100644
--- a/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_zxy.prj
+++ b/Tests/Data/Mechanics/Linear/MFront/cube_1e0_orthotropic_zxy.prj
@@ -72,6 +72,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_45xy_z.prj b/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_45xy_z.prj
index 402a2c51abc6dd8d4804b27220ca6beb42f3e18a..4cd74765a3d0e97a48e224fa5215ededab275137 100644
--- a/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_45xy_z.prj
+++ b/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_45xy_z.prj
@@ -70,6 +70,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_xyz.prj b/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_xyz.prj
index b0daf22fd4c880e0724c98a2de45a363513bf0f1..b9394ee2c66e6a7b50e3abecbb69efc0865060a8 100644
--- a/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_xyz.prj
+++ b/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_xyz.prj
@@ -70,6 +70,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_y-xz.prj b/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_y-xz.prj
index 361dcfb3da37cc12ad063604bafe5e39858015b9..6e972cd7af07caf0b2976726ade63f6654cf2a01 100644
--- a/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_y-xz.prj
+++ b/Tests/Data/Mechanics/Linear/MFront/square_1e0_orthotropic_y-xz.prj
@@ -70,6 +70,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/MaterialForces/bar.prj b/Tests/Data/Mechanics/Linear/MaterialForces/bar.prj
index 88d132d01ba4976737f067d0ad11f10f67407399..cc489786998185f9625047a507d79206417aae57 100644
--- a/Tests/Data/Mechanics/Linear/MaterialForces/bar.prj
+++ b/Tests/Data/Mechanics/Linear/MaterialForces/bar.prj
@@ -63,6 +63,7 @@
                 <variable>free_energy_density</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>MaterialForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/MaterialForces/bar_3D.prj b/Tests/Data/Mechanics/Linear/MaterialForces/bar_3D.prj
index 690a80c56fdc0a503d45ee29e8b33de0b92fcecd..2d94d9d5d62371a3bca82849c3617ee7ed87341f 100644
--- a/Tests/Data/Mechanics/Linear/MaterialForces/bar_3D.prj
+++ b/Tests/Data/Mechanics/Linear/MaterialForces/bar_3D.prj
@@ -63,6 +63,7 @@
                 <variable>free_energy_density</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>MaterialForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_xyz.prj b/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_xyz.prj
index 1cf3dc28832cc3d622bf30195313d718470e9360..59af96049794c878980f25799604d5eb36343cba 100644
--- a/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_xyz.prj
+++ b/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_xyz.prj
@@ -62,6 +62,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_yzx.prj b/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_yzx.prj
index dc1ebd7ae8bd04027a00ec4e90af9ccf1531bc17..56b085b89605af83a73f4185af2ea9c9ba010341 100644
--- a/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_yzx.prj
+++ b/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_yzx.prj
@@ -62,6 +62,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_zxy.prj b/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_zxy.prj
index c7625e32032ca2c8bcd5f5b89c2f95560e97485a..8f1d47624b6986218ee939003edd91c880e87b57 100644
--- a/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_zxy.prj
+++ b/Tests/Data/Mechanics/Linear/Orthotropy/cube_1e0_orthotropic_zxy.prj
@@ -62,6 +62,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_45xy_z.prj b/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_45xy_z.prj
index 4242203c02af8fa7f0368a0e174faf449140e356..52d9348444a98b54a8a0804f901f390a276bb0c0 100644
--- a/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_45xy_z.prj
+++ b/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_45xy_z.prj
@@ -62,6 +62,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_xyz.prj b/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_xyz.prj
index ea5767c859f981c90f8fa62e2e4bd4b8a314b61f..fc4d48c86be9fc19fb6fdef6d3701b423010564c 100644
--- a/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_xyz.prj
+++ b/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_xyz.prj
@@ -62,6 +62,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_y-xz.prj b/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_y-xz.prj
index 8c7250f704020f26456fd5337ad8bdac9b92b5e1..6fd9e38f76fac3baa35b78fb5cbdebd153c9fcfa 100644
--- a/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_y-xz.prj
+++ b/Tests/Data/Mechanics/Linear/Orthotropy/square_1e0_orthotropic_y-xz.prj
@@ -62,6 +62,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/PrincipalStress/sphere.prj b/Tests/Data/Mechanics/Linear/PrincipalStress/sphere.prj
index d658cf04dc625bd68c2670595004432febf1c2a5..b21a564ff70e247c9a2f5825f43afaacda526d12 100644
--- a/Tests/Data/Mechanics/Linear/PrincipalStress/sphere.prj
+++ b/Tests/Data/Mechanics/Linear/PrincipalStress/sphere.prj
@@ -67,6 +67,7 @@
                 <variable>pressure</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>principal_stress_values</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/cube_1e0.prj b/Tests/Data/Mechanics/Linear/cube_1e0.prj
index 5439c3e016f81842ab2db2457eac34e3381070eb..92de2bacb7deb659c95bfd90684d7ecf6edd156c 100644
--- a/Tests/Data/Mechanics/Linear/cube_1e0.prj
+++ b/Tests/Data/Mechanics/Linear/cube_1e0.prj
@@ -61,6 +61,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/cube_1e0_SNES.prj b/Tests/Data/Mechanics/Linear/cube_1e0_SNES.prj
index 649f8cad505362f27c5b036849722ae1f5ce9343..295f7f8bb91ba6c7cbac33b11a91fd84c5375f95 100644
--- a/Tests/Data/Mechanics/Linear/cube_1e0_SNES.prj
+++ b/Tests/Data/Mechanics/Linear/cube_1e0_SNES.prj
@@ -61,6 +61,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/Linear/cube_1e0_simple_shear.prj b/Tests/Data/Mechanics/Linear/cube_1e0_simple_shear.prj
index bf37fdbf731175677cf2070168e52d01fe5e200a..ec811e911a8a508035c108bb0db244341b4087c8 100644
--- a/Tests/Data/Mechanics/Linear/cube_1e0_simple_shear.prj
+++ b/Tests/Data/Mechanics/Linear/cube_1e0_simple_shear.prj
@@ -61,6 +61,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/ModifiedCamClay/model_triaxtest.prj b/Tests/Data/Mechanics/ModifiedCamClay/model_triaxtest.prj
index aca1f614fa8676bced45332a49bfe778b6ad6462..253aeae50a8d53a681d8318a9eac3b89363b0c88 100644
--- a/Tests/Data/Mechanics/ModifiedCamClay/model_triaxtest.prj
+++ b/Tests/Data/Mechanics/ModifiedCamClay/model_triaxtest.prj
@@ -106,7 +106,6 @@
                         </pair>
                     </timesteps>
                 </time_stepping>
-
             </process>
         </processes>
         <output>
@@ -122,6 +121,7 @@
                 <variable>displacement</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>PreConsolidationPressure</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Mechanics/ModifiedCamClay/square_1e0_biax.prj b/Tests/Data/Mechanics/ModifiedCamClay/square_1e0_biax.prj
index 15a09fedabe6ebbbf1f174771eaea5c116ce2419..098f19b084ecfd492aee686fa52c4a5a812166c6 100644
--- a/Tests/Data/Mechanics/ModifiedCamClay/square_1e0_biax.prj
+++ b/Tests/Data/Mechanics/ModifiedCamClay/square_1e0_biax.prj
@@ -74,6 +74,7 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>sigma</variable>
+                <variable>PreConsolidationPressure</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/ModifiedCamClay/square_1e0_shear.prj b/Tests/Data/Mechanics/ModifiedCamClay/square_1e0_shear.prj
index 6ad52e168e32a591e7c1bc50c418790b1ae2b84d..5f7072ae1641a483bb370470d3447b956c1e9457 100644
--- a/Tests/Data/Mechanics/ModifiedCamClay/square_1e0_shear.prj
+++ b/Tests/Data/Mechanics/ModifiedCamClay/square_1e0_shear.prj
@@ -75,6 +75,7 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>sigma</variable>
+                <variable>PreConsolidationPressure</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/MohrCoulombAbboSloan/load_test_mc.prj b/Tests/Data/Mechanics/MohrCoulombAbboSloan/load_test_mc.prj
index 3c094297ac8fd0a5eba050e9b211ea7620885ab3..604c200f2f3cef2a3ad0356b680a2487870c6e21 100644
--- a/Tests/Data/Mechanics/MohrCoulombAbboSloan/load_test_mc.prj
+++ b/Tests/Data/Mechanics/MohrCoulombAbboSloan/load_test_mc.prj
@@ -78,6 +78,7 @@
                 <variable>epsilon</variable>
                 <variable>ElasticStrain</variable>
                 <variable>EquivalentPlasticStrain</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/MohrCoulombAbboSloan/oedometer.prj b/Tests/Data/Mechanics/MohrCoulombAbboSloan/oedometer.prj
index 070e6538feab827c3114045cd22173145d723e13..42f5b33cf449ff191fad10be4276af9f23c3b42b 100644
--- a/Tests/Data/Mechanics/MohrCoulombAbboSloan/oedometer.prj
+++ b/Tests/Data/Mechanics/MohrCoulombAbboSloan/oedometer.prj
@@ -74,6 +74,11 @@
             </timesteps>
             <variables>
                 <variable>displacement</variable>
+                <variable>ElasticStrain</variable>
+                <variable>EquivalentPlasticStrain</variable>
+                <variable>NodalForces</variable>
+                <variable>epsilon</variable>
+                <variable>sigma</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/MohrCoulombAbboSloan/slope.prj b/Tests/Data/Mechanics/MohrCoulombAbboSloan/slope.prj
index b25f8278a0449c1c5b1284ba46438e019215c4a5..faa1fe6856f66ac17279726a564be2b41da7b0c2 100644
--- a/Tests/Data/Mechanics/MohrCoulombAbboSloan/slope.prj
+++ b/Tests/Data/Mechanics/MohrCoulombAbboSloan/slope.prj
@@ -73,6 +73,7 @@
                 <variable>epsilon</variable>
                 <variable>ElasticStrain</variable>
                 <variable>EquivalentPlasticStrain</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Mechanics/PLLC/uniax_compression.prj b/Tests/Data/Mechanics/PLLC/uniax_compression.prj
index 56e422837acb037f69933cb431cf5e34d4d0c57a..44cf17c8750882ff4cd8bd060440533df4628e82 100644
--- a/Tests/Data/Mechanics/PLLC/uniax_compression.prj
+++ b/Tests/Data/Mechanics/PLLC/uniax_compression.prj
@@ -74,6 +74,7 @@
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>ElasticStrain</variable>
+                <variable>NodalForces</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay.prj b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay.prj
index 235b2f14b2e20bcf2d0fc0b13b2bb61d8cadf261..d1afc2d7c55e825813ccf1aa10228b6154d4e8ec 100644
--- a/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/AdvectionDiffusionSorptionDecay/1D_AdvectionDiffusionSorptionDecay.prj
@@ -172,7 +172,8 @@
             </timesteps>
             <variables>
                 <variable>Cs</variable>
-                <variable>pressure</variable>
+                <variable>LiquidMassFlowRate</variable>
+                <variable>CsFlowRate</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ConTracer/ConTracer_1d.prj b/Tests/Data/Parabolic/ComponentTransport/ConTracer/ConTracer_1d.prj
index b4334b8480bf048840f9409d24e15b638de3df19..3284154bcc34fa82c19ad371e54a293c017dc7a5 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ConTracer/ConTracer_1d.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/ConTracer/ConTracer_1d.prj
@@ -165,6 +165,8 @@
             <variables>
                 <variable>Cs</variable>
                 <variable>pressure</variable>
+                <variable>molar_flux</variable>
+                <variable>CsFlowRate</variable>
                 <variable>darcy_velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
diff --git a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay.prj b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay.prj
index 482e6778f7c2bb6148ee62a43e5fdc57c98f1dcd..ea67fec379a6eb8931d0cc63c1063e8c2cd4406c 100644
--- a/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/DiffusionSorptionDecay/1D_DiffusionSorptionDecay.prj
@@ -173,6 +173,7 @@
             <variables>
                 <variable>Cs</variable>
                 <variable>pressure</variable>
+                <variable>CsFlowRate</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/ComponentTransport/MultiLayerDiffusion/1D_MultiLayerDiffusion.prj b/Tests/Data/Parabolic/ComponentTransport/MultiLayerDiffusion/1D_MultiLayerDiffusion.prj
index 773faddb7d06e809cd8934d48b7a45c77648328d..633752332524aebb3347b7a71f6f9ebd00c09fea 100644
--- a/Tests/Data/Parabolic/ComponentTransport/MultiLayerDiffusion/1D_MultiLayerDiffusion.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/MultiLayerDiffusion/1D_MultiLayerDiffusion.prj
@@ -238,7 +238,9 @@
             </timesteps>
             <variables>
                 <variable>HTO</variable>
+                <variable>HTOFlowRate</variable>
                 <variable>pressure</variable>
+                <variable>molar_flux</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/DecayChain/GlobalImplicitApproach/1d_decay_chain_GIA.prj b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/DecayChain/GlobalImplicitApproach/1d_decay_chain_GIA.prj
index 872c37de0947a84e74f159657f2d570878aecb2a..2b06c6c8a4cd6ace85ad4eb30bfc1fe514a442be 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/DecayChain/GlobalImplicitApproach/1d_decay_chain_GIA.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/DecayChain/GlobalImplicitApproach/1d_decay_chain_GIA.prj
@@ -238,6 +238,12 @@
                 <variable>[Pa-231]</variable>
                 <variable>[Ac-227]</variable>
                 <variable>pressure</variable>
+                <variable>LiquidMassFlowRate</variable>
+                <variable>[Cm-247]FlowRate</variable>
+                <variable>[Pu-239]FlowRate</variable>
+                <variable>[U-235]FlowRate</variable>
+                <variable>[Pa-231]FlowRate</variable>
+                <variable>[Ac-227]FlowRate</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcite.prj b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcite.prj
index 1884ba5743a81cf6da5018cfdfe3361e27432046..3585064fdc4fe7b8a96600206fa111500b3bdf46 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcite.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcite.prj
@@ -396,6 +396,8 @@
                 <variable>H</variable>
                 <variable>pressure</variable>
                 <variable>darcy_velocity</variable>
+                <variable>Dolomite(dis)_avg</variable>
+                <variable>Calcite_avg</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calciteDissolvePrecipitateOnly.prj b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calciteDissolvePrecipitateOnly.prj
index 0ea5207012aa679ee67de586597a193b4fb9bc98..8580fa7cc0c6123ca2984b5646f3d40f671847d2 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calciteDissolvePrecipitateOnly.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calciteDissolvePrecipitateOnly.prj
@@ -396,6 +396,8 @@
                 <variable>H</variable>
                 <variable>pressure</variable>
                 <variable>darcy_velocity</variable>
+                <variable>Dolomite(dis)_avg</variable>
+                <variable>Calcite_avg</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange.prj b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange.prj
index 9c6fbbd59215f755881ee4eba71e16bde7d1ddff..7dc9c706fbd8430606651d1ced181546e36d9a4b 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/EquilibriumPhase/calcitePorosityChange.prj
@@ -407,6 +407,9 @@
                 <variable>H</variable>
                 <variable>pressure</variable>
                 <variable>darcy_velocity</variable>
+                <variable>Calcite_avg</variable>
+                <variable>Dolomite(dis)_avg</variable>
+                <variable>porosity_avg</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac.prj b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac.prj
index c590da110b1df0cf4a32083455a9b8c9872a3ba6..acbddd7fd077626cdea051e3b025df923f64133f 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac.prj
@@ -309,6 +309,7 @@
                 <variable>Syntheticb</variable>
                 <variable>pressure</variable>
                 <variable>darcy_velocity</variable>
+                <variable>Productc_avg</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_flag_formula.prj b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_flag_formula.prj
index 3de38238e98f53000589945ad80291873a4976c2..7cfffd3563144764a56ab8b297572edf7d68b267 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_flag_formula.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_flag_formula.prj
@@ -290,6 +290,7 @@
                 <variable>Syntheticb</variable>
                 <variable>pressure</variable>
                 <variable>darcy_velocity</variable>
+                <variable>porosity_avg</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_small_domain.prj b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_small_domain.prj
index 61796036cbfb08526eec7e80d9ad90fc4e2de879..67268df7733bd3dab6fafb9df78f640021b8d3e0 100644
--- a/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_small_domain.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/ReactiveTransport/KineticReactant/1d_isofrac_small_domain.prj
@@ -309,6 +309,7 @@
                 <variable>Syntheticb</variable>
                 <variable>pressure</variable>
                 <variable>darcy_velocity</variable>
+                <variable>Productc_avg</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/surfaceflux_component-transport_cube_1e3.prj b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/surfaceflux_component-transport_cube_1e3.prj
index 78c4a7960b164e65ffb6998bed4943d014ed1b77..b6d4901050dc317f2f6d4cd0add11626b0c658c8 100644
--- a/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/surfaceflux_component-transport_cube_1e3.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/SimpleSynthetics/surfaceflux_component-transport_cube_1e3.prj
@@ -140,6 +140,7 @@
                 <variable>Si</variable>
                 <variable>pressure</variable>
                 <variable>darcy_velocity</variable>
+                <variable>specific_flux</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/ComponentTransport/StaggeredScheme/surfaceflux_component-transport_cube_1e3.prj b/Tests/Data/Parabolic/ComponentTransport/StaggeredScheme/surfaceflux_component-transport_cube_1e3.prj
index aeb269e44bb52b8734d2ba78b45cf82c2ffe01e0..9b602a67817cc0f98c32833ccdd1ee55cb45d5f1 100644
--- a/Tests/Data/Parabolic/ComponentTransport/StaggeredScheme/surfaceflux_component-transport_cube_1e3.prj
+++ b/Tests/Data/Parabolic/ComponentTransport/StaggeredScheme/surfaceflux_component-transport_cube_1e3.prj
@@ -182,6 +182,7 @@
                 <variable>Si</variable>
                 <variable>pressure</variable>
                 <variable>darcy_velocity</variable>
+                <variable>specific_flux</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/calculatesurfaceflux_ht_cube_1e3.prj b/Tests/Data/Parabolic/HT/SimpleSynthetics/calculatesurfaceflux_ht_cube_1e3.prj
index 44f921f8d42be8f4cce729bb27ee8406a789f0d7..58e0f827732cdbb98e45171591f417709d0ed9ec 100644
--- a/Tests/Data/Parabolic/HT/SimpleSynthetics/calculatesurfaceflux_ht_cube_1e3.prj
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/calculatesurfaceflux_ht_cube_1e3.prj
@@ -166,6 +166,7 @@
                 <variable>T</variable>
                 <variable>p</variable>
                 <variable>darcy_velocity</variable>
+                <variable>specific_flux</variable>
             </variables>
             <meshes>
                 <mesh>cube_1x1x1_hex_1e3</mesh>
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/calculatesurfaceflux_ht_cube_1e4.prj b/Tests/Data/Parabolic/HT/SimpleSynthetics/calculatesurfaceflux_ht_cube_1e4.prj
index fdb7a2720035a162b0db8cea5ae7905c4deb13e5..f8dee6dd70d9158381450059403aae19c3851f0d 100644
--- a/Tests/Data/Parabolic/HT/SimpleSynthetics/calculatesurfaceflux_ht_cube_1e4.prj
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/calculatesurfaceflux_ht_cube_1e4.prj
@@ -171,6 +171,7 @@
                 <variable>T</variable>
                 <variable>p</variable>
                 <variable>darcy_velocity</variable>
+                <variable>specific_flux</variable>
             </variables>
             <meshes>
                 <mesh>cube_1x1x1_hex_1e4</mesh>
diff --git a/Tests/Data/Parabolic/LiquidFlow/AxiSymTheis/axisym_theis.prj b/Tests/Data/Parabolic/LiquidFlow/AxiSymTheis/axisym_theis.prj
index 5a12c1785b7a20e615cc928f679341f9da38424c..621bdc6bf407fe19c0edb2de2aa184aabe78039a 100644
--- a/Tests/Data/Parabolic/LiquidFlow/AxiSymTheis/axisym_theis.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/AxiSymTheis/axisym_theis.prj
@@ -120,6 +120,7 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v </variable>
+                <variable> OGS5_pressure </variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1.8e1_calculatesurfaceflux.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1.8e1_calculatesurfaceflux.prj
index 563da1555a10636e2cd89de6e1cf95b0cafad0f8..dc06dc32ea3f0a5ef8e28af74295bd354015ab1b 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1.8e1_calculatesurfaceflux.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1.8e1_calculatesurfaceflux.prj
@@ -62,6 +62,7 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v </variable>
+                <variable> specific_flux </variable>
             </variables>
             <meshes>
                 <mesh>square_1x1_tri_1.8e1</mesh>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1e1_calculatesurfaceflux.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1e1_calculatesurfaceflux.prj
index 0cb01ecdb015d956d662b649985515cd7303c698..3a2093f81827270e572d48710675ffd2fde82b64 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1e1_calculatesurfaceflux.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/2D/square_1e1_calculatesurfaceflux.prj
@@ -62,6 +62,7 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v </variable>
+                <variable> specific_flux </variable>
             </variables>
             <meshes>
                 <mesh>square_1x1_quad_1e1</mesh>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27.vtu
new file mode 120000
index 0000000000000000000000000000000000000000..bc5b276db6282dc50691497465aab3ce1f378a58
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27.vtu
@@ -0,0 +1 @@
+../cuboid_1x1x1_hex_27.vtu
\ No newline at end of file
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
new file mode 120000
index 0000000000000000000000000000000000000000..812ca537ec8f7f950d0c25c41adeec6c8fa48229
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
@@ -0,0 +1 @@
+../cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
\ No newline at end of file
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet_multiple_outputs_different_variables.xml b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet_multiple_outputs_different_variables.xml
new file mode 100644
index 0000000000000000000000000000000000000000..45cb82a5a56689fe65b0812ca1950ef37b0a6673
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet_multiple_outputs_different_variables.xml
@@ -0,0 +1,134 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProjectDiff base_file="cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj">
+    <add sel="/*/time_loop/output" pos="after">
+        <outputs>
+            <output>
+                <type>VTK</type>
+                <prefix>first_output_config_{:meshname}</prefix>
+                <suffix>_ts_{:timestep}_t_{:time}</suffix>
+                <timesteps>
+                    <pair>
+                        <repeat> 1 </repeat>
+                        <each_steps> 1 </each_steps>
+                    </pair>
+                </timesteps>
+                <variables>
+                    <variable> pressure </variable>
+                    <variable> v </variable>
+                    <variable> VolumetricFlowRate </variable>
+                </variables>
+                <meshes>
+                    <mesh>cuboid_1x1x1_hex_27</mesh>
+                    <mesh>cuboid_1x1x1_hex_27_top_boundary</mesh>
+                    <mesh>cuboid_1x1x1_hex_27_bottom_boundary</mesh>
+                </meshes>
+            </output>
+            <output>
+                <type>VTK</type>
+                <prefix>second_output_config_{:meshname}</prefix>
+                <suffix>_ts_{:timestep}_t_{:time}</suffix>
+                <timesteps>
+                    <pair>
+                        <repeat> 1 </repeat>
+                        <each_steps> 2 </each_steps>
+                    </pair>
+                </timesteps>
+                <variables>
+                    <variable> v </variable>
+                </variables>
+                <meshes>
+                    <mesh>cuboid_1x1x1_hex_27</mesh>
+                </meshes>
+            </output>
+            <output>
+                <type>VTK</type>
+                <prefix>third_output_config_{:meshname}</prefix>
+                <suffix>_ts_{:timestep}_t_{:time}</suffix>
+                <timesteps>
+                    <pair>
+                        <repeat> 1 </repeat>
+                        <each_steps> 2 </each_steps>
+                    </pair>
+                </timesteps>
+                <variables>
+                    <variable> VolumetricFlowRate </variable>
+                </variables>
+                <meshes>
+                    <mesh>cuboid_1x1x1_hex_27</mesh>
+                </meshes>
+            </output>
+            <output>
+                <type>VTK</type>
+                <prefix>fourth_output_config_{:meshname}</prefix>
+                <suffix>_ts_{:timestep}_t_{:time}</suffix>
+                <timesteps>
+                    <pair>
+                        <repeat> 1 </repeat>
+                        <each_steps> 2 </each_steps>
+                    </pair>
+                </timesteps>
+                <variables>
+                </variables>
+                <meshes>
+                    <mesh>cuboid_1x1x1_hex_27</mesh>
+                </meshes>
+            </output>
+        </outputs>
+    </add>
+    <remove sel="/*/time_loop/output"/>
+    <add sel="/*/linear_solvers" pos="after">
+        <test_definition>
+            <vtkdiff>
+                <regex>first_output_config_cuboid_1x1x1_hex_27.*.vtu</regex>
+                <field>pressure</field>
+                <absolute_tolerance>1e-10</absolute_tolerance>
+                <relative_tolerance>1e-15</relative_tolerance>
+            </vtkdiff>
+            <vtkdiff>
+                <regex>first_output_config_cuboid_1x1x1_hex_27.*.vtu</regex>
+                <field>v</field>
+                <absolute_tolerance>1e-10</absolute_tolerance>
+                <relative_tolerance>1e-15</relative_tolerance>
+            </vtkdiff>
+            <vtkdiff>
+                <regex>first_output_config_cuboid_1x1x1_hex_27.*.vtu</regex>
+                <field>VolumetricFlowRate</field>
+                <absolute_tolerance>1e-10</absolute_tolerance>
+                <relative_tolerance>1e-15</relative_tolerance>
+            </vtkdiff>
+
+            <vtkdiff>
+                <regex>second_output_config_cuboid_1x1x1_hex_27.*.vtu</regex>
+                <field>v</field>
+                <absolute_tolerance>1e-10</absolute_tolerance>
+                <relative_tolerance>1e-15</relative_tolerance>
+            </vtkdiff>
+
+            <vtkdiff>
+                <regex>third_output_config_cuboid_1x1x1_hex_27.*.vtu</regex>
+                <field>VolumetricFlowRate</field>
+                <absolute_tolerance>1e-10</absolute_tolerance>
+                <relative_tolerance>1e-15</relative_tolerance>
+            </vtkdiff>
+
+            <vtkdiff>
+                <regex>fourth_output_config_cuboid_1x1x1_hex_27.*.vtu</regex>
+                <field>pressure</field>
+                <absolute_tolerance>1e-10</absolute_tolerance>
+                <relative_tolerance>1e-15</relative_tolerance>
+            </vtkdiff>
+            <vtkdiff>
+                <regex>fourth_output_config_cuboid_1x1x1_hex_27.*.vtu</regex>
+                <field>v</field>
+                <absolute_tolerance>1e-10</absolute_tolerance>
+                <relative_tolerance>1e-15</relative_tolerance>
+            </vtkdiff>
+            <vtkdiff>
+                <regex>fourth_output_config_cuboid_1x1x1_hex_27.*.vtu</regex>
+                <field>VolumetricFlowRate</field>
+                <absolute_tolerance>1e-10</absolute_tolerance>
+                <relative_tolerance>1e-15</relative_tolerance>
+            </vtkdiff>
+        </test_definition>
+    </add>
+</OpenGeoSysProjectDiff>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_bottom_boundary.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_bottom_boundary.vtu
new file mode 120000
index 0000000000000000000000000000000000000000..8c8f99b15c87df80ded19cd78b46df48f45c7fd8
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_bottom_boundary.vtu
@@ -0,0 +1 @@
+../cuboid_1x1x1_hex_27_bottom_boundary.vtu
\ No newline at end of file
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_top_boundary.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_top_boundary.vtu
new file mode 120000
index 0000000000000000000000000000000000000000..35723b44db4facca7cd93798ea9fdc96234f0448
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/cuboid_1x1x1_hex_27_top_boundary.vtu
@@ -0,0 +1 @@
+../cuboid_1x1x1_hex_27_top_boundary.vtu
\ No newline at end of file
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_bottom_boundary_ts_0_t_0.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_bottom_boundary_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..076ae4e9f34479692f655e6e14673599c025c631
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_bottom_boundary_ts_0_t_0.000000.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="16"                   NumberOfCells="9"                   >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="0"                    RangeMax="0"                    offset="0"                   />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="5000000"              RangeMax="5000000"              offset="60"                  />
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="6.1270569481e-17"     RangeMax="4.841887057e-16"      offset="132"                 />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.4142135624"         offset="684"                 />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="820"                 />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="956"                 />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="1044"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAEwAAAAAAAAA=eF5jYGBguCAU7MgwQDQAwRsXYQ==AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAfQEAAAAAAAA=eF7z8e6KLP+1xOZR3fPpcepzbRbPzigskV+15/OkpT6iv+bv0THea7pbpcamIGf/i2ciZTZxGZ/ttr6fuGeN+h0ZiczVNnbbj6h2/t5q41qaJn5PvnnPp1Ljqhmbp9gUHrDYbjVryx6NMz/6nCd12ziA1S21+VDb9MzDbJNNW7ztisVvl+05BDY3bc9OhatCEc/X2dzqPeeqNanH5p5Ejt/zvKY9PSKbWvvki21u5jZu3vd0+p4Tt3qmFnut3FMCMd/mi9TECUuYu21cnA8kna/fsSdJ9ta66f+rbQQmaoaYt87fE8vq8+bG0ml7rDob/AQfzbGxOxfGdpFh+Z6TbmdLWiqm7ikDmzNlj5Ni38cHMxfaBNoINj66MdUmEyp+o8zNsev1Ihvl4zVhq77ttykCiy/Zc2euuYKjY8eelYe2JhbqrLRZckrwF6PWPJvZ8Tn/gxmX7Fm06YnGnTPVcHceyY1vYvuYaWP8JXt+4JQNNrPB/m6zAQDL7cWfAQAAAAAAAAAAgAAAAAAAAIABAAAAAAAARQAAAAAAAAA=eF5jYMAOQsHgqj128acY4hDwAYc4pnkwPnbxpxjiEPABhzimu1DNRxd/iiEOAR9wiCPkYSxUd6CLP8UQR+UjxAHSjzsZAQAAAAAAAAAAgAAAAAAAACABAAAAAAAARQAAAAAAAAA=eF5dzTkSgDAQxEBz2GBz+P+/JUBKZpOuUjBbyn87VlxCe8M1tB+4obsndqyhfWAL7Rf6x90bH+yh/cUR2if65wN1qAEPAQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAIAAAAAAAAAA=eF5jYYAADijNA6UFoLQIlJaA0jJQWgFKq0BpABToALU=AQAAAAAAAAAAgAAAAAAAAAkAAAAAAAAACwAAAAAAAAA=eF7j5IQCAAGeAFI=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_bottom_boundary_ts_1_t_43200.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_bottom_boundary_ts_1_t_43200.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..9530497388b580f437fc66eff62befd786f70a0c
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_bottom_boundary_ts_1_t_43200.000000.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="16"                   NumberOfCells="9"                   >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="-0.094723294723"      RangeMax="-0.023680823681"      offset="0"                   />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="1000000"              RangeMax="1000000"              offset="120"                 />
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="0.63938223938"        RangeMax="0.63938223938"        offset="192"                 />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.4142135624"         offset="664"                 />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="800"                 />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="936"                 />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="1024"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAANwAAAAAAAAA=eF7708yoe9J+xv4/YHrF/r9Q+h9U/DeUD5HfAZXfsf8/VPw3ir4dUH0wGqZvBtxcdPMBkhhBzw==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAEwAAAAAAAAA=eF5jZGBgaGjRc2QcIBoAttsXQQ==AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAPwEAAAAAAAA=eF4TCs7OrvHttLneZsy7av18myj3n3svlj7Zn/+w9M3s7Lo9Nc8nXRZXqrZJhoovuzf97Mx/CTbzGU8laT9KhYtnVuxecW/dBBuLVH7hwkuLbAqg4npM3dIBi7P2NBRpxntsnWQTBxWfu5C37tLtiXsMLl5vVrGesCcRKi7rvGPx9v6CPdwNx5waK+vg6nujeiIZtNv2rOTI2Cq7NckmHSp+VGrRxgOJjTYfpn3+sHeDm00CVNzvc8oxxsVJe3pyrl9MsHHdkwQV5zzbX1L9qNHmeGnarTO85Xtg7s/64zdFxDLTJv/xna9rMqbDzZe/H6T73LvKpv5H6MVK5/Y9KVDxH+bb395d7r4nyzaTTcK/E26vtGbG7o31DXsOxXrqvFFth7s/9Z3gwZDJATbXHv/VP/Buzh5YOAMAtMjTpQ==AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAARQAAAAAAAAA=eF5jYMAOQsHgqj128acY4hDwAYc4pnkwPnbxpxjiEPABhzimu1DNRxd/iiEOAR9wiCPkYSxUd6CLP8UQR+UjxAHSjzsZAQAAAAAAAAAAgAAAAAAAACABAAAAAAAARQAAAAAAAAA=eF5dzTkSgDAQxEBz2GBz+P+/JUBKZpOuUjBbyn87VlxCe8M1tB+4obsndqyhfWAL7Rf6x90bH+yh/cUR2if65wN1qAEPAQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAIAAAAAAAAAA=eF5jYYAADijNA6UFoLQIlJaA0jJQWgFKq0BpABToALU=AQAAAAAAAAAAgAAAAAAAAAkAAAAAAAAACwAAAAAAAAA=eF7j5IQCAAGeAFI=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_bottom_boundary_ts_2_t_86400.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_bottom_boundary_ts_2_t_86400.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..4dbdbd30102a99648e8e781cbb822f69b9cb07a0
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_bottom_boundary_ts_2_t_86400.000000.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="16"                   NumberOfCells="9"                   >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="-0.071042471042"      RangeMax="-0.017760617761"      offset="0"                   />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="1000000"              RangeMax="1000000"              offset="120"                 />
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="0.63938223938"        RangeMax="0.63938223938"        offset="192"                 />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.4142135624"         offset="668"                 />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="804"                 />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="940"                 />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="1028"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOQAAAAAAAAA=eF77o3Tw4TX9Sfv/gOlF+/9B6f9o4hB6E1R+E1R+EZo+mPim/YzKyOKT4OYyQMUh8pP2AwCQqU+FAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAEwAAAAAAAAA=eF5jZGBgaGjRc2QcIBoAttsXQQ==AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAARAEAAAAAAAA=eF7jUmioWX+hcs/js5+ebtbQ3RPh/nPvxdIn+3fpmf9L82raw3Bkr7lvxASbZKi4R7Z9w3fXFptnPVdTQxwabVKh4qUVu1fcWzfBJuG77Y2IRXNt8qHiykJbTu7iqrMJsBVLiORy2pMAFT/Iv7ru3qOOPX5Rt/a8uNhsAxM/tuvfaoWu7D0tL5Vm/FCu3RMPFXcS+j5rS3TznpUv91zvz6+wSYGK8zSEXk+wnmDDYNTK9rDqs3UcVNzhWs2pj4XONnbmYUcY3rXuSYKKdzcdPt0yO8rm6bY7fMzKoXsSoeIu7xqPTLpZZmOrom5s5t1hkw4V14v5nnrtg98ejscz6lb9nWiTBhVvsIl4YrczbE+BSE+glUMF3F/npBe09RQl7tkvK9uiH7IY7q8rSdtUJk8L2bOuueSdj+CCPVFQcQCQc9PFAQAAAAAAAAAAgAAAAAAAAIABAAAAAAAARQAAAAAAAAA=eF5jYMAOQsHgqj128acY4hDwAYc4pnkwPnbxpxjiEPABhzimu1DNRxd/iiEOAR9wiCPkYSxUd6CLP8UQR+UjxAHSjzsZAQAAAAAAAAAAgAAAAAAAACABAAAAAAAARQAAAAAAAAA=eF5dzTkSgDAQxEBz2GBz+P+/JUBKZpOuUjBbyn87VlxCe8M1tB+4obsndqyhfWAL7Rf6x90bH+yh/cUR2if65wN1qAEPAQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAIAAAAAAAAAA=eF5jYYAADijNA6UFoLQIlJaA0jJQWgFKq0BpABToALU=AQAAAAAAAAAAgAAAAAAAAAkAAAAAAAAACwAAAAAAAAA=eF7j5IQCAAGeAFI=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_top_boundary_ts_0_t_0.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_top_boundary_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..ac51aa7d5cb16ddd68c7623ff72cc902c4c24a38
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_top_boundary_ts_0_t_0.000000.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="16"                   NumberOfCells="9"                   >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="0"                    RangeMax="0"                    offset="0"                   />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="5000000"              RangeMax="5000000"              offset="60"                  />
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="4.585651416e-17"      RangeMax="6.3143132445e-16"     offset="132"                 />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="1"                    RangeMax="1.7320508076"         offset="680"                 />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="812"                 />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="948"                 />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="1036"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAADAAAAAAAAAA=eF5jYBhYAAAAgAABAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAEwAAAAAAAAA=eF5jYGBguCAU7MgwQDQAwRsXYQ==AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAeAEAAAAAAAA=eF4ruhxc9tRk9h73PEWZ5JL0PS93L1Ey02rYU1F5s0CrKWdP0z+RSs8jk2xmnxL8xajVZ7NOYsMerfVL91R+u6HmvHGRzYKc/S+eiUyz2b6/IuCX56I9Qeszen/83bCn6oDFdqtZR/Z8+Nm86n1Q+Z4g7xa3irMz93wCm1+w54Kkh05gnvSenPcP7U8s6bBxEdnU2idfbBOk7vVy3/tsm/mT/P//Cl1gMxlq74KIbZM/r5qypyHU2bA+vcTm/SWFW3OvL9zzTuJVcsz9eXu8KwPWT1LcsqcAbO8UGwPjd6dDNlXtEdWfz2MhvWaPwN2fMqnps/csa85Wv87danPgkDdL6Uq3PS7bj6h2/l665wTTi8ylPitt7mye6jSrpmZPLtSct0U7C7wjFtjs0VCp58ias2c+2L9te/j+P3fnbpphs6fu+fQ49V6bLqg72WpKQ7ZMnWVzIdZHOVBhEtyc36UzTJXuTtmjH7m0ZJHDBhsGKAAAqka+TQ==AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAQQAAAAAAAAA=eF5jYMAGPtiHgsFVe+ziTzHEUWl0cQYGVPMQ5mMXf4ohjq4Pu/lPcZiPLv4UQxxdHt187OpwuQ/TXFw0AEYRTgk=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAARQAAAAAAAAA=eF5djDkOwCAAw6AXbenx/98yEC/OYsmyUspcDfdwk1/CQx1+DZs6fuiv8JSnv9Xh+e3q+KF/w0ee/lOH5/dXNwBz+AEPAQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAIAAAAAAAAAA=eF5jYYAADijNA6UFoLQIlJaA0jJQWgFKq0BpABToALU=AQAAAAAAAAAAgAAAAAAAAAkAAAAAAAAACwAAAAAAAAA=eF7j5IQCAAGeAFI=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_top_boundary_ts_1_t_43200.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_top_boundary_ts_1_t_43200.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..6d0bea559e0e1dcb5d08029a3a2f161a051ffc98
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_top_boundary_ts_1_t_43200.000000.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="16"                   NumberOfCells="9"                   >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="0.029601029601"       RangeMax="0.1184041184"         offset="0"                   />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="10000000"             RangeMax="10000000"             offset="124"                 />
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="0.63938223938"        RangeMax="0.63938223938"        offset="196"                 />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="1"                    RangeMax="1.7320508076"         offset="676"                 />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="808"                 />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="944"                 />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="1032"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOwAAAAAAAAA=eF778cSxYrf/PPvPYHqd/Vco/QMq/g3K/wKm99l/h9K/oOIw+V9Q8a9QGiYOM+cnCh+mb549AHz+SK4=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAEwAAAAAAAAA=eF5jYGBguCCU7MgwQDQA/xsYYQ==AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAARgEAAAAAAAA=eF7bNeU8+5yozXtmhjob1qdP2VPk/nPvxdIn+3l62E/l/Vu2Z1tp8AGleSttEqDi+66yprKLbd/z7feTNvam9XBxY1mXacl2G/esu+li3Oy5Y08oVFwp/WnY1dzOPWbKrKyzNBftSYGK27hpfmIST9rjcPDlhR+eoXuioeIWVcFsj0Iy9tSJ/Tk/IXwq3PwvbZMMEn2m7WHYs0TJTCtgTyRUvGvPknlx8Rv2bJmSmL5l7t49aVBxh9nyGSFX0m2OLynZzpu3Z08cVFzra/Tdv90LbBJezZnI2pq8Jx4qvtJ/I7u30UabMpcGY1m3FXtCoOLtbEHL6z5ussk4bn/wl/TBPWFQcZ2+hVuaulbZXLWYoyw6YY0NzPwaP7Flm6U22TTt5CvoNp6xJwIqPkV+YtxiZnub6M7k/68/HrJJhooDAPJm0P4=AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAQQAAAAAAAAA=eF5jYMAGPtiHgsFVe+ziTzHEUWl0cQYGVPMQ5mMXf4ohjq4Pu/lPcZiPLv4UQxxdHt187OpwuQ/TXFw0AEYRTgk=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAARQAAAAAAAAA=eF5djDkOwCAAw6AXbenx/98yEC/OYsmyUspcDfdwk1/CQx1+DZs6fuiv8JSnv9Xh+e3q+KF/w0ee/lOH5/dXNwBz+AEPAQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAIAAAAAAAAAA=eF5jYYAADijNA6UFoLQIlJaA0jJQWgFKq0BpABToALU=AQAAAAAAAAAAgAAAAAAAAAkAAAAAAAAACwAAAAAAAAA=eF7j5IQCAAGeAFI=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_top_boundary_ts_2_t_86400.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_top_boundary_ts_2_t_86400.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..c4b381b6d8400795963b8aa78c2686b2ef221f56
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_top_boundary_ts_2_t_86400.000000.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="16"                   NumberOfCells="9"                   >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="0.017760617761"       RangeMax="0.071042471042"       offset="0"                   />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="10000000"             RangeMax="10000000"             offset="120"                 />
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="0.63938223938"        RangeMax="0.63938223938"        offset="192"                 />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="1"                    RangeMax="1.7320508076"         offset="664"                 />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="796"                 />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="932"                 />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="1020"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOAAAAAAAAAA=eF77oXTw4TX9SfbfwPQi+18o9CT7H1D+dzC9yf4nlP6FVf0muDqYeX+g5sDUfYfSX6DiANNWSic=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAEwAAAAAAAAA=eF5jYGBguCCU7MgwQDQA/xsYYQ==AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAPwEAAAAAAAA=eF77YOD2IPj/hj3HJvbrre1vsCly/7n3YumT/e6K79mOPFu2h3P6uhXSF1bZJEDFp5rO/W4aum1Pas7/kt8v58PFL3DFFDMwbtlzsWyG9rol+/eEQMUXNDlyv/Js3VOpUznN+8SCPWlQcQ6+f1fc5OJs3ri/X6Qwq3pPLFTcY0rX1K8iyXtSwObU28RBxd0630/h/TJ1z4JcRZnkknQbmPoDJ5uM2+9v2PPk47Qz7uqH9iRBxRsKHhRqcjjZ6Ln/u3ineO8emDmpy9P0XykusBEQn8pw9X6YTTxUvGXSy7PGiRttTvyvl9Bc1Qx3fydb0PK6j5tsVipcFYp4vm9PBFS8mlFL7uOqtTa3f0Zt+FCyGO6eyVoL1mQIb7Fp2clX0G08A65eWnFi3GJme5vF/312fPE/aJMMFQcAdIjWSw==AQAAAAAAAAAAgAAAAAAAAIABAAAAAAAAQQAAAAAAAAA=eF5jYMAGPtiHgsFVe+ziTzHEUWl0cQYGVPMQ5mMXf4ohjq4Pu/lPcZiPLv4UQxxdHt187OpwuQ/TXFw0AEYRTgk=AQAAAAAAAAAAgAAAAAAAACABAAAAAAAARQAAAAAAAAA=eF5djDkOwCAAw6AXbenx/98yEC/OYsmyUspcDfdwk1/CQx1+DZs6fuiv8JSnv9Xh+e3q+KF/w0ee/lOH5/dXNwBz+AEPAQAAAAAAAAAAgAAAAAAAAEgAAAAAAAAAIAAAAAAAAAA=eF5jYYAADijNA6UFoLQIlJaA0jJQWgFKq0BpABToALU=AQAAAAAAAAAAgAAAAAAAAAkAAAAAAAAACwAAAAAAAAA=eF7j5IQCAAGeAFI=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..f5397a52e14c596c1f22b5f74751dc972eecab5b
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="64"                   NumberOfCells="27"                  >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="0"                    RangeMax="0"                    offset="0"                   />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="5000000"              RangeMax="5000000"              offset="64"                  />
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="1.1277624687e-17"     RangeMax="6.3143132445e-16"     offset="136"                 />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.7320508076"         offset="2020"                />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="2292"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="2712"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="2852"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAAACAAAAAAAADgAAAAAAAAA=eF5jYBgFIxkAAAIAAAE=AQAAAAAAAAAAgAAAAAAAAAACAAAAAAAAFQAAAAAAAAA=eF5jYGBguCAU7MgwSo9IGgAos12BAQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAYwUAAAAAAAA=eF49lHk41Hkcx0e1aqsNpajIkbDDRlHJfgqRcmw5puKJrCS3HTLOWfd9TJrDjHPGmEmjMpOjVJ8tD3rotqsW88iWLbpdKemw+/zM0/fPz/N93p/v5/X+fN+uLvnecTMiGEoe4R4xqoSa8uCoGJ06nGSKXVfO8NHU4o8tVw3oQA2//nxYPRaOBE/ubB49heeMBrQ0Q87CzksdG/I+NcNu2nGNQZ0MnKBZJPIa2RDVanXJuqwJje9MM+yZBWBL3BPD2O/pw3u3NkC2344zNW9OYxuhexwv6z5c7jVSD/Kie7vJzEIY1AzfNxKZjoXqDVkMnRPQH5HWeO0ZF7vkhZwTzhKMmdOHd2tOFYvmF4CDfevR+ykteFRbXs+dTQLVUz9StmXx0fc719d94hK0zkvdpzZUATvvHVT+k1SLNx3vxmTGczCW0GHjLj3G+OPSanAHtbShPg6EKOp9sY52+a+EsL6TfrDu/XWIJuoiHKjcpmtnl4uStmb/KFMJiG6pzSiRq6DcL3zWU0mEwoanxgN3kr69syPCL115PAQs3oXx3dkyKCfmzoYZSYSbXwsPLewPL9Mix0GqgrNBQCd9oZ0MPe+fpXczeDgYyzOpF/Gx61e/zvsSHjo7FPl9aZWCSOFXn1FLD0U7BD4Lzutss0lBb4JbBjp8L2hY7psEKk7qXXoX6fj0L115ZW81VGb1J7McK+BrXHeotzkD+4l6GspfB7FG1PMwcxX4iE1/g9BQWtXTgQzMWG1QHdzDhTGxkrGRPAE0Bthh+2NlIGhqXmgwJcAcglsW7prrC5LzVS5mf9NQcO1F97TTFpggfGdiv+H2546ZSZBRbW3xVoWF7YTvDJxWabe1LwhE7qslWj0kAVQRPJPRChY1/NRTCzKlW0dNhmgQqvBl2ROKdWldDej6q5nHN3Kxc44PmFqv7NXIrQCZL7PS7GUUlq5zATu5EFZQCx5cSQyA4Bt3XGwWe2NLe6hme90hpPocG9Tj+mNqsAllSOSPkXP68GnjjgRSTQ1WqEzd6H7IgxvjJvm3BsqQVsKRWwu5wBpTDrzZmoeGQYnWBlNUJPn4GFebxOMxWdlgPzsBlta7sXt9aiE+lzbcczsCYh7Pt2oOPI9nFPtu7kHxZ2+OxsSZ5CfBUyycOFSUlMbPxivjNbb/XmWiZkTBeMGVClTdZOEn0ijGQtLione/iEG2VZ9qP82EYqK/BJmfLYcsf64G43uxlvlf2LCKmC8BdXPKuzjaAmDNLzCa51ACQesKug5I6XBmsXhqeAEH/c1rh1bsYYF+zvu4mQfFUL/Gkyldn4J9qz9YuNt4AYk4Ydtt1SLVZ1cVg6G5e1p2SBtqE/OmgjrZmaLfFw6mG60P65OrcTOxvzm46LRWQNuGEBD4SE9O7BRCK+HLSVSNWy5lTNbD15cffhBMXVDok0jRPZ6xzyzLcU+knlZATBC+uCrS30pOxfiEfio5PRzTv6onOHUwoZzYBwbUa8qQLBVjwvs+Q/sLQhAQPEvg0vV4txknIXpIg4umv8gwkfCxA8c+ZtSNesShh0umY/zdUpwg9KnYvXqvqXvkWgwffWLTJcoFB0XOeBg5v7g2GgZ85v7ZmQP/81P0FXhdZE3WsTH1gP2mlKAYGJ37R/hW82WAzz9V6JLgJmXqNSFVsT/mFm9vUxoScaUZf6nV2nOo+uijVmBQOZ7OCDPqXZIFrW0uC2gSR3SYy0Xsmvc8ROwqgYFGzq4yOh0jFDpvoi9TXbwEgMYGKYtCK5A/lxu4bHZkz5J0HiCR20WQr3inMp1GaeKUQbev63p3XeY3nU803hb9R2w08xbHCG1l3/j/B0Ip9VY=AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAqwAAAAAAAAA=eF6Fk9sNgDAIRd3MKRzdEfjvCCY2xNyX8IOctlcK5Ti8Xa/dp+clfNsKXPU69ryEb1uBa16oz7yEb1uBf+v9hXkwL+EYz3VC/ek+KY9c739fwlvfc+3bv14JR/25//4c14ffywr8s/R/5PqeWtdz7R/GuS6ok+6p+fG+yaP+PB/bp76pXsees6X+pHNprvK8duy53o/Pef2pT7lOvM76fl/KT3WTfwC9AmKgAQAAAAAAAAAAgAAAAAAAAMAGAAAAAAAAGQEAAAAAAAA=eF6F0rdWAlAURFEMiAlEERMg0YACZkEB+f+/ovCckjWv2dXMfcUUCv9vC4u4i1U8xTqeo7lt3EN7zJ3hBdpjbgdLaI+5Gl6iPf7Tewe4j/7Te9d4hea8d4j2mPPeDdpjzn8foT3m/HcD7fGf3ivjMfpP77Wwiea8V0F7zHnvFu0x579P0B5z/ruN9qR9dLCLdzjAtA9zPbxHe9I+zPXxAe1J+/Cf3hviI6Z9mPPeE9qT9mHOfz+jPWkf/tN7Yxxh2oc5703QnrQPc/77Be3p4KZ9vOIbfuEnpn2Ye8dvtCftw9wHTtGetA//6b0fnGHahznv/aI9aR/m/Pcc7Un78J/eW+IC0z7Mee8P7Un7MOe/V2jPGpr9GpU=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAASAAAAAAAAAA=eF4txRECgEAAALAwDMMwDMPwMAzDMAzDMAwPwzAMw7DnBG2yNPllzl24dOXajYNbd+49ePTk2YtXb47effj05duPX3+ShQvRAQAAAAAAAAAAgAAAAAAAABsAAAAAAAAACwAAAAAAAAA=eF7j4cEJABHTAUU=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_ts_1_t_43200.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_ts_1_t_43200.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..c6098689218d3f3c6a312d22f4aaed134c2f4559
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_ts_1_t_43200.000000.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="64"                   NumberOfCells="27"                  >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="-0.1184041184"        RangeMax="0.1184041184"         offset="0"                   />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="1000000"              RangeMax="10000000"             offset="276"                 />
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="0.63938223938"        RangeMax="0.63938223938"        offset="440"                 />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.7320508076"         offset="2048"                />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="2320"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="2740"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="2880"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAAACAAAAAAAArQAAAAAAAAA=eF5ljj0OgkAQRpH7eIaJJ/A69noIKmor6k1I2Np7SE+iiPybMG82rlYv39uZb3Y67/Y3yfy08epnuOBHsr473p1f8WO059gz2l4Wen/7oXBfkss3s5D1nhP6wzx3ZI68++vpmOvhwN3mfjiVx9xPGwv/hAO+Ic8bK/+GSzRX+BFvc2mtPqm1Z2XOmOK1L5eHemmh+RdZ/1VJB7kb3gd8C81bTx9l28vlA/KuDVg=AQAAAAAAAAAAgAAAAAAAAAACAAAAAAAAWgAAAAAAAAA=eF5jZGBgaGjRc2QcIJoBTPs5/geCeiD9D0rD+H/R+DD1jFCaAY1mhtJMaOYyoeljRtL3YXMUnAapf4+FZkKTZ0Tjw/TDxJnR+OjmwGiQvgtCyQNGAwBrRJyKAQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAlQQAAAAAAAA=eF491Hk03FcUB/BEj7UZSVPB4aiGRkyMClokXkJCqKqlDIJax49Biihpgk4sLQYRI0NKrDNii4axZCZcibGEo0RsQTVyKnFiCdnV3v7zm/fv59xz3/19332//a7h4QkObDSRakyprS9FXrbr7cOxz+9H/hP78kY4CxJe5I6qaMUjBumVM9cHC3b8Uenu/kC9WQJ76MW26pnbOciM2Pvp+REeiiL9S6lMdWd+GCRGU/2+aclFvqQXl1NYI9McODo8kfKFeQ4EkK5hJeILr0bBx4kPTiddYuH6K15Znrv0UqFGjtmi0RKIQkjvUeMJOgKS0Ov8d6/bG2yQP+mO74Ie7OYHQta5iWF/dAYCSZcfvBoTP5uEemOD/xqg/AyS+cO2HLlKx0JR5LO/P9Qxr+P+mk9d9F/Yx6HLa27Dl6zSIIj0NVPh8pMqWwg7ESqj6sTG56pTmW2Cy4nQ6WNHe3koDc9PrHwipl9zRo+fbRt0rBSBJGeL1rQrcvc5qHif0CeooBxJ+k/nTenSRzLASKHLUHQoCQjSLecZOxHUC0jFZqpu2qMIf9fiMtPU8mg24rwNJQ5O1oFknj67YAoDaiC+uJXoFPJwPZ2f3edg4QnswUXhOi0dSdxA+6RZz1ge5M6yPQwF9bj+8fl134U4HhzXHburs1WC+/u10GnHPG8CI6U/M8o5GucwSaS8muXWwc2vjQfU5yrAj/Q2RWY5U4MPZzcWNty3OdhVlJeP0QxYiBg+dZ/j9Ae+l6/CvlUUGQlgrUHw7HwhD+d5kGi/Rbhx4fHS5E+jKXfAh/QZ83m79JRfYVKdsufPR5W4/txhUNdUbQZ2eGHX3FoMnnM+Rv59l1MzOtU029RM70GSPXz13jFDr5qPKIlyMj1J3TiHjFKpfn15IUrYf8H7w3Y1drVdS2M59fnogHLpzMOIDpxPjf7JKiVaD6xGxoz0GRXgcy0nFM3nxVy4iJQu2o4m4v2f/D5MJag7D62urPfspFZhd/eIjpXZuQVxBuOHzRVFuP++uA33lbI0lP1UR/V0/TXsZwhR0BC1GKxM67Z20tvxuZ1ptTWfld8D4RvT1ZicJpw/S1RWEv4EoEQhWfEX5zL4gXQT6gl5K04r0qFF7m3sFIMkT5dAGOkeqgYtwsBzE6qQpN6WN9Bnsi2Gh2rVDntfMfG9hOiOhB4PLkP/Kth/3r7Zhudp5Q7JFnk1QYGbleHlEC5Ek74nS7Y/YqcS7sS6dmiV1OD6e+PShKyyEFY3n6fKJtdjN9awzmecFMDtKWvjFDsRuJGuFTLnPv4jG0y0paULqTz8fpEN9a2USiBYihcerdm5gTfpZnGuMrN0JrCUt4ZyPPJw//epuUcDvsuHXVChZXLEGTwl+wAVJb5+DdDMDQhpLm6HYMn93tBk0sdCUG9FjJASATi3Ix+8n2xnliH/xSKO9G8MnH+Nk0DW3kiALlgnGmvYVAOd9DQZlyrWm0bE7LUQb6iLwZ10WnZ5c3JGLRo3K9I+kFOH9zzBUbmySa0RJd9VjMo0/h3Oks7V5PjyP7JA3mzGf0tvOvF/+399QFqzAQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAqwAAAAAAAAA=eF6Fk9sNgDAIRd3MKRzdEfjvCCY2xNyX8IOctlcK5Ti8Xa/dp+clfNsKXPU69ryEb1uBa16oz7yEb1uBf+v9hXkwL+EYz3VC/ek+KY9c739fwlvfc+3bv14JR/25//4c14ffywr8s/R/5PqeWtdz7R/GuS6ok+6p+fG+yaP+PB/bp76pXsees6X+pHNprvK8duy53o/Pef2pT7lOvM76fl/KT3WTfwC9AmKgAQAAAAAAAAAAgAAAAAAAAMAGAAAAAAAAGQEAAAAAAAA=eF6F0rdWAlAURFEMiAlEERMg0YACZkEB+f+/ovCckjWv2dXMfcUUCv9vC4u4i1U8xTqeo7lt3EN7zJ3hBdpjbgdLaI+5Gl6iPf7Tewe4j/7Te9d4hea8d4j2mPPeDdpjzn8foT3m/HcD7fGf3ivjMfpP77Wwiea8V0F7zHnvFu0x579P0B5z/ruN9qR9dLCLdzjAtA9zPbxHe9I+zPXxAe1J+/Cf3hviI6Z9mPPeE9qT9mHOfz+jPWkf/tN7Yxxh2oc5703QnrQPc/77Be3p4KZ9vOIbfuEnpn2Ye8dvtCftw9wHTtGetA//6b0fnGHahznv/aI9aR/m/Pcc7Un78J/eW+IC0z7Mee8P7Un7MOe/V2jPGpr9GpU=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAASAAAAAAAAAA=eF4txRECgEAAALAwDMMwDMPwMAzDMAzDMAwPwzAMw7DnBG2yNPllzl24dOXajYNbd+49ePTk2YtXb47effj05duPX3+ShQvRAQAAAAAAAAAAgAAAAAAAABsAAAAAAAAACwAAAAAAAAA=eF7j4cEJABHTAUU=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..c66e07cee9811919ee33ef7c5b85f5753cf2ed44
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/first_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="64"                   NumberOfCells="27"                  >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="-0.071042471042"      RangeMax="0.071042471042"       offset="0"                   />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="1000000"              RangeMax="10000000"             offset="264"                 />
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="0.63938223938"        RangeMax="0.63938223938"        offset="424"                 />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.7320508076"         offset="2032"                />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="2304"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="2724"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="2864"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAAACAAAAAAAAowAAAAAAAAA=eF77o3Tw4TX9Sfv/gOlF+/9B6f9o4hB6E1R+E1R+EZo+mPim/YzKyOKT4OYyQMUh8pP2M4BBgg2EbtjDgAIKoPwONBqmrgGqr8AGlT8BjZ+wB00duj1Q8Slo5qODBTD70OyHAZj75kDFZ6DJT0BzBwPDD0j42H+DhI/9LxR6kv0PKP87JFztf0LpX1jVb4Krg5n3B2oOTN13KP0FKg4A0XqyBw==AQAAAAAAAAAAgAAAAAAAAAACAAAAAAAAVQAAAAAAAAA=eF5jZGBgaGjRc2QcIPo/ENS3+DkygPl+cD6M/ofGh6lDp2HyjFA+E5o4M5o4C5L+D5ujwPpANIwP0vceiUaXR1eHbg4zDnOZsIhfEEoeMBoAHZSUlA==AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAkwQAAAAAAAA=eF49lHkw3GcYx2kwjtJkK6JVV6QRJCgJTfNKZLXuEA2KbJTIYoq6rwg2rBDMWCw2jnZlWTfbEAkerMQ4WgwmREQddSQr6o6NdtCZzm/ffz/zzvN+j+d9ZTVodxqGY2F+cHOx8dQZcLPabR+JXOhs1Tfd97VNArHudtMrbgzkQ3DrgEs04Xd0tJQ5RnU2v4uoBI+MaaucrmcgL6HZhNvDEhRMcC1SU3+rbAK6aqbk5S5LBi+Cd31SkzD9Vxo4eEzC25FkJOI9rfs1GukBQBccZ33QiocfCU4mCQubridDlQBeZgXHoFsE/5jm8tLrAgOJGaVIzd3euuBJcPPxO79vhFqgi6au3WKrKXCT4BlJz/+gF3mgxeYphUNaLuBN8G9X73bnvIpCZie0jU3s0pAfwfUpQur4ugNIz7MSqveykS/Bacht4WKLK4QoZjp9Yx6DfQ2psO9lhnlDp6oq3cCZg329uNl8IjffGeqTI1btj7DBg+CMoqB/ZBxykePpoUW6Ow/5E3zvvERv4mIxDAmsvArKmCDyS9NhiiufsUXqBsfi8mYLMc8ku0VvKXEQc45ROTSSjkR+ixto9hBdCgF7fs/Kg2jYr8LTnUe1kyFokJwlGU2nIREPdCqvCnyXBOz3mSa9j7Pw+bRM9+e2htFoPiDhRgqJA6Kcl8VrRyfIbCjO0Iu/V8rGOagNkxyFEr9Akow3W2mHhXtMC5VsSJ7gAGOKF/Xn5xw8X3XUVJZBrQVD1tKbnuuxuK+v9K5tmjKboHJSa7q8KBaJ5pi7LBeZrOTATpi98vjRVrhBcEWLwzKpRwOBnzPdYLr9GCgEP/XWt7vipxrI8tuW5jdxsM5kC135kk0u7AwsDcgVAO5rK+MHzhH1QCTl/6WvMi8D71sN2aZaPoqP+BVXIiVDufj8h79JbWUdT5BkFfXnkyQ+zsfXwDOS6dwCta2uo3FnH+B7T1YLnT+L4AJ3ea/KjlqOfSl7TuRRPIoQQ0Vo7HTJGs8xfFXRk11Ug9SmjHq52nwQvcfx16VyHR/lIyOt+HDhKhv3OH42IWO/ngVLOl0kP7tirHOXQ8muM+QDA8qOm+hexXp43r2pgzkANoIwsZgdPt5Pun0ddbGMiwT5qaOUyS6c50Tsg9AWZh009tdar3zNwflPUoJvC9qegIT/zLvzS3FIxA1jGmtsTJqR+pzFpJNnHxLNWTe0nL12wIOe7Cz9uiwaCiO4leaaVPcSF2QK6itVhqux/rxzJcJzLs1ADTyI+FfwK+bDspRwMfEmGIli6dWXdYIzwdlJl+WWbVIg9nRsvl0fG0TvV1ph/4WlmidasVp7qFEYh/VbM9Pz3iv6wK3/5yQiUf6W99eY8tt5wA7S/MInwg/74vcnGafO8GBhI3/ASvsZ3ltayGyojjQZ6Vvtj0yFt+MeqRW+BsuabHT4WJ7Y2Iwr7p2eIxg09v4N9R0kKutUJ2P996W+r0jYeISqNMZIbm868P8cJ66rtlFdh17vevDWIzhYT64uu9b/0yZEb1EIyTBm4fMqmtmenEOXEOfA/um2Yxf+z/8DoPxfvg==AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAqwAAAAAAAAA=eF6Fk9sNgDAIRd3MKRzdEfjvCCY2xNyX8IOctlcK5Ti8Xa/dp+clfNsKXPU69ryEb1uBa16oz7yEb1uBf+v9hXkwL+EYz3VC/ek+KY9c739fwlvfc+3bv14JR/25//4c14ffywr8s/R/5PqeWtdz7R/GuS6ok+6p+fG+yaP+PB/bp76pXsees6X+pHNprvK8duy53o/Pef2pT7lOvM76fl/KT3WTfwC9AmKgAQAAAAAAAAAAgAAAAAAAAMAGAAAAAAAAGQEAAAAAAAA=eF6F0rdWAlAURFEMiAlEERMg0YACZkEB+f+/ovCckjWv2dXMfcUUCv9vC4u4i1U8xTqeo7lt3EN7zJ3hBdpjbgdLaI+5Gl6iPf7Tewe4j/7Te9d4hea8d4j2mPPeDdpjzn8foT3m/HcD7fGf3ivjMfpP77Wwiea8V0F7zHnvFu0x579P0B5z/ruN9qR9dLCLdzjAtA9zPbxHe9I+zPXxAe1J+/Cf3hviI6Z9mPPeE9qT9mHOfz+jPWkf/tN7Yxxh2oc5703QnrQPc/77Be3p4KZ9vOIbfuEnpn2Ye8dvtCftw9wHTtGetA//6b0fnGHahznv/aI9aR/m/Pcc7Un78J/eW+IC0z7Mee8P7Un7MOe/V2jPGpr9GpU=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAASAAAAAAAAAA=eF4txRECgEAAALAwDMMwDMPwMAzDMAzDMAwPwzAMw7DnBG2yNPllzl24dOXajYNbd+49ePTk2YtXb47effj05duPX3+ShQvRAQAAAAAAAAAAgAAAAAAAABsAAAAAAAAACwAAAAAAAAA=eF7j4cEJABHTAUU=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/fourth_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/fourth_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..ab596d5b086fe32dc8e7ee0ab7755a3ddb38c157
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/fourth_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu
@@ -0,0 +1,28 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <FieldData>
+      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="25" format="appended" RangeMin="45"                   RangeMax="121"                  offset="0"                   />
+    </FieldData>
+    <Piece NumberOfPoints="64"                   NumberOfCells="27"                  >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="0"                    RangeMax="0"                    offset="88"                  />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="5000000"              RangeMax="5000000"              offset="152"                 />
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="1.1277624687e-17"     RangeMax="6.3143132445e-16"     offset="224"                 />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.7320508076"         offset="2108"                />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="2380"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="2800"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="2940"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAABkAAAAAAAAAIQAAAAAAAAA=eF4z0zPRM9E1MjXRTU9LtTBPSjYy1UvJLCqpBABPOwcXAQAAAAAAAAAAgAAAAAAAAAACAAAAAAAADgAAAAAAAAA=eF5jYBgFIxkAAAIAAAE=AQAAAAAAAAAAgAAAAAAAAAACAAAAAAAAFQAAAAAAAAA=eF5jYGBguCAU7MgwSo9IGgAos12BAQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAYwUAAAAAAAA=eF49lHk41Hkcx0e1aqsNpajIkbDDRlHJfgqRcmw5puKJrCS3HTLOWfd9TJrDjHPGmEmjMpOjVJ8tD3rotqsW88iWLbpdKemw+/zM0/fPz/N93p/v5/X+fN+uLvnecTMiGEoe4R4xqoSa8uCoGJ06nGSKXVfO8NHU4o8tVw3oQA2//nxYPRaOBE/ubB49heeMBrQ0Q87CzksdG/I+NcNu2nGNQZ0MnKBZJPIa2RDVanXJuqwJje9MM+yZBWBL3BPD2O/pw3u3NkC2344zNW9OYxuhexwv6z5c7jVSD/Kie7vJzEIY1AzfNxKZjoXqDVkMnRPQH5HWeO0ZF7vkhZwTzhKMmdOHd2tOFYvmF4CDfevR+ykteFRbXs+dTQLVUz9StmXx0fc719d94hK0zkvdpzZUATvvHVT+k1SLNx3vxmTGczCW0GHjLj3G+OPSanAHtbShPg6EKOp9sY52+a+EsL6TfrDu/XWIJuoiHKjcpmtnl4uStmb/KFMJiG6pzSiRq6DcL3zWU0mEwoanxgN3kr69syPCL115PAQs3oXx3dkyKCfmzoYZSYSbXwsPLewPL9Mix0GqgrNBQCd9oZ0MPe+fpXczeDgYyzOpF/Gx61e/zvsSHjo7FPl9aZWCSOFXn1FLD0U7BD4Lzutss0lBb4JbBjp8L2hY7psEKk7qXXoX6fj0L115ZW81VGb1J7McK+BrXHeotzkD+4l6GspfB7FG1PMwcxX4iE1/g9BQWtXTgQzMWG1QHdzDhTGxkrGRPAE0Bthh+2NlIGhqXmgwJcAcglsW7prrC5LzVS5mf9NQcO1F97TTFpggfGdiv+H2546ZSZBRbW3xVoWF7YTvDJxWabe1LwhE7qslWj0kAVQRPJPRChY1/NRTCzKlW0dNhmgQqvBl2ROKdWldDej6q5nHN3Kxc44PmFqv7NXIrQCZL7PS7GUUlq5zATu5EFZQCx5cSQyA4Bt3XGwWe2NLe6hme90hpPocG9Tj+mNqsAllSOSPkXP68GnjjgRSTQ1WqEzd6H7IgxvjJvm3BsqQVsKRWwu5wBpTDrzZmoeGQYnWBlNUJPn4GFebxOMxWdlgPzsBlta7sXt9aiE+lzbcczsCYh7Pt2oOPI9nFPtu7kHxZ2+OxsSZ5CfBUyycOFSUlMbPxivjNbb/XmWiZkTBeMGVClTdZOEn0ijGQtLione/iEG2VZ9qP82EYqK/BJmfLYcsf64G43uxlvlf2LCKmC8BdXPKuzjaAmDNLzCa51ACQesKug5I6XBmsXhqeAEH/c1rh1bsYYF+zvu4mQfFUL/Gkyldn4J9qz9YuNt4AYk4Ydtt1SLVZ1cVg6G5e1p2SBtqE/OmgjrZmaLfFw6mG60P65OrcTOxvzm46LRWQNuGEBD4SE9O7BRCK+HLSVSNWy5lTNbD15cffhBMXVDok0jRPZ6xzyzLcU+knlZATBC+uCrS30pOxfiEfio5PRzTv6onOHUwoZzYBwbUa8qQLBVjwvs+Q/sLQhAQPEvg0vV4txknIXpIg4umv8gwkfCxA8c+ZtSNesShh0umY/zdUpwg9KnYvXqvqXvkWgwffWLTJcoFB0XOeBg5v7g2GgZ85v7ZmQP/81P0FXhdZE3WsTH1gP2mlKAYGJ37R/hW82WAzz9V6JLgJmXqNSFVsT/mFm9vUxoScaUZf6nV2nOo+uijVmBQOZ7OCDPqXZIFrW0uC2gSR3SYy0Xsmvc8ROwqgYFGzq4yOh0jFDpvoi9TXbwEgMYGKYtCK5A/lxu4bHZkz5J0HiCR20WQr3inMp1GaeKUQbev63p3XeY3nU803hb9R2w08xbHCG1l3/j/B0Ip9VY=AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAqwAAAAAAAAA=eF6Fk9sNgDAIRd3MKRzdEfjvCCY2xNyX8IOctlcK5Ti8Xa/dp+clfNsKXPU69ryEb1uBa16oz7yEb1uBf+v9hXkwL+EYz3VC/ek+KY9c739fwlvfc+3bv14JR/25//4c14ffywr8s/R/5PqeWtdz7R/GuS6ok+6p+fG+yaP+PB/bp76pXsees6X+pHNprvK8duy53o/Pef2pT7lOvM76fl/KT3WTfwC9AmKgAQAAAAAAAAAAgAAAAAAAAMAGAAAAAAAAGQEAAAAAAAA=eF6F0rdWAlAURFEMiAlEERMg0YACZkEB+f+/ovCckjWv2dXMfcUUCv9vC4u4i1U8xTqeo7lt3EN7zJ3hBdpjbgdLaI+5Gl6iPf7Tewe4j/7Te9d4hea8d4j2mPPeDdpjzn8foT3m/HcD7fGf3ivjMfpP77Wwiea8V0F7zHnvFu0x579P0B5z/ruN9qR9dLCLdzjAtA9zPbxHe9I+zPXxAe1J+/Cf3hviI6Z9mPPeE9qT9mHOfz+jPWkf/tN7Yxxh2oc5703QnrQPc/77Be3p4KZ9vOIbfuEnpn2Ye8dvtCftw9wHTtGetA//6b0fnGHahznv/aI9aR/m/Pcc7Un78J/eW+IC0z7Mee8P7Un7MOe/V2jPGpr9GpU=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAASAAAAAAAAAA=eF4txRECgEAAALAwDMMwDMPwMAzDMAzDMAwPwzAMw7DnBG2yNPllzl24dOXajYNbd+49ePTk2YtXb47effj05duPX3+ShQvRAQAAAAAAAAAAgAAAAAAAABsAAAAAAAAACwAAAAAAAAA=eF7j4cEJABHTAUU=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/fourth_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/fourth_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..9c53612fc2143b59ade27db2d8302371606c5083
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/fourth_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu
@@ -0,0 +1,28 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <FieldData>
+      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="25" format="appended" RangeMin="45"                   RangeMax="121"                  offset="0"                   />
+    </FieldData>
+    <Piece NumberOfPoints="64"                   NumberOfCells="27"                  >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="-0.071042471042"      RangeMax="0.071042471042"       offset="88"                  />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="1000000"              RangeMax="10000000"             offset="352"                 />
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="0.63938223938"        RangeMax="0.63938223938"        offset="512"                 />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.7320508076"         offset="2120"                />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="2392"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="2812"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="2952"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAABkAAAAAAAAAIQAAAAAAAAA=eF4z0zPRM9E1MjXRTU9LtTBPSjYy1UvJLCqpBABPOwcXAQAAAAAAAAAAgAAAAAAAAAACAAAAAAAAowAAAAAAAAA=eF77o3Tw4TX9Sfv/gOlF+/9B6f9o4hB6E1R+E1R+EZo+mPim/YzKyOKT4OYyQMUh8pP2M4BBgg2EbtjDgAIKoPwONBqmrgGqr8AGlT8BjZ+wB00duj1Q8Slo5qODBTD70OyHAZj75kDFZ6DJT0BzBwPDD0j42H+DhI/9LxR6kv0PKP87JFztf0LpX1jVb4Krg5n3B2oOTN13KP0FKg4A0XqyBw==AQAAAAAAAAAAgAAAAAAAAAACAAAAAAAAVQAAAAAAAAA=eF5jZGBgaGjRc2QcIPo/ENS3+DkygPl+cD6M/ofGh6lDp2HyjFA+E5o4M5o4C5L+D5ujwPpANIwP0vceiUaXR1eHbg4zDnOZsIhfEEoeMBoAHZSUlA==AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAkwQAAAAAAAA=eF49lHkw3GcYx2kwjtJkK6JVV6QRJCgJTfNKZLXuEA2KbJTIYoq6rwg2rBDMWCw2jnZlWTfbEAkerMQ4WgwmREQddSQr6o6NdtCZzm/ffz/zzvN+j+d9ZTVodxqGY2F+cHOx8dQZcLPabR+JXOhs1Tfd97VNArHudtMrbgzkQ3DrgEs04Xd0tJQ5RnU2v4uoBI+MaaucrmcgL6HZhNvDEhRMcC1SU3+rbAK6aqbk5S5LBi+Cd31SkzD9Vxo4eEzC25FkJOI9rfs1GukBQBccZ33QiocfCU4mCQubridDlQBeZgXHoFsE/5jm8tLrAgOJGaVIzd3euuBJcPPxO79vhFqgi6au3WKrKXCT4BlJz/+gF3mgxeYphUNaLuBN8G9X73bnvIpCZie0jU3s0pAfwfUpQur4ugNIz7MSqveykS/Bacht4WKLK4QoZjp9Yx6DfQ2psO9lhnlDp6oq3cCZg329uNl8IjffGeqTI1btj7DBg+CMoqB/ZBxykePpoUW6Ow/5E3zvvERv4mIxDAmsvArKmCDyS9NhiiufsUXqBsfi8mYLMc8ku0VvKXEQc45ROTSSjkR+ixto9hBdCgF7fs/Kg2jYr8LTnUe1kyFokJwlGU2nIREPdCqvCnyXBOz3mSa9j7Pw+bRM9+e2htFoPiDhRgqJA6Kcl8VrRyfIbCjO0Iu/V8rGOagNkxyFEr9Akow3W2mHhXtMC5VsSJ7gAGOKF/Xn5xw8X3XUVJZBrQVD1tKbnuuxuK+v9K5tmjKboHJSa7q8KBaJ5pi7LBeZrOTATpi98vjRVrhBcEWLwzKpRwOBnzPdYLr9GCgEP/XWt7vipxrI8tuW5jdxsM5kC135kk0u7AwsDcgVAO5rK+MHzhH1QCTl/6WvMi8D71sN2aZaPoqP+BVXIiVDufj8h79JbWUdT5BkFfXnkyQ+zsfXwDOS6dwCta2uo3FnH+B7T1YLnT+L4AJ3ea/KjlqOfSl7TuRRPIoQQ0Vo7HTJGs8xfFXRk11Ug9SmjHq52nwQvcfx16VyHR/lIyOt+HDhKhv3OH42IWO/ngVLOl0kP7tirHOXQ8muM+QDA8qOm+hexXp43r2pgzkANoIwsZgdPt5Pun0ddbGMiwT5qaOUyS6c50Tsg9AWZh009tdar3zNwflPUoJvC9qegIT/zLvzS3FIxA1jGmtsTJqR+pzFpJNnHxLNWTe0nL12wIOe7Cz9uiwaCiO4leaaVPcSF2QK6itVhqux/rxzJcJzLs1ADTyI+FfwK+bDspRwMfEmGIli6dWXdYIzwdlJl+WWbVIg9nRsvl0fG0TvV1ph/4WlmidasVp7qFEYh/VbM9Pz3iv6wK3/5yQiUf6W99eY8tt5wA7S/MInwg/74vcnGafO8GBhI3/ASvsZ3ltayGyojjQZ6Vvtj0yFt+MeqRW+BsuabHT4WJ7Y2Iwr7p2eIxg09v4N9R0kKutUJ2P996W+r0jYeISqNMZIbm868P8cJ66rtlFdh17vevDWIzhYT64uu9b/0yZEb1EIyTBm4fMqmtmenEOXEOfA/um2Yxf+z/8DoPxfvg==AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAqwAAAAAAAAA=eF6Fk9sNgDAIRd3MKRzdEfjvCCY2xNyX8IOctlcK5Ti8Xa/dp+clfNsKXPU69ryEb1uBa16oz7yEb1uBf+v9hXkwL+EYz3VC/ek+KY9c739fwlvfc+3bv14JR/25//4c14ffywr8s/R/5PqeWtdz7R/GuS6ok+6p+fG+yaP+PB/bp76pXsees6X+pHNprvK8duy53o/Pef2pT7lOvM76fl/KT3WTfwC9AmKgAQAAAAAAAAAAgAAAAAAAAMAGAAAAAAAAGQEAAAAAAAA=eF6F0rdWAlAURFEMiAlEERMg0YACZkEB+f+/ovCckjWv2dXMfcUUCv9vC4u4i1U8xTqeo7lt3EN7zJ3hBdpjbgdLaI+5Gl6iPf7Tewe4j/7Te9d4hea8d4j2mPPeDdpjzn8foT3m/HcD7fGf3ivjMfpP77Wwiea8V0F7zHnvFu0x579P0B5z/ruN9qR9dLCLdzjAtA9zPbxHe9I+zPXxAe1J+/Cf3hviI6Z9mPPeE9qT9mHOfz+jPWkf/tN7Yxxh2oc5703QnrQPc/77Be3p4KZ9vOIbfuEnpn2Ye8dvtCftw9wHTtGetA//6b0fnGHahznv/aI9aR/m/Pcc7Un78J/eW+IC0z7Mee8P7Un7MOe/V2jPGpr9GpU=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAASAAAAAAAAAA=eF4txRECgEAAALAwDMMwDMPwMAzDMAzDMAwPwzAMw7DnBG2yNPllzl24dOXajYNbd+49ePTk2YtXb47effj05duPX3+ShQvRAQAAAAAAAAAAgAAAAAAAABsAAAAAAAAACwAAAAAAAAA=eF7j4cEJABHTAUU=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/second_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/second_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..278c1e6ee70f91fbf8b0517f60ab9e560fd69ee5
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/second_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu
@@ -0,0 +1,23 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="64"                   NumberOfCells="27"                  >
+      <PointData>
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="1.1277624687e-17"     RangeMax="6.3143132445e-16"     offset="0"                   />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.7320508076"         offset="1884"                />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="2156"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="2576"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="2716"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAYwUAAAAAAAA=eF49lHk41Hkcx0e1aqsNpajIkbDDRlHJfgqRcmw5puKJrCS3HTLOWfd9TJrDjHPGmEmjMpOjVJ8tD3rotqsW88iWLbpdKemw+/zM0/fPz/N93p/v5/X+fN+uLvnecTMiGEoe4R4xqoSa8uCoGJ06nGSKXVfO8NHU4o8tVw3oQA2//nxYPRaOBE/ubB49heeMBrQ0Q87CzksdG/I+NcNu2nGNQZ0MnKBZJPIa2RDVanXJuqwJje9MM+yZBWBL3BPD2O/pw3u3NkC2344zNW9OYxuhexwv6z5c7jVSD/Kie7vJzEIY1AzfNxKZjoXqDVkMnRPQH5HWeO0ZF7vkhZwTzhKMmdOHd2tOFYvmF4CDfevR+ykteFRbXs+dTQLVUz9StmXx0fc719d94hK0zkvdpzZUATvvHVT+k1SLNx3vxmTGczCW0GHjLj3G+OPSanAHtbShPg6EKOp9sY52+a+EsL6TfrDu/XWIJuoiHKjcpmtnl4uStmb/KFMJiG6pzSiRq6DcL3zWU0mEwoanxgN3kr69syPCL115PAQs3oXx3dkyKCfmzoYZSYSbXwsPLewPL9Mix0GqgrNBQCd9oZ0MPe+fpXczeDgYyzOpF/Gx61e/zvsSHjo7FPl9aZWCSOFXn1FLD0U7BD4Lzutss0lBb4JbBjp8L2hY7psEKk7qXXoX6fj0L115ZW81VGb1J7McK+BrXHeotzkD+4l6GspfB7FG1PMwcxX4iE1/g9BQWtXTgQzMWG1QHdzDhTGxkrGRPAE0Bthh+2NlIGhqXmgwJcAcglsW7prrC5LzVS5mf9NQcO1F97TTFpggfGdiv+H2546ZSZBRbW3xVoWF7YTvDJxWabe1LwhE7qslWj0kAVQRPJPRChY1/NRTCzKlW0dNhmgQqvBl2ROKdWldDej6q5nHN3Kxc44PmFqv7NXIrQCZL7PS7GUUlq5zATu5EFZQCx5cSQyA4Bt3XGwWe2NLe6hme90hpPocG9Tj+mNqsAllSOSPkXP68GnjjgRSTQ1WqEzd6H7IgxvjJvm3BsqQVsKRWwu5wBpTDrzZmoeGQYnWBlNUJPn4GFebxOMxWdlgPzsBlta7sXt9aiE+lzbcczsCYh7Pt2oOPI9nFPtu7kHxZ2+OxsSZ5CfBUyycOFSUlMbPxivjNbb/XmWiZkTBeMGVClTdZOEn0ijGQtLione/iEG2VZ9qP82EYqK/BJmfLYcsf64G43uxlvlf2LCKmC8BdXPKuzjaAmDNLzCa51ACQesKug5I6XBmsXhqeAEH/c1rh1bsYYF+zvu4mQfFUL/Gkyldn4J9qz9YuNt4AYk4Ydtt1SLVZ1cVg6G5e1p2SBtqE/OmgjrZmaLfFw6mG60P65OrcTOxvzm46LRWQNuGEBD4SE9O7BRCK+HLSVSNWy5lTNbD15cffhBMXVDok0jRPZ6xzyzLcU+knlZATBC+uCrS30pOxfiEfio5PRzTv6onOHUwoZzYBwbUa8qQLBVjwvs+Q/sLQhAQPEvg0vV4txknIXpIg4umv8gwkfCxA8c+ZtSNesShh0umY/zdUpwg9KnYvXqvqXvkWgwffWLTJcoFB0XOeBg5v7g2GgZ85v7ZmQP/81P0FXhdZE3WsTH1gP2mlKAYGJ37R/hW82WAzz9V6JLgJmXqNSFVsT/mFm9vUxoScaUZf6nV2nOo+uijVmBQOZ7OCDPqXZIFrW0uC2gSR3SYy0Xsmvc8ROwqgYFGzq4yOh0jFDpvoi9TXbwEgMYGKYtCK5A/lxu4bHZkz5J0HiCR20WQr3inMp1GaeKUQbev63p3XeY3nU803hb9R2w08xbHCG1l3/j/B0Ip9VY=AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAqwAAAAAAAAA=eF6Fk9sNgDAIRd3MKRzdEfjvCCY2xNyX8IOctlcK5Ti8Xa/dp+clfNsKXPU69ryEb1uBa16oz7yEb1uBf+v9hXkwL+EYz3VC/ek+KY9c739fwlvfc+3bv14JR/25//4c14ffywr8s/R/5PqeWtdz7R/GuS6ok+6p+fG+yaP+PB/bp76pXsees6X+pHNprvK8duy53o/Pef2pT7lOvM76fl/KT3WTfwC9AmKgAQAAAAAAAAAAgAAAAAAAAMAGAAAAAAAAGQEAAAAAAAA=eF6F0rdWAlAURFEMiAlEERMg0YACZkEB+f+/ovCckjWv2dXMfcUUCv9vC4u4i1U8xTqeo7lt3EN7zJ3hBdpjbgdLaI+5Gl6iPf7Tewe4j/7Te9d4hea8d4j2mPPeDdpjzn8foT3m/HcD7fGf3ivjMfpP77Wwiea8V0F7zHnvFu0x579P0B5z/ruN9qR9dLCLdzjAtA9zPbxHe9I+zPXxAe1J+/Cf3hviI6Z9mPPeE9qT9mHOfz+jPWkf/tN7Yxxh2oc5703QnrQPc/77Be3p4KZ9vOIbfuEnpn2Ye8dvtCftw9wHTtGetA//6b0fnGHahznv/aI9aR/m/Pcc7Un78J/eW+IC0z7Mee8P7Un7MOe/V2jPGpr9GpU=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAASAAAAAAAAAA=eF4txRECgEAAALAwDMMwDMPwMAzDMAzDMAwPwzAMw7DnBG2yNPllzl24dOXajYNbd+49ePTk2YtXb47effj05duPX3+ShQvRAQAAAAAAAAAAgAAAAAAAABsAAAAAAAAACwAAAAAAAAA=eF7j4cEJABHTAUU=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/second_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/second_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..734d88e76abd5f83ce1f5f68fdf815bac0637ea0
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/second_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu
@@ -0,0 +1,23 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="64"                   NumberOfCells="27"                  >
+      <PointData>
+        <DataArray type="Float64" Name="v" NumberOfComponents="3" format="appended" RangeMin="0.63938223938"        RangeMax="0.63938223938"        offset="0"                   />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.7320508076"         offset="1608"                />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="1880"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="2300"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="2440"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAkwQAAAAAAAA=eF49lHkw3GcYx2kwjtJkK6JVV6QRJCgJTfNKZLXuEA2KbJTIYoq6rwg2rBDMWCw2jnZlWTfbEAkerMQ4WgwmREQddSQr6o6NdtCZzm/ffz/zzvN+j+d9ZTVodxqGY2F+cHOx8dQZcLPabR+JXOhs1Tfd97VNArHudtMrbgzkQ3DrgEs04Xd0tJQ5RnU2v4uoBI+MaaucrmcgL6HZhNvDEhRMcC1SU3+rbAK6aqbk5S5LBi+Cd31SkzD9Vxo4eEzC25FkJOI9rfs1GukBQBccZ33QiocfCU4mCQubridDlQBeZgXHoFsE/5jm8tLrAgOJGaVIzd3euuBJcPPxO79vhFqgi6au3WKrKXCT4BlJz/+gF3mgxeYphUNaLuBN8G9X73bnvIpCZie0jU3s0pAfwfUpQur4ugNIz7MSqveykS/Bacht4WKLK4QoZjp9Yx6DfQ2psO9lhnlDp6oq3cCZg329uNl8IjffGeqTI1btj7DBg+CMoqB/ZBxykePpoUW6Ow/5E3zvvERv4mIxDAmsvArKmCDyS9NhiiufsUXqBsfi8mYLMc8ku0VvKXEQc45ROTSSjkR+ixto9hBdCgF7fs/Kg2jYr8LTnUe1kyFokJwlGU2nIREPdCqvCnyXBOz3mSa9j7Pw+bRM9+e2htFoPiDhRgqJA6Kcl8VrRyfIbCjO0Iu/V8rGOagNkxyFEr9Akow3W2mHhXtMC5VsSJ7gAGOKF/Xn5xw8X3XUVJZBrQVD1tKbnuuxuK+v9K5tmjKboHJSa7q8KBaJ5pi7LBeZrOTATpi98vjRVrhBcEWLwzKpRwOBnzPdYLr9GCgEP/XWt7vipxrI8tuW5jdxsM5kC135kk0u7AwsDcgVAO5rK+MHzhH1QCTl/6WvMi8D71sN2aZaPoqP+BVXIiVDufj8h79JbWUdT5BkFfXnkyQ+zsfXwDOS6dwCta2uo3FnH+B7T1YLnT+L4AJ3ea/KjlqOfSl7TuRRPIoQQ0Vo7HTJGs8xfFXRk11Ug9SmjHq52nwQvcfx16VyHR/lIyOt+HDhKhv3OH42IWO/ngVLOl0kP7tirHOXQ8muM+QDA8qOm+hexXp43r2pgzkANoIwsZgdPt5Pun0ddbGMiwT5qaOUyS6c50Tsg9AWZh009tdar3zNwflPUoJvC9qegIT/zLvzS3FIxA1jGmtsTJqR+pzFpJNnHxLNWTe0nL12wIOe7Cz9uiwaCiO4leaaVPcSF2QK6itVhqux/rxzJcJzLs1ADTyI+FfwK+bDspRwMfEmGIli6dWXdYIzwdlJl+WWbVIg9nRsvl0fG0TvV1ph/4WlmidasVp7qFEYh/VbM9Pz3iv6wK3/5yQiUf6W99eY8tt5wA7S/MInwg/74vcnGafO8GBhI3/ASvsZ3ltayGyojjQZ6Vvtj0yFt+MeqRW+BsuabHT4WJ7Y2Iwr7p2eIxg09v4N9R0kKutUJ2P996W+r0jYeISqNMZIbm868P8cJ66rtlFdh17vevDWIzhYT64uu9b/0yZEb1EIyTBm4fMqmtmenEOXEOfA/um2Yxf+z/8DoPxfvg==AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAqwAAAAAAAAA=eF6Fk9sNgDAIRd3MKRzdEfjvCCY2xNyX8IOctlcK5Ti8Xa/dp+clfNsKXPU69ryEb1uBa16oz7yEb1uBf+v9hXkwL+EYz3VC/ek+KY9c739fwlvfc+3bv14JR/25//4c14ffywr8s/R/5PqeWtdz7R/GuS6ok+6p+fG+yaP+PB/bp76pXsees6X+pHNprvK8duy53o/Pef2pT7lOvM76fl/KT3WTfwC9AmKgAQAAAAAAAAAAgAAAAAAAAMAGAAAAAAAAGQEAAAAAAAA=eF6F0rdWAlAURFEMiAlEERMg0YACZkEB+f+/ovCckjWv2dXMfcUUCv9vC4u4i1U8xTqeo7lt3EN7zJ3hBdpjbgdLaI+5Gl6iPf7Tewe4j/7Te9d4hea8d4j2mPPeDdpjzn8foT3m/HcD7fGf3ivjMfpP77Wwiea8V0F7zHnvFu0x579P0B5z/ruN9qR9dLCLdzjAtA9zPbxHe9I+zPXxAe1J+/Cf3hviI6Z9mPPeE9qT9mHOfz+jPWkf/tN7Yxxh2oc5703QnrQPc/77Be3p4KZ9vOIbfuEnpn2Ye8dvtCftw9wHTtGetA//6b0fnGHahznv/aI9aR/m/Pcc7Un78J/eW+IC0z7Mee8P7Un7MOe/V2jPGpr9GpU=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAASAAAAAAAAAA=eF4txRECgEAAALAwDMMwDMPwMAzDMAzDMAwPwzAMw7DnBG2yNPllzl24dOXajYNbd+49ePTk2YtXb47effj05duPX3+ShQvRAQAAAAAAAAAAgAAAAAAAABsAAAAAAAAACwAAAAAAAAA=eF7j4cEJABHTAUU=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/third_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/third_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..6560d5b266b3137649445db9daf35c194e604e2e
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/third_output_config_cuboid_1x1x1_hex_27_ts_0_t_0.000000.vtu
@@ -0,0 +1,23 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="64"                   NumberOfCells="27"                  >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="0"                    RangeMax="0"                    offset="0"                   />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.7320508076"         offset="64"                  />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="336"                 />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="756"                 />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="896"                 />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAAACAAAAAAAADgAAAAAAAAA=eF5jYBgFIxkAAAIAAAE=AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAqwAAAAAAAAA=eF6Fk9sNgDAIRd3MKRzdEfjvCCY2xNyX8IOctlcK5Ti8Xa/dp+clfNsKXPU69ryEb1uBa16oz7yEb1uBf+v9hXkwL+EYz3VC/ek+KY9c739fwlvfc+3bv14JR/25//4c14ffywr8s/R/5PqeWtdz7R/GuS6ok+6p+fG+yaP+PB/bp76pXsees6X+pHNprvK8duy53o/Pef2pT7lOvM76fl/KT3WTfwC9AmKgAQAAAAAAAAAAgAAAAAAAAMAGAAAAAAAAGQEAAAAAAAA=eF6F0rdWAlAURFEMiAlEERMg0YACZkEB+f+/ovCckjWv2dXMfcUUCv9vC4u4i1U8xTqeo7lt3EN7zJ3hBdpjbgdLaI+5Gl6iPf7Tewe4j/7Te9d4hea8d4j2mPPeDdpjzn8foT3m/HcD7fGf3ivjMfpP77Wwiea8V0F7zHnvFu0x579P0B5z/ruN9qR9dLCLdzjAtA9zPbxHe9I+zPXxAe1J+/Cf3hviI6Z9mPPeE9qT9mHOfz+jPWkf/tN7Yxxh2oc5703QnrQPc/77Be3p4KZ9vOIbfuEnpn2Ye8dvtCftw9wHTtGetA//6b0fnGHahznv/aI9aR/m/Pcc7Un78J/eW+IC0z7Mee8P7Un7MOe/V2jPGpr9GpU=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAASAAAAAAAAAA=eF4txRECgEAAALAwDMMwDMPwMAzDMAzDMAwPwzAMw7DnBG2yNPllzl24dOXajYNbd+49ePTk2YtXb47effj05duPX3+ShQvRAQAAAAAAAAAAgAAAAAAAABsAAAAAAAAACwAAAAAAAAA=eF7j4cEJABHTAUU=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/third_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/third_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..e57d177e27b9de0ceac45fc1a2ff6aaf417d5f94
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/MultipleOutputsDifferentVariablesSections/third_output_config_cuboid_1x1x1_hex_27_ts_2_t_86400.000000.vtu
@@ -0,0 +1,23 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="64"                   NumberOfCells="27"                  >
+      <PointData>
+        <DataArray type="Float64" Name="VolumetricFlowRate" format="appended" RangeMin="-0.071042471042"      RangeMax="0.071042471042"       offset="0"                   />
+      </PointData>
+      <CellData>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.7320508076"         offset="264"                 />
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="536"                 />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="956"                 />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="1096"                />
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+  <AppendedData encoding="base64">
+   _AQAAAAAAAAAAgAAAAAAAAAACAAAAAAAAowAAAAAAAAA=eF77o3Tw4TX9Sfv/gOlF+/9B6f9o4hB6E1R+E1R+EZo+mPim/YzKyOKT4OYyQMUh8pP2M4BBgg2EbtjDgAIKoPwONBqmrgGqr8AGlT8BjZ+wB00duj1Q8Slo5qODBTD70OyHAZj75kDFZ6DJT0BzBwPDD0j42H+DhI/9LxR6kv0PKP87JFztf0LpX1jVb4Krg5n3B2oOTN13KP0FKg4A0XqyBw==AQAAAAAAAAAAgAAAAAAAAAAGAAAAAAAAqwAAAAAAAAA=eF6Fk9sNgDAIRd3MKRzdEfjvCCY2xNyX8IOctlcK5Ti8Xa/dp+clfNsKXPU69ryEb1uBa16oz7yEb1uBf+v9hXkwL+EYz3VC/ek+KY9c739fwlvfc+3bv14JR/25//4c14ffywr8s/R/5PqeWtdz7R/GuS6ok+6p+fG+yaP+PB/bp76pXsees6X+pHNprvK8duy53o/Pef2pT7lOvM76fl/KT3WTfwC9AmKgAQAAAAAAAAAAgAAAAAAAAMAGAAAAAAAAGQEAAAAAAAA=eF6F0rdWAlAURFEMiAlEERMg0YACZkEB+f+/ovCckjWv2dXMfcUUCv9vC4u4i1U8xTqeo7lt3EN7zJ3hBdpjbgdLaI+5Gl6iPf7Tewe4j/7Te9d4hea8d4j2mPPeDdpjzn8foT3m/HcD7fGf3ivjMfpP77Wwiea8V0F7zHnvFu0x579P0B5z/ruN9qR9dLCLdzjAtA9zPbxHe9I+zPXxAe1J+/Cf3hviI6Z9mPPeE9qT9mHOfz+jPWkf/tN7Yxxh2oc5703QnrQPc/77Be3p4KZ9vOIbfuEnpn2Ye8dvtCftw9wHTtGetA//6b0fnGHahznv/aI9aR/m/Pcc7Un78J/eW+IC0z7Mee8P7Un7MOe/V2jPGpr9GpU=AQAAAAAAAAAAgAAAAAAAANgAAAAAAAAASAAAAAAAAAA=eF4txRECgEAAALAwDMMwDMPwMAzDMAzDMAwPwzAMw7DnBG2yNPllzl24dOXajYNbd+49ePTk2YtXb47effj05duPX3+ShQvRAQAAAAAAAAAAgAAAAAAAABsAAAAAAAAACwAAAAAAAAA=eF7j4cEJABHTAUU=
+  </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/Parallel/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/Parallel/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
index 829af844a0e7cf24e15672fedfae490996b5985e..3df5dd23281a99adb119153414f8c67dd61e406e 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/Parallel/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Hex/Parallel/cuboid_1x1x1_hex_27_Dirichlet_Dirichlet.prj
@@ -62,6 +62,7 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v </variable>
+                <variable> specific_flux </variable>
             </variables>
             <meshes>
                 <mesh>cuboid_1x1x1_hex_27</mesh>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Pyramid/cuboid_1x1x1_pyramid_6000_calculatesurfaceflux.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Pyramid/cuboid_1x1x1_pyramid_6000_calculatesurfaceflux.prj
index b8024d2b4a2f312e97acb2dbe2c5abf35c6d3f17..700a326f97e23c809684608733ee2e90dd40f66e 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Pyramid/cuboid_1x1x1_pyramid_6000_calculatesurfaceflux.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/3D/Pyramid/cuboid_1x1x1_pyramid_6000_calculatesurfaceflux.prj
@@ -62,6 +62,7 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v </variable>
+                <variable> specific_flux </variable>
             </variables>
             <meshes>
                 <mesh>cuboid_1x1x1_pyramid_6000</mesh>
diff --git a/Tests/Data/Parabolic/LiquidFlow/Flux/cube_1e3_calculatesurfaceflux.prj b/Tests/Data/Parabolic/LiquidFlow/Flux/cube_1e3_calculatesurfaceflux.prj
index 44dce9e3156bd05a2504cec8a268785225de60e2..20c7eec8f32b08080198caceb233543b64dedacc 100644
--- a/Tests/Data/Parabolic/LiquidFlow/Flux/cube_1e3_calculatesurfaceflux.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/Flux/cube_1e3_calculatesurfaceflux.prj
@@ -62,6 +62,8 @@
             <variables>
                 <variable> pressure </variable>
                 <variable> v </variable>
+                <variable> specific_flux </variable>
+                <variable> VolumetricFlowRate </variable>
             </variables>
             <meshes>
                 <mesh>cube_1x1x1_hex_1e3</mesh>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/FunctionParameterTest_XDMF.prj b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/FunctionParameterTest_XDMF.prj
index 8adbfbf6b659eefbe88242b81b6143bcdd78e6fc..7753f7d2465a4dc136e10d026001cf907c427e72 100644
--- a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/FunctionParameterTest_XDMF.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/FunctionParameterTest_XDMF.prj
@@ -145,6 +145,10 @@ boundary meshes will be written as well -->
             <variables>
                 <variable> pressure </variable>
                 <variable> v </variable>
+                <variable> bulk_node_ids </variable>
+                <variable> bulk_element_ids </variable>
+                <variable> MaterialIDs </variable>
+                <variable> VolumetricFlowRate </variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF_compression_off/FunctionParameterTest_XDMF.prj b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF_compression_off/FunctionParameterTest_XDMF.prj
index a8f186dbc50bd5b6e31b1f063f4b71a02006bd92..773a7fe487038ef4dfa9df18b2624d59581df956 100644
--- a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF_compression_off/FunctionParameterTest_XDMF.prj
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF_compression_off/FunctionParameterTest_XDMF.prj
@@ -146,6 +146,10 @@ boundary meshes will be written as well -->
             <variables>
                 <variable> pressure </variable>
                 <variable> v </variable>
+                <variable> bulk_node_ids </variable>
+                <variable> bulk_element_ids </variable>
+                <variable> MaterialIDs </variable>
+                <variable> VolumetricFlowRate </variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet/line_60_heat.prj b/Tests/Data/Parabolic/T/1D_dirichlet/line_60_heat.prj
index ea4547158b0c8e12a34ace95957ef205e165a5b9..006833f04df0b00b9f4249db63e0c008b4dde9c9 100644
--- a/Tests/Data/Parabolic/T/1D_dirichlet/line_60_heat.prj
+++ b/Tests/Data/Parabolic/T/1D_dirichlet/line_60_heat.prj
@@ -78,6 +78,8 @@
             <variables>
                 <variable> temperature </variable>
                 <variable> heat_flux </variable>
+                <variable> HeatFlowRate </variable>
+                <variable> MaterialIDs </variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Parabolic/T/1D_neumann/newton.prj b/Tests/Data/Parabolic/T/1D_neumann/newton.prj
index 21a10fd2fe3b486df46641ca64b3b3eeaca16c4f..1b4a9fbe1ca11f796274d7e33fd384969f13a504 100644
--- a/Tests/Data/Parabolic/T/1D_neumann/newton.prj
+++ b/Tests/Data/Parabolic/T/1D_neumann/newton.prj
@@ -86,6 +86,7 @@
             <variables>
                 <variable> temperature </variable>
                 <variable> heat_flux </variable>
+                <variable> HeatFlowRate </variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/Parabolic/T/1D_neumann/picard.prj b/Tests/Data/Parabolic/T/1D_neumann/picard.prj
index a6eecf0f316c60ee28163abdd62527956daf6bda..b80f034c53027fe1f1ced7436ad1daa0dc397a5c 100644
--- a/Tests/Data/Parabolic/T/1D_neumann/picard.prj
+++ b/Tests/Data/Parabolic/T/1D_neumann/picard.prj
@@ -86,6 +86,7 @@
             <variables>
                 <variable> temperature </variable>
                 <variable> heat_flux </variable>
+                <variable> HeatFlowRate </variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/PhaseField/PForthotropy_jupyter_notebook/shear.prj b/Tests/Data/PhaseField/PForthotropy_jupyter_notebook/shear.prj
index 5d00757b932048978e3e362057295c8e5d6f74fa..8511788f0e8d0232a80f3fe09a53c05c905f405d 100644
--- a/Tests/Data/PhaseField/PForthotropy_jupyter_notebook/shear.prj
+++ b/Tests/Data/PhaseField/PForthotropy_jupyter_notebook/shear.prj
@@ -124,6 +124,9 @@
                 <variable>sigma_tensile</variable>
                 <variable>epsilon</variable>
                 <variable>epsilon_tensile</variable>
+                <variable>NodalForces</variable>
+                <variable>bulk_node_ids</variable>
+                <variable>bulk_element_ids</variable>
             </variables>
             <type>VTK</type>
             <prefix>shear_deg_0</prefix>
diff --git a/Tests/Data/PhaseField/beam_jupyter_notebook/beam.prj b/Tests/Data/PhaseField/beam_jupyter_notebook/beam.prj
index 484c0e7e06adf25d69634bd46c06ca3686091de1..5b87db37e0783baab8ced7d5b127d4698f3e19e4 100644
--- a/Tests/Data/PhaseField/beam_jupyter_notebook/beam.prj
+++ b/Tests/Data/PhaseField/beam_jupyter_notebook/beam.prj
@@ -108,6 +108,7 @@
                 <variable>phasefield</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <type>VTK</type>
             <prefix>AT1_iso_compressive</prefix>
diff --git a/Tests/Data/PhaseField/tpb_jupyter_notebook/TPB.prj b/Tests/Data/PhaseField/tpb_jupyter_notebook/TPB.prj
index 07fe9d9f9f87cbdcce6321c3ebffb9e181ff2bd9..798c3b8f951dec069ca50b34e2efbfa9d37e3db9 100644
--- a/Tests/Data/PhaseField/tpb_jupyter_notebook/TPB.prj
+++ b/Tests/Data/PhaseField/tpb_jupyter_notebook/TPB.prj
@@ -111,6 +111,7 @@
                 <variable>phasefield</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
+                <variable>NodalForces</variable>
             </variables>
             <type>VTK</type>
             <prefix>TPB_</prefix>
diff --git a/Tests/Data/RichardsMechanics/A2.prj b/Tests/Data/RichardsMechanics/A2.prj
index 93c4b23e5a7f9122806fca746f6a5746e047afaf..4d43d1148b4b18d0b3f80d9dd17fff7e12a67398 100644
--- a/Tests/Data/RichardsMechanics/A2.prj
+++ b/Tests/Data/RichardsMechanics/A2.prj
@@ -245,6 +245,8 @@
                 <variable>epsilon</variable>
                 <variable>saturation</variable>
                 <variable>velocity</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos.prj b/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos.prj
index 157877133958766b239c9303fadc4acf68b02d65..da44a66d250cd4868cedcc3003376306b0979cfa 100644
--- a/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos.prj
+++ b/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos.prj
@@ -175,6 +175,9 @@
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
                 <variable>saturation</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>pressure_interpolated</variable>
             </variables>
             <fixed_output_times>
                 0.06
diff --git a/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos_QN.prj b/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos_QN.prj
index cb29f5c74a0fd82980c94076f56730b8372b4345..61d1486686973e4b652324ed158276d31a7fbdc6 100644
--- a/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos_QN.prj
+++ b/Tests/Data/RichardsMechanics/LiakopoulosHM/liakopoulos_QN.prj
@@ -195,6 +195,9 @@
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
                 <variable>saturation</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>pressure_interpolated</variable>
             </variables>
             <fixed_output_times>
                 0.06
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_quasinewton.prj b/Tests/Data/RichardsMechanics/RichardsFlow_2d_quasinewton.prj
index c40b6ec7842a92fb3285e522bc983825ba516c6a..f6558166ef5d4bbd920e9615ca9bf6bcb2085ff8 100644
--- a/Tests/Data/RichardsMechanics/RichardsFlow_2d_quasinewton.prj
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_quasinewton.prj
@@ -179,6 +179,9 @@
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
                 <variable>saturation</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>saturation_avg</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small.prj b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small.prj
index d5e6a0c33ad16f76626a1a97ad9435b723597dbc..e6c7038efccc4fda6d2688893be29bdac55bc42b 100644
--- a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small.prj
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small.prj
@@ -179,6 +179,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/RichardsMechanics/RichardsFlow_2d_small_masslumping.prj b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping.prj
index d115b1c322804ef33cbc4baeb0d51636dbe02fb5..dab3766064b675e176c59b29e4199f6c49fabebb 100644
--- a/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping.prj
+++ b/Tests/Data/RichardsMechanics/RichardsFlow_2d_small_masslumping.prj
@@ -175,6 +175,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/RichardsMechanics/alternative_mass_balance_anzInterval_10.prj b/Tests/Data/RichardsMechanics/alternative_mass_balance_anzInterval_10.prj
index 301c1be4a6559d81510aefa7afdf4d8b67580a12..463e805f0cf242178efdc3fd0e7789c2b1ea43d1 100644
--- a/Tests/Data/RichardsMechanics/alternative_mass_balance_anzInterval_10.prj
+++ b/Tests/Data/RichardsMechanics/alternative_mass_balance_anzInterval_10.prj
@@ -148,6 +148,9 @@
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
                 <variable>saturation</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>pressure_interpolated</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law.prj b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law.prj
index 1bbfae01321949d7785cb8b9ee92dc3be465e812..c46de104167bc822891913643c0a1942c87b2ebb 100644
--- a/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law.prj
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_power_law.prj
@@ -162,6 +162,7 @@
                 <variable>velocity</variable>
                 <variable>saturation</variable>
                 <variable>porosity</variable>
+                <variable>porosity_avg</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff.prj b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff.prj
index d9eb518062e8559e18ef66272b4b2800d547b1ec..f7ac54f318ffc365f0008ab94071cc3918f65012 100644
--- a/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff.prj
+++ b/Tests/Data/RichardsMechanics/bishops_effective_stress_saturation_cutoff.prj
@@ -162,6 +162,7 @@
                 <variable>velocity</variable>
                 <variable>saturation</variable>
                 <variable>porosity</variable>
+                <variable>porosity_avg</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/RichardsMechanics/confined_compression_fully_saturated.prj b/Tests/Data/RichardsMechanics/confined_compression_fully_saturated.prj
index fc16d2711675a1db75182feb2ac99ddc979ef2bb..64c184ebb301e337db7cd742d3d092e3a9653c29 100644
--- a/Tests/Data/RichardsMechanics/confined_compression_fully_saturated.prj
+++ b/Tests/Data/RichardsMechanics/confined_compression_fully_saturated.prj
@@ -170,6 +170,8 @@
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
                 <variable>saturation</variable>
+                <variable>NodalForces</variable>
+                <variable>pressure_interpolated</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/RichardsMechanics/deformation_dependent_porosity.prj b/Tests/Data/RichardsMechanics/deformation_dependent_porosity.prj
index d3edf8356999b20873ce08840f2618f36e9eb440..339e02108d3e89e2980206cab056b1686edf8cf2 100644
--- a/Tests/Data/RichardsMechanics/deformation_dependent_porosity.prj
+++ b/Tests/Data/RichardsMechanics/deformation_dependent_porosity.prj
@@ -179,6 +179,7 @@
                 <variable>velocity</variable>
                 <variable>saturation</variable>
                 <variable>porosity</variable>
+                <variable>porosity_avg</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/RichardsMechanics/deformation_dependent_porosity_swelling.prj b/Tests/Data/RichardsMechanics/deformation_dependent_porosity_swelling.prj
index 8e1cca41544274a2a65d054bde84efe7983c8509..ddc948c1989f655bc2dd82f1f040816879d62820 100644
--- a/Tests/Data/RichardsMechanics/deformation_dependent_porosity_swelling.prj
+++ b/Tests/Data/RichardsMechanics/deformation_dependent_porosity_swelling.prj
@@ -208,6 +208,7 @@
                 <variable>porosity</variable>
                 <variable>transport_porosity</variable>
                 <variable>dry_density_solid</variable>
+                <variable>porosity_avg</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/RichardsMechanics/flow_fully_saturated.prj b/Tests/Data/RichardsMechanics/flow_fully_saturated.prj
index 10fc7b1301842d43588a226b750a1de9e115cb1a..05fc576b85ef8ceb394bd49acd579669a95dd2ae 100644
--- a/Tests/Data/RichardsMechanics/flow_fully_saturated.prj
+++ b/Tests/Data/RichardsMechanics/flow_fully_saturated.prj
@@ -147,6 +147,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/RichardsMechanics/flow_fully_saturated_anisotropic.prj b/Tests/Data/RichardsMechanics/flow_fully_saturated_anisotropic.prj
index 181d0f9424beaf10c5ce1ac7d2c7d3c589f9c476..62adf715a5e0b4042cccff2a861ddbd507221220 100644
--- a/Tests/Data/RichardsMechanics/flow_fully_saturated_anisotropic.prj
+++ b/Tests/Data/RichardsMechanics/flow_fully_saturated_anisotropic.prj
@@ -147,6 +147,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/RichardsMechanics/flow_fully_saturated_coordinate_system.prj b/Tests/Data/RichardsMechanics/flow_fully_saturated_coordinate_system.prj
index 59e29426d16be397d60c178febebb550821a14d5..6517f8e7c743c48e42abba3e286729b881129b41 100644
--- a/Tests/Data/RichardsMechanics/flow_fully_saturated_coordinate_system.prj
+++ b/Tests/Data/RichardsMechanics/flow_fully_saturated_coordinate_system.prj
@@ -147,6 +147,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/RichardsMechanics/flow_fully_saturated_linear.prj b/Tests/Data/RichardsMechanics/flow_fully_saturated_linear.prj
index 81e69bf88307826d43a2804ca687771d447d321e..ce270f2ff4176f2d2d68be87d2a3570324715101 100644
--- a/Tests/Data/RichardsMechanics/flow_fully_saturated_linear.prj
+++ b/Tests/Data/RichardsMechanics/flow_fully_saturated_linear.prj
@@ -147,6 +147,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/RichardsMechanics/gravity.prj b/Tests/Data/RichardsMechanics/gravity.prj
index d1322cda7a8d20f5d3794892a06b3eb44363f803..d25eb7c26644f07b42ea574771ff123aa5e53b97 100644
--- a/Tests/Data/RichardsMechanics/gravity.prj
+++ b/Tests/Data/RichardsMechanics/gravity.prj
@@ -181,6 +181,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/RichardsMechanics/mechanics_linear.prj b/Tests/Data/RichardsMechanics/mechanics_linear.prj
index 738583f76b402b7217173df053547c4f1540bb8c..b4b600895eebe98ca51feb566b85b368b45cd021 100644
--- a/Tests/Data/RichardsMechanics/mechanics_linear.prj
+++ b/Tests/Data/RichardsMechanics/mechanics_linear.prj
@@ -155,6 +155,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/RichardsMechanics/mfront_restart_part1.prj b/Tests/Data/RichardsMechanics/mfront_restart_part1.prj
index 7a79360c41c46a4553c123e5637df45cb066574e..bac3da050a75f01e3c39222f41c6c380d09e47b6 100644
--- a/Tests/Data/RichardsMechanics/mfront_restart_part1.prj
+++ b/Tests/Data/RichardsMechanics/mfront_restart_part1.prj
@@ -169,6 +169,8 @@
                 <variable>velocity</variable>
                 <variable>saturation</variable>
                 <variable>porosity</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/RichardsMechanics/mfront_restart_part2.prj b/Tests/Data/RichardsMechanics/mfront_restart_part2.prj
index 79bba885bd6c4bd3a2fa68355cdb1d10f96f7b6c..676d1ed501bf72e0d5730573c63c0b5ce865ad1c 100644
--- a/Tests/Data/RichardsMechanics/mfront_restart_part2.prj
+++ b/Tests/Data/RichardsMechanics/mfront_restart_part2.prj
@@ -168,6 +168,8 @@
                 <variable>velocity</variable>
                 <variable>saturation</variable>
                 <variable>porosity</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/RichardsMechanics/orthotropic_power_law_permeability_xyz.prj b/Tests/Data/RichardsMechanics/orthotropic_power_law_permeability_xyz.prj
index 50e2c9567dc6b92715844db2f810270acd108052..ed93a51b8f4f6c235f7727e259ed035a51cf8893 100644
--- a/Tests/Data/RichardsMechanics/orthotropic_power_law_permeability_xyz.prj
+++ b/Tests/Data/RichardsMechanics/orthotropic_power_law_permeability_xyz.prj
@@ -176,6 +176,7 @@
                 <variable>velocity</variable>
                 <variable>saturation</variable>
                 <variable>porosity</variable>
+                <variable>porosity_avg</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/RichardsMechanics/orthotropic_swelling_xy.prj b/Tests/Data/RichardsMechanics/orthotropic_swelling_xy.prj
index e25d8d3de14bc424ebc70ccbfc2cdca61e0e74a5..3ea953afdb1406d4820210a859cd040a2967fe24 100644
--- a/Tests/Data/RichardsMechanics/orthotropic_swelling_xy.prj
+++ b/Tests/Data/RichardsMechanics/orthotropic_swelling_xy.prj
@@ -172,6 +172,7 @@
                 <variable>velocity</variable>
                 <variable>saturation</variable>
                 <variable>porosity</variable>
+                <variable>porosity_avg</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/RichardsMechanics/orthotropic_swelling_xyz.prj b/Tests/Data/RichardsMechanics/orthotropic_swelling_xyz.prj
index f2693dc7b9f72fb9c15fa9420ad0b1ef9ed2eb8a..c5e61fe66f582a5c10c1f736ed5cbcc74fd61647 100644
--- a/Tests/Data/RichardsMechanics/orthotropic_swelling_xyz.prj
+++ b/Tests/Data/RichardsMechanics/orthotropic_swelling_xyz.prj
@@ -175,6 +175,7 @@
                 <variable>velocity</variable>
                 <variable>saturation</variable>
                 <variable>porosity</variable>
+                <variable>porosity_avg</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/StokesFlow/2DPlanarPermeableFracture/2DPlanarPermeableFracture.prj b/Tests/Data/StokesFlow/2DPlanarPermeableFracture/2DPlanarPermeableFracture.prj
index 32c1a54a4334812c218605374f8745f729761528..adf48b2e9c1a6738c33bfdef0fb4f1c45c823b9c 100644
--- a/Tests/Data/StokesFlow/2DPlanarPermeableFracture/2DPlanarPermeableFracture.prj
+++ b/Tests/Data/StokesFlow/2DPlanarPermeableFracture/2DPlanarPermeableFracture.prj
@@ -102,6 +102,7 @@
             <variables>
                 <variable>liquid_velocity</variable>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/StokesFlow/ParallelPlate.prj b/Tests/Data/StokesFlow/ParallelPlate.prj
index 00b177e33ca83d139b4b7eb24f5017bd4172b260..b032bbb414371b2615e8e38498fc418db07428ce 100644
--- a/Tests/Data/StokesFlow/ParallelPlate.prj
+++ b/Tests/Data/StokesFlow/ParallelPlate.prj
@@ -73,6 +73,7 @@
             <variables>
                 <variable>liquid_velocity</variable>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/StokesFlow/ParallelPlateWithBodyForce.prj b/Tests/Data/StokesFlow/ParallelPlateWithBodyForce.prj
index 0dd416ae970b69e40ba0e2402fcc7ce243ea3cd7..0d464f3adc0cba4cf73559dd29265e2821320356 100644
--- a/Tests/Data/StokesFlow/ParallelPlateWithBodyForce.prj
+++ b/Tests/Data/StokesFlow/ParallelPlateWithBodyForce.prj
@@ -73,6 +73,7 @@
             <variables>
                 <variable>liquid_velocity</variable>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/TH2M/H2/dissolution_diffusion/bourgeat.prj b/Tests/Data/TH2M/H2/dissolution_diffusion/bourgeat.prj
index 382c7afa5e52780153f927a9420fbbd889665a49..d27e26c1793df4f32dd826cbd0a9c7f0bc3e6e49 100644
--- a/Tests/Data/TH2M/H2/dissolution_diffusion/bourgeat.prj
+++ b/Tests/Data/TH2M/H2/dissolution_diffusion/bourgeat.prj
@@ -328,6 +328,10 @@
                 <variable>xnCG</variable>
                 <variable>xmCG</variable>
                 <variable>xmWL</variable>
+                <variable>d_CG</variable>
+                <variable>d_CL</variable>
+                <variable>d_WG</variable>
+                <variable>d_WL</variable>
             </variables>
             <fixed_output_times>
                 1e10
diff --git a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square.prj b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square.prj
index 513444bbb8ee16d86c352d4bf005737039c25d20..1b1f09e5b902d76cb23919a1229ff4a8a8a07339 100644
--- a/Tests/Data/TH2M/H2M/OrthotropicSwelling/square.prj
+++ b/Tests/Data/TH2M/H2M/OrthotropicSwelling/square.prj
@@ -269,6 +269,7 @@
                 <variable>xmCG</variable>
                 <variable>xmWL</variable>
                 <variable>xnCG</variable>
+                <variable>liquid_pressure_interpolated</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
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/TH2M/HM/Confined_Compression/HM_confined_compression_gas.prj b/Tests/Data/TH2M/HM/Confined_Compression/HM_confined_compression_gas.prj
index 584bab055f1841735ddbe9b504c5d922701e2336..89a26eb1d50b72d11b89307a2373bce1026f017f 100644
--- a/Tests/Data/TH2M/HM/Confined_Compression/HM_confined_compression_gas.prj
+++ b/Tests/Data/TH2M/HM/Confined_Compression/HM_confined_compression_gas.prj
@@ -274,6 +274,10 @@
                 <variable>gas_density</variable>
                 <variable>porosity</variable>
                 <variable>saturation</variable>
+                <variable>GasMassFlowRate</variable>
+                <variable>HeatFlowRate</variable>
+                <variable>LiquidMassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/TH2M/HM/Confined_Compression/HM_confined_compression_liquid.prj b/Tests/Data/TH2M/HM/Confined_Compression/HM_confined_compression_liquid.prj
index 55e7aa778ee1ba49e4e0ed71b8a559d29b0b4d6f..a49b3214827361a1c51d2b940529da0471dd5f0f 100644
--- a/Tests/Data/TH2M/HM/Confined_Compression/HM_confined_compression_liquid.prj
+++ b/Tests/Data/TH2M/HM/Confined_Compression/HM_confined_compression_liquid.prj
@@ -281,6 +281,10 @@
                 <variable>gas_density</variable>
                 <variable>porosity</variable>
                 <variable>saturation</variable>
+                <variable>GasMassFlowRate</variable>
+                <variable>HeatFlowRate</variable>
+                <variable>LiquidMassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/TH2M/HM/flow_fully_saturated.prj b/Tests/Data/TH2M/HM/flow_fully_saturated.prj
index 2305a8d20c40e7296282ad0f0c648526b96e8f1b..cbc93139e2add5dc901f575789ad5f175ef59479 100644
--- a/Tests/Data/TH2M/HM/flow_fully_saturated.prj
+++ b/Tests/Data/TH2M/HM/flow_fully_saturated.prj
@@ -253,6 +253,11 @@
                 <variable>gas_density</variable>
                 <variable>porosity</variable>
                 <variable>saturation</variable>
+                <variable>GasMassFlowRate</variable>
+                <variable>HeatFlowRate</variable>
+                <variable>LiquidMassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>liquid_pressure_interpolated</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/TH2M/HM/flow_fully_saturated_gas.prj b/Tests/Data/TH2M/HM/flow_fully_saturated_gas.prj
index b1b1f8144290a31ba2100ada2a86d8375ddb5721..b5dbc9a9266ec53f93d1fae7e342281348d841d7 100644
--- a/Tests/Data/TH2M/HM/flow_fully_saturated_gas.prj
+++ b/Tests/Data/TH2M/HM/flow_fully_saturated_gas.prj
@@ -253,6 +253,11 @@
                 <variable>gas_density</variable>
                 <variable>porosity</variable>
                 <variable>saturation</variable>
+                <variable>GasMassFlowRate</variable>
+                <variable>HeatFlowRate</variable>
+                <variable>LiquidMassFlowRate</variable>
+                <variable>NodalForces</variable>
+                <variable>liquid_pressure_interpolated</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/TH2M/TH/idealGasLaw/compression_gas.prj b/Tests/Data/TH2M/TH/idealGasLaw/compression_gas.prj
index f5dbf9ddb68851994639a5260bf53f5c07d51f0d..da12c2468e5c91d7ae81f2c8e816bf3bee52f129 100644
--- a/Tests/Data/TH2M/TH/idealGasLaw/compression_gas.prj
+++ b/Tests/Data/TH2M/TH/idealGasLaw/compression_gas.prj
@@ -272,6 +272,8 @@
                 <variable>rhoL_hL_dot</variable>
                 <variable>k_rel_G</variable>
                 <variable>k_rel_L</variable>
+                <variable>xmCG</variable>
+                <variable>xnCG</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/TH2M/TH2/heatpipe/heat_pipe_rough.prj b/Tests/Data/TH2M/TH2/heatpipe/heat_pipe_rough.prj
index 6298b1b6c033eb0cb5ecf0a5b4f9b732f466c187..227cda034bcf420b56d3f909db784def9a3d3e13 100644
--- a/Tests/Data/TH2M/TH2/heatpipe/heat_pipe_rough.prj
+++ b/Tests/Data/TH2M/TH2/heatpipe/heat_pipe_rough.prj
@@ -378,6 +378,7 @@
                 <variable>rhoL_hL_dot</variable>
                 <variable>k_rel_G</variable>
                 <variable>k_rel_L</variable>
+                <variable>xnCG</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/TH2M/TH2/heatpipe/heat_pipe_strict.prj b/Tests/Data/TH2M/TH2/heatpipe/heat_pipe_strict.prj
index ce3aeb3e9632709c17765cb65f19574af2ea3ffc..0acd9a74a25d35ab0db531fe2a119f2d7f12dbe4 100644
--- a/Tests/Data/TH2M/TH2/heatpipe/heat_pipe_strict.prj
+++ b/Tests/Data/TH2M/TH2/heatpipe/heat_pipe_strict.prj
@@ -373,6 +373,7 @@
                 <variable>rhoL_hL_dot</variable>
                 <variable>k_rel_G</variable>
                 <variable>k_rel_L</variable>
+                <variable>xnCG</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/TH2M/TH2/unicube/unicube.prj b/Tests/Data/TH2M/TH2/unicube/unicube.prj
index 4b69dd3bd5755bd35b4c8799ebaa2da1b0e41e58..0bec798a7ef9ee6034eb31f06cb0f2dc430308c8 100644
--- a/Tests/Data/TH2M/TH2/unicube/unicube.prj
+++ b/Tests/Data/TH2M/TH2/unicube/unicube.prj
@@ -344,6 +344,7 @@
                 <variable>rhoL_hL_dot</variable>
                 <variable>k_rel_G</variable>
                 <variable>k_rel_L</variable>
+                <variable>xnCG</variable>
             </variables>
             <fixed_output_times>
                 6500
diff --git a/Tests/Data/TH2M/THM/Confined_Compression/THM_confined_compression_gas.prj b/Tests/Data/TH2M/THM/Confined_Compression/THM_confined_compression_gas.prj
index 213ff4ff85853f5889ed1905dab2b86ddae02f35..7295de98bcf54369aa66ad705962bedba06136ea 100644
--- a/Tests/Data/TH2M/THM/Confined_Compression/THM_confined_compression_gas.prj
+++ b/Tests/Data/TH2M/THM/Confined_Compression/THM_confined_compression_gas.prj
@@ -280,6 +280,10 @@
                 <variable>gas_density</variable>
                 <variable>porosity</variable>
                 <variable>saturation</variable>
+                <variable>GasMassFlowRate</variable>
+                <variable>HeatFlowRate</variable>
+                <variable>LiquidMassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/TH2M/THM/Confined_Compression/THM_confined_compression_liquid.prj b/Tests/Data/TH2M/THM/Confined_Compression/THM_confined_compression_liquid.prj
index 3b6c121bd357ed4b20f6d55f06fdb9250bcfa81e..47df0cf8dcde4087e11204d554445a8bd5753337 100644
--- a/Tests/Data/TH2M/THM/Confined_Compression/THM_confined_compression_liquid.prj
+++ b/Tests/Data/TH2M/THM/Confined_Compression/THM_confined_compression_liquid.prj
@@ -280,6 +280,10 @@
                 <variable>gas_density</variable>
                 <variable>porosity</variable>
                 <variable>saturation</variable>
+                <variable>GasMassFlowRate</variable>
+                <variable>HeatFlowRate</variable>
+                <variable>LiquidMassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/TH2M/submesh_residuum_assembly/T.xml b/Tests/Data/TH2M/submesh_residuum_assembly/T.xml
index e7dfe781ef1f5976345097121a6b8496e31ebd1f..d209d325f1032d5d189f3fedb34135d01612ce61 100644
--- a/Tests/Data/TH2M/submesh_residuum_assembly/T.xml
+++ b/Tests/Data/TH2M/submesh_residuum_assembly/T.xml
@@ -47,7 +47,20 @@
     <add sel="/*/time_loop/submesh_residuum_output">
         <prefix>T_{:meshname}</prefix>
     </add>
-    
+
+    <add sel="/*/time_loop/output/variables">
+        <variable>HeatFlowRate</variable>
+        <variable>GasMassFlowRate</variable>
+        <variable>LiquidMassFlowRate</variable>
+        <variable>NodalForces</variable>
+    </add>
+    <add sel="/*/time_loop/submesh_residuum_output/variables">
+        <variable>HeatFlowRate</variable>
+        <variable>GasMassFlowRate</variable>
+        <variable>LiquidMassFlowRate</variable>
+        <variable>NodalForces</variable>
+    </add>
+
     <add sel="/*">
         <test_definition>
             <vtkdiff>
diff --git a/Tests/Data/TH2M/submesh_residuum_assembly/p_G.xml b/Tests/Data/TH2M/submesh_residuum_assembly/p_G.xml
index 6635e23a25690ea6e1a7beabd9d7b360312c9dba..c9c204aab9b1fc165894dfd03d1f462ea68a85bc 100644
--- a/Tests/Data/TH2M/submesh_residuum_assembly/p_G.xml
+++ b/Tests/Data/TH2M/submesh_residuum_assembly/p_G.xml
@@ -81,6 +81,19 @@
         <prefix>p_G_{:meshname}</prefix>
     </add>
     
+    <add sel="/*/time_loop/output/variables">
+        <variable>HeatFlowRate</variable>
+        <variable>GasMassFlowRate</variable>
+        <variable>LiquidMassFlowRate</variable>
+        <variable>NodalForces</variable>
+    </add>
+    <add sel="/*/time_loop/submesh_residuum_output/variables">
+        <variable>HeatFlowRate</variable>
+        <variable>GasMassFlowRate</variable>
+        <variable>LiquidMassFlowRate</variable>
+        <variable>NodalForces</variable>
+    </add>
+
     <add sel="/*">
         <test_definition>
             <vtkdiff>
diff --git a/Tests/Data/TH2M/submesh_residuum_assembly/p_cap.xml b/Tests/Data/TH2M/submesh_residuum_assembly/p_cap.xml
index b2393d35734205e777829c41dda9f68a83a1f96e..bb662ee5e18199ecf339fcb7858de52a7f43a087 100644
--- a/Tests/Data/TH2M/submesh_residuum_assembly/p_cap.xml
+++ b/Tests/Data/TH2M/submesh_residuum_assembly/p_cap.xml
@@ -61,6 +61,19 @@
         <prefix>p_cap_{:meshname}</prefix>
     </add>
     
+    <add sel="/*/time_loop/output/variables">
+        <variable>HeatFlowRate</variable>
+        <variable>GasMassFlowRate</variable>
+        <variable>LiquidMassFlowRate</variable>
+        <variable>NodalForces</variable>
+    </add>
+    <add sel="/*/time_loop/submesh_residuum_output/variables">
+        <variable>HeatFlowRate</variable>
+        <variable>GasMassFlowRate</variable>
+        <variable>LiquidMassFlowRate</variable>
+        <variable>NodalForces</variable>
+    </add>
+
     <add sel="/*">
         <test_definition>
             <vtkdiff>
diff --git a/Tests/Data/TH2M/submesh_residuum_assembly/u.xml b/Tests/Data/TH2M/submesh_residuum_assembly/u.xml
index 7667a9c3cc9c7fbb8d3d7200694d6b8ef05164ad..1c98385ec8c288ad3b96d1060b6b9120b4fbc4ce 100644
--- a/Tests/Data/TH2M/submesh_residuum_assembly/u.xml
+++ b/Tests/Data/TH2M/submesh_residuum_assembly/u.xml
@@ -67,6 +67,19 @@
         <prefix>u_{:meshname}</prefix>
     </add>
     
+    <add sel="/*/time_loop/output/variables">
+        <variable>HeatFlowRate</variable>
+        <variable>GasMassFlowRate</variable>
+        <variable>LiquidMassFlowRate</variable>
+        <variable>NodalForces</variable>
+    </add>
+    <add sel="/*/time_loop/submesh_residuum_output/variables">
+        <variable>HeatFlowRate</variable>
+        <variable>GasMassFlowRate</variable>
+        <variable>LiquidMassFlowRate</variable>
+        <variable>NodalForces</variable>
+    </add>
+
     <add sel="/*">
         <test_definition>
             <vtkdiff>
diff --git a/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem.prj b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem.prj
index dfeccc53982bba5a447c89a80335b3d0a705c75b..e33691fcdf4b9a70da4da68fd00ba8a2ef80c96b 100644
--- a/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem.prj
+++ b/Tests/Data/ThermoHydroMechanics/1D_freezing_column_Stefan/Stefan_problem.prj
@@ -272,6 +272,7 @@
             <variables>
                 <variable> temperature </variable>
                 <variable> heat_flux </variable>
+                <variable>HeatFlowRate</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
         </output>
diff --git a/Tests/Data/ThermoHydroMechanics/A2/A2.prj b/Tests/Data/ThermoHydroMechanics/A2/A2.prj
index 9536ebba2a1b2d98ed4b2fde039568bd919ecf0d..d2950d807fa8c1bce65c7a9eac722a51478aa114 100644
--- a/Tests/Data/ThermoHydroMechanics/A2/A2.prj
+++ b/Tests/Data/ThermoHydroMechanics/A2/A2.prj
@@ -219,6 +219,9 @@
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
+                <variable>HeatFlowRate</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/ThermoHydroMechanics/A2/A2_heating.prj b/Tests/Data/ThermoHydroMechanics/A2/A2_heating.prj
index 9a14db99eabcb173dcdf4c9c803a98b4573a789d..7ad7e818a10884065837a619e7cf07be21029630 100644
--- a/Tests/Data/ThermoHydroMechanics/A2/A2_heating.prj
+++ b/Tests/Data/ThermoHydroMechanics/A2/A2_heating.prj
@@ -178,6 +178,9 @@
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
                 <variable>velocity</variable>
+                <variable>HeatFlowRate</variable>
+                <variable>MassFlowRate</variable>
+                <variable>NodalForces</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/ThermoHydroMechanics/Linear/HeatTransportInStationaryFlow/HeatTransportInStationaryFlow.prj b/Tests/Data/ThermoHydroMechanics/Linear/HeatTransportInStationaryFlow/HeatTransportInStationaryFlow.prj
index 66e91e05a2336fb5633e3a8da313eba0bd5b7272..3fff86479df22cbb291d3f07cfa1c81aa568e994 100644
--- a/Tests/Data/ThermoHydroMechanics/Linear/HeatTransportInStationaryFlow/HeatTransportInStationaryFlow.prj
+++ b/Tests/Data/ThermoHydroMechanics/Linear/HeatTransportInStationaryFlow/HeatTransportInStationaryFlow.prj
@@ -149,7 +149,9 @@
             </timesteps>
             <variables>
                 <variable>temperature</variable>
+                <variable>temperature_interpolated</variable>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
                 <variable>velocity</variable>
             </variables>
             <suffix>_ts_{:timestep}_t_{:time}</suffix>
diff --git a/Tests/Data/ThermoHydroMechanics/Linear/Point_injection/pointheatsource_quadratic-mesh.prj b/Tests/Data/ThermoHydroMechanics/Linear/Point_injection/pointheatsource_quadratic-mesh.prj
index 136a5aa1dcf7406dc056ffa3979d11b31c3e2a4a..a773c71c7c3a247370524bd82f3da7fb6d1eee94 100644
--- a/Tests/Data/ThermoHydroMechanics/Linear/Point_injection/pointheatsource_quadratic-mesh.prj
+++ b/Tests/Data/ThermoHydroMechanics/Linear/Point_injection/pointheatsource_quadratic-mesh.prj
@@ -154,7 +154,9 @@
             <variables>
                 <variable>displacement</variable>
                 <variable>pressure</variable>
+                <variable>pressure_interpolated</variable>
                 <variable>temperature</variable>
+                <variable>temperature_interpolated</variable>
                 <variable>sigma</variable>
                 <variable>epsilon</variable>
             </variables>
diff --git a/Tests/Data/ThermoRichardsFlow/SimplifiedMechanics/TRhyd_unsaturated_bishopstest.prj b/Tests/Data/ThermoRichardsFlow/SimplifiedMechanics/TRhyd_unsaturated_bishopstest.prj
index 442f3cc073e7d1d24a1b15d84abb73910e1970dc..60c3aea8b8f4c8fb043fcf7e0d648fab55a70ccc 100644
--- a/Tests/Data/ThermoRichardsFlow/SimplifiedMechanics/TRhyd_unsaturated_bishopstest.prj
+++ b/Tests/Data/ThermoRichardsFlow/SimplifiedMechanics/TRhyd_unsaturated_bishopstest.prj
@@ -188,6 +188,7 @@
                 <variable>temperature</variable>
                 <variable>velocity</variable>
                 <variable>saturation</variable>
+                <variable>MassFlowRate</variable>
             </variables>
         </output>
     </time_loop>
diff --git a/Tests/Data/ThermoRichardsFlow/SimplifiedMechanics/TRuni_unsaturated_bishopstest.prj b/Tests/Data/ThermoRichardsFlow/SimplifiedMechanics/TRuni_unsaturated_bishopstest.prj
index 5d651ea7bab3e22b8513daa2ab878af79a167d99..d5920915fc705cdecffd68b204ea0a45efc58a40 100644
--- a/Tests/Data/ThermoRichardsFlow/SimplifiedMechanics/TRuni_unsaturated_bishopstest.prj
+++ b/Tests/Data/ThermoRichardsFlow/SimplifiedMechanics/TRuni_unsaturated_bishopstest.prj
@@ -188,6 +188,7 @@
                 <variable>temperature</variable>
                 <variable>velocity</variable>
                 <variable>saturation</variable>
+                <variable>MassFlowRate</variable>
             </variables>
         </output>
     </time_loop>
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>