From 3d1ed8274fb45cbb812430d132e85e270711f576 Mon Sep 17 00:00:00 2001
From: Christian Silbermann <christian.silbermann@ifgt.tu-freiberg.de>
Date: Tue, 23 Apr 2024 13:46:36 +0200
Subject: [PATCH] [T/THM] Test setup, data and time stepping modified

---
 .../RestartMCC/mfront_restart_part1.prj       | 46 ++++++------
 .../RestartMCC/mfront_restart_part1.pvd       |  8 ++-
 .../mfront_restart_part1_t_1000.000000.vtu    | 70 +++++++++---------
 .../RestartMCC/mfront_restart_part2.pvd       | 11 +--
 .../RestartMCC/mfront_restart_part2.xml       |  7 +-
 .../mfront_restart_part2_t_1000.000000.vtu    | 70 +++++++++---------
 .../mfront_restart_part2_t_2000.000000.vtu    | 72 +++++++++----------
 7 files changed, 143 insertions(+), 141 deletions(-)

diff --git a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1.prj b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1.prj
index a124c3079e7..d9a785078df 100644
--- a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1.prj
+++ b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1.prj
@@ -128,20 +128,21 @@
                 <convergence_criterion>
                     <type>PerComponentDeltaX</type>
                     <norm_type>NORM2</norm_type>
-                    <abstols>1e-15 1e-9 1e-14 1e-14</abstols>
+                    <abstols>1e-15 1e-10 1e-14 1e-14</abstols>
                 </convergence_criterion>
                 <time_discretization>
                     <type>BackwardEuler</type>
                 </time_discretization>
                 <time_stepping>
-                    <type>IterationNumberBasedTimeStepping</type>
+                    <type>FixedTimeStepping</type>
                     <t_initial>0</t_initial>
                     <t_end>1000</t_end>
-                    <initial_dt>100</initial_dt>
-                    <minimum_dt>10</minimum_dt>
-                    <maximum_dt>1000</maximum_dt>
-                    <number_iterations>1   5   10   15</number_iterations>
-                    <multiplier>1.2 1.0 0.8 0.6</multiplier>
+                    <timesteps>
+                        <pair>
+                            <repeat>10</repeat>
+                            <delta_t>100</delta_t>
+                        </pair>
+                    </timesteps>
                 </time_stepping>
             </process>
         </processes>
@@ -196,15 +197,16 @@
             <type>Constant</type>
             <value>1.76</value>
         </parameter>
-        <!-- Initial and boundary values -->
+        <!-- Initial effective stress field -->
         <parameter>
             <name>sigma0</name>
             <type>Function</type>
-            <expression>-5e3</expression>
-            <expression>-5e3</expression>
-            <expression>-5e3</expression>
+            <expression>-5e4</expression>
+            <expression>-5e4</expression>
+            <expression>-5e4</expression>
             <expression>0</expression>
         </parameter>
+        <!-- Initial and boundary values -->
         <parameter>
             <name>displacement0</name>
             <type>Constant</type>
@@ -221,14 +223,14 @@
             <value>293.15</value>
         </parameter>
         <parameter>
-            <name>top_pressure</name>
+            <name>total_axial_stress</name>
             <type>Constant</type>
-            <value>-2.7e5</value>
+            <value>-3.35e5</value>
         </parameter>
         <parameter>
             <name>pressure_ic</name>
             <type>Constant</type>
-            <value>-5e3</value>
+            <value>5e4</value>
         </parameter>
     </parameters>
     <process_variables>
@@ -274,6 +276,13 @@
                     <component>0</component>
                     <parameter>zero</parameter>
                 </boundary_condition>
+                <!--fix right in radial direction-->
+                <boundary_condition>
+                    <mesh>square_1x1_quad_1e0_right</mesh>
+                    <type>Dirichlet</type>
+                    <component>0</component>
+                    <parameter>zero</parameter>
+                </boundary_condition>
                 <!--fix bottom in axial direction-->
                 <boundary_condition>
                     <mesh>square_1x1_quad_1e0_bottom</mesh>
@@ -286,14 +295,7 @@
                     <mesh>square_1x1_quad_1e0_top</mesh>
                     <type>Neumann</type>
                     <component>1</component>
-                    <parameter>top_pressure</parameter>
-                </boundary_condition>
-                <!--compression in -radial direction-->
-                <boundary_condition>
-                    <mesh>square_1x1_quad_1e0_right</mesh>
-                    <type>Dirichlet</type>
-                    <component>0</component>
-                    <parameter>zero</parameter>
+                    <parameter>total_axial_stress</parameter>
                 </boundary_condition>
             </boundary_conditions>
         </process_variable>
diff --git a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1.pvd b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1.pvd
index b1b1f8a15f1..6e8cc7fd77b 100644
--- a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1.pvd
+++ b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1.pvd
@@ -6,9 +6,11 @@
     <DataSet timestep="200" group="" part="0" file="mfront_restart_part1_t_200.000000.vtu"/>
     <DataSet timestep="300" group="" part="0" file="mfront_restart_part1_t_300.000000.vtu"/>
     <DataSet timestep="400" group="" part="0" file="mfront_restart_part1_t_400.000000.vtu"/>
-    <DataSet timestep="520" group="" part="0" file="mfront_restart_part1_t_520.000000.vtu"/>
-    <DataSet timestep="664" group="" part="0" file="mfront_restart_part1_t_664.000000.vtu"/>
-    <DataSet timestep="836.8" group="" part="0" file="mfront_restart_part1_t_836.800000.vtu"/>
+    <DataSet timestep="500" group="" part="0" file="mfront_restart_part1_t_500.000000.vtu"/>
+    <DataSet timestep="600" group="" part="0" file="mfront_restart_part1_t_600.000000.vtu"/>
+    <DataSet timestep="700" group="" part="0" file="mfront_restart_part1_t_700.000000.vtu"/>
+    <DataSet timestep="800" group="" part="0" file="mfront_restart_part1_t_800.000000.vtu"/>
+    <DataSet timestep="900" group="" part="0" file="mfront_restart_part1_t_900.000000.vtu"/>
     <DataSet timestep="1000" group="" part="0" file="mfront_restart_part1_t_1000.000000.vtu"/>
   </Collection>
 </VTKFile>
diff --git a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1_t_1000.000000.vtu b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1_t_1000.000000.vtu
index ec13ecf3e1e..751948e2027 100644
--- a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1_t_1000.000000.vtu
+++ b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part1_t_1000.000000.vtu
@@ -3,54 +3,54 @@
   <UnstructuredGrid>
     <FieldData>
       <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="761" format="appended" RangeMin="34"                   RangeMax="125"                  offset="0"                   />
-      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="25" format="appended" RangeMin="45"                   RangeMax="121"                  offset="284"                 />
-      <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.012104037981"       RangeMax="0.012104037981"       offset="372"                 />
-      <DataArray type="Float64" Name="epsilon_m_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.012104037981"       RangeMax="0.012104037981"       offset="484"                 />
-      <DataArray type="Float64" Name="material_state_variable_ElasticStrain_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.0034383732392"      RangeMax="0.0034383732392"      offset="596"                 />
-      <DataArray type="Float64" Name="material_state_variable_EquivalentPlasticStrain_ip" NumberOfTuples="4" format="appended" RangeMin="0.0089577474413"      RangeMax="0.0089577474413"      offset="728"                 />
-      <DataArray type="Float64" Name="material_state_variable_PlasticVolumetricStrain_ip" NumberOfTuples="4" format="appended" RangeMin="-0.0073892070871"     RangeMax="-0.0073892070871"     offset="800"                 />
-      <DataArray type="Float64" Name="material_state_variable_PreConsolidationPressure_ip" NumberOfTuples="4" format="appended" RangeMin="242368.59489"         RangeMax="242368.59489"         offset="872"                 />
-      <DataArray type="Float64" Name="material_state_variable_VolumeRatio_ip" NumberOfTuples="4" format="appended" RangeMin="1.7388253014"         RangeMax="1.7388253014"         offset="944"                 />
-      <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="371354.0957"          RangeMax="371354.0957"          offset="1016"                />
+      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="19" format="appended" RangeMin="45"                   RangeMax="103"                  offset="284"                 />
+      <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.012518851362"       RangeMax="0.012518851362"       offset="364"                 />
+      <DataArray type="Float64" Name="epsilon_m_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.012518851362"       RangeMax="0.012518851362"       offset="472"                 />
+      <DataArray type="Float64" Name="material_state_variable_ElasticStrain_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.0035517417046"      RangeMax="0.0035517417046"      offset="588"                 />
+      <DataArray type="Float64" Name="material_state_variable_EquivalentPlasticStrain_ip" NumberOfTuples="4" format="appended" RangeMin="0.0095861092158"      RangeMax="0.0095861092158"      offset="716"                 />
+      <DataArray type="Float64" Name="material_state_variable_PlasticVolumetricStrain_ip" NumberOfTuples="4" format="appended" RangeMin="-0.0086274863678"     RangeMax="-0.0086274863678"     offset="792"                 />
+      <DataArray type="Float64" Name="material_state_variable_PreConsolidationPressure_ip" NumberOfTuples="4" format="appended" RangeMin="250968.14442"         RangeMax="250968.14442"         offset="864"                 />
+      <DataArray type="Float64" Name="material_state_variable_VolumeRatio_ip" NumberOfTuples="4" format="appended" RangeMin="1.7381041629"         RangeMax="1.7381041629"         offset="940"                 />
+      <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="387298.31369"         RangeMax="387298.31369"         offset="1012"                />
     </FieldData>
     <Piece NumberOfPoints="4"                    NumberOfCells="1"                   >
       <PointData>
