diff --git a/ProcessLib/HeatConduction/Tests.cmake b/ProcessLib/HeatConduction/Tests.cmake
index ecd6d10fe47cca8998b9b52fa3c6c8f2fdb4db07..2f51573f5b56c1d15bea50225aa47871428a0729 100644
--- a/ProcessLib/HeatConduction/Tests.cmake
+++ b/ProcessLib/HeatConduction/Tests.cmake
@@ -4,6 +4,7 @@ if (NOT OGS_USE_MPI)
OgsTest(PROJECTFILE Parabolic/T/2D_freezing_disk/circle_disk.prj RUNTIME 2)
OgsTest(PROJECTFILE Parabolic/T/2D_Robin/square_1e4_robin.prj RUNTIME 1)
OgsTest(PROJECTFILE Parabolic/T/2D_Robin/square_1e4_robin_newton.xml RUNTIME 1)
+ OgsTest(PROJECTFILE Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_IceWaterMix_Scaled.prj RUNTIME 5)
endif()
AddTest(
diff --git a/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_IceWaterMix_Scaled.prj b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_IceWaterMix_Scaled.prj
new file mode 100644
index 0000000000000000000000000000000000000000..6a15a9da558e14315056c5c1e7afbbf25289daef
--- /dev/null
+++ b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_IceWaterMix_Scaled.prj
@@ -0,0 +1,281 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<!--units: kg, m, s, K and derived units: N, Pa=N/m^2 -->
+<OpenGeoSysProject>
+ <meshes>
+ <!-- meshes with labels 0,1,2,3 are bottom, right, top and left, resp. -->
+ <mesh>ManSol3_mesh_unit_square_QU_10x10.vtu</mesh>
+ <mesh>ManSol3_mesh_unit_square_QU_10x10_0.vtu</mesh>
+ <mesh>ManSol3_mesh_unit_square_QU_10x10_1.vtu</mesh>
+ <mesh>ManSol3_mesh_unit_square_QU_10x10_2.vtu</mesh>
+ <mesh>ManSol3_mesh_unit_square_QU_10x10_3.vtu</mesh>
+ </meshes>
+ <processes>
+ <process>
+ <name>HeatConduction</name>
+ <type>HEAT_CONDUCTION</type>
+ <integration_order>2</integration_order>
+ <process_variables>
+ <process_variable>temperature</process_variable>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="heat_flux" output_name="heat_flux"/>
+ </secondary_variables>
+ </process>
+ </processes>
+ <media>
+ <medium id="4">
+ <phases>
+ <phase>
+ <!-- Pore water phase -->
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>0.58</value> <!-- W/m/K-->
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>4.190</value> <!-- J/kg/K, this is 10^(-3) rescaled -->
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>1.000</value> <!--kg/m^3, this is 10^(-3) rescaled -->
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <!-- Pore ice phase -->
+ <type>FrozenLiquid</type>
+ <properties>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>2.2</value> <!-- W/m/K-->
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>2.090</value> <!-- J/kg/K, this is 10^(-3) rescaled -->
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>0.920</value> <!--kg/m^3, this is 10^(-3) rescaled -->
+ </property>
+ </properties>
+ </phase>
+ <phase>
+ <!-- Porous solid phase -->
+ <type>Solid</type>
+ <properties>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>1.1</value> <!-- W/m/K-->
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>Constant</type>
+ <value>0.900</value> <!-- J/kg/K, this is 10^(-3) rescaled -->
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>2.000</value> <!--kg/m^3, this is 10^(-3) rescaled -->
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.5</value>
+ </property>
+ <property>
+ <name>volume_fraction</name>
+ <type>TemperatureDependentFraction</type>
+ <steepness>2</steepness> <!-- k in the Sigmoid function -->
+ <characteristic_temperature>273.15</characteristic_temperature> <!-- melting temperature T_m -->
+ </property>
+ <property>
+ <name>density</name>
+ <type>VolumeFractionAverage</type>
+ </property>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>VolumeFractionAverage</type>
+ </property>
+ <property>
+ <name>specific_heat_capacity</name>
+ <type>SpecificHeatCapacityWithLatentHeat</type>
+ <specific_latent_heat>334.000</specific_latent_heat> <!-- J/kg, this is 10^(-3) rescaled -->
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <parameters>
+ <parameter>
+ <name>HeatSource</name> <!-- this where we insert Q_T(x,y,t) -->
+ <type>Function</type>
+ <expression>
+ (0.5*2.000*0.900
+ +0.5*exp(2*(17*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))-7))/(1+exp(2*(17*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))-7)))*1.000*4.190
+ +0.5/(1+exp(2*(17*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))-7)))*0.920*2.090
+ +3.34e2*0.920*0.5*2*exp(2*(17*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))-7))/(1+exp(2*(17*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))-7)))^2 )
+ *0.5*pi*17*y*((1-x)*cos(0.5*pi*t)-x*sin(0.5*pi*t))
+
+ +(2.2-0.58)*0.5*2*exp(2*(17*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))-7))/(1+exp(2*(17*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))-7)))^2
+ *17^2*(y^2*(1-sin(pi*t))+x*(1-x)*sin(pi*t)+(x-0.5)*cos(pi*t)+x^2-x+0.5)
+ </expression>
+ <!-- Here is the expression for Q_T in general terms:
+ ( (1-phi)*rho_SR*c_pS
+ +phi*exp(k*(b*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))+c-Tm))/(1+exp(k*(b*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))+c-Tm)))*rho_FR*c_pF
+ +phi/(1+exp(k*(b*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))+c-Tm)))*rho_IR*c_pI
+ +ell*rho_IR*phi*k*exp(k*(b*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))+c-Tm))/(1+exp(k*(b*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))+c-Tm)))^2 )
+ *0.5*pi*b*y*((1-x)*cos(0.5*pi*t)-x*sin(0.5*pi*t))
+
+ +(lambda_IR-lambda_FR)*phi*k*exp(k*(b*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))+c-Tm))/(1+exp(k*(b*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))+c-Tm)))^2
+ *b^2*(y^2*(1-sin(pi*t))+x*(1-x)*sin(pi*t)+(x-0.5)*cos(pi*t)+x^2-x+0.5)
+ -->
+ </parameter>
+ <parameter>
+ <name>T0</name>
+ <type>Function</type>
+ <expression> 17*x*y+266.15 </expression> <!-- the IC function T_0 -->
+ </parameter>
+ <parameter>
+ <name>T1b</name>
+ <type>Function</type>
+ <expression> 266.15 </expression> <!-- this is T1 on bottom side -->
+ </parameter>
+ <parameter>
+ <name>T1r</name>
+ <type>Function</type>
+ <expression> 17*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))+266.15 </expression> <!-- this is T1 on right side -->
+ </parameter>
+ <parameter>
+ <name>T1t</name>
+ <type>Function</type>
+ <expression> 17*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))+266.15 </expression> <!-- this is T1 on top side -->
+ </parameter>
+ <parameter>
+ <name>T1l</name>
+ <type>Function</type>
+ <expression> 17*y*(x*cos(0.5*pi*t)+(1-x)*sin(0.5*pi*t))+266.15 </expression> <!-- this is T1 on left side -->
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>temperature</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>T0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <mesh>ManSol3_mesh_unit_square_QU_10x10_0</mesh>
+ <type>Dirichlet</type>
+ <parameter>T1b</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>ManSol3_mesh_unit_square_QU_10x10_1</mesh>
+ <type>Dirichlet</type>
+ <parameter>T1r</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>ManSol3_mesh_unit_square_QU_10x10_2</mesh>
+ <type>Dirichlet</type>
+ <parameter>T1t</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>ManSol3_mesh_unit_square_QU_10x10_3</mesh>
+ <type>Dirichlet</type>
+ <parameter>T1l</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ <source_terms>
+ <source_term>
+ <mesh>ManSol3_mesh_unit_square_QU_10x10</mesh>
+ <type>Volumetric</type>
+ <parameter>HeatSource</parameter>
+ </source_term>
+ </source_terms>
+ </process_variable>
+ </process_variables>
+ <time_loop>
+ <processes>
+ <process ref="HeatConduction">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>2e-12</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial>0.0</t_initial>
+ <t_end>1.0</t_end>
+ <timesteps>
+ <pair>
+ <repeat>200</repeat>
+ <delta_t>0.005</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>ManSol_IceWaterMix_T-plot</prefix>
+ <timesteps>
+ <pair>
+ <repeat>10000</repeat>
+ <each_steps>10</each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable> temperature </variable>
+ <variable> heat_flux </variable>
+ </variables>
+ <suffix>_ts_{:timestep}_t_{:gtime}</suffix>
+ </output>
+ </time_loop>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_picard</name>
+ <type>Picard</type>
+ <max_iter>200</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <eigen>
+ <solver_type>SparseLU</solver_type>
+ <scaling>true</scaling>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <!--primary field-->
+ <vtkdiff>
+ <file>ManSol_IceWaterMix_T-plot_ts_50_t_0.25.vtu</file>
+ <field>temperature</field>
+ <absolute_tolerance>2e-6</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ <vtkdiff>
+ <file>ManSol_IceWaterMix_T-plot_ts_100_t_0.5.vtu</file>
+ <field>temperature</field>
+ <absolute_tolerance>5e-6</absolute_tolerance>
+ <relative_tolerance>0</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10.vtu b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..301092706e95424c83117b20942f49576012963c
--- /dev/null
+++ b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10.vtu
@@ -0,0 +1,25 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="121" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="120" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="4" RangeMax="4" offset="1304" />
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="99" offset="1848" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="2928" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="6812" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="11092" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="12172" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _yAMAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAEAAAAAAAAAAQQAAAAAAAABCAAAAAAAAAEMAAAAAAAAARAAAAAAAAABFAAAAAAAAAEYAAAAAAAAARwAAAAAAAABIAAAAAAAAAEkAAAAAAAAASgAAAAAAAABLAAAAAAAAAEwAAAAAAAAATQAAAAAAAABOAAAAAAAAAE8AAAAAAAAAUAAAAAAAAABRAAAAAAAAAFIAAAAAAAAAUwAAAAAAAABUAAAAAAAAAFUAAAAAAAAAVgAAAAAAAABXAAAAAAAAAFgAAAAAAAAAWQAAAAAAAABaAAAAAAAAAFsAAAAAAAAAXAAAAAAAAABdAAAAAAAAAF4AAAAAAAAAXwAAAAAAAABgAAAAAAAAAGEAAAAAAAAAYgAAAAAAAABjAAAAAAAAAGQAAAAAAAAAZQAAAAAAAABmAAAAAAAAAGcAAAAAAAAAaAAAAAAAAABpAAAAAAAAAGoAAAAAAAAAawAAAAAAAABsAAAAAAAAAG0AAAAAAAAAbgAAAAAAAABvAAAAAAAAAHAAAAAAAAAAcQAAAAAAAAByAAAAAAAAAHMAAAAAAAAAdAAAAAAAAAB1AAAAAAAAAHYAAAAAAAAAdwAAAAAAAAB4AAAAAAAAAA==kAEAAAAAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAABAAAAAQAAAAEAAAAIAMAAAAAAAAAAAAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAEAAAAAAAAAAQQAAAAAAAABCAAAAAAAAAEMAAAAAAAAARAAAAAAAAABFAAAAAAAAAEYAAAAAAAAARwAAAAAAAABIAAAAAAAAAEkAAAAAAAAASgAAAAAAAABLAAAAAAAAAEwAAAAAAAAATQAAAAAAAABOAAAAAAAAAE8AAAAAAAAAUAAAAAAAAABRAAAAAAAAAFIAAAAAAAAAUwAAAAAAAABUAAAAAAAAAFUAAAAAAAAAVgAAAAAAAABXAAAAAAAAAFgAAAAAAAAAWQAAAAAAAABaAAAAAAAAAFsAAAAAAAAAXAAAAAAAAABdAAAAAAAAAF4AAAAAAAAAXwAAAAAAAABgAAAAAAAAAGEAAAAAAAAAYgAAAAAAAABjAAAAAAAAAA==WAsAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADwPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADwPwAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAABvZZmZmZm5PwAAAAAAAAAAAAAAAAAAAABhW5mZmZnJPwAAAAAAAAAAAAAAAAAAAACU/zIzMzPTPwAAAAAAAAAAAAAAAAAAAAB3UZmZmZnZPwAAAAAAAAAAAAAAAAAAAAAfpP/////fPwAAAAAAAAAAAAAAAAAAAAAlDjMzMzPjPwAAAAAAAAAAAAAAAAAAAACbSmZmZmbmPwAAAAAAAAAAAAAAAAAAAAASh5mZmZnpPwAAAAAAAAAAAAAAAAAAAACKw8zMzMzsPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADwP29lmZmZmbk/AAAAAAAAAAAAAAAAAADwP2FbmZmZmck/AAAAAAAAAAAAAAAAAADwP5T/MjMzM9M/AAAAAAAAAAAAAAAAAADwP3dRmZmZmdk/AAAAAAAAAAAAAAAAAADwPx+k/////98/AAAAAAAAAAAAAAAAAADwPyUOMzMzM+M/AAAAAAAAAAAAAAAAAADwP5tKZmZmZuY/AAAAAAAAAAAAAAAAAADwPxKHmZmZmek/AAAAAAAAAAAAAAAAAADwP4rDzMzMzOw/AAAAAAAAAAAmvszMzMzsPwAAAAAAAPA/AAAAAAAAAACYmZmZmZnpPwAAAAAAAPA/AAAAAAAAAADMfmZmZmbmPwAAAAAAAPA/AAAAAAAAAAD/YzMzMzPjPwAAAAAAAPA/AAAAAAAAAABvSAAAAADgPwAAAAAAAPA/AAAAAAAAAAC5DpqZmZnZPwAAAAAAAPA/AAAAAAAAAAALizMzMzPTPwAAAAAAAPA/AAAAAAAAAAC5DpqZmZnJPwAAAAAAAPA/AAAAAAAAAACyDpqZmZm5PwAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAACa+zMzMzOw/AAAAAAAAAAAAAAAAAAAAAJiZmZmZmek/AAAAAAAAAAAAAAAAAAAAAMx+ZmZmZuY/AAAAAAAAAAAAAAAAAAAAAP9jMzMzM+M/AAAAAAAAAAAAAAAAAAAAAG9IAAAAAOA/AAAAAAAAAAAAAAAAAAAAALkOmpmZmdk/AAAAAAAAAAAAAAAAAAAAAAuLMzMzM9M/AAAAAAAAAAAAAAAAAAAAALkOmpmZmck/AAAAAAAAAAAAAAAAAAAAALIOmpmZmbk/AAAAAAAAAAD1a5mZmZm5Pw7zmZmZmbk/AAAAAAAAAABqeJmZmZm5P8nxmZmZmck/AAAAAAAAAAArhpmZmZm5P1tzMzMzM9M/AAAAAAAAAAAflJmZmZm5P/TpmZmZmdk/AAAAAAAAAACGopmZmZm5P20rAAAAAOA/AAAAAAAAAABJspmZmZm5P+RUMzMzM+M/AAAAAAAAAABxxJmZmZm5P0h6ZmZmZuY/AAAAAAAAAADy2ZmZmZm5PwOgmZmZmek/AAAAAAAAAAAO85mZmZm5P0DLzMzMzOw/AAAAAAAAAAD9aJmZmZnJP/LZmZmZmbk/AAAAAAAAAADBd5mZmZnJPyjYmZmZmck/AAAAAAAAAAAUhpmZmZnJP4hfMzMzM9M/AAAAAAAAAADjk5mZmZnJP97OmZmZmdk/AAAAAAAAAADcoZmZmZnJP2obAAAAAOA/AAAAAAAAAAATsZmZmZnJP9ZJMzMzM+M/AAAAAAAAAADIwpmZmZnJP2x1ZmZmZuY/AAAAAAAAAAAl2JmZmZnJP+2gmZmZmek/AAAAAAAAAADG8ZmZmZnJP9XOzMzMzOw/AAAAAAAAAABxDTMzMzPTP3HEmZmZmbk/AAAAAAAAAADgGTMzMzPTP8jCmZmZmck/AAAAAAAAAADMJDMzMzPTP9dPMzMzM9M/AAAAAAAAAADELjMzMzPTPyy7mZmZmdk/AAAAAAAAAACbODMzMzPTP/UQAAAAAOA/AAAAAAAAAABEQzMzMzPTP6BBMzMzM+M/AAAAAAAAAADVTzMzMzPTP6dwZmZmZuY/AAAAAAAAAACIXzMzMzPTP3CfmZmZmek/AAAAAAAAAABbczMzMzPTPyTPzMzMzOw/AAAAAAAAAADJZ5mZmZnZP0mymZmZmbk/AAAAAAAAAABjeZmZmZnZPxOxmZmZmck/AAAAAAAAAACAh5mZmZnZP0RDMzMzM9M/AAAAAAAAAAC5k5mZmZnZPw+smZmZmdk/AAAAAAAAAABzn5mZmZnZPy4JAAAAAOA/AAAAAAAAAAARrJmZmZnZPw87MzMzM+M/AAAAAAAAAAAqu5mZmZnZPx5sZmZmZuY/AAAAAAAAAADezpmZmZnZPwedmZmZmek/AAAAAAAAAADy6ZmZmZnZP03OzMzMzOw/AAAAAAAAAADexv/////fP4aimZmZmbk/AAAAAAAAAABh3P/////fP9mhmZmZmck/AAAAAAAAAAAd7P/////fP5s4MzMzM9M/AAAAAAAAAAAt+f/////fP3OfmZmZmdk/AAAAAAAAAACzAgAAAADgP7MCAAAAAOA/AAAAAAAAAAAwCQAAAADgP041MzMzM+M/AAAAAAAAAAD0EAAAAADgP75nZmZmZuY/AAAAAAAAAABqGwAAAADgP0OamZmZmek/AAAAAAAAAABuKwAAAADgPwfNzMzMzOw/AAAAAAAAAACZGTMzMzPjPyGUmZmZmbk/AAAAAAAAAACVIjMzMzPjP+OTmZmZmck/AAAAAAAAAACzKTMzMzPjP8QuMzMzM9M/AAAAAAAAAAC5LzMzMzPjP7mTmZmZmdk/AAAAAAAAAABPNTMzMzPjPyv5/////98/AAAAAAAAAAAQOzMzMzPjP7cvMzMzM+M/AAAAAAAAAACgQTMzMzPjP0pjZmZmZuY/AAAAAAAAAADYSTMzMzPjP0GXmZmZmek/AAAAAAAAAADjVDMzMzPjP4jLzMzMzOw/AAAAAAAAAAA4UmZmZmbmPyyGmZmZmbk/AAAAAAAAAADYWGZmZmbmPxSGmZmZmck/AAAAAAAAAAByXmZmZmbmP80kMzMzM9M/AAAAAAAAAABKY2ZmZmbmP4SHmZmZmdk/AAAAAAAAAAC/Z2ZmZmbmPx3s/////98/AAAAAAAAAAAfbGZmZmbmP7MpMzMzM+M/AAAAAAAAAACncGZmZmbmP29eZmZmZuY/AAAAAAAAAABudWZmZmbmP+uTmZmZmek/AAAAAAAAAABJemZmZmbmP9zJzMzMzOw/AAAAAAAAAADPi5mZmZnpP2l4mZmZmbk/AAAAAAAAAAAjkJmZmZnpP8R3mZmZmck/AAAAAAAAAADsk5mZmZnpP+IZMzMzM9M/AAAAAAAAAABAl5mZmZnpP2N5mZmZmdk/AAAAAAAAAABDmpmZmZnpP2Pc/////98/AAAAAAAAAAAInZmZmZnpP5UiMzMzM+M/AAAAAAAAAABwn5mZmZnpP9dYZmZmZuY/AAAAAAAAAADsoJmZmZnpPyCQmZmZmek/AAAAAAAAAAAEoJmZmZnpP/PHzMzMzOw/AAAAAAAAAADSxczMzMzsP/ZrmZmZmbk/AAAAAAAAAAD0x8zMzMzsP/1omZmZmck/AAAAAAAAAADcyczMzMzsP3ENMzMzM9M/AAAAAAAAAACHy8zMzMzsP8lnmZmZmdk/AAAAAAAAAAAHzczMzMzsP97G/////98/AAAAAAAAAABOzszMzMzsP5kZMzMzM+M/AAAAAAAAAAAmz8zMzMzsPzdSZmZmZuY/AAAAAAAAAADUzszMzMzsP8+LmZmZmek/AAAAAAAAAAA+y8zMzMzsP9LFzMzMzOw/AAAAAAAAAAA=gAwAAAAAAAAAAAAAAAAAAAQAAAAAAAAAKAAAAAAAAAAnAAAAAAAAACcAAAAAAAAAKAAAAAAAAAApAAAAAAAAACYAAAAAAAAAJgAAAAAAAAApAAAAAAAAACoAAAAAAAAAJQAAAAAAAAAlAAAAAAAAACoAAAAAAAAAKwAAAAAAAAAkAAAAAAAAACQAAAAAAAAAKwAAAAAAAAAsAAAAAAAAACMAAAAAAAAAIwAAAAAAAAAsAAAAAAAAAC0AAAAAAAAAIgAAAAAAAAAiAAAAAAAAAC0AAAAAAAAALgAAAAAAAAAhAAAAAAAAACEAAAAAAAAALgAAAAAAAAAvAAAAAAAAACAAAAAAAAAAIAAAAAAAAAAvAAAAAAAAADAAAAAAAAAAHwAAAAAAAAAfAAAAAAAAADAAAAAAAAAAHgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAxAAAAAAAAACgAAAAAAAAAKAAAAAAAAAAxAAAAAAAAADIAAAAAAAAAKQAAAAAAAAApAAAAAAAAADIAAAAAAAAAMwAAAAAAAAAqAAAAAAAAACoAAAAAAAAAMwAAAAAAAAA0AAAAAAAAACsAAAAAAAAAKwAAAAAAAAA0AAAAAAAAADUAAAAAAAAALAAAAAAAAAAsAAAAAAAAADUAAAAAAAAANgAAAAAAAAAtAAAAAAAAAC0AAAAAAAAANgAAAAAAAAA3AAAAAAAAAC4AAAAAAAAALgAAAAAAAAA3AAAAAAAAADgAAAAAAAAALwAAAAAAAAAvAAAAAAAAADgAAAAAAAAAOQAAAAAAAAAwAAAAAAAAADAAAAAAAAAAOQAAAAAAAAAdAAAAAAAAAB4AAAAAAAAABQAAAAAAAAAGAAAAAAAAADoAAAAAAAAAMQAAAAAAAAAxAAAAAAAAADoAAAAAAAAAOwAAAAAAAAAyAAAAAAAAADIAAAAAAAAAOwAAAAAAAAA8AAAAAAAAADMAAAAAAAAAMwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAANAAAAAAAAAA0AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA1AAAAAAAAADUAAAAAAAAAPgAAAAAAAAA/AAAAAAAAADYAAAAAAAAANgAAAAAAAAA/AAAAAAAAAEAAAAAAAAAANwAAAAAAAAA3AAAAAAAAAEAAAAAAAAAAQQAAAAAAAAA4AAAAAAAAADgAAAAAAAAAQQAAAAAAAABCAAAAAAAAADkAAAAAAAAAOQAAAAAAAABCAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAAQwAAAAAAAAA6AAAAAAAAADoAAAAAAAAAQwAAAAAAAABEAAAAAAAAADsAAAAAAAAAOwAAAAAAAABEAAAAAAAAAEUAAAAAAAAAPAAAAAAAAAA8AAAAAAAAAEUAAAAAAAAARgAAAAAAAAA9AAAAAAAAAD0AAAAAAAAARgAAAAAAAABHAAAAAAAAAD4AAAAAAAAAPgAAAAAAAABHAAAAAAAAAEgAAAAAAAAAPwAAAAAAAAA/AAAAAAAAAEgAAAAAAAAASQAAAAAAAABAAAAAAAAAAEAAAAAAAAAASQAAAAAAAABKAAAAAAAAAEEAAAAAAAAAQQAAAAAAAABKAAAAAAAAAEsAAAAAAAAAQgAAAAAAAABCAAAAAAAAAEsAAAAAAAAAGwAAAAAAAAAcAAAAAAAAAAcAAAAAAAAACAAAAAAAAABMAAAAAAAAAEMAAAAAAAAAQwAAAAAAAABMAAAAAAAAAE0AAAAAAAAARAAAAAAAAABEAAAAAAAAAE0AAAAAAAAATgAAAAAAAABFAAAAAAAAAEUAAAAAAAAATgAAAAAAAABPAAAAAAAAAEYAAAAAAAAARgAAAAAAAABPAAAAAAAAAFAAAAAAAAAARwAAAAAAAABHAAAAAAAAAFAAAAAAAAAAUQAAAAAAAABIAAAAAAAAAEgAAAAAAAAAUQAAAAAAAABSAAAAAAAAAEkAAAAAAAAASQAAAAAAAABSAAAAAAAAAFMAAAAAAAAASgAAAAAAAABKAAAAAAAAAFMAAAAAAAAAVAAAAAAAAABLAAAAAAAAAEsAAAAAAAAAVAAAAAAAAAAaAAAAAAAAABsAAAAAAAAACAAAAAAAAAAJAAAAAAAAAFUAAAAAAAAATAAAAAAAAABMAAAAAAAAAFUAAAAAAAAAVgAAAAAAAABNAAAAAAAAAE0AAAAAAAAAVgAAAAAAAABXAAAAAAAAAE4AAAAAAAAATgAAAAAAAABXAAAAAAAAAFgAAAAAAAAATwAAAAAAAABPAAAAAAAAAFgAAAAAAAAAWQAAAAAAAABQAAAAAAAAAFAAAAAAAAAAWQAAAAAAAABaAAAAAAAAAFEAAAAAAAAAUQAAAAAAAABaAAAAAAAAAFsAAAAAAAAAUgAAAAAAAABSAAAAAAAAAFsAAAAAAAAAXAAAAAAAAABTAAAAAAAAAFMAAAAAAAAAXAAAAAAAAABdAAAAAAAAAFQAAAAAAAAAVAAAAAAAAABdAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAJAAAAAAAAAAoAAAAAAAAAXgAAAAAAAABVAAAAAAAAAFUAAAAAAAAAXgAAAAAAAABfAAAAAAAAAFYAAAAAAAAAVgAAAAAAAABfAAAAAAAAAGAAAAAAAAAAVwAAAAAAAABXAAAAAAAAAGAAAAAAAAAAYQAAAAAAAABYAAAAAAAAAFgAAAAAAAAAYQAAAAAAAABiAAAAAAAAAFkAAAAAAAAAWQAAAAAAAABiAAAAAAAAAGMAAAAAAAAAWgAAAAAAAABaAAAAAAAAAGMAAAAAAAAAZAAAAAAAAABbAAAAAAAAAFsAAAAAAAAAZAAAAAAAAABlAAAAAAAAAFwAAAAAAAAAXAAAAAAAAABlAAAAAAAAAGYAAAAAAAAAXQAAAAAAAABdAAAAAAAAAGYAAAAAAAAAGAAAAAAAAAAZAAAAAAAAAAoAAAAAAAAACwAAAAAAAABnAAAAAAAAAF4AAAAAAAAAXgAAAAAAAABnAAAAAAAAAGgAAAAAAAAAXwAAAAAAAABfAAAAAAAAAGgAAAAAAAAAaQAAAAAAAABgAAAAAAAAAGAAAAAAAAAAaQAAAAAAAABqAAAAAAAAAGEAAAAAAAAAYQAAAAAAAABqAAAAAAAAAGsAAAAAAAAAYgAAAAAAAABiAAAAAAAAAGsAAAAAAAAAbAAAAAAAAABjAAAAAAAAAGMAAAAAAAAAbAAAAAAAAABtAAAAAAAAAGQAAAAAAAAAZAAAAAAAAABtAAAAAAAAAG4AAAAAAAAAZQAAAAAAAABlAAAAAAAAAG4AAAAAAAAAbwAAAAAAAABmAAAAAAAAAGYAAAAAAAAAbwAAAAAAAAAXAAAAAAAAABgAAAAAAAAACwAAAAAAAAAMAAAAAAAAAHAAAAAAAAAAZwAAAAAAAABnAAAAAAAAAHAAAAAAAAAAcQAAAAAAAABoAAAAAAAAAGgAAAAAAAAAcQAAAAAAAAByAAAAAAAAAGkAAAAAAAAAaQAAAAAAAAByAAAAAAAAAHMAAAAAAAAAagAAAAAAAABqAAAAAAAAAHMAAAAAAAAAdAAAAAAAAABrAAAAAAAAAGsAAAAAAAAAdAAAAAAAAAB1AAAAAAAAAGwAAAAAAAAAbAAAAAAAAAB1AAAAAAAAAHYAAAAAAAAAbQAAAAAAAABtAAAAAAAAAHYAAAAAAAAAdwAAAAAAAABuAAAAAAAAAG4AAAAAAAAAdwAAAAAAAAB4AAAAAAAAAG8AAAAAAAAAbwAAAAAAAAB4AAAAAAAAABYAAAAAAAAAFwAAAAAAAAAMAAAAAAAAAAEAAAAAAAAADQAAAAAAAABwAAAAAAAAAHAAAAAAAAAADQAAAAAAAAAOAAAAAAAAAHEAAAAAAAAAcQAAAAAAAAAOAAAAAAAAAA8AAAAAAAAAcgAAAAAAAAByAAAAAAAAAA8AAAAAAAAAEAAAAAAAAABzAAAAAAAAAHMAAAAAAAAAEAAAAAAAAAARAAAAAAAAAHQAAAAAAAAAdAAAAAAAAAARAAAAAAAAABIAAAAAAAAAdQAAAAAAAAB1AAAAAAAAABIAAAAAAAAAEwAAAAAAAAB2AAAAAAAAAHYAAAAAAAAAEwAAAAAAAAAUAAAAAAAAAHcAAAAAAAAAdwAAAAAAAAAUAAAAAAAAABUAAAAAAAAAeAAAAAAAAAB4AAAAAAAAABUAAAAAAAAAAgAAAAAAAAAWAAAAAAAAAA==IAMAAAAAAAAEAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAQAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAcAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAoAAAAAAAAACwAAAAAAAAAMAAAAAAAAAA0AAAAAAAAADgAAAAAAAAAPAAAAAAAAABAAAAAAAAAAEQAAAAAAAAASAAAAAAAAABMAAAAAAAAAFAAAAAAAAAAVAAAAAAAAABYAAAAAAAAAFwAAAAAAAAAYAAAAAAAAABkAAAAAAAAAGgAAAAAAAAAbAAAAAAAAABwAAAAAAAAAHQAAAAAAAAAeAAAAAAAAAB8AAAAAAAAAIAAAAAAAAAAhAAAAAAAAACIAAAAAAAAAIwAAAAAAAAAkAAAAAAAAACUAAAAAAAAAJgAAAAAAAAAnAAAAAAAAACgAAAAAAAAAKQAAAAAAAAAqAAAAAAAAACsAAAAAAAAALAAAAAAAAAAtAAAAAAAAAC4AAAAAAAAALwAAAAAAAAAwAAAAAAAAADEAAAAAAAAAMgAAAAAAAAAzAAAAAAAAADQAAAAAAAAANQAAAAAAAAA2AAAAAAAAADcAAAAAAAAAOAAAAAAAAAA5AAAAAAAAADoAAAAAAAAAOwAAAAAAAAA8AAAAAAAAAD0AAAAAAAAAPgAAAAAAAAA/AAAAAAAAAAAAQAAAAAAAAQBAAAAAAAACAEAAAAAAAAMAQAAAAAAABABAAAAAAAAFAEAAAAAAAAYAQAAAAAAABwBAAAAAAAAIAEAAAAAAAAkAQAAAAAAACgBAAAAAAAALAEAAAAAAAAwAQAAAAAAADQBAAAAAAAAOAEAAAAAAAA8AQAAAAAAAEABAAAAAAAARAEAAAAAAABIAQAAAAAAAEwBAAAAAAAAUAEAAAAAAABUAQAAAAAAAFgBAAAAAAAAXAEAAAAAAABgAQAAAAAAAGQBAAAAAAAAaAEAAAAAAABsAQAAAAAAAHABAAAAAAAAdAEAAAAAAAB4AQAAAAAAAHwBAAAAAAAAgAEAAAAAAACEAQAAAAAAAIgBAAAAAAAAjAEAAAAAAACQAQAAAAAAAA==ZAAAAAAAAAAJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJCQkJ
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_0.vtu b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_0.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..f3253d9e549809ae30efcc554132456364de468d
--- /dev/null
+++ b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_0.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="11" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="12" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="90" offset="128" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="248" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="612" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="836" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="956" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _WAAAAAAAAAAAAAAAAAAAAAEAAAAAAAAABAAAAAAAAAAFAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAIAAAAAAAAAAkAAAAAAAAACgAAAAAAAAALAAAAAAAAAAwAAAAAAAAAUAAAAAAAAAAAAAAAAAAAAAoAAAAAAAAAFAAAAAAAAAAeAAAAAAAAACgAAAAAAAAAMgAAAAAAAAA8AAAAAAAAAEYAAAAAAAAAUAAAAAAAAABaAAAAAAAAAA==CAEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAADwPwAAAAAAAAAAAAAAAAAAAABvZZmZmZm5PwAAAAAAAAAAAAAAAAAAAABhW5mZmZnJPwAAAAAAAAAAAAAAAAAAAACU/zIzMzPTPwAAAAAAAAAAAAAAAAAAAAB3UZmZmZnZPwAAAAAAAAAAAAAAAAAAAAAfpP/////fPwAAAAAAAAAAAAAAAAAAAAAlDjMzMzPjPwAAAAAAAAAAAAAAAAAAAACbSmZmZmbmPwAAAAAAAAAAAAAAAAAAAAASh5mZmZnpPwAAAAAAAAAAAAAAAAAAAACKw8zMzMzsPwAAAAAAAAAAAAAAAAAAAAA=oAAAAAAAAAAAAAAAAAAAAAIAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAHAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAkAAAAAAAAACgAAAAAAAAAKAAAAAAAAAAEAAAAAAAAAUAAAAAAAAAACAAAAAAAAAAQAAAAAAAAABgAAAAAAAAAIAAAAAAAAAAoAAAAAAAAADAAAAAAAAAAOAAAAAAAAABAAAAAAAAAAEgAAAAAAAAAUAAAAAAAAAA==CgAAAAAAAAADAwMDAwMDAwMD
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_1.vtu b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_1.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..1ab139b0db28849453822458b63377d723cd1fc5
--- /dev/null
+++ b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_1.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="11" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="1" RangeMax="21" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="90" RangeMax="99" offset="128" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1.4142135624" offset="248" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="612" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="836" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="956" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _WAAAAAAAAAABAAAAAAAAAAIAAAAAAAAADQAAAAAAAAAOAAAAAAAAAA8AAAAAAAAAEAAAAAAAAAARAAAAAAAAABIAAAAAAAAAEwAAAAAAAAAUAAAAAAAAABUAAAAAAAAAUAAAAAAAAABaAAAAAAAAAFsAAAAAAAAAXAAAAAAAAABdAAAAAAAAAF4AAAAAAAAAXwAAAAAAAABgAAAAAAAAAGEAAAAAAAAAYgAAAAAAAABjAAAAAAAAAA==CAEAAAAAAAAAAAAAAADwPwAAAAAAAAAAAAAAAAAAAAAAAAAAAADwPwAAAAAAAPA/AAAAAAAAAAAAAAAAAADwP29lmZmZmbk/AAAAAAAAAAAAAAAAAADwP2FbmZmZmck/AAAAAAAAAAAAAAAAAADwP5T/MjMzM9M/AAAAAAAAAAAAAAAAAADwP3dRmZmZmdk/AAAAAAAAAAAAAAAAAADwPx+k/////98/AAAAAAAAAAAAAAAAAADwPyUOMzMzM+M/AAAAAAAAAAAAAAAAAADwP5tKZmZmZuY/AAAAAAAAAAAAAAAAAADwPxKHmZmZmek/AAAAAAAAAAAAAAAAAADwP4rDzMzMzOw/AAAAAAAAAAA=oAAAAAAAAAAAAAAAAAAAAAIAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAHAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAkAAAAAAAAACgAAAAAAAAAKAAAAAAAAAAEAAAAAAAAAUAAAAAAAAAACAAAAAAAAAAQAAAAAAAAABgAAAAAAAAAIAAAAAAAAAAoAAAAAAAAADAAAAAAAAAAOAAAAAAAAABAAAAAAAAAAEgAAAAAAAAAUAAAAAAAAAA==CgAAAAAAAAADAwMDAwMDAwMD
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_2.vtu b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_2.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..92a9c96c4d029c9bd74dd8705ea0fb2b01bf8441
--- /dev/null
+++ b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_2.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="11" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="2" RangeMax="30" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="9" RangeMax="99" offset="128" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="1" RangeMax="1.4142135624" offset="248" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="612" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="836" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="956" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _WAAAAAAAAAACAAAAAAAAAAMAAAAAAAAAFgAAAAAAAAAXAAAAAAAAABgAAAAAAAAAGQAAAAAAAAAaAAAAAAAAABsAAAAAAAAAHAAAAAAAAAAdAAAAAAAAAB4AAAAAAAAAUAAAAAAAAABjAAAAAAAAAFkAAAAAAAAATwAAAAAAAABFAAAAAAAAADsAAAAAAAAAMQAAAAAAAAAnAAAAAAAAAB0AAAAAAAAAEwAAAAAAAAAJAAAAAAAAAA==CAEAAAAAAAAAAAAAAADwPwAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAAmvszMzMzsPwAAAAAAAPA/AAAAAAAAAACYmZmZmZnpPwAAAAAAAPA/AAAAAAAAAADMfmZmZmbmPwAAAAAAAPA/AAAAAAAAAAD/YzMzMzPjPwAAAAAAAPA/AAAAAAAAAABvSAAAAADgPwAAAAAAAPA/AAAAAAAAAAC5DpqZmZnZPwAAAAAAAPA/AAAAAAAAAAALizMzMzPTPwAAAAAAAPA/AAAAAAAAAAC5DpqZmZnJPwAAAAAAAPA/AAAAAAAAAACyDpqZmZm5PwAAAAAAAPA/AAAAAAAAAAA=oAAAAAAAAAAAAAAAAAAAAAIAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAHAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAkAAAAAAAAACgAAAAAAAAAKAAAAAAAAAAEAAAAAAAAAUAAAAAAAAAACAAAAAAAAAAQAAAAAAAAABgAAAAAAAAAIAAAAAAAAAAoAAAAAAAAADAAAAAAAAAAOAAAAAAAAABAAAAAAAAAAEgAAAAAAAAAUAAAAAAAAAA==CgAAAAAAAAADAwMDAwMDAwMD
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_3.vtu b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_3.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..88f94aa981ec49ebe7f5d7b2e63cb7659d067d51
--- /dev/null
+++ b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_mesh_unit_square_QU_10x10_3.vtu
@@ -0,0 +1,24 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="11" NumberOfCells="10" >
+ <PointData>
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="39" offset="0" />
+ </PointData>
+ <CellData>
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="9" offset="128" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1" offset="248" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="612" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="836" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="956" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _WAAAAAAAAAAAAAAAAAAAAAMAAAAAAAAAHwAAAAAAAAAgAAAAAAAAACEAAAAAAAAAIgAAAAAAAAAjAAAAAAAAACQAAAAAAAAAJQAAAAAAAAAmAAAAAAAAACcAAAAAAAAAUAAAAAAAAAAJAAAAAAAAAAgAAAAAAAAABwAAAAAAAAAGAAAAAAAAAAUAAAAAAAAABAAAAAAAAAADAAAAAAAAAAIAAAAAAAAAAQAAAAAAAAAAAAAAAAAAAA==CAEAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAAPA/AAAAAAAAAAAAAAAAAAAAACa+zMzMzOw/AAAAAAAAAAAAAAAAAAAAAJiZmZmZmek/AAAAAAAAAAAAAAAAAAAAAMx+ZmZmZuY/AAAAAAAAAAAAAAAAAAAAAP9jMzMzM+M/AAAAAAAAAAAAAAAAAAAAAG9IAAAAAOA/AAAAAAAAAAAAAAAAAAAAALkOmpmZmdk/AAAAAAAAAAAAAAAAAAAAAAuLMzMzM9M/AAAAAAAAAAAAAAAAAAAAALkOmpmZmck/AAAAAAAAAAAAAAAAAAAAALIOmpmZmbk/AAAAAAAAAAA=oAAAAAAAAAABAAAAAAAAAAIAAAAAAAAAAgAAAAAAAAADAAAAAAAAAAMAAAAAAAAABAAAAAAAAAAEAAAAAAAAAAUAAAAAAAAABQAAAAAAAAAGAAAAAAAAAAYAAAAAAAAABwAAAAAAAAAHAAAAAAAAAAgAAAAAAAAACAAAAAAAAAAJAAAAAAAAAAkAAAAAAAAACgAAAAAAAAAKAAAAAAAAAAAAAAAAAAAAUAAAAAAAAAACAAAAAAAAAAQAAAAAAAAABgAAAAAAAAAIAAAAAAAAAAoAAAAAAAAADAAAAAAAAAAOAAAAAAAAABAAAAAAAAAAEgAAAAAAAAAUAAAAAAAAAA==CgAAAAAAAAADAwMDAwMDAwMD
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol_IceWaterMix_T-plot_ts_100_t_0.5.vtu b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol_IceWaterMix_T-plot_ts_100_t_0.5.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..39fb9214d445bc1d604bae843e488752da04951b
--- /dev/null
+++ b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol_IceWaterMix_T-plot_ts_100_t_0.5.vtu
@@ -0,0 +1,31 @@
+<?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="19" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="121" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-5.3624597615" RangeMax="6.3126874387" offset="80" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="120" offset="1432" />
+ <DataArray type="Float64" Name="heat_flux" NumberOfComponents="2" format="appended" RangeMin="10.097414142" RangeMax="19.834346567" offset="1744" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="266.15" RangeMax="278.17081528" offset="4288" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="4" RangeMax="4" offset="5296" />
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="99" offset="5364" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="5644" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="6776" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="7468" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="7776" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9Y1M9VNNzRNSk40TjMwBQArIASwAQAAAAAAAAAAgAAAAAAAAMgDAAAAAAAA0wMAAAAAAAA=eF4ByAM3/PTrZlKjo++/CX9zQW3N77+7MGXzhCbYP8wWhUroa+Q/ZD77uvCl/7+5IngI9Kf/vz+QAGUOqv+/32eNjDas/78XfAGuq67/v8RndgEEsv+/hYzm+A63/7+5UQ6cfb7/vy0x9y1xyP+/Pj13vHlojb98ZqI5doadv2hkCt4lr6u/kJ1rMyoxz78kHnFwCQ75v+SZsn9YLgnABFycnd8M9b+MMXMLoHrXv3woml1oFtC/1PcNIhRa6T/09WS+sjHrP0MMt5KDCO0/t0Fss/Ph7j9W5OYy+F7wP2f/zEtsTPE/UFpX4RM18j9qUI8wIRLzP0D45JlL3fM/BmMIpzM10T8Q4MXKv5vaP4Ae3WHKQPg/m3ntJa5IDUBcp8pjtxn+Px8TVbnA09M/fvPJ67/zsT90I9k9yjaiPzCeg82Cq5E/sCdKhUTtnD9onfJu+OytPxY4LUjk9r0/8EFTQXoA4T9ATVD8VvkJQGs3zCIxQBlAnN6weimFBEDEAW6XeInlP+mF3kAVJtw/EsAgKDq0lT+AdRR+UX6mP27l0RdVvrY/J4BPxK4x2j81aq8TawYEQGQo6LjGaRNAZAZpX5Rp/z8Yj3Rnz07gP9Q68fD8I9U/SvgDll/2jD9SkD0hVR+ePzIDE0eMC68/EhOxiWhi0j+igP2j/Cb8P4A8iVO3JgtAnP/bY8/I9T/UiapHSijWP3mMfmjJQ8w/miOvt5UIfT/GAaWIDoSOP+vsYm9tmqA/7s3xjj4mxT/WZXrlH0HwP/TBoDi/8/4/jLEikv9P6D/gLBun3GXHP+Jkvw4xf7w/kFhJkiE2Aj8IPG4cWS45Pwk1WUtxSnE/jx/6uKgepj+p2cqKCW3RP2Atwas8aN4/yMKiSW05xD/h+ViZBNmTP27Zvhm/s10/zT5iMb2/fL96MNb4KPGMvynbp2Cij5i/b8Ab0dQttL/gtMmxbCrev3hIyrZBf++/bO/hdJJm3L/kxIPim2/Cv94vnh+Tkbu/yHz8UvPRjL82kiJd4lWdv4Xcs33wuKq/DLtrjX61yb8W7r9eeHDzvyDOrv6njAPA6KH8M3W68L8YWrJBK63Tvw0dipH6zMu/eVafBgSilb+voqwemBmmv36rJM4HlbS/nMxPSQiq1L/qYMzEU1b/vx4Jg2B9OQ/AAOfVtDxb+r/wZqrzgyLevwae8LOV6NS/5FbHYw7bnL/H7hgPP4itv9xi0F+Xzbu/HqTixFB53L95CPf/Fp4FwAhM1qYocxXAXw7Buf79AcBW1Tf47Evkv6951CGu6tu/bK/2vQ==AQAAAAAAAAAAgAAAAAAAAMgDAAAAAAAAyAAAAAAAAAA=eF4txdNCAwAAAMDFZdtatu1aWrZt+/N76O7lAoF/UY52jGMd53gHneBEJznZKU51mtOd4UxnOds5znWe813gQhe52CUudZnLXeFKV7naNQ651nWud4Mb3eRmt7jVbW53hzvd5W73uNd97veABz3kYY941GMe94QnPeVpz3jWYc953gte9JIjXvaKV73mdW9401ve9o53ved9H/jQRz72iU995nNf+NJXvvaNb33nez/40U9+9otf/eZ3f/jTX/72j3/9BxAkHF0=AQAAAAAAAAAAgAAAAAAAAJAHAAAAAAAAUwcAAAAAAAA=eF4t1Xs41GkbB3Ar0bZObcohOfV7d9Ok1UyE1Tw5LJKENeXUTLSKEpt3QxExuLQOMUZpqIyzSSjjMGXneUbIsQNZkV7nU1mnphLFPtP1/vm5rt/hfu77/v5+UlKZFxLjydQbRal117spSEpqkCPtYitcyf6/c4Y4Gg0PqSEMivXgbgJJJaTtcXvTJ/QxV7L5aqli6daS3dTVX7/rremS3O+rOl5mSG2PmRlY0ynxpP3SfUPqYkvuAO2ZxKOOSy6G1JYAZ5OcJxKnMR417qDWOd8hr+3AnssRD7E2U9Pl3Bcn2rAv66nIvdISBuup+AxILGWd75ltIrRJqrx4QXK9lJ2KprS1MIGT6zqAny/+g9bKHWLBgMa54Qxcf1TXeoX3hQXwbmLg1ljsaWijbNbFh+1KZIWHLyjIqksYWh3YAKtXXK/pt1JQspPzvdlPz2G7Cduwd5GM/hc8FZS+LIC7ruxy7igmodX/jGdvCSuEpYyjWUkdBIo9s84p83ge9NCv9qg0I5BjjLfXweUsWBxdLMoj4/40tDgHvq2h2q2ln30NsMVHqybP8qllvI+7ug5gu5lsd7KqoX5mbF9r5opdetFwQaaOGjb6k3D8OPaLHlnrTCGVv20xuT8Iu8IuPrZGQH0ZstFkIBI71f+p+ZlM6l9txIehBOxq890vhpCwVJ/b3xyLPS3jemFTt9D+DIWj5UcgcvaMIFHrIWANP+6R1LcqngeKcu1AEFRsIKl/6Z+iZI2ATqAVuN4uGp9Pl29ED9K+C0JGbNcN4/MPnrp688NCHtC5ZeL3Avdnpa/x6YYEDtBISuyQ9E/KWUe5LoQNtgcwlCX9DXO/FUbMJgHFwuiv/bd3ijxgunIJMA7Pf52P22jhp7FbEeDJ2Nm5Jrw/V1NsWOncRJDtd3v3GLZmdaGVBj0dzCTQTU2x5X8xfpQdex04bI44fvMxnpflbxvzTnLAcqO3zGExGYX6zWh3bygABgEOL0u8SUiowNS0pzUDG9GWKJ0KAk2Kt4YtxjWAYy5RJlO+BLptvqmgSsAH1uq+srmnCWTTkW/TrhQOKtTmynXx/vqUx+z/IepPYM5fDgrB3hvZKP+NgAWsmdIODthmRtZZstYZIKFd+a5TMwWFj2asz89nATmXCOH6d2QkjRr3mvpzgbuXKIvuQEI3LsWb7JwqAXcS+lVFtwnkG9k09GaCB7TjnhzcHEgg80O3855554CxWjV2Fp6fqTbyp24NB4aL+t+F4v228Z7n+z9IBNqZ3U8nsFtgzxdZlQzwSdC33+g5BRVP/q5Wr5EFDh/tulCM+6O/IrA658YFZpqTzI4ZMnqkZDYg1ikCvLCB8E4aCcmMq/av6KbBwZ8tIzYWESg9GpVVemXAxAiFnSm4Hm7Y7Nvw3FBosaVCNzyOQDKpsfKn4y6CJsclHx7Or70SvHdfLxm8pZ8yHMF2/KJ+7ORqJpBXqmpdxh47PaKhSM8D1wSDAaiBgn7boeDMTCwFGk51vQ5vyUhdcW0OQ68AiJzETxr+IKFyflm1eLAMvhF2nS6tJZCrvPv05do78Mhnmh/Vn0CMviEXmswNmP+t6PjHCwTSiA5plq4/DwKWWNWd7RQ0pTzfN8v6EyRVfePphH1/f3aBDZsD7s0xTuzA34/+ZkulvOhiUD7sF6YCsYOSaWTFStAq6B+ymSYjMk+krOB6HbAd+1e0uCRU2DxioKnJh4NznJSSVgKFvPpsV99WDhe2mZWo0wnEf+9Nj/x4E2bvUVJ9juujpd3RVzrmA37mpKxzw+8bibllbDxzCci+Nw3ltFDQg6BgnpYRG9SNyWcrNFFQn+KS/QlREZird1+l1VDQ8InX7+4erQT792lenF8go7+f8SJLjIKh1fvPzKZ2EuqVT9efCroPX4rLTjb1EGhneu4/LNVSyGPK2Tq4ECiJuTw2zM2CnqnPxMYeBKrzHMjapuYPRW4LPe9wPhPYroUzWefgXLD2sb/xfmzf3H+qU8cC/GIXtEFfREEfZhxbFFsywLe7Dv03qRz3U63nx22RReCDS0bH5S9kZHgg48pjag5s3hU6Pf6WhALT3Mdj2aUwJm+ToudrAlUU5b570VkAg51/eDJhR6CZ2Zezv268Dg+csLgpOoTz1sGNbFtlwmFZi7ZzeD5XjQ0OOg6xoQv58EQ5zs/8976kfRY5UFYYWsXE85m6Il/6MD4P2l/j+fB5OG9765NdtG5Cg5wN8HtZCtrjnGIr1VkM2W6yWk4fSehQwr0Ru9p8eCifHfF7P4Hs5pvaLc9zYbTgRr2/JYH0LqnkVHuxYaqi3h1/awKBsB1H4sxT4U8dr3bsxHlSZWzx2GeXC/dEsMih+P/i5Nkr1ztSAdGR2qcPGnG+7J8ZZ3gIoQe3Rj6/jILIulbNTIXHUMbcdvLUGgrSaWxwTOishGuWA44/nSAhrwHi/OirPPi8e7ZUu5dAFzXpDFobF648kKdPUgl0KV496sq76zBze1OPszmB/gUNFs+nAQAAAAAAAAAAgAAAAAAAAMgDAAAAAAAA0QIAAAAAAAA=eF6d0mtIU1EAB3BTQ6FIhWlmTEuXpTJ0WZKP8qjTrKkZpoiC+ShN7LE7tFGSlo8eph80MVNs4iMX5YRZSV3xFcIUlzYhmxojX0hao3wOWes/rvShvnm//OGec3/nf8656el4GoQkfTOnJEopnUr9zfR/xrea1cXH9zW8FhK6iTX+RSEkY+8XHu2ZFJJfU32RsVohsTatsSg3pYhkP9U7ZEsRz6DwXEs3ivQkOx3lB1Dk7O3VH/lnqP/6bTV5m55XCuNnBjPrNbow6382Z/p4zDL9bvYzfRXNTP/dd5n9zC7NRNUjHdmB8eN4/8JRU8fCPOXFY5oofJfoaudwD06rfN6EhmuVksspPESR7Em7O2p/iqy8UoVEoAfX3ln6FE50U9jBMTgrBSOe1nA22GZTp+BkfeBpC+FcCWlWyeHoHqoX4uFsM5sxlcJxn9GseMGxfOtmUQdHFateGIVjqO1+sAOOOEd2nw+nqnjv8i04LcKk7kY4st6drHNwBtsMCc/gjF8q38WFw5P3J9XCEbjmXB2BU5RNrW2H01rV+i4Azo2w6nYRnAFBc0wFnITA61QGnBxNXmQnnD7y8rEPHDF9ZOQJHK28fWIITmXp4SzDhJAcaFvr4cERLeY5psApsHGdK4JD0YaOEjj2mV2sCTiSip/TfDgGLwldDWfNL0I/AIcl3mDp4DStLuo4cOpEOloA53S4ajQfTr1eG9QGh6vOGFmH82Z+1SUGjkOJoN3oePv+9jM6+s5lyQocmYm60wFOpyzHwxfOt1BSVgxH0WBmq4RzQawos8Z/Evz1uc95ON+lq6QGzqfQLifjvvyiy0fX4ahZKpHRGfStsvKG87HyREwlHPe5ODIN55reJtEZDs/b3ikTTnJcdqnxvqpmGxOM96W2uJxmPOfUDj7NhcPxN88/CYdw7Pxb4LC7JkKX4LAdh/t5cJRpw7UiOH8AbwjE8A==AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAAEQAAAAAAAAA=eF5jYWBgYBnFgwYDAD0/AZE=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAArwAAAAAAAAA=eF4txdciAgAAAMCQWSFbJCuzsjNDiuyMyo7y///gwd3LBQL/WtzqNgfd7g53usvd7nHIYUfc6z73O+oBD3rIwx7xqMc87pgnPOm4p5zwtGc86znPO+kFL3rJy17xqlNOO+M1r3vDm97ytnec9a73vO8DH/rIOR/7xKfO+8wFF33uC5d86Stf+8a3vvO9y37wo5/87IqrrvnFr37zuz/86S/X/e0fN9z0r/8AXBwTVw==AQAAAAAAAAAAgAAAAAAAAFgLAAAAAAAAMAMAAAAAAAA=eF6FlU1IVFEUx4cKLCK1j41E2YdKpWk0RELQU2yhQkrSQhAMGcpS8yMJrUW0iBnJjaBRKqFG0KZNuoghMCqzhvyAynRgsCEjkCZ46Ri4CJvnuf/e+8972VnMcH9z3333d+65Z1yutULX1ub/+93lajjfHQu/bV5VhcEDNt61ctTtdn+w8ZYzxvygjac+XInFZxtPT4wt456z8d7iC7H4ZuPbfMb68zZ++9VYLCL/9HT20zVnP11z9tM1Zz9dc/bTNWc/XXP20zVnP12L98sYto7Nc7zXbX3e5GM3rO8z+YrHuj9LPRQan2Eb9yf2WPxNvrndbclX/PyAjQ+ucr+NI9jPDPYzg/3MYD8z2M8M9jOD/eLn2+8H+7lcS7VSf4mLXIc1rTIO/OQ6zPQKr2jm96Z2CY/Ocx16Hwi/nGmMTK+iQeFfz3IemkaEF17nvC0Eha/v5zxj33nv+Fx+X5R94znwFy3CD8yw1w6vjNvOsddcp/DZCfYKDQiv2WmMTK/tT4R/KmKvty+F111lr3S1jx9xXqMq71MT7NW0RfaHPIGHU4TjPeBjacKnS9lrJFvGWc/YqzdX+FLy6qp/+akC4f357DWl1n3UyF7IY2Mfe6Fu0ibZK1At+0BdgHuuCUdewW/6hGNf4P5O4UmP2au5T8bZm4yR6bVVzUs6wV6HVF5217EX6iDhPnstqLo/HXdes6PSf3EPwKtCwoMD7LUrIhznAH54WTg8wIfWGZ9hDd/gR5RnyTH2iqpzHa5mL9RxQQ971at7mzDOXt0psu7eLva6s0847g340EHhqDtwf45wnBt4qdp35jL/fyWrc3qaw16oy2IPe82oe5h/l73mVN9pi+sbuWXyfJaXvWbKhaNPgF+pFD6exl7Yxy0fez1Xecc5g6eqOkOewHGvGirZq171ke+d7FWk+mYowF6T7TLvUit77e8QPtLCXhG17pcU9spTecS9BEfdeELstVHdE9QFOPrCdDl7RVQf3NPBXhsUX3zDXu9fy/hXLXtF1Tz8D4AjL+ij4D5VB+hD4Kh73GPwEnXPcQ/AM1RfO17GXh/VfJwD+En1XniA/wHeUUBhAQAAAAAAAAAAgAAAAAAAAIAMAAAAAAAA5AEAAAAAAAA=eF5d0EVXEAAABGHs7u5WUOnuMmilQVA6pBt+vgdnLruX7705bkHB/53FQvwQ2ovwfWj/iO9C+yd8G9o/45vQXoyvQ3sJvgrtpfgytJfhi9D+HM+gv53DciwM7RVYFNor8WNor8JPob0aP4f2GiwO7bVYEtrrsDS012NZaH+G/uhv57EBy0N7I1aE9iasDO3NWBXaW7A6tLdiTWhvw9rQ3o51ob0D60P7U/RHf7uAndgQ2r9gY2j/ik2h/Rs2h/bv2BLau7A1tHdjW2jvwfbQ3osdof0J+qO/XcQ+7Azt/fgltA/g19D+A7+F9p/4PbQPYldoH8Lu0D6MPaF9BHtD+2P0R3+7hKPYF9rHsD+0j+NAaJ/AH6F9En+G9l84GNqncCi0T+NwaP+NI6H9Efqjv13GPzga2mdwLLTP4nhon8OJ0D6Pk6F9AX+F9kWcCu1LOB3al/F3aH+I/uhvV3AF/4T2VZwJ7X9xNrSv4VxoX8f50L6BC6F9ExdD+xYuhfZtXA7tD9Af/e0q7uBKaN/F1dC+h39D+z6uhfYDXA/th7gR2o9wM7Qf41ZoP8Ht0H4f/dHfTuE13Ant13E3tN/AvdB+E/dD+y08CO238TC038Gj0H4Xj0P7PTwJ7afRH/8B/hprUw==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAxgAAAAAAAAA=eF4txRF0QgEAAMB6LwiCIAiCIAiCIPgQBEEQBEEQBEEQBEEQBEEQBINBEARBEARBMBgEQRAMBoNBEARBEATdyUVCb1HHHHfCSaecdsZZ55x3wYGLLrnsiquuue6Gm2657Y677rnvgYceeeyJp575w5+ee+GlV15746133vvL3z746JPP/vGv//zvi6+++e6Hnw6F30UcdcxxJ5x0ymlnnHXOeRccuOiSy6646prrbrjpltvuuOue+x546JHHnnjqmV/7BioOAQAAAAAAAAAAgAAAAAAAAGQAAAAAAAAADAAAAAAAAAA=eF7j5KQ9AACx7gOF
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol_IceWaterMix_T-plot_ts_50_t_0.25.vtu b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol_IceWaterMix_T-plot_ts_50_t_0.25.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..d1690e6e96f1cc560650b87ea523ef391884609e
--- /dev/null
+++ b/Tests/Data/Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol_IceWaterMix_T-plot_ts_50_t_0.25.vtu
@@ -0,0 +1,31 @@
+<?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="19" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="121" NumberOfCells="100" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlowRate" format="appended" RangeMin="-2.4407851286" RangeMax="12.43635259" offset="80" />
+ <DataArray type="UInt64" Name="bulk_node_ids" format="appended" RangeMin="0" RangeMax="120" offset="1432" />
+ <DataArray type="Float64" Name="heat_flux" NumberOfComponents="2" format="appended" RangeMin="10.734266857" RangeMax="26.362006114" offset="1744" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="266.15" RangeMax="281.85595205" offset="4264" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="4" RangeMax="4" offset="5488" />
+ <DataArray type="UInt64" Name="bulk_element_ids" format="appended" RangeMin="0" RangeMax="99" offset="5556" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="1.4142135624" offset="5836" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="6968" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="7660" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="7968" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABMAAAAAAAAAGwAAAAAAAAA=eF4z0zPRM9Y1M9VNNzRNSk40TjMwBQArIASwAQAAAAAAAAAAgAAAAAAAAMgDAAAAAAAA0wMAAAAAAAA=eF4ByAM3/C6myWJMrOK/JRr4WUvw87+r1Cy2IpvuP1+dnPbzcwxAp/jiGg9+879uqEPp3+71v+jX7cDNX/i/D/KbosnQ+r8MXsdq2EH9vwIBp1Qss/+/Wyme/JoSAcCD2eFBQ0wCwAxNf1q6hgPAtFVSwS6Ewj/2wLtX7HTSP0FEjMtpnNg/AJc7g3KRdD8InyziZNLQv+rmzFgT39E/SB8UdOp63D+/604K0njgP7om1zqgeOI/o53nzi5g8j/Dd/1Iy/rxP+TH5QeUlPE/B3k564s08T83/Ln40/7wP7n89JXm6PE/AOtY/Obj+T9hh1zdiLEPQEu8qHrWJyBA3FHk9i+dCkCo7XU1337bPwA3PRbmSeG/5WW5Nv095r+QZf239Crkv6TFSmzDd+C/kvBj5Z3T2L/7Wg7OJ5LQv4cWo5o9k8C/iJOTgiM1pD8jRoem+UO0P+OxqI5Ex74/n9gHqa2ixT/3yrQyeJ7QPxrZkWBFHOI/MEf7P+WJ/T/urR5+CwgYQEA6Uptp3yhAnA6McpXioD+zK4FXJPewP3kxUL50Ero/7W6QJ2WQwz/UcdL0y2LTP5R/CLflXPA/ithKFl+uD0C6YzXI8KMiQDiBKcEyQSNAunVlJA0gmz96M2GyRFerP8BC+/7AfrU//ADxcY1uwj+hXTR3jBPbP7QsMx+V5f4/AtOBo9kdGEBcdnhKzmIdQIQNXriN1gxAOWz7tct6lD8zx/YM/8GkP5VIK4a3ELE/JChRJgaowj8sa18wV+3kP1zhVzDbyQdAP4EFdx29FUAlkdiiOvUIQG/fd/BtJe8/yDdxFVuqiz8qQut+RVecPxMOseggd6k/1CJtyuY4xD+aL4kVpP/tP2+rSAtsEAlA6+fwevKUBkDyjXhwoMPsPwn87ZOyjdM/GtYimFm7fD9cXvy5xjWOP34MF+SHfaA/lAiCPZsfxT8QNRGLNx3wP3Z6WujY+P8/sOpnbfa27D/62nFX2gDMP/QlLfMZi70/RULRvJrmQD/+DjzwGyJbP8BU2Peqr4Y/fMQu05luvD/uWwe0at7hP/yg7h79AOM/LoaIJGrewT/I53B1ub2UP7ZDwnhYW4A/Nh32gWGKeL9eWuaC9DiIvzQ5kWnWvZW/bHTbzo7YvL8wfwjVVa3Vv+wfmfgot8a/56pxrbjstL8GtQdFySq0v7Z2cH3bYra/nIGreNKfib8aY6EZhb+av5Hqd9DP2LS/WCrMUYiO478jiezrubbyv7ChnImo+Ni/bMG90KbtxL91STeKX9PEv/70YypQO8e/FVX+dg==AQAAAAAAAAAAgAAAAAAAAMgDAAAAAAAAyAAAAAAAAAA=eF4txdNCAwAAAMDFZdtatu1aWrZt+/N76O7lAoF/UY52jGMd53gHneBEJznZKU51mtOd4UxnOds5znWe813gQhe52CUudZnLXeFKV7naNQ651nWud4Mb3eRmt7jVbW53hzvd5W73uNd97veABz3kYY941GMe94QnPeVpz3jWYc953gte9JIjXvaKV73mdW9401ve9o53ved9H/jQRz72iU995nNf+NJXvvaNb33nez/40U9+9otf/eZ3f/jTX/72j3/9BxAkHF0=AQAAAAAAAAAAgAAAAAAAAJAHAAAAAAAAPwcAAAAAAAA=eF4V1Xs4E2ocB3C6TmrVxCa22dra2hAWoeSSdH3spCzRaVIqJYdYV4XuyiOklSN0ug3npJTUKbZvHUooxYqwjusoTpejyyrpvOev9/k8z/vH9/m93/d9jYx2ZNYr09XWn5Y2vN/Oh5FRk1/pRLVa5itrc3rthkesaz73XrFQJDN9/CRMhHm2ZWdUdSLkZsbHbTKdQvaHGj8QKdWTD76OdEkQEpfUBh/JVrPWMtW3D4mJJ69V9JxVB1jQekJT7YmpO6vSL6ibbqbMNalxIL6kOvIiX309Xr6v29+JeJWLpOWyeoz/fV1LvYT4WMSfumL10OBAu8sKZ2KZRhJwQ736vb/DplYXYvp3S/mfauvWfFF1mCscx60MPL/NoJly4US5G8mfdld9ZiiBgvO/HNvSRdy9PvBnfSIN7365WP1xwA39s6cVlsvp2JRyJ3GcmStGp7T4S/dMxDePeTMZg7ZwCKdjawUNvcn9i3RLRBgeHv1R1GyOYyOUE6+uESG1K+uvL70M9K0ZdruCzKdtXfHG0e+tMGZse8j/89rff9Lk//n9a39neaxeiB22jry7xFuNN0R8WivAPv8ZA+XEqw6co+5s40OkZZSWEvs05MFJzoPKj1NeQxz9yitD387FXJNXc3L6WHhJ895eEMFBTWjt154hFvw9fELyn7IxaqPx4UntbPhlV2S4mzFx+qmvNP7CFIiaSxu7a9hYpdH6usWL8EMqz/nRycPR7M6d36hC3A9Jc7rmz4f+IHt+gpiPyLYQBmszH+sTlM9W23JRGJoBehwfBcqqGPYSFp5l+j+s38bHDauZhwNnM3AieJb9WdKf0Kai1BRbGq5SXO1aiGXuZTm57hQUrQt60Uc8VOa+a4q7QfOHuiTxHbFjwoZ53FkGTdBrnyxf0p/z2tvuWbMoCKa/O7uF2O6g2XpvOxp8o+c/zyGePtZSQfNg4IzpSPfdxLoTD7RJi1k42ucxLn6PEJ84irfRAi7skyPll2KFmD16+jNLUz7Mezj7fl0oRO7lthu6lQKsbzDNjcsg67ICL6PhAnxlT47l3uEhXvPv5vzZBo2VftmdRtLfiTcjF1I9KMjned6wOyzGjyvlQz72NMgH9n7IIk5a4GJpPIcBPs38ojexH31plNlCFi517A+tTxJjbJ74BovBhYgaLLgbKkblmsQ5wiweJqzI9hp9TYT+5cq6gwk88Fgrxu33Jn1YYGFi8cYGX7fTZq3x4gAHXpb+7GHQHDu0bk8yuT8xTwcKlSTPT/rYtuHH7aGXFr9aR/K0TrXrySL+/e+4a1Qyn8Sf0stlZH9i/8uwlZ4snPz2MThirz0mHAn5HtrLQZeh5HTORHs46L1VqaRXzSXh1Ey2GHEJWWd0dDbqFp6bxSS9iVHXBJcXWWPDyMNLvrrYoLuZGaIieYpU/Wlscn/PPE8o3EHy1CpsN4YRn6JeVteT8xKYTf2UR7xLXNnZ78bAMImIWfvQAaUjEv0r+CxQP/qkCEsc0Hva50DQfA4aX0qLL2mno62u++uHGBZUeSUeZ02ECK+8EOM4YIkglr682JiLdJXXx1FUa2jboxuTwjiICfwekznHoPHxLWv4QN6P2mE6vwCSJ3koOUUodcKiDZrHj0gfqxK3yCyJo/52SbrnzMBKPj1OssgJnlFdNTP0TPgZzpkbj3eCx9i9F4f6WZDvMRr39rM9XK97yWhhDMx27Um2HsbHpUnccxhvCVvV1LuNjVwEHL+Sxv5iBQrrtcKtlYvIqh1Dx3wMml3CUbeqyfulzyh9OMOTAvoI011lxIlU4e4kMQ1LNhYIFxPnXrV4M9mBgaD7QTbsBxKM+aI8f8iRCQ67qn5gtASMlyX5c1dbItZuYcrDVDFMn9V2ZMy3wAfFWNemizxkRg3Z1L1hIOe5CUUi5UF7XRs41GcFG2+6oieIhwWVvp8VSwwapVHx/haZM+RNNZXzFlCw7ktlxD+Bzqhvq9bpyXldD5q2tWy5M3I7213NOAy03dfc6fJyhjRz5ciOuVYoeiufohopwTSLAI8xZuYwYqa489KEsDl6K7C33Rwee34rMK7ng9lwPkClY6ArQTq1+gkfj0N0AcxuKwjmHM/ktvAh1d5M/22VQXMu0yj4aYsLBrusTjWHUPD+uCo8/IULgp51Tsj2poGRvqxiq9YFfMtN25+PZ2Bv3j/erwpc8HTxtbrdoEP29tGYJo0j/Idfri70NcNJJpd5t1WAIE/JqaIac5j6mnC8lwlQFfnY5F4zA/VN0188kQlwa2m4sqbTCuLYmotsuQCDV7OHLYsyaKKuuv8hJf8H4y/LUw2xFPh1xax2DXVFVWFBhT6chpGKiHtXZK4ozdnBjLZmYJvTYHBH2UwoJyf1KjiTwEvtmjD4ZDqsnco3L9LSsCZW4bmgUYi0nRLj2HpzDHSZm15+IURKwMxYWQcDvqf6THg6IbxyP4/ofW2FXa0N8/I6hPgPwRCr5g==AQAAAAAAAAAAgAAAAAAAAMgDAAAAAAAAcwMAAAAAAAA=eF6lkv9PFGQcx+NmAYJTthMEppNvHVbaQZ5A88YbBDFQjpNQC+HK6trFFz8fSYEKBoLDpYKXJPQFqWOhIHBABAVmSMmXSxOyTjAOlYhDEHA5NoQBPT+0/oGeXz7bs+fZ5/1+v95arTgGgvbf6VLcNddQxmi+Vq3LWs7/3f/fmZhvLixoI4Qatrdc+I2wtr35Ts8kYdbytN3EUwzz/Fm543pGhJvtKxuDGG2B6bnRuxnP77VWHUxiVBze++vpfEb9zqFjjwsZvU/OrArJYaz9weHz48x4J8NT1nuA0eQX1OjyMiNmXKXUhDPGDdruyi2MvPis2CmZeC8tHlK4MrJnvJqlTzAq15VJamYIkxHOqm0TBAUXfXb7LuG9T+3G2EwI+zF3s/01wooHczlfXCH8Ln33ekAL4Zxy0vVGDcHs47GwvY6QdTmhQtYqfCpUOrtOQok6O/F+H8Ez/f1TJguhNPLw8MX7hG2GlLBTYt+Bgg7vVBuhq7SLd61gvDirWLpsJNhPcbPL94T4sOA6MhHWn3R1MN0itEs6rV5/EYwutersvwlvloxE9C8RRqNWz/uL/57SaW22G0NZWucU0CD0hm+x2rYTznQccRv4hcCqZ6aqhY7Oh8NFWQ8IyYHOeTHzhKp2jY234JxhLDeOrmFkWvp6j4uc5IH1G4yNhEPjX7nPdhD8X7+ZEXKT0N279MmJPwlxiSct5keEmEmHXI9lIs9X5claKSO/aIcl2ZvhJTGnmTYz3qbFMFkToc/sGD52lVBTued8lfDjnm6oShJ5JEfFSTcJHXKr8/RDR8YHi76BBesEt2HnWD+54NgzNuYeyjCFfJxZ/g0By5472N8t+Jzx+dBpkGAzFO8XNU0IykzwOCphpK1a6GpazZiyRA/85Cv6W+0UsHIrw7fkSLBOJXjHXlq5RvDTS8qX9/xMeHb3hCRT8M6P26jeILik5t2t7rdn5GhqHx0SOuCpmfF5gZF6+16wNYKh/jr6RuN+RutL+nT9t4T9+5pfaxP5JvhXtN4bIei0mjSJ8FOoN7UonMTeP2TfRYk8i/ec33dMyUjZ1HDuUizjykW9ck7HSJKMSBxFbxYa36qtF30ptP3yllrkcmFX4tU7op8DW1MGTwguRy2R+iCRR2T849hRoeP0tGbxIw3jjeuDZaHpjH8AL0G3vw==AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAAEQAAAAAAAAA=eF5jYWBgYBnFgwYDAD0/AZE=AQAAAAAAAAAAgAAAAAAAACADAAAAAAAArwAAAAAAAAA=eF4txdciAgAAAMCQWSFbJCuzsjNDiuyMyo7y///gwd3LBQL/WtzqNgfd7g53usvd7nHIYUfc6z73O+oBD3rIwx7xqMc87pgnPOm4p5zwtGc86znPO+kFL3rJy17xqlNOO+M1r3vDm97ytnec9a73vO8DH/rIOR/7xKfO+8wFF33uC5d86Stf+8a3vvO9y37wo5/87IqrrvnFr37zuz/86S/X/e0fN9z0r/8AXBwTVw==AQAAAAAAAAAAgAAAAAAAAFgLAAAAAAAAMAMAAAAAAAA=eF6FlU1IVFEUx4cKLCK1j41E2YdKpWk0RELQU2yhQkrSQhAMGcpS8yMJrUW0iBnJjaBRKqFG0KZNuoghMCqzhvyAynRgsCEjkCZ46Ri4CJvnuf/e+8972VnMcH9z3333d+65Z1yutULX1ub/+93lajjfHQu/bV5VhcEDNt61ctTtdn+w8ZYzxvygjac+XInFZxtPT4wt456z8d7iC7H4ZuPbfMb68zZ++9VYLCL/9HT20zVnP11z9tM1Zz9dc/bTNWc/XXP20zVnP12L98sYto7Nc7zXbX3e5GM3rO8z+YrHuj9LPRQan2Eb9yf2WPxNvrndbclX/PyAjQ+ucr+NI9jPDPYzg/3MYD8z2M8M9jOD/eLn2+8H+7lcS7VSf4mLXIc1rTIO/OQ6zPQKr2jm96Z2CY/Ocx16Hwi/nGmMTK+iQeFfz3IemkaEF17nvC0Eha/v5zxj33nv+Fx+X5R94znwFy3CD8yw1w6vjNvOsddcp/DZCfYKDQiv2WmMTK/tT4R/KmKvty+F111lr3S1jx9xXqMq71MT7NW0RfaHPIGHU4TjPeBjacKnS9lrJFvGWc/YqzdX+FLy6qp/+akC4f357DWl1n3UyF7IY2Mfe6Fu0ibZK1At+0BdgHuuCUdewW/6hGNf4P5O4UmP2au5T8bZm4yR6bVVzUs6wV6HVF5217EX6iDhPnstqLo/HXdes6PSf3EPwKtCwoMD7LUrIhznAH54WTg8wIfWGZ9hDd/gR5RnyTH2iqpzHa5mL9RxQQ971at7mzDOXt0psu7eLva6s0847g340EHhqDtwf45wnBt4qdp35jL/fyWrc3qaw16oy2IPe82oe5h/l73mVN9pi+sbuWXyfJaXvWbKhaNPgF+pFD6exl7Yxy0fez1Xecc5g6eqOkOewHGvGirZq171ke+d7FWk+mYowF6T7TLvUit77e8QPtLCXhG17pcU9spTecS9BEfdeELstVHdE9QFOPrCdDl7RVQf3NPBXhsUX3zDXu9fy/hXLXtF1Tz8D4AjL+ij4D5VB+hD4Kh73GPwEnXPcQ/AM1RfO17GXh/VfJwD+En1XniA/wHeUUBhAQAAAAAAAAAAgAAAAAAAAIAMAAAAAAAA5AEAAAAAAAA=eF5d0EVXEAAABGHs7u5WUOnuMmilQVA6pBt+vgdnLruX7705bkHB/53FQvwQ2ovwfWj/iO9C+yd8G9o/45vQXoyvQ3sJvgrtpfgytJfhi9D+HM+gv53DciwM7RVYFNor8WNor8JPob0aP4f2GiwO7bVYEtrrsDS012NZaH+G/uhv57EBy0N7I1aE9iasDO3NWBXaW7A6tLdiTWhvw9rQ3o51ob0D60P7U/RHf7uAndgQ2r9gY2j/ik2h/Rs2h/bv2BLau7A1tHdjW2jvwfbQ3osdof0J+qO/XcQ+7Azt/fgltA/g19D+A7+F9p/4PbQPYldoH8Lu0D6MPaF9BHtD+2P0R3+7hKPYF9rHsD+0j+NAaJ/AH6F9En+G9l84GNqncCi0T+NwaP+NI6H9Efqjv13GPzga2mdwLLTP4nhon8OJ0D6Pk6F9AX+F9kWcCu1LOB3al/F3aH+I/uhvV3AF/4T2VZwJ7X9xNrSv4VxoX8f50L6BC6F9ExdD+xYuhfZtXA7tD9Af/e0q7uBKaN/F1dC+h39D+z6uhfYDXA/th7gR2o9wM7Qf41ZoP8Ht0H4f/dHfTuE13Ant13E3tN/AvdB+E/dD+y08CO238TC038Gj0H4Xj0P7PTwJ7afRH/8B/hprUw==AQAAAAAAAAAAgAAAAAAAACADAAAAAAAAxgAAAAAAAAA=eF4txRF0QgEAAMB6LwiCIAiCIAiCIPgQBEEQBEEQBEEQBEEQBEEQBINBEARBEARBMBgEQRAMBoNBEARBEATdyUVCb1HHHHfCSaecdsZZ55x3wYGLLrnsiquuue6Gm2657Y677rnvgYceeeyJp575w5+ee+GlV15746133vvL3z746JPP/vGv//zvi6+++e6Hnw6F30UcdcxxJ5x0ymlnnHXOeRccuOiSy6646prrbrjpltvuuOue+x546JHHnnjqmV/7BioOAQAAAAAAAAAAgAAAAAAAAGQAAAAAAAAADAAAAAAAAAA=eF7j5KQ9AACx7gOF
+ </AppendedData>
+</VTKFile>
diff --git a/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/Heat_conduction_phase_change_(manufactured_solution).pdf b/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/Heat_conduction_phase_change_(manufactured_solution).pdf
new file mode 100644
index 0000000000000000000000000000000000000000..3c6a36d41ed858b01ded64d9fbbecffcf0db0b7f
Binary files /dev/null and b/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/Heat_conduction_phase_change_(manufactured_solution).pdf differ
diff --git a/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/ManSol_IceWaterMix_(comparison).png b/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/ManSol_IceWaterMix_(comparison).png
new file mode 100644
index 0000000000000000000000000000000000000000..b455d6bd748540bdd66d1360338481a2e8aadb1c
Binary files /dev/null and b/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/ManSol_IceWaterMix_(comparison).png differ
diff --git a/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/ManSol_IceWaterMix_(maple_plots).png b/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/ManSol_IceWaterMix_(maple_plots).png
new file mode 100644
index 0000000000000000000000000000000000000000..144f7fa383fb3bbdfe042a76bac22194a5623b79
Binary files /dev/null and b/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/ManSol_IceWaterMix_(maple_plots).png differ
diff --git a/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/index.md b/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/index.md
new file mode 100644
index 0000000000000000000000000000000000000000..7e4a9aac57aed7a084d96f697348474aff701dc7
--- /dev/null
+++ b/web/content/docs/benchmarks/heatconduction/heatconduction-freezing-melting/index.md
@@ -0,0 +1,31 @@
++++
+author = "Tymofiy Gerasimov"
+date = "2022-11-03"
+title = "Heat conduction with phase change (manufactured solution)"
+project = ["Parabolic/T/2D_Ice_melting-forming_manuf_solution/ManSol3_IceWaterMix_Scaled.prj"]
+image = "ManSol_IceWaterMix_(comparison).png"
++++
+
+## Problem description
+
+This benchmark is testing the extended version of a classical heat conduction equation which is capable of modeling ice formation and melting in water-saturated porous medium. It is sometimes termed 'heat conduction equation with phase change', or, with a slight abuse of notations, simply 'T+freezing' equation.
+
+Since the equation is strongly non-linear in the temperature variable $T$, a carefully designed code verification is a must and is performed. To this end, we use the concept of a manufactured solution, when a prescribed function (in this case, it is $T$) is plugged in the corresponding boundary-value problem yielding the heat source-term, as well as it also provides the initial and boundary conditions. This recovered/extracted data-set is then used in the OGS code as the corresponding input, and the obtained numerical solution can be compared with the prescribed (manufactured) counterpart.
+
+See [this PDF](Heat_conduction_phase_change_(manufactured_solution).pdf) for the detailed description.
+
+## Test case
+
+In the benchmark test, we restrict ourselves to the problem formulated in the unit square and on the unit time-interval. The manufactured solution is constructed in such a way that its evolution in time mimics the simultaneous ice melting-forming process in the domain, as in the following figure:
+
+{{< img src="ManSol_IceWaterMix_(maple_plots).png" >}}
+
+The plane in white color in figure represents zero (melting) temperature in Celsius, such that the zero-level set of $T$ mimics the interface between ice and water fractions which moves in time.
+
+Next figure depicts comparison of the manufactured solution and the numerical one at different time steps. The temperature given in Kelvins. Note that we have made the vertical range of the OGS solutions in the ParaView plots $10^{-1}$ rescaled, to make the comparison feasible. Also, we have tuned the color legend in the ParaView plots such that the ice and water fractions can be visible/identified.
+
+{{< img src="ManSol_IceWaterMix_(comparison).png" >}}
+
+The test is presented in:
+
+{{< data-link "2D_Ice_melting-forming_manuf_solution/ManSol3_IceWaterMix_Scaled.prj" >}},