diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/FunctionParameterTest_XDMF.prj b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/FunctionParameterTest_XDMF.prj
new file mode 100644
index 0000000000000000000000000000000000000000..929e7f6369b3f990e83112fcf4a8fff9c9d195bb
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/FunctionParameterTest_XDMF.prj
@@ -0,0 +1,186 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<!-- Almost same as FunctionParameterTest - changed XDMF output and
+boundary meshes will be written as well -->
+<OpenGeoSysProject>
+ <meshes>
+ <mesh>square_5x5_tris_32.vtu</mesh>
+ <mesh>square_5x5_tris_32_left_boundary.vtu</mesh>
+ <mesh>square_5x5_tris_32_right_boundary.vtu</mesh>
+ </meshes>
+ <processes>
+ <process>
+ <name>LiquidFlow</name>
+ <type>LIQUID_FLOW</type>
+ <integration_order>2</integration_order>
+ <darcy_gravity>
+ <axis_id>1</axis_id>
+ <g>0</g>
+ </darcy_gravity>
+ <process_variables>
+ <process_variable>pressure</process_variable>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="darcy_velocity" output_name="v"/>
+ </secondary_variables>
+ </process>
+ </processes>
+ <media>
+ <medium id="0">
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value> 0.0011373 </value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value> 1000 </value>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>3.2439e-12</value>
+ </property>
+ <property>
+ <name>reference_temperature</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.17</value>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value> 0.0 </value>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <parameters>
+ <parameter>
+ <name>p_0</name>
+ <type>Function</type>
+ <expression>z</expression>
+ </parameter>
+ <parameter>
+ <name>p_polyline_left</name>
+ <type>Function</type>
+ <expression>x+1</expression>
+ </parameter>
+ <parameter>
+ <name>p_polyline_right</name>
+ <type>Function</type>
+ <expression>x-5</expression>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>pressure</name>
+ <order>1</order>
+ <components>1</components>
+ <initial_condition>p_0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <type>Dirichlet</type>
+ <mesh>square_5x5_tris_32_left_boundary</mesh>
+ <parameter>p_polyline_left</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <type>Dirichlet</type>
+ <mesh>square_5x5_tris_32_right_boundary</mesh>
+ <parameter>p_polyline_right</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <time_loop>
+ <processes>
+ <process ref="LiquidFlow">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1.e-6</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial> 0 </t_initial>
+ <t_end> 1 </t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>XDMF</type>
+ <prefix>{:meshname}</prefix>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ <meshes>
+ <mesh>square_5x5_tris_32</mesh>
+ <mesh>square_5x5_tris_32_left_boundary</mesh>
+ <mesh>square_5x5_tris_32_right_boundary</mesh>
+ </meshes>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ </timesteps>
+ <output_iteration_results>false</output_iteration_results>
+ <variables>
+ <variable> pressure </variable>
+ <variable> v </variable>
+ </variables>
+ </output>
+ </time_loop>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_picard</name>
+ <type>Picard</type>
+ <max_iter>10</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <lis>-i cg -p jacobi -tol 1e-15 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>DIAGONAL</precon_type>
+ <max_iteration_step>1000</max_iteration_step>
+ <error_tolerance>1e-16</error_tolerance>
+ </eigen>
+ <petsc>
+ <prefix>lf</prefix>
+ <parameters>-lf_ksp_type cg -lf_pc_type bjacobi -lf_ksp_rtol 1e-16 -lf_ksp_max_it 10000</parameters>
+ </petsc>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <regex>square_5x5_tris_32_ts_1_t_1.000000.vtu</regex>
+ <field>pressure</field>
+ <absolute_tolerance>1e-15</absolute_tolerance>
+ <relative_tolerance>1e-14</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/FunctionParameterTest_XDMF.prj b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/FunctionParameterTest_XDMF.prj
new file mode 100644
index 0000000000000000000000000000000000000000..929e7f6369b3f990e83112fcf4a8fff9c9d195bb
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/FunctionParameterTest_XDMF.prj
@@ -0,0 +1,186 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<!-- Almost same as FunctionParameterTest - changed XDMF output and
+boundary meshes will be written as well -->
+<OpenGeoSysProject>
+ <meshes>
+ <mesh>square_5x5_tris_32.vtu</mesh>
+ <mesh>square_5x5_tris_32_left_boundary.vtu</mesh>
+ <mesh>square_5x5_tris_32_right_boundary.vtu</mesh>
+ </meshes>
+ <processes>
+ <process>
+ <name>LiquidFlow</name>
+ <type>LIQUID_FLOW</type>
+ <integration_order>2</integration_order>
+ <darcy_gravity>
+ <axis_id>1</axis_id>
+ <g>0</g>
+ </darcy_gravity>
+ <process_variables>
+ <process_variable>pressure</process_variable>
+ </process_variables>
+ <secondary_variables>
+ <secondary_variable internal_name="darcy_velocity" output_name="v"/>
+ </secondary_variables>
+ </process>
+ </processes>
+ <media>
+ <medium id="0">
+ <phases>
+ <phase>
+ <type>AqueousLiquid</type>
+ <properties>
+ <property>
+ <name>viscosity</name>
+ <type>Constant</type>
+ <value> 0.0011373 </value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value> 1000 </value>
+ </property>
+ </properties>
+ </phase>
+ </phases>
+ <properties>
+ <property>
+ <name>permeability</name>
+ <type>Constant</type>
+ <value>3.2439e-12</value>
+ </property>
+ <property>
+ <name>reference_temperature</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </property>
+ <property>
+ <name>porosity</name>
+ <type>Constant</type>
+ <value>0.17</value>
+ </property>
+ <property>
+ <name>storage</name>
+ <type>Constant</type>
+ <value> 0.0 </value>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <parameters>
+ <parameter>
+ <name>p_0</name>
+ <type>Function</type>
+ <expression>z</expression>
+ </parameter>
+ <parameter>
+ <name>p_polyline_left</name>
+ <type>Function</type>
+ <expression>x+1</expression>
+ </parameter>
+ <parameter>
+ <name>p_polyline_right</name>
+ <type>Function</type>
+ <expression>x-5</expression>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>pressure</name>
+ <order>1</order>
+ <components>1</components>
+ <initial_condition>p_0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <type>Dirichlet</type>
+ <mesh>square_5x5_tris_32_left_boundary</mesh>
+ <parameter>p_polyline_left</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <type>Dirichlet</type>
+ <mesh>square_5x5_tris_32_right_boundary</mesh>
+ <parameter>p_polyline_right</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <time_loop>
+ <processes>
+ <process ref="LiquidFlow">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1.e-6</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial> 0 </t_initial>
+ <t_end> 1 </t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>XDMF</type>
+ <prefix>{:meshname}</prefix>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ <meshes>
+ <mesh>square_5x5_tris_32</mesh>
+ <mesh>square_5x5_tris_32_left_boundary</mesh>
+ <mesh>square_5x5_tris_32_right_boundary</mesh>
+ </meshes>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <each_steps>1</each_steps>
+ </pair>
+ </timesteps>
+ <output_iteration_results>false</output_iteration_results>
+ <variables>
+ <variable> pressure </variable>
+ <variable> v </variable>
+ </variables>
+ </output>
+ </time_loop>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_picard</name>
+ <type>Picard</type>
+ <max_iter>10</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <lis>-i cg -p jacobi -tol 1e-15 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>DIAGONAL</precon_type>
+ <max_iteration_step>1000</max_iteration_step>
+ <error_tolerance>1e-16</error_tolerance>
+ </eigen>
+ <petsc>
+ <prefix>lf</prefix>
+ <parameters>-lf_ksp_type cg -lf_pc_type bjacobi -lf_ksp_rtol 1e-16 -lf_ksp_max_it 10000</parameters>
+ </petsc>
+ </linear_solver>
+ </linear_solvers>
+ <test_definition>
+ <vtkdiff>
+ <regex>square_5x5_tris_32_ts_1_t_1.000000.vtu</regex>
+ <field>pressure</field>
+ <absolute_tolerance>1e-15</absolute_tolerance>
+ <relative_tolerance>1e-14</relative_tolerance>
+ </vtkdiff>
+ </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32.vtu b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..af7bdb6b4b1f69d6f3aa8d7328b55cc447297416
Binary files /dev/null and b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32.vtu differ
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary.vtu b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..52b20854c8fa372bdc09ddcf39b91e323c1982ae
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary.vtu
@@ -0,0 +1,51 @@
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="5" NumberOfCells="4">
+ <PointData Scalars="bulk_node_ids">
+ <DataArray type="UInt64" IdType="1" Name="bulk_node_ids" format="ascii" RangeMin="0" RangeMax="20">
+ 5 0 10 15 20
+ </DataArray>
+ </PointData>
+ <CellData Scalars="bulk_element_ids">
+ <DataArray type="UInt64" IdType="1" Name="bulk_element_ids" format="ascii" RangeMin="0" RangeMax="24">
+ 0 8 16 24
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="ascii" RangeMin="0" RangeMax="5">
+ 0 1.25 3 0 0 3
+ 0 2.5 3 0 3.75 3
+ 0 5 3
+ <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+ <Value index="0">
+ 0
+ </Value>
+ <Value index="1">
+ 5
+ </Value>
+ </InformationKey>
+ <InformationKey name="L2_NORM_FINITE_RANGE" location="vtkDataArray" length="2">
+ <Value index="0">
+ 0
+ </Value>
+ <Value index="1">
+ 5
+ </Value>
+ </InformationKey>
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="ascii" RangeMin="0" RangeMax="4">
+ 0 1 2 0 3 2
+ 4 3
+ </DataArray>
+ <DataArray type="Int64" Name="offsets" format="ascii" RangeMin="2" RangeMax="8">
+ 2 4 6 8
+ </DataArray>
+ <DataArray type="UInt8" Name="types" format="ascii" RangeMin="3" RangeMax="3">
+ 3 3 3 3
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_0_t_0.000000.h5 b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_0_t_0.000000.h5
new file mode 100644
index 0000000000000000000000000000000000000000..82daee09f5d9bab4d4ef6cf1c24d6067e7457f8a
Binary files /dev/null and b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_0_t_0.000000.h5 differ
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_0_t_0.000000.xdmf b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_0_t_0.000000.xdmf
new file mode 100644
index 0000000000000000000000000000000000000000..4b083795314d3b95df55654400ae6c373befc54d
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_0_t_0.000000.xdmf
@@ -0,0 +1,19 @@
+<?xml version="1.0" encoding="utf-8"?>
+<Xdmf xmlns:xi="http://www.w3.org/2001/XInclude" Version="3.0">
+ <Domain>
+ <Grid Name="Grid">
+ <Geometry Origin="" Type="XYZ">
+ <DataItem DataType="Float" Dimensions="5 3" Format="XML" Precision="8">0 1.25 3 0 0 3 0 2.5 3 0 3.75 3 0 5 3</DataItem>
+ </Geometry>
+ <Topology Dimensions="4" Type="Mixed">
+ <DataItem DataType="Int" Dimensions="16" Format="XML" Precision="8">2 2 0 1 2 2 2 0 2 2 3 2 2 2 4 3</DataItem>
+ </Topology>
+ <Attribute Center="Grid" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="OGS_VERSION" Type="None">
+ <DataItem DataType="Char" Dimensions="41" Format="XML" Precision="1">54 46 51 46 50 45 52 50 56 45 103 57 98 49 100 97 48 98 56 57 46 100 105 114 116 121 46 50 48 50 48 49 48 49 50 48 55 52 56 52 52</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="pressure" Type="None">
+ <DataItem DataType="Float" Dimensions="5" Format="XML" Precision="8">3 3 3 3 3</DataItem>
+ </Attribute>
+ </Grid>
+ </Domain>
+</Xdmf>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_1_t_1.000000.h5 b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_1_t_1.000000.h5
new file mode 100644
index 0000000000000000000000000000000000000000..82daee09f5d9bab4d4ef6cf1c24d6067e7457f8a
Binary files /dev/null and b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_1_t_1.000000.h5 differ
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_1_t_1.000000.xdmf b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_1_t_1.000000.xdmf
new file mode 100644
index 0000000000000000000000000000000000000000..16cb9288e16d27dea51a00203ea718a4add28889
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_left_boundary_ts_1_t_1.000000.xdmf
@@ -0,0 +1,19 @@
+<?xml version="1.0" encoding="utf-8"?>
+<Xdmf xmlns:xi="http://www.w3.org/2001/XInclude" Version="3.0">
+ <Domain>
+ <Grid Name="Grid">
+ <Geometry Origin="" Type="XYZ">
+ <DataItem DataType="Float" Dimensions="5 3" Format="XML" Precision="8">0 1.25 3 0 0 3 0 2.5 3 0 3.75 3 0 5 3</DataItem>
+ </Geometry>
+ <Topology Dimensions="4" Type="Mixed">
+ <DataItem DataType="Int" Dimensions="16" Format="XML" Precision="8">2 2 0 1 2 2 2 0 2 2 3 2 2 2 4 3</DataItem>
+ </Topology>
+ <Attribute Center="Grid" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="OGS_VERSION" Type="None">
+ <DataItem DataType="Char" Dimensions="41" Format="XML" Precision="1">54 46 51 46 50 45 52 50 56 45 103 57 98 49 100 97 48 98 56 57 46 100 105 114 116 121 46 50 48 50 48 49 48 49 50 48 55 52 56 52 52</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="pressure" Type="None">
+ <DataItem DataType="Float" Dimensions="5" Format="XML" Precision="8">1 1 1 1 1</DataItem>
+ </Attribute>
+ </Grid>
+ </Domain>
+</Xdmf>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary.vtu b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..8df915d4a239106d1144e779cc796a10dd270e87
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary.vtu
@@ -0,0 +1,51 @@
+<VTKFile type="UnstructuredGrid" version="1.0" byte_order="LittleEndian" header_type="UInt64">
+ <UnstructuredGrid>
+ <Piece NumberOfPoints="5" NumberOfCells="4">
+ <PointData Scalars="bulk_node_ids">
+ <DataArray type="UInt64" IdType="1" Name="bulk_node_ids" format="ascii" RangeMin="4" RangeMax="24">
+ 4 9 14 19 24
+ </DataArray>
+ </PointData>
+ <CellData Scalars="bulk_element_ids">
+ <DataArray type="UInt64" IdType="1" Name="bulk_element_ids" format="ascii" RangeMin="7" RangeMax="31">
+ 7 15 23 31
+ </DataArray>
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="ascii" RangeMin="5" RangeMax="7.0710678118654755">
+ 5 0 3 5 1.25 3
+ 5 2.5 3 5 3.75 3
+ 5 5 3
+ <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+ <Value index="0">
+ 5
+ </Value>
+ <Value index="1">
+ 7.0710678119
+ </Value>
+ </InformationKey>
+ <InformationKey name="L2_NORM_FINITE_RANGE" location="vtkDataArray" length="2">
+ <Value index="0">
+ 5
+ </Value>
+ <Value index="1">
+ 7.0710678119
+ </Value>
+ </InformationKey>
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="ascii" RangeMin="0" RangeMax="4">
+ 0 1 1 2 2 3
+ 3 4
+ </DataArray>
+ <DataArray type="Int64" Name="offsets" format="ascii" RangeMin="2" RangeMax="8">
+ 2 4 6 8
+ </DataArray>
+ <DataArray type="UInt8" Name="types" format="ascii" RangeMin="3" RangeMax="3">
+ 3 3 3 3
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_0_t_0.000000.h5 b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_0_t_0.000000.h5
new file mode 100644
index 0000000000000000000000000000000000000000..82daee09f5d9bab4d4ef6cf1c24d6067e7457f8a
Binary files /dev/null and b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_0_t_0.000000.h5 differ
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_0_t_0.000000.xdmf b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_0_t_0.000000.xdmf
new file mode 100644
index 0000000000000000000000000000000000000000..19e83b7ee362229dcdc347096cf372490bbc6eaf
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_0_t_0.000000.xdmf
@@ -0,0 +1,19 @@
+<?xml version="1.0" encoding="utf-8"?>
+<Xdmf xmlns:xi="http://www.w3.org/2001/XInclude" Version="3.0">
+ <Domain>
+ <Grid Name="Grid">
+ <Geometry Origin="" Type="XYZ">
+ <DataItem DataType="Float" Dimensions="5 3" Format="XML" Precision="8">5 0 3 5 1.25 3 5 2.5 3 5 3.75 3 5 5 3</DataItem>
+ </Geometry>
+ <Topology Dimensions="4" Type="Mixed">
+ <DataItem DataType="Int" Dimensions="16" Format="XML" Precision="8">2 2 0 1 2 2 1 2 2 2 2 3 2 2 3 4</DataItem>
+ </Topology>
+ <Attribute Center="Grid" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="OGS_VERSION" Type="None">
+ <DataItem DataType="Char" Dimensions="41" Format="XML" Precision="1">54 46 51 46 50 45 52 50 56 45 103 57 98 49 100 97 48 98 56 57 46 100 105 114 116 121 46 50 48 50 48 49 48 49 50 48 55 52 56 52 52</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="pressure" Type="None">
+ <DataItem DataType="Float" Dimensions="5" Format="XML" Precision="8">3 3 3 3 3</DataItem>
+ </Attribute>
+ </Grid>
+ </Domain>
+</Xdmf>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_1_t_1.000000.h5 b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_1_t_1.000000.h5
new file mode 100644
index 0000000000000000000000000000000000000000..82daee09f5d9bab4d4ef6cf1c24d6067e7457f8a
Binary files /dev/null and b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_1_t_1.000000.h5 differ
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_1_t_1.000000.xdmf b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_1_t_1.000000.xdmf
new file mode 100644
index 0000000000000000000000000000000000000000..76fd582bb3203c3963ddc6972b8a0b4efc2864a6
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_right_boundary_ts_1_t_1.000000.xdmf
@@ -0,0 +1,19 @@
+<?xml version="1.0" encoding="utf-8"?>
+<Xdmf xmlns:xi="http://www.w3.org/2001/XInclude" Version="3.0">
+ <Domain>
+ <Grid Name="Grid">
+ <Geometry Origin="" Type="XYZ">
+ <DataItem DataType="Float" Dimensions="5 3" Format="XML" Precision="8">5 0 3 5 1.25 3 5 2.5 3 5 3.75 3 5 5 3</DataItem>
+ </Geometry>
+ <Topology Dimensions="4" Type="Mixed">
+ <DataItem DataType="Int" Dimensions="16" Format="XML" Precision="8">2 2 0 1 2 2 1 2 2 2 2 3 2 2 3 4</DataItem>
+ </Topology>
+ <Attribute Center="Grid" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="OGS_VERSION" Type="None">
+ <DataItem DataType="Char" Dimensions="41" Format="XML" Precision="1">54 46 51 46 50 45 52 50 56 45 103 57 98 49 100 97 48 98 56 57 46 100 105 114 116 121 46 50 48 50 48 49 48 49 50 48 55 52 56 52 52</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="pressure" Type="None">
+ <DataItem DataType="Float" Dimensions="5" Format="XML" Precision="8">0 0 0 0 0</DataItem>
+ </Attribute>
+ </Grid>
+ </Domain>
+</Xdmf>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_ts_0_t_0.000000.xdmf b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_ts_0_t_0.000000.xdmf
new file mode 100644
index 0000000000000000000000000000000000000000..489226749d6631690c1168a744567424c84ac5b6
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_ts_0_t_0.000000.xdmf
@@ -0,0 +1,28 @@
+<?xml version="1.0" encoding="utf-8"?>
+<Xdmf xmlns:xi="http://www.w3.org/2001/XInclude" Version="3.0">
+ <Domain>
+ <Grid Name="Grid">
+ <Geometry Origin="" Type="XYZ">
+ <DataItem DataType="Float" Dimensions="25 3" Format="XML" Precision="8">0 0 3 1.25 0 3 2.5 0 3 3.75 0 3 5 0 3 0 1.25 3 1.25 1.25 3 2.5 1.25 3 3.75 1.25 3 5 1.25 3 0 2.5 3 1.25 2.5 3 2.5 2.5 3 3.75 2.5 3 5 2.5 3 0 3.75 3 1.25 3.75 3 2.5 3.75 3 3.75 3.75 3 5 3.75 3 0 5 3 1.25 5 3 2.5 5 3 3.75 5 3 5 5 3</DataItem>
+ </Geometry>
+ <Topology Dimensions="32" Type="Mixed">
+ <DataItem DataType="Int" Dimensions="128" Format="HDF" Precision="8">square_5x5_tris_32_ts_0_t_0.000000.h5:Data7</DataItem>
+ </Topology>
+ <Attribute Center="Grid" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="OGS_VERSION" Type="None">
+ <DataItem DataType="Char" Dimensions="41" Format="XML" Precision="1">54 46 51 46 50 45 53 56 54 45 103 56 102 56 53 56 53 54 49 97 46 100 105 114 116 121 46 50 48 50 48 49 48 50 54 49 51 50 52 51 53</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="HydraulicFlow" Type="None">
+ <DataItem DataType="Float" Dimensions="25" Format="XML" Precision="8">0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="MaterialIDs" Type="None">
+ <DataItem DataType="Int" Dimensions="25" Format="XML" Precision="4">0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="pressure" Type="None">
+ <DataItem DataType="Float" Dimensions="25" Format="XML" Precision="8">3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3 3</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="v" Type="None">
+ <DataItem DataType="Float" Dimensions="25 2" Format="XML" Precision="8">0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0</DataItem>
+ </Attribute>
+ </Grid>
+ </Domain>
+</Xdmf>
diff --git a/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_ts_1_t_1.000000.xdmf b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_ts_1_t_1.000000.xdmf
new file mode 100644
index 0000000000000000000000000000000000000000..e9a508ac22ff4ffbca83868dde7b34bf96e2a8ca
--- /dev/null
+++ b/Tests/Data/Parabolic/LiquidFlow/SimpleSynthetics/XDMF/square_5x5_tris_32_ts_1_t_1.000000.xdmf
@@ -0,0 +1,28 @@
+<?xml version="1.0" encoding="utf-8"?>
+<Xdmf xmlns:xi="http://www.w3.org/2001/XInclude" Version="3.0">
+ <Domain>
+ <Grid Name="Grid">
+ <Geometry Origin="" Type="XYZ">
+ <DataItem DataType="Float" Dimensions="25 3" Format="XML" Precision="8">0 0 3 1.25 0 3 2.5 0 3 3.75 0 3 5 0 3 0 1.25 3 1.25 1.25 3 2.5 1.25 3 3.75 1.25 3 5 1.25 3 0 2.5 3 1.25 2.5 3 2.5 2.5 3 3.75 2.5 3 5 2.5 3 0 3.75 3 1.25 3.75 3 2.5 3.75 3 3.75 3.75 3 5 3.75 3 0 5 3 1.25 5 3 2.5 5 3 3.75 5 3 5 5 3</DataItem>
+ </Geometry>
+ <Topology Dimensions="32" Type="Mixed">
+ <DataItem DataType="Int" Dimensions="128" Format="HDF" Precision="8">square_5x5_tris_32_ts_1_t_1.000000.h5:Data7</DataItem>
+ </Topology>
+ <Attribute Center="Grid" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="OGS_VERSION" Type="None">
+ <DataItem DataType="Char" Dimensions="41" Format="XML" Precision="1">54 46 51 46 50 45 53 56 54 45 103 56 102 56 53 56 53 54 49 97 46 100 105 114 116 121 46 50 48 50 48 49 48 50 54 49 51 50 52 51 53</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="HydraulicFlow" Type="None">
+ <DataItem DataType="Float" Dimensions="25" Format="XML" Precision="8">-2.8522817198628322e-09 2.8522817198628317e-09 0 4.2784225797942485e-09 -4.2784225797942485e-09 -5.7045634397256652e-09 5.7045634397256635e-09 0 8.5568451595884969e-09 -8.5568451595884969e-09 -5.7045634397256652e-09 5.7045634397256635e-09 0 8.5568451595884969e-09 -8.5568451595884969e-09 -5.7045634397256652e-09 5.7045634397256635e-09 0 8.5568451595884969e-09 -8.5568451595884969e-09 -2.8522817198628322e-09 2.8522817198628317e-09 0 4.2784225797942485e-09 -4.2784225797942485e-09</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="MaterialIDs" Type="None">
+ <DataItem DataType="Int" Dimensions="25" Format="XML" Precision="4">0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="pressure" Type="None">
+ <DataItem DataType="Float" Dimensions="25" Format="XML" Precision="8">1 0.75000000000000011 0.50000000000000067 0.25000000000000044 0 1 0.75000000000000033 0.50000000000000067 0.25000000000000044 0 1 0.75 0.50000000000000067 0.2500000000000005 0 1 0.75000000000000033 0.50000000000000078 0.2500000000000005 0 1 0.75000000000000011 0.50000000000000078 0.2500000000000005 0</DataItem>
+ </Attribute>
+ <Attribute Center="Node" ElementCell="" ElementDegree="0" ElementFamily="" ItemType="" Name="v" Type="None">
+ <DataItem DataType="Float" Dimensions="25 2" Format="XML" Precision="8">5.7045634397256606e-10 -1.7756826847496619e-25 5.7045634397256585e-10 -2.3675769129995494e-25 5.7045634397256668e-10 0 5.7045634397256761e-10 0 5.7045634397256771e-10 0 5.7045634397256616e-10 2.3675769129995494e-25 5.7045634397256596e-10 1.1837884564997749e-25 5.7045634397256637e-10 -8.8784134237483108e-26 5.7045634397256751e-10 -5.9189422824988735e-26 5.7045634397256771e-10 0 5.7045634397256637e-10 -2.3675769129995494e-25 5.7045634397256585e-10 -2.9594711412494396e-26 5.7045634397256616e-10 5.9189422824988746e-26 5.7045634397256751e-10 -5.9189422824988735e-26 5.7045634397256792e-10 -5.9189422824988746e-26 5.7045634397256596e-10 1.1837884564997747e-25 5.7045634397256596e-10 -1.1837884564997742e-25 5.7045634397256637e-10 -1.7756826847496622e-25 5.7045634397256761e-10 -2.9594711412494373e-26 5.7045634397256792e-10 0 5.7045634397256637e-10 0 5.7045634397256585e-10 2.3675769129995494e-25 5.7045634397256596e-10 1.1837884564997749e-25 5.7045634397256751e-10 0 5.7045634397256792e-10 0</DataItem>
+ </Attribute>
+ </Grid>
+ </Domain>
+</Xdmf>