-        <DataArray type="Float64" Name="ElasticStrain" NumberOfComponents="4" format="appended" RangeMin="0.0034383732392"      RangeMax="0.0034383732392"      offset="1144"                />
-        <DataArray type="Float64" Name="EquivalentPlasticStrain" format="appended" RangeMin="0.0089577474413"      RangeMax="0.0089577474413"      offset="1288"                />
+        <DataArray type="Float64" Name="ElasticStrain" NumberOfComponents="4" format="appended" RangeMin="0.0035517417046"      RangeMax="0.0035517417046"      offset="1144"                />
+        <DataArray type="Float64" Name="EquivalentPlasticStrain" format="appended" RangeMin="0.0095861092158"      RangeMax="0.0095861092158"      offset="1288"                />
         <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0"                   RangeMax="5.6843418861e-14"     offset="1364"                />
-        <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-9.689962817e-09"     RangeMax="-3.1871269002e-21"    offset="1436"                />
-        <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="282743.33882"         RangeMax="780991.55066"         offset="1524"                />
-        <DataArray type="Float64" Name="PlasticVolumetricStrain" format="appended" RangeMin="-0.0073892070871"     RangeMax="-0.0073892070871"     offset="1652"                />
-        <DataArray type="Float64" Name="PreConsolidationPressure" format="appended" RangeMin="242368.59489"         RangeMax="242368.59489"         offset="1728"                />
-        <DataArray type="Float64" Name="VolumeRatio" format="appended" RangeMin="1.7388253014"         RangeMax="1.7388253014"         offset="1804"                />
-        <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="0.012104037981"       offset="1880"                />
-        <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.012104037981"       RangeMax="0.012104037981"       offset="1956"                />
-        <DataArray type="Float64" Name="epsilon_m" NumberOfComponents="4" format="appended" RangeMin="0.012104037981"       RangeMax="0.012104037981"       offset="2076"                />
-        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="-5000"                RangeMax="-4997.6866918"        offset="2196"                />
-        <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="-5000"                RangeMax="-4997.6866918"        offset="2272"                />
-        <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="371354.0957"          RangeMax="371354.0957"          offset="2348"                />
-        <DataArray type="Float64" Name="temperature" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2480"                />
-        <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2548"                />
-        <DataArray type="Float64" Name="velocity" NumberOfComponents="2" format="appended" RangeMin="2.313308221e-09"      RangeMax="2.313308221e-09"      offset="2616"                />
+        <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-2.4136393841e-09"    RangeMax="6.4085817958e-21"     offset="1436"                />
+        <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="350811.17965"         RangeMax="1019491.2018"         offset="1532"                />
+        <DataArray type="Float64" Name="PlasticVolumetricStrain" format="appended" RangeMin="-0.0086274863678"     RangeMax="-0.0086274863678"     offset="1656"                />
+        <DataArray type="Float64" Name="PreConsolidationPressure" format="appended" RangeMin="250968.14442"         RangeMax="250968.14442"         offset="1728"                />
+        <DataArray type="Float64" Name="VolumeRatio" format="appended" RangeMin="1.7381041629"         RangeMax="1.7381041629"         offset="1804"                />
+        <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="0.012518851362"       offset="1880"                />
+        <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.012518851362"       RangeMax="0.012518851362"       offset="1956"                />
+        <DataArray type="Float64" Name="epsilon_m" NumberOfComponents="4" format="appended" RangeMin="0.012518851362"       RangeMax="0.012518851362"       offset="2080"                />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="50000"                RangeMax="50000.576214"         offset="2204"                />
+        <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="50000"                RangeMax="50000.576214"         offset="2280"                />
+        <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="387298.31369"         RangeMax="387298.31369"         offset="2356"                />
+        <DataArray type="Float64" Name="temperature" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2488"                />
+        <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2556"                />
+        <DataArray type="Float64" Name="velocity" NumberOfComponents="2" format="appended" RangeMin="5.7621395825e-10"     RangeMax="5.7621395826e-10"     offset="2624"                />
       </PointData>
       <CellData>
-        <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2724"                />
-        <DataArray type="Float64" Name="fluid_density_avg" format="appended" RangeMin="1000"                 RangeMax="1000"                 offset="2784"                />
-        <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="371354.0957"          RangeMax="371354.0957"          offset="2848"                />
-        <DataArray type="UInt8" Name="temperature_active" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2940"                />
-        <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001"                RangeMax="0.001"                offset="2996"                />
+        <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2732"                />
+        <DataArray type="Float64" Name="fluid_density_avg" format="appended" RangeMin="1000"                 RangeMax="1000"                 offset="2792"                />
+        <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="387298.31369"         RangeMax="387298.31369"         offset="2856"                />
+        <DataArray type="UInt8" Name="temperature_active" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2948"                />
+        <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001"                RangeMax="0.001"                offset="3004"                />
       </CellData>
       <Points>
-        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.4142135624"         offset="3064"                />
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.4142135624"         offset="3072"                />
       </Points>
       <Cells>
-        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="3136"                />
-        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="3208"                />
-        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="3268"                />
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="3144"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="3216"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="3276"                />
       </Cells>
     </Piece>
   </UnstructuredGrid>
   <AppendedData encoding="base64">
