diff --git a/ProcessLib/HT/Tests.cmake b/ProcessLib/HT/Tests.cmake
index 2dcd18844fa10401f0e30160a2258be972bcd199..1834925d80797f47dc328017241d320b55ec705c 100644
--- a/ProcessLib/HT/Tests.cmake
+++ b/ProcessLib/HT/Tests.cmake
@@ -244,3 +244,37 @@ AddTest(
CoupledPressureParabolicTemperatureParabolic_ts_10_expected.vtu CoupledPressureParabolicTemperatureParabolic_pcs_0_ts_10_t_1.000000.vtu p p 1e-10 1e-16
CoupledPressureParabolicTemperatureParabolic_ts_10_expected.vtu CoupledPressureParabolicTemperatureParabolic_pcs_0_ts_10_t_1.000000.vtu darcy_velocity darcy_velocity 1e-10 1e-16
)
+
+AddTest(
+ NAME HT_Balance
+ PATH Parabolic/HT/SimpleSynthetics
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS balance_ht_cube_1e3.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT OGS_USE_MPI
+ DIFF_DATA
+ flux_1e3_t_0.000000_expected.vtu flux_1e3_t_0.000000.vtu flux flux 1e-10 1e-16
+ flux_1e3_t_0.000010_expected.vtu flux_1e3_t_0.000010.vtu flux flux 1e-10 1e-16
+ flux_1e3_t_0.001010_expected.vtu flux_1e3_t_0.001010.vtu flux flux 1e-10 1e-16
+ flux_1e3_t_0.101010_expected.vtu flux_1e3_t_0.101010.vtu flux flux 1e-10 1e-16
+ flux_1e3_t_1.101010_expected.vtu flux_1e3_t_1.101010.vtu flux flux 1e-10 1e-16
+ flux_1e3_t_10.000000_expected.vtu flux_1e3_t_10.000000.vtu flux flux 1e-10 1e-16
+)
+
+AddTest(
+ NAME LARGE_HT_Balance
+ PATH Parabolic/HT/SimpleSynthetics
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS balance_ht_cube_1e4.prj
+ WRAPPER time
+ TESTER vtkdiff
+ REQUIREMENTS NOT OGS_USE_MPI
+ DIFF_DATA
+ flux_1e4_t_0.000000_expected.vtu flux_1e4_t_0.000000.vtu flux flux 1e-10 1e-16
+ flux_1e4_t_0.000010_expected.vtu flux_1e4_t_0.000010.vtu flux flux 1e-10 1e-16
+ flux_1e4_t_0.001010_expected.vtu flux_1e4_t_0.001010.vtu flux flux 1e-10 1e-16
+ flux_1e4_t_0.101010_expected.vtu flux_1e4_t_0.101010.vtu flux flux 1e-10 1e-16
+ flux_1e4_t_1.101010_expected.vtu flux_1e4_t_1.101010.vtu flux flux 1e-10 1e-16
+ flux_1e4_t_10.000000_expected.vtu flux_1e4_t_10.000000.vtu flux flux 1e-10 1e-16
+)
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/balance_ht_cube_1e3.prj b/Tests/Data/Parabolic/HT/SimpleSynthetics/balance_ht_cube_1e3.prj
new file mode 100644
index 0000000000000000000000000000000000000000..9cf77be7cd49f350c119087940581bf5c5d0e0bc
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/balance_ht_cube_1e3.prj
@@ -0,0 +1,271 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>cube_1x1x1_hex_1e3.vtu</mesh>
+ <geometry>cube_1x1x1.gml</geometry>
+ <processes>
+ <process>
+ <name>SimpleThermalConvection</name>
+ <type>HT</type>
+ <integration_order>2</integration_order>
+ <process_variables>
+ <temperature>T</temperature>
+ <pressure>p</pressure>
+ </process_variables>
+ <fluid>
+ <density>
+ <type>TemperatureDependent</type>
+ <rho0>1000</rho0>
+ <temperature0>20</temperature0>
+ <beta>4.3e-4</beta>
+ </density>
+ <viscosity>
+ <type>Constant</type>
+ <value>1.0e-3</value>
+ </viscosity>
+ <specific_heat_capacity>
+ <type>Constant</type>
+ <value>4200</value>
+ </specific_heat_capacity>
+ </fluid>
+ <porous_medium>
+ <porous_medium id="0">
+ <permeability>
+ <permeability_tensor_entries>kappa1</permeability_tensor_entries>
+ <type>Constant</type>
+ </permeability>
+ <porosity>
+ <type>Constant</type>
+ <porosity_parameter>constant_porosity_parameter</porosity_parameter>
+ </porosity>
+ <storage>
+ <type>Constant</type>
+ <value>0.0</value>
+ </storage>
+ </porous_medium>
+ </porous_medium>
+ <density_solid>rho_solid</density_solid>
+ <specific_heat_capacity_solid>c_solid</specific_heat_capacity_solid>
+ <thermal_conductivity_solid>lambda_solid</thermal_conductivity_solid>
+ <thermal_conductivity_fluid>lambda_fluid</thermal_conductivity_fluid>
+ <thermal_dispersivity>
+ <longitudinal>alpha_l</longitudinal>
+ <transversal>alpha_t</transversal>
+ </thermal_dispersivity>
+ <specific_body_force>0 0 0</specific_body_force>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="darcy_velocity" output_name="darcy_velocity"/>
+ </secondary_variables>
+ <calculatesurfaceflux>
+ <mesh>cube_1x1x1_hex_1e3_complete_surface.vtu</mesh>
+ <property_name>flux</property_name>
+ <output_mesh>flux_1e3.vtu</output_mesh>
+ </calculatesurfaceflux>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="SimpleThermalConvection">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1.e-3</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <output>
+ <variables>
+ <variable>T</variable>
+ <variable>p</variable>
+ <variable>darcy_velocity</variable>
+ </variables>
+ </output>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial> 0.0 </t_initial>
+ <t_end> 10 </t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1e-7</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1e-5</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1e-3</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1e-1</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>10</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>ConstViscosityThermalConvection</prefix>
+ <timesteps>
+ <pair>
+ <repeat> 10 </repeat>
+ <each_steps> 1 </each_steps>
+ </pair>
+ </timesteps>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>rho_solid</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <parameter>
+ <name>rho_fluid</name>
+ <type>Constant</type>
+ <value>1000</value>
+ </parameter>
+ <parameter>
+ <name>lambda_solid</name>
+ <type>Constant</type>
+ <value>3.0</value>
+ </parameter>
+ <parameter>
+ <name>lambda_fluid</name>
+ <type>Constant</type>
+ <value>0.65</value>
+ </parameter>
+ <parameter>
+ <name>alpha_l</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <parameter>
+ <name>alpha_t</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <parameter>
+ <name>c_solid</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>T0</name>
+ <type>Constant</type>
+ <value>20</value>
+ </parameter>
+ <parameter>
+ <name>P0</name>
+ <type>Constant</type>
+ <value>2.99e7</value>
+ </parameter>
+ <parameter>
+ <name>p_Dirichlet_left</name>
+ <type>Constant</type>
+ <value>3e7</value>
+ </parameter>
+ <parameter>
+ <name>p_Dirichlet_right</name>
+ <type>Constant</type>
+ <value>2.99e7</value>
+ </parameter>
+ <parameter>
+ <name>t_Dirichlet_left</name>
+ <type>Constant</type>
+ <value>20</value>
+ </parameter>
+ <parameter>
+ <name>t_Dirichlet_right</name>
+ <type>Constant</type>
+ <value>10</value>
+ </parameter>
+ <parameter>
+ <name>constant_porosity_parameter</name>
+ <type>Constant</type>
+ <value>0.001</value>
+ </parameter>
+ <parameter>
+ <name>kappa1</name>
+ <type>Constant</type>
+ <values>1.e-14 0 0 0 1.e-14 0 0 0 1.e-14</values>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>T</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>T0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <parameter>t_Dirichlet_left</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <parameter>t_Dirichlet_right</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>p</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>P0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <parameter>p_Dirichlet_left</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <parameter>p_Dirichlet_right</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_picard</name>
+ <type>Picard</type>
+ <max_iter>1000</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <lis>-i bicgstab -p jacobi -tol 1e-12 -maxiter 10000</lis>
+ <petsc>
+ <prefix>T</prefix>
+ <parameters> -T_ksp_type bcgs -T_pc_type bjacobi -T_ksp_rtol 1e-16 -T_ksp_max_it 4000</parameters>
+ </petsc>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-12</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/balance_ht_cube_1e4.prj b/Tests/Data/Parabolic/HT/SimpleSynthetics/balance_ht_cube_1e4.prj
new file mode 100644
index 0000000000000000000000000000000000000000..dcc3175bdb87af4a7465a3bb5d9042afb5b68a48
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/balance_ht_cube_1e4.prj
@@ -0,0 +1,271 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <mesh>cube_1x1x1_hex_1e4.vtu</mesh>
+ <geometry>cube_1x1x1.gml</geometry>
+ <processes>
+ <process>
+ <name>SimpleThermalConvection</name>
+ <type>HT</type>
+ <integration_order>2</integration_order>
+ <process_variables>
+ <temperature>T</temperature>
+ <pressure>p</pressure>
+ </process_variables>
+ <fluid>
+ <density>
+ <type>TemperatureDependent</type>
+ <rho0>1000</rho0>
+ <temperature0>20</temperature0>
+ <beta>4.3e-4</beta>
+ </density>
+ <viscosity>
+ <type>Constant</type>
+ <value>1.0e-3</value>
+ </viscosity>
+ <specific_heat_capacity>
+ <type>Constant</type>
+ <value>4200</value>
+ </specific_heat_capacity>
+ </fluid>
+ <porous_medium>
+ <porous_medium id="0">
+ <permeability>
+ <permeability_tensor_entries>kappa1</permeability_tensor_entries>
+ <type>Constant</type>
+ </permeability>
+ <porosity>
+ <type>Constant</type>
+ <porosity_parameter>constant_porosity_parameter</porosity_parameter>
+ </porosity>
+ <storage>
+ <type>Constant</type>
+ <value>0.0</value>
+ </storage>
+ </porous_medium>
+ </porous_medium>
+ <density_solid>rho_solid</density_solid>
+ <specific_heat_capacity_solid>c_solid</specific_heat_capacity_solid>
+ <thermal_conductivity_solid>lambda_solid</thermal_conductivity_solid>
+ <thermal_conductivity_fluid>lambda_fluid</thermal_conductivity_fluid>
+ <thermal_dispersivity>
+ <longitudinal>alpha_l</longitudinal>
+ <transversal>alpha_t</transversal>
+ </thermal_dispersivity>
+ <specific_body_force>0 0 0</specific_body_force>
+ <secondary_variables>
+ <secondary_variable type="static" internal_name="darcy_velocity" output_name="darcy_velocity"/>
+ </secondary_variables>
+ <calculatesurfaceflux>
+ <mesh>cube_1x1x1_hex_1e4_complete_surface.vtu</mesh>
+ <property_name>flux</property_name>
+ <output_mesh>flux_1e4.vtu</output_mesh>
+ </calculatesurfaceflux>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="SimpleThermalConvection">
+ <nonlinear_solver>basic_picard</nonlinear_solver>
+ <convergence_criterion>
+ <type>DeltaX</type>
+ <norm_type>NORM2</norm_type>
+ <abstol>1.e-3</abstol>
+ </convergence_criterion>
+ <time_discretization>
+ <type>BackwardEuler</type>
+ </time_discretization>
+ <output>
+ <variables>
+ <variable>T</variable>
+ <variable>p</variable>
+ <variable>darcy_velocity</variable>
+ </variables>
+ </output>
+ <time_stepping>
+ <type>FixedTimeStepping</type>
+ <t_initial> 0.0 </t_initial>
+ <t_end> 10 </t_end>
+ <timesteps>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1e-7</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1e-5</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1e-3</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1e-1</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>1</delta_t>
+ </pair>
+ <pair>
+ <repeat>1</repeat>
+ <delta_t>10</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>ConstViscosityThermalConvection</prefix>
+ <timesteps>
+ <pair>
+ <repeat> 10 </repeat>
+ <each_steps> 1 </each_steps>
+ </pair>
+ </timesteps>
+ </output>
+ </time_loop>
+ <parameters>
+ <parameter>
+ <name>rho_solid</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <parameter>
+ <name>rho_fluid</name>
+ <type>Constant</type>
+ <value>1000</value>
+ </parameter>
+ <parameter>
+ <name>lambda_solid</name>
+ <type>Constant</type>
+ <value>3.0</value>
+ </parameter>
+ <parameter>
+ <name>lambda_fluid</name>
+ <type>Constant</type>
+ <value>0.65</value>
+ </parameter>
+ <parameter>
+ <name>alpha_l</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <parameter>
+ <name>alpha_t</name>
+ <type>Constant</type>
+ <value>0.0</value>
+ </parameter>
+ <parameter>
+ <name>c_solid</name>
+ <type>Constant</type>
+ <value>0</value>
+ </parameter>
+ <parameter>
+ <name>T0</name>
+ <type>Constant</type>
+ <value>20</value>
+ </parameter>
+ <parameter>
+ <name>P0</name>
+ <type>Constant</type>
+ <value>2.99e7</value>
+ </parameter>
+ <parameter>
+ <name>p_Dirichlet_left</name>
+ <type>Constant</type>
+ <value>3e7</value>
+ </parameter>
+ <parameter>
+ <name>p_Dirichlet_right</name>
+ <type>Constant</type>
+ <value>2.99e7</value>
+ </parameter>
+ <parameter>
+ <name>t_Dirichlet_left</name>
+ <type>Constant</type>
+ <value>20</value>
+ </parameter>
+ <parameter>
+ <name>t_Dirichlet_right</name>
+ <type>Constant</type>
+ <value>10</value>
+ </parameter>
+ <parameter>
+ <name>constant_porosity_parameter</name>
+ <type>Constant</type>
+ <value>0.001</value>
+ </parameter>
+ <parameter>
+ <name>kappa1</name>
+ <type>Constant</type>
+ <values>1.e-14 0 0 0 1.e-14 0 0 0 1.e-14</values>
+ </parameter>
+ </parameters>
+ <process_variables>
+ <process_variable>
+ <name>T</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>T0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <parameter>t_Dirichlet_left</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <parameter>t_Dirichlet_right</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ <process_variable>
+ <name>p</name>
+ <components>1</components>
+ <order>1</order>
+ <initial_condition>P0</initial_condition>
+ <boundary_conditions>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>left</geometry>
+ <type>Dirichlet</type>
+ <parameter>p_Dirichlet_left</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <geometrical_set>cube_1x1x1_geometry</geometrical_set>
+ <geometry>right</geometry>
+ <type>Dirichlet</type>
+ <parameter>p_Dirichlet_right</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ </process_variable>
+ </process_variables>
+ <nonlinear_solvers>
+ <nonlinear_solver>
+ <name>basic_picard</name>
+ <type>Picard</type>
+ <max_iter>1000</max_iter>
+ <linear_solver>general_linear_solver</linear_solver>
+ </nonlinear_solver>
+ </nonlinear_solvers>
+ <linear_solvers>
+ <linear_solver>
+ <name>general_linear_solver</name>
+ <lis>-i bicgstab -p jacobi -tol 1e-12 -maxiter 10000</lis>
+ <petsc>
+ <prefix>T</prefix>
+ <parameters> -T_ksp_type bcgs -T_pc_type bjacobi -T_ksp_rtol 1e-16 -T_ksp_max_it 4000</parameters>
+ </petsc>
+ <eigen>
+ <solver_type>BiCGSTAB</solver_type>
+ <precon_type>ILUT</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-12</error_tolerance>
+ </eigen>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1.gml b/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1.gml
new file mode 100644
index 0000000000000000000000000000000000000000..e9c428f64b8c241cb7cb158c371d33ded17803d8
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1.gml
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:6fe4e4940b7d080889c63573d45036e0e10c9edc265b18c041e2a153c872e3b6
+size 1689
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e3.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e3.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..500d9f815a5726e303aa4f8c220220c9bb517c7d
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e3.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:f3e77b980dc0a4cefc19ba02eca4ad07ded3dc97497a01f169c719fc738a6bc9
+size 23197
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e3_complete_surface.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e3_complete_surface.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..e4ea771a0a2a6a76ac3039d26a375d858648d326
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e3_complete_surface.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:27c8e6e020a45ab15dd370c841ac7342573805646ff76aedd0c66398a5338295
+size 19396
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e4.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e4.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..af2f291bbbdc13093f328e6ad3f1883ec1ab5056
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e4.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:234ccbf944fd4b45fa4c00a8f4e2b453a452bd308b8a7fb24c85575f6d012ff4
+size 253962
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e4_complete_surface.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e4_complete_surface.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..3ece58fa5a8c08fbc42745518cd0a4d54eda8e9d
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/cube_1x1x1_hex_1e4_complete_surface.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:d95c0af23a1b5304b980c5030f71e3f072d05477e89cd6c8109151f46938e645
+size 346076
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.000000_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.000000_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..b286e8f7a75e4e3942f9720f66c433e6bb5e2923
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.000000_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:9aac0c1ed26ab48425d4a711e46dfdecd86773c3243d2419a7317e091582c6b6
+size 92355
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.000010_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.000010_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..6d562c4301a7aff1ed54da82981040e002f19c8e
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.000010_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:bcfd2ee73234feec6f451a9cbfba14ed4c812d1aca14c0dd3afcc0229e6b55e1
+size 92355
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.001010_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.001010_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..0dd61188b34f443047cafb1d4cccd51dd7d2e2b1
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.001010_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:3548f243709a6509360e09b190db5e7390b2e724e459685c25349df797fab755
+size 92355
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.101010_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.101010_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..c1d133d81ebe6a787752d7bc954780dac1207803
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_0.101010_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:0e2c0817047e18278ef2b257520bb5ae54126e1740aa70f80ab517dd3adf07f8
+size 92355
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_1.101010_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_1.101010_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..b464da689450142130f2ca9d18d8eea3b92c7404
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_1.101010_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:12ff850d259f4dcb2c2710db899906fa5ad82bc1adba4b121d5b8b7fa413a2ce
+size 92355
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_10.000000_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_10.000000_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..d7e6fe969f42a23c4abb837fe94d353bd40a5552
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e3_t_10.000000_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:2f853726b6aa195c0205049f5082363d1d01c270e80badcc772a2cd0bddb460f
+size 92355
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.000000_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.000000_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..583df7e8ce9e5e3824901206c76848ba7b5937fc
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.000000_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:1ab62ef860e28f2a6bf802cce20fa9a51195fd9ea81595961f142e8497a16ddc
+size 377211
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.000010_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.000010_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..d4ac15597025fbfe29fa879fa730755c623618f1
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.000010_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:dd214808868a38242181838b3415121a1671801076848d2d5abef89871da1e8f
+size 377212
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.001010_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.001010_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..9ecd56d416109674504dd7f9365e083947c80535
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.001010_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:18a134b4466a80827f5758ddc36c4aaee275dda4a6b68bc4e405384a92c03c1a
+size 377212
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.101010_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.101010_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..48e909f010548b91c276f832573a841abe38256d
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_0.101010_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:b212f7adec5b069ca07e14f978ee98ba5757bfcad085eca97980b939fa354a31
+size 377212
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_1.101010_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_1.101010_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..17a1492ea6ae63be330efa1109d87eb1105432be
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_1.101010_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:c2db5e3b1290821952e56dbd118222bcf7426e6b6f6ace2155fd5f1a45455738
+size 377212
diff --git a/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_10.000000_expected.vtu b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_10.000000_expected.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..0f7e5b86a6125541cd405db2e79e19ab753de336
--- /dev/null
+++ b/Tests/Data/Parabolic/HT/SimpleSynthetics/flux_1e4_t_10.000000_expected.vtu
@@ -0,0 +1,3 @@
+version https://git-lfs.github.com/spec/v1
+oid sha256:48b1460f8f619d8ce5131f7d5519ed7cfe22d4a6a04b894636d085d83275c115
+size 377212