-   _AQAAAAAAAAAAgAAAAAAAAPkCAAAAAAAAtAAAAAAAAAA=eF61jz0LwkAMhv/LzV0Up67iLgouIkfaxhK4u9RcriCl/93zY3FRUbvlzUuehwyGgmIroMTBdpyTBRE4R1Puh6eSpUEx5bwwATya0kRqPVjqTN4kX6FYPtqafccBg2bAYizeILCL5HLh/4P5CeJBUQicjZon20MOlUO7chCV6q0K0FSKU6IeXMasP5XNvpY9FDt2yaPKxDLBJYfIjprbZc4xJsFpbPefNter14LDeAEk3BwxAQAAAAAAAAAAgAAAAAAAABkAAAAAAAAAIQAAAAAAAAA=eF4z0zPVM9Q1NDHUTU9JTE41NEox0kvJLCqpBABQQQcwAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAMwAAAAAAAAA=eF5jYICA6pc2VQynOvZDuQxrd7/T2cCnvAfGR5d32x49SX2fI1y+hoB+dHmYfgDBtSCFAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAMwAAAAAAAAA=eF5jYICAqpc2VQynOvZDuQxrd7/T2cCnvAfGR5d32x49SX2fI1y+hoB+dHmYfgDA5SCDAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAQQAAAAAAAAA=eF6zcV67y/yW5/5H02LCpy3J2m8D5R9rPSGy5BnTHhsc8i+cY1/PChXaYw/lP4DKw/jXoPpxyTO4QPQDAPP7QqA=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF5L7S9uKohosk+F0uloNADnEA3NAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF4z2tmio+pSt98ISpui0QDTUA0LAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF5j+x1yiGUaryMblGZHowHL+Av/AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5zNV9TZXXpt70rDhoA+jgPoQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPwAAAAAAAAA=eF5LCBbc2NLFenCWwLT42ycEDiZA+ZwWOpvD7vbvTcAhv2W1yNzLL+btTYTyZ0LlYXwVqH5c8neh+gG64jxlAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAASQAAAAAAAAA=eF6zcV67y/yW5/7H02LCpy3J2m8D5W+7lFBbf/rHbico/xFUHsY/AZW3gPKfQOVhfKGgTDHeaaJ77KD8e1B5GN8FKg8AU2FB7g==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF5L7S9uKohoss+E0ilQOg1KAwDnOA3NAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYACDBgYU2sAWxgcAGaUB7g==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAIQAAAAAAAAA=eF5jYAACqXW7QRSD7LbdvMdL2GXPm+zjA9Mu+wB1MQl9AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAPQAAAAAAAAA=eF6b1MZ42O7SVVuxgi9Rdx0FHa2Sr5VmFigclADzFR1trawN/JVm2YpC5A8uWN86NwUoLwmRPwgAiLQcJQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF4z2tmio+pSt98SShtBaRMoDQDTwA0OAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF5j/R1yiGUaryMXlGaB0qxQGgDL+Av9AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF5zNV9TZXXpt707lHaB0k5QGgD6QA+fAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGAAAAAAAAAA=eF5jYCAOVL20qWI41bEfxq+F8gFnRAhYAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOAAAAAAAAAA=eF5jYICAqpc2VQynOvYzwMBcSe0XFXN2w7gN6PIrUOUr0OT/Mk1qcKlz2wPjo5v/DyoPAG6kH5M=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOAAAAAAAAAA=eF5jYICAipc2VQynOvYzwMBcSe0XFXN2w7gN6PIrUOXR9f9lmtTgUue2B8avRpP/B5UHAG3MH5I=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF7j8OU4vb518wEuKM0AAh0IGgDJ4guXAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF7j8OU4vb518wEuKM0AAh0IGgDJ4guXAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAQgAAAAAAAAA=eF6LDxbc2NLFenCWwLT42ycEDsZD+S3NV3K7nrXsTUWTh/EnQ+UT0ORh/E9OD+7LGC3aCzNvOpr5ss4QeQBeUj0FAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAALwAAAAAAAAA=eF7bxsK/xuvvdOt3XVEHWe4r222F8kVrIXy+r2vfXPiUsxsmLwflC0DlAeqvH8E=AQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADgAAAAAAAAA=eF5jYAACh34HAAImARA=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF5LDhbc2NLFenC2wLT42ycEDiZB+YkacjYGEw/uBQDiMA3OAQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF5jAAAAAQABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF77s/LjJd+kAHsAGYAEpw==AQAAAAAAAAAAgAAAAAAAAGAAAAAAAAAAFQAAAAAAAAA=eF5jYMAHPtjjlcaQh/ER4gCW5AS9AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYIAARijNDKWZoDQAAHgABw==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYYAAAAAoAAU=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jBAAACgAK
+   _AQAAAAAAAAAAgAAAAAAAAPkCAAAAAAAAtAAAAAAAAAA=eF61jz0LwkAMhv/LzV0Up67iLgouIkfaxhK4u9RcriCl/93zY3FRUbvlzUuehwyGgmIroMTBdpyTBRE4R1Puh6eSpUEx5bwwATya0kRqPVjqTN4kX6FYPtqafccBg2bAYizeILCL5HLh/4P5CeJBUQicjZon20MOlUO7chCV6q0K0FSKU6IeXMasP5XNvpY9FDt2yaPKxDLBJYfIjprbZc4xJsFpbPefNter14LDeAEk3BwxAQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPVM9Q1NDHQTbc0M09NMU1MBQApRAS5AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAALwAAAAAAAAA=eF5jYICAhRrFJ6sXd+6Hchmkprzc25/81RrGR5f3c2kt+yDBaYNLnlj9AOozIE0=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAANAAAAAAAAAA=eF5jYICARRrFJ6sXd+6Hchmkprzc25/81RrGR5f3c2kt+yDBaQPjLyagH10eph8A66MgUw==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPwAAAAAAAAA=eF7LtC+SX/Daez+vVwHP4bjs/ZlQPr8Ms5vDLEabTBzyGc1pVfyHZG3yoHw+qHwGlK8A1Y9LfgVUPwAQNTREAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF4rYJx/5uHCZvsCKF2ERgMADgMQhQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF4rioqf57m6cX8RlC5EowEAdw/zAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF7L23BM/eAiPsc8KJ2PRgMA/2AO9w==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF7Lsuv1dzv/2z4LBw0A9gQPTQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAQAAAAAAAAAA=eF7LF/j96dsitoOvq3/dkEsVPJgP5bf9aln8Y02vbT4O+doY74pu0ZW2hVD+G6h8AZS/CKofl/xmqH4A0sJJVQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAASwAAAAAAAAA=eF7LsC+SX/Daez+vVwHP4bjs/dlQ/guT3bePiF2yLoLy+aDyWVD+AWOIfCqUzw2Vz4DyWfSvLX+4TckmD8rngcrD1PtA5QH5FzuJAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF4rYJx/5uHCZvsSKJ0PpQuhNAAOKxCFAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFAAAAAAAAAA=eF5jYAABA1sGCGhgQOMDABf6Ads=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJQAAAAAAAAA=eF5jYGBgONByYjeQYmhw3metduS07rLdIvsqnUG0yj4ArXQMkA==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAPAAAAAAAAAA=eF5r+eUd3eR5wVY59dv2nExRR49JeiU9k9QOqoP5qo7CBh69Sab3bFUh8gcX5dvYdwPl1SDyBwGZWRyBAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF4rioqf57m6cX8RGp0PpQEAdw/yAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF4r2HBM/eAiPsciKJ2DRgMA/7gO9w==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF7Lsuv1dzv/2z4XSmdC6QwoDQD2LA9NAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAFgAAAAAAAAA=eF5jYCAOLNQoPlm9uHM/Oh8Aa7oI1w==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOgAAAAAAAAA=eF5jYICAhRrFJ6sXd+6Hchm+nK/q2L/hqTWMvxhN/hea/AI0+dio6OKOOD4bGH8emnw8VB4AsPohxQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOgAAAAAAAAA=eF5jYICAxRrFJ6sXd+6Hchm+nK/q2L/hqTWMvwxN/hea/AI0+dio6OKOOD4bGH8hmnw8VB4AszohzQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFwAAAAAAAAA=eF5zd4woEMp64eAGpRlAAIkGAKD+CQw=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFwAAAAAAAAA=eF5zd4woEMp64eAGpRlAAIkGAKD+CQw=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAQgAAAAAAAAA=eF7LE/j96dsitoOvqn/dkEsVPJgH5e+Ivdv1pC7OthjKfwuVL4LyN0RD5PNw6JfacLXmUs9623wc8tOg8gDYCkq8AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAALwAAAAAAAAA=eF4Turf1+AOug9ZbYu9V551hthOC8v0W3QHzK+3b3E/eWAuXL4HyvaDyACrAIAY=AQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADgAAAAAAAAA=eF5jYAACh34HAAImARA=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF4rEPj96dsitoNvqn/dkEsVPFgA5Ts3Jyjve33EFgA0MhIaAQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF5jAAAAAQABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF77s/LjJd+kAHsAGYAEpw==AQAAAAAAAAAAgAAAAAAAAGAAAAAAAAAAFQAAAAAAAAA=eF5jYMAHPtjjlcaQh/ER4gCW5AS9AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYIAARijNDKWZoDQAAHgABw==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYYAAAAAoAAU=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jBAAACgAK
   </AppendedData>
 </VTKFile>
diff --git a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2.pvd b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2.pvd
index 6bcf394236d..f40de943638 100644
--- a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2.pvd
+++ b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2.pvd
@@ -4,10 +4,13 @@
     <DataSet timestep="1000" group="" part="0" file="mfront_restart_part2_t_1000.000000.vtu"/>
     <DataSet timestep="1100" group="" part="0" file="mfront_restart_part2_t_1100.000000.vtu"/>
     <DataSet timestep="1200" group="" part="0" file="mfront_restart_part2_t_1200.000000.vtu"/>
-    <DataSet timestep="1320" group="" part="0" file="mfront_restart_part2_t_1320.000000.vtu"/>
-    <DataSet timestep="1464" group="" part="0" file="mfront_restart_part2_t_1464.000000.vtu"/>
-    <DataSet timestep="1636.8" group="" part="0" file="mfront_restart_part2_t_1636.800000.vtu"/>
-    <DataSet timestep="1844.16" group="" part="0" file="mfront_restart_part2_t_1844.160000.vtu"/>
+    <DataSet timestep="1300" group="" part="0" file="mfront_restart_part2_t_1300.000000.vtu"/>
+    <DataSet timestep="1400" group="" part="0" file="mfront_restart_part2_t_1400.000000.vtu"/>
+    <DataSet timestep="1500" group="" part="0" file="mfront_restart_part2_t_1500.000000.vtu"/>
+    <DataSet timestep="1600" group="" part="0" file="mfront_restart_part2_t_1600.000000.vtu"/>
+    <DataSet timestep="1700" group="" part="0" file="mfront_restart_part2_t_1700.000000.vtu"/>
+    <DataSet timestep="1800" group="" part="0" file="mfront_restart_part2_t_1800.000000.vtu"/>
+    <DataSet timestep="1900" group="" part="0" file="mfront_restart_part2_t_1900.000000.vtu"/>
     <DataSet timestep="2000" group="" part="0" file="mfront_restart_part2_t_2000.000000.vtu"/>
   </Collection>
 </VTKFile>
diff --git a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2.xml b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2.xml
index 7756d863ab5..68f33fa990d 100644
--- a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2.xml
+++ b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2.xml
@@ -27,18 +27,13 @@
         <parameter>
             <name>pressure_bc</name>
             <type>Constant</type>
-            <value>-5e4</value>
+            <value>1e3</value>
         </parameter>
     </add>
 
     <remove sel="/*/process_variables/process_variable[name='pressure']/boundary_conditions"/>
     <add sel="/*/process_variables/process_variable[name='pressure']">
         <boundary_conditions>
-            <boundary_condition>
-                <mesh>square_1x1_quad_1e0_top</mesh>
-                <type>Dirichlet</type>
-                <parameter>pressure_bc</parameter>
-            </boundary_condition>
             <boundary_condition>
                 <mesh>square_1x1_quad_1e0_bottom</mesh>
                 <type>Dirichlet</type>
diff --git a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2_t_1000.000000.vtu b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2_t_1000.000000.vtu
index 6d7abb6a0dd..87ebf8775a3 100644
--- a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2_t_1000.000000.vtu
+++ b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2_t_1000.000000.vtu
@@ -3,54 +3,54 @@
   <UnstructuredGrid>
     <FieldData>
       <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="761" format="appended" RangeMin="34"                   RangeMax="125"                  offset="0"                   />
-      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="25" format="appended" RangeMin="45"                   RangeMax="121"                  offset="284"                 />
-      <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.012104037981"       RangeMax="0.012104037981"       offset="372"                 />
-      <DataArray type="Float64" Name="epsilon_m_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.012104037981"       RangeMax="0.012104037981"       offset="484"                 />
-      <DataArray type="Float64" Name="material_state_variable_ElasticStrain_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.0034383732392"      RangeMax="0.0034383732392"      offset="596"                 />
-      <DataArray type="Float64" Name="material_state_variable_EquivalentPlasticStrain_ip" NumberOfTuples="4" format="appended" RangeMin="0.0089577474413"      RangeMax="0.0089577474413"      offset="728"                 />
-      <DataArray type="Float64" Name="material_state_variable_PlasticVolumetricStrain_ip" NumberOfTuples="4" format="appended" RangeMin="-0.0073892070871"     RangeMax="-0.0073892070871"     offset="800"                 />
-      <DataArray type="Float64" Name="material_state_variable_PreConsolidationPressure_ip" NumberOfTuples="4" format="appended" RangeMin="242368.59489"         RangeMax="242368.59489"         offset="872"                 />
-      <DataArray type="Float64" Name="material_state_variable_VolumeRatio_ip" NumberOfTuples="4" format="appended" RangeMin="1.7388253014"         RangeMax="1.7388253014"         offset="944"                 />
-      <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="371354.0957"          RangeMax="371354.0957"          offset="1016"                />
+      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="19" format="appended" RangeMin="45"                   RangeMax="103"                  offset="284"                 />
+      <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.012518851362"       RangeMax="0.012518851362"       offset="364"                 />
+      <DataArray type="Float64" Name="epsilon_m_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.012518851362"       RangeMax="0.012518851362"       offset="472"                 />
+      <DataArray type="Float64" Name="material_state_variable_ElasticStrain_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.0035517417046"      RangeMax="0.0035517417046"      offset="588"                 />
+      <DataArray type="Float64" Name="material_state_variable_EquivalentPlasticStrain_ip" NumberOfTuples="4" format="appended" RangeMin="0.0095861092158"      RangeMax="0.0095861092158"      offset="716"                 />
+      <DataArray type="Float64" Name="material_state_variable_PlasticVolumetricStrain_ip" NumberOfTuples="4" format="appended" RangeMin="-0.0086274863678"     RangeMax="-0.0086274863678"     offset="792"                 />
+      <DataArray type="Float64" Name="material_state_variable_PreConsolidationPressure_ip" NumberOfTuples="4" format="appended" RangeMin="250968.14442"         RangeMax="250968.14442"         offset="864"                 />
+      <DataArray type="Float64" Name="material_state_variable_VolumeRatio_ip" NumberOfTuples="4" format="appended" RangeMin="1.7381041629"         RangeMax="1.7381041629"         offset="940"                 />
+      <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="387298.31369"         RangeMax="387298.31369"         offset="1012"                />
     </FieldData>
     <Piece NumberOfPoints="4"                    NumberOfCells="1"                   >
       <PointData>
-        <DataArray type="Float64" Name="ElasticStrain" NumberOfComponents="4" format="appended" RangeMin="0.0034383732392"      RangeMax="0.0034383732392"      offset="1144"                />
-        <DataArray type="Float64" Name="EquivalentPlasticStrain" format="appended" RangeMin="0.0089577474413"      RangeMax="0.0089577474413"      offset="1288"                />
+        <DataArray type="Float64" Name="ElasticStrain" NumberOfComponents="4" format="appended" RangeMin="0.0035517417046"      RangeMax="0.0035517417046"      offset="1144"                />
+        <DataArray type="Float64" Name="EquivalentPlasticStrain" format="appended" RangeMin="0.0095861092158"      RangeMax="0.0095861092158"      offset="1288"                />
         <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0"                   RangeMax="5.6843418861e-14"     offset="1364"                />
-        <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-9.689962817e-09"     RangeMax="-3.1871269002e-21"    offset="1436"                />
-        <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="282743.33882"         RangeMax="780991.55066"         offset="1524"                />
-        <DataArray type="Float64" Name="PlasticVolumetricStrain" format="appended" RangeMin="-0.0073892070871"     RangeMax="-0.0073892070871"     offset="1652"                />
-        <DataArray type="Float64" Name="PreConsolidationPressure" format="appended" RangeMin="242368.59489"         RangeMax="242368.59489"         offset="1728"                />
-        <DataArray type="Float64" Name="VolumeRatio" format="appended" RangeMin="1.7388253014"         RangeMax="1.7388253014"         offset="1804"                />
-        <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="0.012104037981"       offset="1880"                />
-        <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.012104037981"       RangeMax="0.012104037981"       offset="1956"                />
-        <DataArray type="Float64" Name="epsilon_m" NumberOfComponents="4" format="appended" RangeMin="0.012104037981"       RangeMax="0.012104037981"       offset="2076"                />
-        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="-5000"                RangeMax="-4997.6866918"        offset="2196"                />
-        <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="-5000"                RangeMax="-4997.6866918"        offset="2272"                />
-        <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="371354.0957"          RangeMax="371354.0957"          offset="2348"                />
-        <DataArray type="Float64" Name="temperature" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2480"                />
-        <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2548"                />
-        <DataArray type="Float64" Name="velocity" NumberOfComponents="2" format="appended" RangeMin="3.1622776602e+149"    RangeMax="-nan"                 offset="2616"                />
+        <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-2.4136393841e-09"    RangeMax="6.4085817958e-21"     offset="1436"                />
+        <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="350811.17965"         RangeMax="1019491.2018"         offset="1532"                />
+        <DataArray type="Float64" Name="PlasticVolumetricStrain" format="appended" RangeMin="-0.0086274863678"     RangeMax="-0.0086274863678"     offset="1656"                />
+        <DataArray type="Float64" Name="PreConsolidationPressure" format="appended" RangeMin="250968.14442"         RangeMax="250968.14442"         offset="1728"                />
+        <DataArray type="Float64" Name="VolumeRatio" format="appended" RangeMin="1.7381041629"         RangeMax="1.7381041629"         offset="1804"                />
+        <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="0.012518851362"       offset="1880"                />
+        <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.012518851362"       RangeMax="0.012518851362"       offset="1956"                />
+        <DataArray type="Float64" Name="epsilon_m" NumberOfComponents="4" format="appended" RangeMin="0.012518851362"       RangeMax="0.012518851362"       offset="2080"                />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="50000"                RangeMax="50000.576214"         offset="2204"                />
+        <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="50000"                RangeMax="50000.576214"         offset="2280"                />
+        <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="387298.31369"         RangeMax="387298.31369"         offset="2356"                />
+        <DataArray type="Float64" Name="temperature" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2488"                />
+        <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2556"                />
+        <DataArray type="Float64" Name="velocity" NumberOfComponents="2" format="appended" RangeMin="3.1622776602e+149"    RangeMax="-nan"                 offset="2624"                />
       </PointData>
       <CellData>
-        <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2684"                />
-        <DataArray type="Float64" Name="fluid_density_avg" format="appended" RangeMin="1e+299"               RangeMax="-1e+299"              offset="2744"                />
-        <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="371354.0957"          RangeMax="371354.0957"          offset="2808"                />
-        <DataArray type="UInt8" Name="temperature_active" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2900"                />
-        <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="1e+299"               RangeMax="-1e+299"              offset="2956"                />
+        <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2692"                />
+        <DataArray type="Float64" Name="fluid_density_avg" format="appended" RangeMin="1e+299"               RangeMax="-1e+299"              offset="2752"                />
+        <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="387298.31369"         RangeMax="387298.31369"         offset="2816"                />
+        <DataArray type="UInt8" Name="temperature_active" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2908"                />
+        <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="1e+299"               RangeMax="-1e+299"              offset="2964"                />
       </CellData>
       <Points>
-        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.4142135624"         offset="3020"                />
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.4142135624"         offset="3028"                />
       </Points>
       <Cells>
-        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="3092"                />
-        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="3164"                />
-        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="3224"                />
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="3100"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="3172"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="3232"                />
       </Cells>
     </Piece>
   </UnstructuredGrid>
   <AppendedData encoding="base64">
-   _AQAAAAAAAAAAgAAAAAAAAPkCAAAAAAAAtAAAAAAAAAA=eF61jz0LwkAMhv/LzV0Up67iLgouIkfaxhK4u9RcriCl/93zY3FRUbvlzUuehwyGgmIroMTBdpyTBRE4R1Puh6eSpUEx5bwwATya0kRqPVjqTN4kX6FYPtqafccBg2bAYizeILCL5HLh/4P5CeJBUQicjZon20MOlUO7chCV6q0K0FSKU6IeXMasP5XNvpY9FDt2yaPKxDLBJYfIjprbZc4xJsFpbPefNter14LDeAEk3BwxAQAAAAAAAAAAgAAAAAAAABkAAAAAAAAAIQAAAAAAAAA=eF4z0zPVM9Q1NDHUTU9JTE41NEox0kvJLCqpBABQQQcwAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAMwAAAAAAAAA=eF5jYICA6pc2VQynOvZDuQxrd7/T2cCnvAfGR5d32x49SX2fI1y+hoB+dHmYfgDBtSCFAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAMwAAAAAAAAA=eF5jYICAqpc2VQynOvZDuQxrd7/T2cCnvAfGR5d32x49SX2fI1y+hoB+dHmYfgDA5SCDAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAQQAAAAAAAAA=eF6zcV67y/yW5/5H02LCpy3J2m8D5R9rPSGy5BnTHhsc8i+cY1/PChXaYw/lP4DKw/jXoPpxyTO4QPQDAPP7QqA=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF5L7S9uKohosk+F0uloNADnEA3NAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF4z2tmio+pSt98ISpui0QDTUA0LAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF5j+x1yiGUaryMblGZHowHL+Av/AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5zNV9TZXXpt70rDhoA+jgPoQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPwAAAAAAAAA=eF5LCBbc2NLFenCWwLT42ycEDiZA+ZwWOpvD7vbvTcAhv2W1yNzLL+btTYTyZ0LlYXwVqH5c8neh+gG64jxlAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAASQAAAAAAAAA=eF6zcV67y/yW5/7H02LCpy3J2m8D5W+7lFBbf/rHbico/xFUHsY/AZW3gPKfQOVhfKGgTDHeaaJ77KD8e1B5GN8FKg8AU2FB7g==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF5L7S9uKohoss+E0ilQOg1KAwDnOA3NAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYACDBgYU2sAWxgcAGaUB7g==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAIQAAAAAAAAA=eF5jYAACqXW7QRSD7LbdvMdL2GXPm+zjA9Mu+wB1MQl9AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAPQAAAAAAAAA=eF6b1MZ42O7SVVuxgi9Rdx0FHa2Sr5VmFigclADzFR1trawN/JVm2YpC5A8uWN86NwUoLwmRPwgAiLQcJQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF4z2tmio+pSt98SShtBaRMoDQDTwA0OAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF5j/R1yiGUaryMXlGaB0qxQGgDL+Av9AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF5zNV9TZXXpt707lHaB0k5QGgD6QA+fAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGAAAAAAAAAA=eF5jYCAOVL20qWI41bEfxq+F8gFnRAhYAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOAAAAAAAAAA=eF5jYICAqpc2VQynOvYzwMBcSe0XFXN2w7gN6PIrUOUr0OT/Mk1qcKlz2wPjo5v/DyoPAG6kH5M=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOAAAAAAAAAA=eF5jYICAipc2VQynOvYzwMBcSe0XFXN2w7gN6PIrUOXR9f9lmtTgUue2B8avRpP/B5UHAG3MH5I=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF7j8OU4vb518wEuKM0AAh0IGgDJ4guXAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF7j8OU4vb518wEuKM0AAh0IGgDJ4guXAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAQgAAAAAAAAA=eF6LDxbc2NLFenCWwLT42ycEDsZD+S3NV3K7nrXsTUWTh/EnQ+UT0ORh/E9OD+7LGC3aCzNvOpr5ss4QeQBeUj0FAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEAAAAAAAAAA=eF5jYACBH/UMZNIAW+cLuQ==AQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBH/UAAncBeA==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF5LCBbc2NLFenCmwLT42ycEDiZA+a+SQ6v3Z0/fCwDrYA9AAQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF5jAAAAAQABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBH/UAAncBeA==AQAAAAAAAAAAgAAAAAAAAGAAAAAAAAAAFQAAAAAAAAA=eF5jYMAHPtjjlcaQh/ER4gCW5AS9AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYIAARijNDKWZoDQAAHgABw==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYYAAAAAoAAU=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jBAAACgAK
+   _AQAAAAAAAAAAgAAAAAAAAPkCAAAAAAAAtAAAAAAAAAA=eF61jz0LwkAMhv/LzV0Up67iLgouIkfaxhK4u9RcriCl/93zY3FRUbvlzUuehwyGgmIroMTBdpyTBRE4R1Puh6eSpUEx5bwwATya0kRqPVjqTN4kX6FYPtqafccBg2bAYizeILCL5HLh/4P5CeJBUQicjZon20MOlUO7chCV6q0K0FSKU6IeXMasP5XNvpY9FDt2yaPKxDLBJYfIjprbZc4xJsFpbPefNter14LDeAEk3BwxAQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPVM9Q1NDHQTbc0M09NMU1MBQApRAS5AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAALwAAAAAAAAA=eF5jYICAhRrFJ6sXd+6Hchmkprzc25/81RrGR5f3c2kt+yDBaYNLnlj9AOozIE0=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAANAAAAAAAAAA=eF5jYICARRrFJ6sXd+6Hchmkprzc25/81RrGR5f3c2kt+yDBaQPjLyagH10eph8A66MgUw==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPwAAAAAAAAA=eF7LtC+SX/Daez+vVwHP4bjs/ZlQPr8Ms5vDLEabTBzyGc1pVfyHZG3yoHw+qHwGlK8A1Y9LfgVUPwAQNTREAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF4rYJx/5uHCZvsCKF2ERgMADgMQhQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF4rioqf57m6cX8RlC5EowEAdw/zAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF7L23BM/eAiPsc8KJ2PRgMA/2AO9w==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF7Lsuv1dzv/2z4LBw0A9gQPTQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAQAAAAAAAAAA=eF7LF/j96dsitoOvq3/dkEsVPJgP5bf9aln8Y02vbT4O+doY74pu0ZW2hVD+G6h8AZS/CKofl/xmqH4A0sJJVQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAASwAAAAAAAAA=eF7LsC+SX/Daez+vVwHP4bjs/dlQ/guT3bePiF2yLoLy+aDyWVD+AWOIfCqUzw2Vz4DyWfSvLX+4TckmD8rngcrD1PtA5QH5FzuJAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF4rYJx/5uHCZvsSKJ0PpQuhNAAOKxCFAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFAAAAAAAAAA=eF5jYAABA1sGCGhgQOMDABf6Ads=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJQAAAAAAAAA=eF5jYGBgONByYjeQYmhw3metduS07rLdIvsqnUG0yj4ArXQMkA==AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAPAAAAAAAAAA=eF5r+eUd3eR5wVY59dv2nExRR49JeiU9k9QOqoP5qo7CBh69Sab3bFUh8gcX5dvYdwPl1SDyBwGZWRyBAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF4rioqf57m6cX8RGp0PpQEAdw/yAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF4r2HBM/eAiPsciKJ2DRgMA/7gO9w==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF7Lsuv1dzv/2z4XSmdC6QwoDQD2LA9NAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAFgAAAAAAAAA=eF5jYCAOLNQoPlm9uHM/Oh8Aa7oI1w==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOgAAAAAAAAA=eF5jYICAhRrFJ6sXd+6Hchm+nK/q2L/hqTWMvxhN/hea/AI0+dio6OKOOD4bGH8emnw8VB4AsPohxQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOgAAAAAAAAA=eF5jYICAxRrFJ6sXd+6Hchm+nK/q2L/hqTWMvwxN/hea/AI0+dio6OKOOD4bGH8hmnw8VB4AszohzQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFwAAAAAAAAA=eF5zd4woEMp64eAGpRlAAIkGAKD+CQw=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFwAAAAAAAAA=eF5zd4woEMp64eAGpRlAAIkGAKD+CQw=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAQgAAAAAAAAA=eF7LE/j96dsitoOvqn/dkEsVPJgH5e+Ivdv1pC7OthjKfwuVL4LyN0RD5PNw6JfacLXmUs9623wc8tOg8gDYCkq8AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAEAAAAAAAAAA=eF5jYACBH/UMZNIAW+cLuQ==AQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBH/UAAncBeA==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF4rEPj96dsitoNvqn/dkEsVPFgA5c/znXYhxWGBLQA2hBH5AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF5jAAAAAQABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADQAAAAAAAAA=eF5jYACBH/UAAncBeA==AQAAAAAAAAAAgAAAAAAAAGAAAAAAAAAAFQAAAAAAAAA=eF5jYMAHPtjjlcaQh/ER4gCW5AS9AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYIAARijNDKWZoDQAAHgABw==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYYAAAAAoAAU=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jBAAACgAK
   </AppendedData>
 </VTKFile>
diff --git a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2_t_2000.000000.vtu b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2_t_2000.000000.vtu
index 23a6468f9aa..dff22b3d478 100644
--- a/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2_t_2000.000000.vtu
+++ b/Tests/Data/ThermoHydroMechanics/RestartMCC/mfront_restart_part2_t_2000.000000.vtu
@@ -3,54 +3,54 @@
   <UnstructuredGrid>
     <FieldData>
       <DataArray type="Int8" Name="IntegrationPointMetaData" NumberOfTuples="761" format="appended" RangeMin="34"                   RangeMax="125"                  offset="0"                   />
-      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="25" format="appended" RangeMin="45"                   RangeMax="121"                  offset="284"                 />
-      <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.018594134882"       RangeMax="0.018594134882"       offset="372"                 />
-      <DataArray type="Float64" Name="epsilon_m_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.018594134882"       RangeMax="0.018594134882"       offset="484"                 />
-      <DataArray type="Float64" Name="material_state_variable_ElasticStrain_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.0039855371996"      RangeMax="0.0039855371996"      offset="596"                 />
-      <DataArray type="Float64" Name="material_state_variable_EquivalentPlasticStrain_ip" NumberOfTuples="4" format="appended" RangeMin="0.014960182305"       RangeMax="0.014960182305"       offset="728"                 />
-      <DataArray type="Float64" Name="material_state_variable_PlasticVolumetricStrain_ip" NumberOfTuples="4" format="appended" RangeMin="-0.013036296298"      RangeMax="-0.013036296298"      offset="804"                 />
-      <DataArray type="Float64" Name="material_state_variable_PreConsolidationPressure_ip" NumberOfTuples="4" format="appended" RangeMin="281653.22505"         RangeMax="281653.22505"         offset="880"                 />
-      <DataArray type="Float64" Name="material_state_variable_VolumeRatio_ip" NumberOfTuples="4" format="appended" RangeMin="1.7275766984"         RangeMax="1.7275766984"         offset="956"                 />
-      <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="435386.10946"         RangeMax="435386.10946"         offset="1028"                />
+      <DataArray type="Int8" Name="OGS_VERSION" NumberOfTuples="19" format="appended" RangeMin="45"                   RangeMax="103"                  offset="284"                 />
+      <DataArray type="Float64" Name="epsilon_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.019520876205"       RangeMax="0.019520876205"       offset="364"                 />
+      <DataArray type="Float64" Name="epsilon_m_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.019520876205"       RangeMax="0.019520876205"       offset="472"                 />
+      <DataArray type="Float64" Name="material_state_variable_ElasticStrain_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="0.0041561497023"      RangeMax="0.0041561497023"      offset="580"                 />
+      <DataArray type="Float64" Name="material_state_variable_EquivalentPlasticStrain_ip" NumberOfTuples="4" format="appended" RangeMin="0.016034685565"       RangeMax="0.016034685565"       offset="708"                 />
+      <DataArray type="Float64" Name="material_state_variable_PlasticVolumetricStrain_ip" NumberOfTuples="4" format="appended" RangeMin="-0.01473980181"       RangeMax="-0.01473980181"       offset="780"                 />
+      <DataArray type="Float64" Name="material_state_variable_PreConsolidationPressure_ip" NumberOfTuples="4" format="appended" RangeMin="293894.56446"         RangeMax="293894.56446"         offset="856"                 />
+      <DataArray type="Float64" Name="material_state_variable_VolumeRatio_ip" NumberOfTuples="4" format="appended" RangeMin="1.7259764233"         RangeMax="1.7259764233"         offset="928"                 />
+      <DataArray type="Float64" Name="sigma_ip" NumberOfComponents="4" NumberOfTuples="4" format="appended" RangeMin="455276.2528"          RangeMax="455276.2528"          offset="1004"                />
     </FieldData>
     <Piece NumberOfPoints="4"                    NumberOfCells="1"                   >
       <PointData>
-        <DataArray type="Float64" Name="ElasticStrain" NumberOfComponents="4" format="appended" RangeMin="0.0039855371996"      RangeMax="0.0039855371996"      offset="1156"                />
-        <DataArray type="Float64" Name="EquivalentPlasticStrain" format="appended" RangeMin="0.014960182305"       RangeMax="0.014960182305"       offset="1300"                />
-        <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0"                   RangeMax="5.6843418861e-14"     offset="1376"                />
-        <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="6.776263578e-21"      RangeMax="1.3552527156e-20"     offset="1448"                />
-        <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="282743.33882"         RangeMax="754012.2212"          offset="1520"                />
-        <DataArray type="Float64" Name="PlasticVolumetricStrain" format="appended" RangeMin="-0.013036296298"      RangeMax="-0.013036296298"      offset="1636"                />
-        <DataArray type="Float64" Name="PreConsolidationPressure" format="appended" RangeMin="281653.22505"         RangeMax="281653.22505"         offset="1716"                />
-        <DataArray type="Float64" Name="VolumeRatio" format="appended" RangeMin="1.7275766984"         RangeMax="1.7275766984"         offset="1792"                />
-        <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="0.018594134882"       offset="1868"                />
-        <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.018594134882"       RangeMax="0.018594134882"       offset="1948"                />
-        <DataArray type="Float64" Name="epsilon_m" NumberOfComponents="4" format="appended" RangeMin="0.018594134882"       RangeMax="0.018594134882"       offset="2068"                />
-        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="-50000"               RangeMax="-50000"               offset="2188"                />
-        <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="-50000"               RangeMax="-50000"               offset="2256"                />
-        <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="435386.10946"         RangeMax="435386.10946"         offset="2324"                />
-        <DataArray type="Float64" Name="temperature" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2468"                />
-        <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2536"                />
-        <DataArray type="Float64" Name="velocity" NumberOfComponents="2" format="appended" RangeMin="3.9929874293e-21"     RangeMax="3.9929874293e-21"     offset="2604"                />
+        <DataArray type="Float64" Name="ElasticStrain" NumberOfComponents="4" format="appended" RangeMin="0.0041561497023"      RangeMax="0.0041561497023"      offset="1132"                />
+        <DataArray type="Float64" Name="EquivalentPlasticStrain" format="appended" RangeMin="0.016034685565"       RangeMax="0.016034685565"       offset="1280"                />
+        <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-0"                   RangeMax="5.6843418861e-14"     offset="1356"                />
+        <DataArray type="Float64" Name="MassFlowRate" format="appended" RangeMin="-5.688163163e-10"     RangeMax="4.539980792e-21"      offset="1428"                />
+        <DataArray type="Float64" Name="NodalForces" NumberOfComponents="2" format="appended" RangeMin="350811.17965"         RangeMax="984359.85971"         offset="1520"                />
+        <DataArray type="Float64" Name="PlasticVolumetricStrain" format="appended" RangeMin="-0.01473980181"       RangeMax="-0.01473980181"       offset="1644"                />
+        <DataArray type="Float64" Name="PreConsolidationPressure" format="appended" RangeMin="293894.56446"         RangeMax="293894.56446"         offset="1720"                />
+        <DataArray type="Float64" Name="VolumeRatio" format="appended" RangeMin="1.7259764233"         RangeMax="1.7259764233"         offset="1796"                />
+        <DataArray type="Float64" Name="displacement" NumberOfComponents="2" format="appended" RangeMin="0"                    RangeMax="0.019520876205"       offset="1872"                />
+        <DataArray type="Float64" Name="epsilon" NumberOfComponents="4" format="appended" RangeMin="0.019520876205"       RangeMax="0.019520876205"       offset="1948"                />
+        <DataArray type="Float64" Name="epsilon_m" NumberOfComponents="4" format="appended" RangeMin="0.019520876205"       RangeMax="0.019520876205"       offset="2068"                />
+        <DataArray type="Float64" Name="pressure" format="appended" RangeMin="1000"                 RangeMax="1000.1357949"         offset="2188"                />
+        <DataArray type="Float64" Name="pressure_interpolated" format="appended" RangeMin="1000"                 RangeMax="1000.1357949"         offset="2272"                />
+        <DataArray type="Float64" Name="sigma" NumberOfComponents="4" format="appended" RangeMin="455276.2528"          RangeMax="455276.2528"          offset="2356"                />
+        <DataArray type="Float64" Name="temperature" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2500"                />
+        <DataArray type="Float64" Name="temperature_interpolated" format="appended" RangeMin="293.15"               RangeMax="293.15"               offset="2568"                />
+        <DataArray type="Float64" Name="velocity" NumberOfComponents="2" format="appended" RangeMin="1.3579489269e-10"     RangeMax="1.3579489269e-10"     offset="2636"                />
       </PointData>
       <CellData>
-        <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2688"                />
-        <DataArray type="Float64" Name="fluid_density_avg" format="appended" RangeMin="1000"                 RangeMax="1000"                 offset="2748"                />
-        <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="435386.10946"         RangeMax="435386.10946"         offset="2812"                />
-        <DataArray type="UInt8" Name="temperature_active" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2904"                />
-        <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001"                RangeMax="0.001"                offset="2960"                />
+        <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2748"                />
+        <DataArray type="Float64" Name="fluid_density_avg" format="appended" RangeMin="1000"                 RangeMax="1000"                 offset="2808"                />
+        <DataArray type="Float64" Name="stress_avg" NumberOfComponents="4" format="appended" RangeMin="455276.2528"          RangeMax="455276.2528"          offset="2872"                />
+        <DataArray type="UInt8" Name="temperature_active" format="appended" RangeMin="0"                    RangeMax="0"                    offset="2964"                />
+        <DataArray type="Float64" Name="viscosity_avg" format="appended" RangeMin="0.001"                RangeMax="0.001"                offset="3020"                />
       </CellData>
       <Points>
-        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.4142135624"         offset="3028"                />
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0"                    RangeMax="1.4142135624"         offset="3088"                />
       </Points>
       <Cells>
-        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="3100"                />
-        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="3172"                />
-        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="3232"                />
+        <DataArray type="Int64" Name="connectivity" format="appended" RangeMin=""                     RangeMax=""                     offset="3160"                />
+        <DataArray type="Int64" Name="offsets" format="appended" RangeMin=""                     RangeMax=""                     offset="3232"                />
+        <DataArray type="UInt8" Name="types" format="appended" RangeMin=""                     RangeMax=""                     offset="3292"                />
       </Cells>
     </Piece>
   </UnstructuredGrid>
   <AppendedData encoding="base64">
-   _AQAAAAAAAAAAgAAAAAAAAPkCAAAAAAAAtAAAAAAAAAA=eF61jz0LwkAMhv/LzV0Up67iLgouIkfaxhK4u9RcriCl/93zY3FRUbvlzUuehwyGgmIroMTBdpyTBRE4R1Puh6eSpUEx5bwwATya0kRqPVjqTN4kX6FYPtqafccBg2bAYizeILCL5HLh/4P5CeJBUQicjZon20MOlUO7chCV6q0K0FSKU6IeXMasP5XNvpY9FDt2yaPKxDLBJYfIjprbZc4xJsFpbPefNter14LDeAEk3BwxAQAAAAAAAAAAgAAAAAAAABkAAAAAAAAAIQAAAAAAAAA=eF4z0zPVM9Q1NDHUTU9JTE41NEox0kvJLCqpBABQQQcwAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAMwAAAAAAAAA=eF5jYICAB723XcK5Ju+Hchn+vCsSbE2y2gPjo8tP/rr3/k6FcLj8IwL60eVh+gFZdiO1AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAMwAAAAAAAAA=eF5jYICAh723XcK5Ju+Hchn+vCsSbE2y2gPjo8tP/7r3/k6FcLj8YwL60eVh+gFb1iPBAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAQAAAAAAAAAA=eF57w//d2HS6//7pd/xeMq7L2/8Gyve/fHmR9Oxf1m9wyMdVH2PUWiVv8wnKnw2Vh/EdoPpxyftD9QMAqW5AtQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF57EftubeHiPvsXUPoNGg0AQ98TDQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF57dFhM6Nimrv2PoPQzNBoALD8SQQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF57kl707IqBoOMTKP0UjQYAFG8P5w==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF4rdT1fpr7kt30pDhoABq8QEQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPgAAAAAAAAA=eF47+NlRYE8150E2BiDoED54EMpvymdfH5rbZnsQh7ytAbvznw2rbI9A+RxQ+aNQfilUPy55Q6h+ALYoMp4=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAASQAAAAAAAAA=eF57xf/d2HS6//6pd/xeMq7L2/8Kyv/kVPHBmOPJ7q9Q/jyoPIx/Byr/EsqfApWH8b1sXLf8zlCx+QDlz4LKw/hWUHkA4vZBcA==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF57FvtubeHiPvsPUPoZlH4NpQFD9xMMAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFAAAAAAAAAA=eF5jYAABA1sGCGhgQOMDABf6Ads=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFAAAAAAAAAA=eF5jYACBA9ZgiuEClEbwATdcBA0=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAANAAAAAAAAAA=eF67tv+6UuSkA7ayBV+i7joKOmo83n5BrFDuoAKYr+i4XnRRRLbqRZj8QTVU+YMA3EgeEw==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGQAAAAAAAAA=eF57dFhM6Nimrv0vofRDKP0USgMALI8SQg==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF57kl707IqBoONzKP0ISj+B0gAUfw/mAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF4rdT1fpr7kt305lC6B0kVQGgAGtxAPAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAGQAAAAAAAAA=eF5jYCAO3O+97RLONXk/jP8YygcAbZAIgQ==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOAAAAAAAAAA=eF5jYICAB723XcK5Ju+HchlOHNdrzihi3APjP0GT/4Amfw9NXsSa8cFx3Vi4PLr5olB5AKppH1A=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOAAAAAAAAAA=eF5jYICAB723XcK5Ju+HchlmHNdrzihi3APjP0OTP4Amj65fyprxwXHdWLj8QzR5Sag8AIuxHv8=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5jYACCrBcHGHDQAHK4CEk=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5jYACCrBcHGHDQAHK4CEk=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAASQAAAAAAAAA=eF47+NlRYE8150FWBiDoED54AMrPqfMpO/CoYO8JKJ8LKn8cyk+Ayu+D8lmg8jD+D9+T90JMNtoegfLZoPKHofw3UHkAwgk3Zw==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAHAAAAAAAAAA=eF5bOeX7Acvbm6wZoGAldv5uGH81Gh8A68YU1w==AQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADgAAAAAAAAA=eF5jYAACh34HAAImARA=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAIgAAAAAAAAA=eF47/NlRYE8150F2BiDoED54GMrnbku5rqG2wBYAzrEMaw==AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF5jAAAAAQABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF77s/LjJd+kAHsAGYAEpw==AQAAAAAAAAAAgAAAAAAAAGAAAAAAAAAAFQAAAAAAAAA=eF5jYMAHPtjjlcaQh/ER4gCW5AS9AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYIAARijNDKWZoDQAAHgABw==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYYAAAAAoAAU=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jBAAACgAK
+   _AQAAAAAAAAAAgAAAAAAAAPkCAAAAAAAAtAAAAAAAAAA=eF61jz0LwkAMhv/LzV0Up67iLgouIkfaxhK4u9RcriCl/93zY3FRUbvlzUuehwyGgmIroMTBdpyTBRE4R1Puh6eSpUEx5bwwATya0kRqPVjqTN4kX6FYPtqafccBg2bAYizeILCL5HLh/4P5CeJBUQicjZon20MOlUO7chCV6q0K0FSKU6IeXMasP5XNvpY9FDt2yaPKxDLBJYfIjprbZc4xJsFpbPefNter14LDeAEk3BwxAQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPVM9Q1NDHQTbc0M09NMU1MBQApRAS5AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAALwAAAAAAAAA=eF5jYICAd/MmXHb/O3k/lMvg+6815cJGVhsYH11+1WsV1bd/RHHKE6sfADytJjE=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAALwAAAAAAAAA=eF5jYICAd/MmXHb/O3k/lMvg+6815cJGVhsYH11+1WsV1bd/RHHKE6sfADytJjE=AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPwAAAAAAAAA=eF6LfPh4Bee6gP18db49KR/y94dC+Z5Bf4qbPZj3ROKQX9dxufJKgcIeHyifHyoP45tA9eOSb4bqBwAQr0PnAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5LfpfdkZo9wT4ZBw0ABp8PjQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF4LKjq99Zl+3/4gKB2ARgMALwMSlQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFQAAAAAAAAA=eF5jnMDkJPVB0JERSjOg0QCM3AjDAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF67KmBbMnPub/urUPoaGg0ABk8QGw==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAPwAAAAAAAAA=eF5rsmH62reF62CzeOWu/UkiBxth/JL/RdqOE/Y24ZCv2HS4s2/76r1VUH4jVL4Syi+C6sclbwDVDwCmuz8kAQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAATAAAAAAAAAA=eF6Le/h4Bee6gP18db49KR/y90dC+QH/rjD/0by02xPKF4LKe0H5BgxXwfKxUD4PVD4Cyo/7urls6UaVPW5o8u5Q/t8vEHkAfzpGTA==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF5LfpfdkZo9wT4VSidB6UQoDQAGrw+MAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYACDBgYU2sAWxgcAGaUB7g==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAIwAAAAAAAAA=eF6rujK3Mqb7896aLyCaeR8DEETs2GINohsKtloDAP5nDQk=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAOgAAAAAAAAA=eF67uirCu+PU472qqd+252SKOs6+vVLyu6DqQTUwX9XxIbtt/dqIZ3tVIPIHt12pSEKSPwgAL+Mggg==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF4LLjq99Zl+3/4gNNoXSgMALysSlQ==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAFgAAAAAAAAA=eF5jmsDkJPVB0JERjf7XD6EBlPUJwg==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAGAAAAAAAAAA=eF67ImBbMnPub/ubUPoqlL4CpQEGbxAbAQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAFgAAAAAAAAA=eF5jYCAOvJs34bL738n70fkAiwILCw==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOQAAAAAAAAA=eF5jYICAN/MmXHb/O3k/lMtgcHRJ4m/h39Yw/gc0eTM0+bdo8s07/9VdE5CygfHRzYfJAwBiMSZ8AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAAOQAAAAAAAAA=eF5jYICAN/MmXHb/O3k/lMtgcHRJ4m/h39Yw/gc0eTM0+bdo8s07/9VdE5CygfHRzYfJAwBiMSZ8AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAHAAAAAAAAAA=eF5jYAACh34HBii94sdCaTHHfofNUBoAcLgJLg==AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAHAAAAAAAAAA=eF5jYAACh34HBii94sdCaTHHfofNUBoAcLgJLg==AQAAAAAAAAAAgAAAAAAAAIAAAAAAAAAASwAAAAAAAAA=eF5rsWH62reF62CzeOWu/UkiB1ugfJ9TV+Z5zo/bW4EmXw7lbzgLkW+G8lug8k1QvlOy2u8phRv3lkH5dVD5Eij/QxJEHgDJDz+LAQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEQAAAAAAAAA=eF5LSwOCoCKHNBw0AMi4DAk=AQAAAAAAAAAAgAAAAAAAAEAAAAAAAAAAMgAAAAAAAAA=eF7rEdL5wDbHwVrjI8OK2ysf7S2A8m+fYwTzp3ArGT3tmLFbDSo/Hcq/AZUHAPUxIak=AQAAAAAAAAAAgAAAAAAAAAQAAAAAAAAADAAAAAAAAAA=eF5jYGBgAAAABAABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAADgAAAAAAAAA=eF5jYAACh34HAAImARA=AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAJAAAAAAAAAA=eF6rs2H62reF62CTeOWu/UkiB2uh/AObd+//cXzVXgD+kBFRAQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF5jAAAAAQABAQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAAEAAAAAAAAAA=eF77s/LjJd+kAHsAGYAEpw==AQAAAAAAAAAAgAAAAAAAAGAAAAAAAAAAFQAAAAAAAAA=eF5jYMAHPtjjlcaQh/ER4gCW5AS9AQAAAAAAAAAAgAAAAAAAACAAAAAAAAAAEwAAAAAAAAA=eF5jYIAARijNDKWZoDQAAHgABw==AQAAAAAAAAAAgAAAAAAAAAgAAAAAAAAACwAAAAAAAAA=eF5jYYAAAAAoAAU=AQAAAAAAAAAAgAAAAAAAAAEAAAAAAAAACQAAAAAAAAA=eF7jBAAACgAK
   </AppendedData>
 </VTKFile>
-- 
GitLab