diff --git a/ProcessLib/HeatConduction/Tests.cmake b/ProcessLib/HeatConduction/Tests.cmake
index 4195d171ec4993a568034730350af021f297dd49..897bc591ada751395f9dee9e07599e5386c4a2f2 100644
--- a/ProcessLib/HeatConduction/Tests.cmake
+++ b/ProcessLib/HeatConduction/Tests.cmake
@@ -166,6 +166,23 @@ AddTest(
REQUIREMENTS NOT OGS_USE_MPI
)
+# test the source term on a subdomain - parallel version
+AddTest(
+ NAME 1D_HeatConduction_dirichlet_SourceTerm_Parallel2
+ PATH Parabolic/T/1D_dirichlet_source-term/2
+ EXECUTABLE ogs
+ EXECUTABLE_ARGS line_1_line_1e2_source_term.prj
+ WRAPPER mpirun
+ WRAPPER_ARGS -np 2
+ REQUIREMENTS OGS_USE_MPI
+ TESTER vtkdiff
+ DIFF_DATA
+ line_1_line_1e2_ts_0_t_0_000000_0.vtu line_1_line_1e2_ts_0_t_0_000000_0.vtu temperature temperature 1e-15 0.0
+ line_1_line_1e2_ts_0_t_0_000000_1.vtu line_1_line_1e2_ts_0_t_0_000000_1.vtu temperature temperature 1e-15 0.0
+ line_1_line_1e2_ts_500_t_39062500_000000_0.vtu line_1_line_1e2_ts_500_t_39062500_000000_0.vtu temperature temperature 1.4e-9 0.0
+ line_1_line_1e2_ts_500_t_39062500_000000_1.vtu line_1_line_1e2_ts_500_t_39062500_000000_1.vtu temperature temperature 1.4e-9 0.0
+)
+
# test the source term on a subdomain with the PETSc embedded executable file
AddTest(
NAME 1D_HeatConduction_dirichlet_SourceTerm_PETSc
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_cell_properties_cfg2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_cell_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..7773aafce42483feac0ad6661ca0ab9859a2976f
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_cell_properties_cfg2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_cell_properties_val2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_cell_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..e2fb55e80edba6100b355483b57a41321d6ebf75
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_cell_properties_val2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_msh_cfg2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_msh_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..9903ca1abe50deb4864cee34e170631c96618d2e
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_msh_cfg2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_msh_ele2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_msh_ele2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..20a04caaf704582c59fdf0b3ea1d84314de85590
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_msh_ele2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_msh_ele_g2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_msh_ele_g2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_msh_nod2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_msh_nod2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..7fe63ffe4138e436b401983188cbdbbb4a207554
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_msh_nod2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_node_properties_cfg2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_node_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..02dba888ca799ac0ebd79efd5517f11823592525
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_node_properties_cfg2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_node_properties_val2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_node_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..51b239409bc815c19b00cb97f62996fb782f24c3
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/SourceTermDomain_partitioned_node_properties_val2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_cell_properties_cfg2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_cell_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..0c3079b923ece3fc90f2c7a6564e96c3c37c83ab
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_cell_properties_cfg2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_cell_properties_val2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_cell_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..bc8840b22b1bcf5594717cc39cf105953e00c24d
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_cell_properties_val2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_msh_cfg2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_msh_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..724ae227a4a24ec8d88bc53c6c558c7afac4329d
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_msh_cfg2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_msh_ele2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_msh_ele2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..34653705c212f15869dae0a5c97dad2e4223862a
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_msh_ele2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_msh_ele_g2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_msh_ele_g2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_msh_nod2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_msh_nod2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..406d5d676d398dfca67b603aa1056d5e5ecf80c6
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_msh_nod2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_node_properties_cfg2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_node_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..e2db701d7737d75f011c619e28288e40b75561fc
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_node_properties_cfg2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_node_properties_val2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_node_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..1b1cb4d44c57c2d7a5122870fa6ac3e62ff7e94e
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_left_partitioned_node_properties_val2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_cell_properties_cfg2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_cell_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..94b67968c900cb2326a54578b95831a36c93fd16
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_cell_properties_cfg2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_cell_properties_val2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_cell_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..e3e90851634f2abe55109932c015129113254e84
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_cell_properties_val2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_msh_cfg2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_msh_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..b3057981a6de8d13768b152200711503c7da04bb
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_msh_cfg2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_msh_ele2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_msh_ele2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..34653705c212f15869dae0a5c97dad2e4223862a
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_msh_ele2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_msh_ele_g2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_msh_ele_g2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..e69de29bb2d1d6434b8b29ae775ad8c2e48c5391
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_msh_nod2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_msh_nod2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..ce7ca9db5516ea79d37cb93a0e46377a8ab33828
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_msh_nod2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_node_properties_cfg2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_node_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..73b6a7585759c84b25ae52ba59df29faae49127c
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_node_properties_cfg2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_node_properties_val2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_node_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..100bd18bdc9e7cdd7a744fa564f5d22b9c629811
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_geometry_right_partitioned_node_properties_val2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_source_term.prj b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_source_term.prj
new file mode 100644
index 0000000000000000000000000000000000000000..e3250a5034eeecb179c9d185eb405d5f850f3e2e
--- /dev/null
+++ b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_source_term.prj
@@ -0,0 +1,158 @@
+<?xml version="1.0" encoding="ISO-8859-1"?>
+<OpenGeoSysProject>
+ <meshes>
+ <mesh>line_1_lines_1e2.vtu</mesh>
+ <mesh>SourceTermDomain.vtu</mesh>
+ <mesh>line_1_geometry_left.vtu</mesh>
+ <mesh>line_1_geometry_right.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_x" output_name="heat_flux_x"/>
+ </secondary_variables>
+ </process>
+ </processes>
+ <time_loop>
+ <processes>
+ <process ref="HeatConduction">
+ <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.0 </t_initial>
+ <t_end> 39062500 </t_end>
+ <timesteps>
+ <pair>
+ <repeat>500</repeat>
+ <delta_t>78125</delta_t>
+ </pair>
+ </timesteps>
+ </time_stepping>
+ </process>
+ </processes>
+ <output>
+ <type>VTK</type>
+ <prefix>line_1_line_1e2</prefix>
+ <timesteps>
+ <pair>
+ <repeat> 1 </repeat>
+ <each_steps> 500 </each_steps>
+ </pair>
+ </timesteps>
+ <variables>
+ <variable> temperature </variable>
+ <variable> heat_flux_x </variable>
+ </variables>
+ <suffix>_ts_{:timestep}_t_{:time}</suffix>
+ </output>
+ </time_loop>
+ <media>
+ <medium id="0">
+ <phases/>
+ <properties>
+ <property>
+ <name>thermal_conductivity</name>
+ <type>Constant</type>
+ <value>5</value>
+ </property>
+ <property>
+ <name>heat_capacity</name>
+ <type>Constant</type>
+ <value>1000</value>
+ </property>
+ <property>
+ <name>density</name>
+ <type>Constant</type>
+ <value>2500</value>
+ </property>
+ <property>
+ <name>reference_temperature</name>
+ <type>Constant</type>
+ <value>293.15</value>
+ </property>
+ </properties>
+ </medium>
+ </media>
+ <parameters>
+ <parameter>
+ <name>T0</name>
+ <type>Constant</type>
+ <value>273.15</value>
+ </parameter>
+ <parameter>
+ <name>T1</name>
+ <type>Constant</type>
+ <value>274.15</value>
+ </parameter>
+ <parameter>
+ <name>st_param</name>
+ <type>Constant</type>
+ <value>5</value>
+ </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>line_1_geometry_left</mesh>
+ <type>Dirichlet</type>
+ <parameter>T1</parameter>
+ </boundary_condition>
+ <boundary_condition>
+ <mesh>line_1_geometry_right</mesh>
+ <type>Dirichlet</type>
+ <parameter>T1</parameter>
+ </boundary_condition>
+ </boundary_conditions>
+ <source_terms>
+ <source_term>
+ <mesh>SourceTermDomain</mesh>
+ <type>Volumetric</type>
+ <parameter>st_param</parameter>
+ </source_term>
+ </source_terms>
+ </process_variable>
+ </process_variables>
+ <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-16 -maxiter 10000</lis>
+ <eigen>
+ <solver_type>CG</solver_type>
+ <precon_type>DIAGONAL</precon_type>
+ <max_iteration_step>10000</max_iteration_step>
+ <error_tolerance>1e-16</error_tolerance>
+ </eigen>
+ <petsc>
+ <prefix>gw</prefix>
+ <parameters>-gw_ksp_type cg -gw_pc_type bjacobi -gw_ksp_rtol 1e-16 -gw_ksp_max_it 10000</parameters>
+ </petsc>
+ </linear_solver>
+ </linear_solvers>
+</OpenGeoSysProject>
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_0_t_0_000000_0.vtu b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_0_t_0_000000_0.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..e9ec30c7f5320180010692b389df1fe2250f14ad
--- /dev/null
+++ b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_0_t_0_000000_0.vtu
@@ -0,0 +1,30 @@
+<?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="21" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="53" NumberOfCells="52" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlux" format="appended" RangeMin="0" RangeMax="0" offset="84" />
+ <DataArray type="Float64" Name="heat_flux_x" format="appended" RangeMin="0" RangeMax="0" offset="148" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="273.15" RangeMax="273.15" offset="212" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="284" />
+ <DataArray type="UInt8" Name="vtkGhostType" format="appended" RangeMin="0" RangeMax="0" offset="344" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0.48" RangeMax="1" offset="404" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="784" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="972" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1164" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAHQAAAAAAAAA=eF4z0zPRM9A1NDK30E03NDU1NTZISTUGADAoBLk=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAADgAAAAAAAAA=eF5jYBgFQwEAAAGoAAE=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAADgAAAAAAAAA=eF5jYBgFQwEAAAGoAAE=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAEwAAAAAAAAA=eF5LSwMCoUKHtFF6UNMAYIuR9g==AQAAAAAAAAAAgAAAAAAAANAAAAAAAAAADAAAAAAAAAA=eF5jYBheAAAA0AABAQAAAAAAAAAAgAAAAAAAADQAAAAAAAAADAAAAAAAAAA=eF5jYCAdAAAANAABAQAAAAAAAAAAgAAAAAAAAPgEAAAAAAAA+gAAAAAAAAA=eF510i1YwlAYhmHzMlSpM5ORqHXkk6HOOjN5Vq2al8lQMdP4E1TYBIpk9p03Ppwv3uW9zva4fPwfu0X35uot4b3R7bBM6B9pJ5q908+xqxfoj6usUTys4G+vdvQqsWH6fb2adtZw/6ycrll6298XvH5sPUy31WhGv7Pnxhv4s81m9M+pHb3lh7fwJ/vMKX1izx3Tmzbb+IYP+nZ0/3tHdP+Zox+4f66j+9mCrvsNdEVXV3R1RVdXu0BXdHVFV1d7+Ivviq6u6OqqDHRFV1d0dVUFuqKrK7q6+gt0RVdXdHV1CHRFV1d0dXUMdEVXV3TdKbA7h18ASgfS7Q==AQAAAAAAAAAAgAAAAAAAAEADAAAAAAAAagAAAAAAAAA=eF5dxTUCwkAAALDDrbgXd4cW+P/bGGAhWRLCV4aznOM8F7jIJS5zhatc44jr3OAmt7jNHe5yj/s84CGPeMwxT3jKM57zgpe84jVveMs73vOBj3ziM1/4yje+84MTTvnJr9/v8O8DnaUKXg==AQAAAAAAAAAAgAAAAAAAAKABAAAAAAAAbQAAAAAAAAA=eF4txdEGwgAAAMAkk2SSSZKZZJIkSZJkJkkmSZIkSdL/f0EP3b1cufRXceCqa647dMNNR2657Y67jp24575TDzz0yGNPPPXMcy+89MprZ8698dY771344KNPPvviq2++++GnX377469/EDoKxQ==AQAAAAAAAAAAgAAAAAAAADQAAAAAAAAADAAAAAAAAAA=eF5jZiYdAAAQWgCd
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_0_t_0_000000_1.vtu b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_0_t_0_000000_1.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..1bc8ddca48b2e060fb1fcde8ece13dd8ca934ace
--- /dev/null
+++ b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_0_t_0_000000_1.vtu
@@ -0,0 +1,30 @@
+<?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="21" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="50" NumberOfCells="49" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlux" format="appended" RangeMin="0" RangeMax="0" offset="84" />
+ <DataArray type="Float64" Name="heat_flux_x" format="appended" RangeMin="0" RangeMax="0" offset="148" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="273.15" RangeMax="273.15" offset="212" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="284" />
+ <DataArray type="UInt8" Name="vtkGhostType" format="appended" RangeMin="0" RangeMax="1" offset="344" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.49" offset="408" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="812" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="992" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1176" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAHQAAAAAAAAA=eF4z0zPRM9A1NDK30E03NDU1NTZISTUGADAoBLk=AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAADgAAAAAAAAA=eF5jYBgFgwkAAAGQAAE=AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAADgAAAAAAAAA=eF5jYBgFgwkAAAGQAAE=AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAAEwAAAAAAAAA=eF5LSwMCoUKHtFF6UNAADQSJsw==AQAAAAAAAAAAgAAAAAAAAMQAAAAAAAAADAAAAAAAAAA=eF5jYBj6AAAAxAABAQAAAAAAAAAAgAAAAAAAADEAAAAAAAAADQAAAAAAAAA=eF5jYCANMAIAADIAAg==AQAAAAAAAAAAgAAAAAAAALAEAAAAAAAADQEAAAAAAAA=eF510itswlAUBmB0NdjOFo3ekMyCrh6WWbCrrsFsAkMNpgnJkpWELGOM7JGWsAeUxwap3mzR67lX/v2P/Nqc5y0U8qNT9GtJ+6qa713wwHR+60GP/O+B31xLDEgeH/yvLl+GJP8tuBH3W6fGHak7Ai/b7mNqjUk/9+Belj3uP5A+J+AVFVPS/xP4NLWyjmZkrmfwhmr0hcz7Cp60pfAb2UMIfinjtiKynzm4jgXZG/q5auid7BNdurTdD7LnT3C1thhd7/8L/ETGMpfkLujNC4kVuVcMLlWLPrq+4xr8KGVTdH3fDfiZKrwld0dX4zo78h6+waNQAl2/kx/wkhp4T94PuiprH8D/Ad6FwlM=AQAAAAAAAAAAgAAAAAAAABADAAAAAAAAZQAAAAAAAAA=eF5dxUkCgQAAAMDsCSFE9qVsWf//OQedZi4TBH81rnODm9ziNnc45C5H3OM+DzjmIY94zAlPeMozTnnOC854ySte84a3vOM9H/jIJ8654DNf+Mo3vnPJD37yi9/84W/1D7vfCWI=AQAAAAAAAAAAgAAAAAAAAIgBAAAAAAAAaQAAAAAAAAA=eF4txdEGwgAAAMDJJEmSJEmSZGaSSZKZzCRJMkmSZP//ET1093KN4C900y233XHXPfc98NAjjz3x1DPPvfDSkWMnXnnt1BtvvfPemXMfXLj00SefffHVN1e+++GnX377469r/wAXYwmTAQAAAAAAAAAAgAAAAAAAADEAAAAAAAAADAAAAAAAAAA=eF5jZiYRAAAOjACU
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_500_t_39062500_000000_0.vtu b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_500_t_39062500_000000_0.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..00bd115c9e63908c044ddcff109d3a58ebe04e46
--- /dev/null
+++ b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_500_t_39062500_000000_0.vtu
@@ -0,0 +1,30 @@
+<?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="21" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="53" NumberOfCells="52" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlux" format="appended" RangeMin="-0.049999999261" RangeMax="0.049999999988" offset="84" />
+ <DataArray type="Float64" Name="heat_flux_x" format="appended" RangeMin="-0.049999999901" RangeMax="0.050000000047" offset="252" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="274.15" RangeMax="274.15495" offset="480" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="964" />
+ <DataArray type="UInt8" Name="vtkGhostType" format="appended" RangeMin="0" RangeMax="0" offset="1024" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0.48" RangeMax="1" offset="1084" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1464" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="1652" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1844" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAHQAAAAAAAAA=eF4z0zPRM9A1NDK30E03NDU1NTZISTUGADAoBLk=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAXAAAAAAAAAA=eF5jYGCYCQT2DAwNQGollAbxQeDCXgh9wJYBBcD4B/aiil+Ail9AE4epQzcPlzoM9YRoqHkP0Pjo7kRXj+F+mL3o/sWhD6f5BOzBpQ45fBwmA+NjP0wdAGHXRJk=AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAAiAAAAAAAAAA=eF510LENwjAUhOFjE8+A6ygexaNkBQIBTGWJiiniPWg8RrpgyO8CC6pPdzo9S5bc3tr7LGlZ12cnZVtyX/IthEcxhU2P+rhTvNKjMF82Ddk11j7WHaYz9xqHxtqr9hP30aAwn9j90aNDYTyyQ48GhWnknfE7RxzQo0E15sPvPr3/a34B2xeh0w==AQAAAAAAAAAAgAAAAAAAAKgBAAAAAAAASgEAAAAAAAA=eF5TlM+pqFIqdPjy5OYaEC0F5edMnsAPovexb11WCaSFam7aguj0D3+vVIDEU5RyQLTYTTdmEJ3vmz2rHEifOthvCKJVzLacKAPStatuxIPoa3J/v5UCaf3Jin0gupPdTRVEP6/O2lMCpF0/9AWD6KUpm18VA2m2m9cbQXS67x8JEH38oMKGIiCtYebqDqK7VmXeKwTSb+X6SkF00ORNPCB6B/v1xQVAWq7mtxWIbvkgfykfSL9JcckE0aE3MxhA9F7f3ul5QFrt0EY9EN1vdu1oLpD+uepXDIhOlpf/kgOkL0x27gbRNhwZSiB6VU3PzmwgLfFxQwCIbk+9+jwLSH+/+bMORKf5yYmB6GuHnNZmAmk383QXEL1jdfftDJD7gd4A0bOnXOEE0bycPxekA+mGWlkLEP35o+P5NJA5IACkK6HxAAAi1qyjAQAAAAAAAAAAgAAAAAAAANAAAAAAAAAADAAAAAAAAAA=eF5jYBheAAAA0AABAQAAAAAAAAAAgAAAAAAAADQAAAAAAAAADAAAAAAAAAA=eF5jYCAdAAAANAABAQAAAAAAAAAAgAAAAAAAAPgEAAAAAAAA+gAAAAAAAAA=eF510i1YwlAYhmHzMlSpM5ORqHXkk6HOOjN5Vq2al8lQMdP4E1TYBIpk9p03Ppwv3uW9zva4fPwfu0X35uot4b3R7bBM6B9pJ5q908+xqxfoj6usUTys4G+vdvQqsWH6fb2adtZw/6ycrll6298XvH5sPUy31WhGv7Pnxhv4s81m9M+pHb3lh7fwJ/vMKX1izx3Tmzbb+IYP+nZ0/3tHdP+Zox+4f66j+9mCrvsNdEVXV3R1RVdXu0BXdHVFV1d7+Ivviq6u6OqqDHRFV1d0dVUFuqKrK7q6+gt0RVdXdHV1CHRFV1d0dXUMdEVXV3TdKbA7h18ASgfS7Q==AQAAAAAAAAAAgAAAAAAAAEADAAAAAAAAagAAAAAAAAA=eF5dxTUCwkAAALDDrbgXd4cW+P/bGGAhWRLCV4aznOM8F7jIJS5zhatc44jr3OAmt7jNHe5yj/s84CGPeMwxT3jKM57zgpe84jVveMs73vOBj3ziM1/4yje+84MTTvnJr9/v8O8DnaUKXg==AQAAAAAAAAAAgAAAAAAAAKABAAAAAAAAbQAAAAAAAAA=eF4txdEGwgAAAMAkk2SSSZKZZJIkSZJkJkkmSZIkSdL/f0EP3b1cufRXceCqa647dMNNR2657Y67jp24575TDzz0yGNPPPXMcy+89MprZ8698dY771344KNPPvviq2++++GnX377469/EDoKxQ==AQAAAAAAAAAAgAAAAAAAADQAAAAAAAAADAAAAAAAAAA=eF5jZiYdAAAQWgCd
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_500_t_39062500_000000_1.vtu b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_500_t_39062500_000000_1.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..133c90c3c7b1ca9ac88490459b9776fd0c39a498
--- /dev/null
+++ b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_line_1e2_ts_500_t_39062500_000000_1.vtu
@@ -0,0 +1,30 @@
+<?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="21" format="appended" RangeMin="45" RangeMax="103" offset="0" />
+ </FieldData>
+ <Piece NumberOfPoints="50" NumberOfCells="49" >
+ <PointData>
+ <DataArray type="Float64" Name="HeatFlux" format="appended" RangeMin="-0.049999999232" RangeMax="0.024999999965" offset="84" />
+ <DataArray type="Float64" Name="heat_flux_x" format="appended" RangeMin="-0.049999999916" RangeMax="-0.037499999889" offset="244" />
+ <DataArray type="Float64" Name="temperature" format="appended" RangeMin="274.15" RangeMax="274.1549" offset="444" />
+ </PointData>
+ <CellData>
+ <DataArray type="Int32" Name="MaterialIDs" format="appended" RangeMin="0" RangeMax="0" offset="912" />
+ <DataArray type="UInt8" Name="vtkGhostType" format="appended" RangeMin="0" RangeMax="1" offset="972" />
+ </CellData>
+ <Points>
+ <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="appended" RangeMin="0" RangeMax="0.49" offset="1036" />
+ </Points>
+ <Cells>
+ <DataArray type="Int64" Name="connectivity" format="appended" RangeMin="" RangeMax="" offset="1440" />
+ <DataArray type="Int64" Name="offsets" format="appended" RangeMin="" RangeMax="" offset="1620" />
+ <DataArray type="UInt8" Name="types" format="appended" RangeMin="" RangeMax="" offset="1804" />
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+ <AppendedData encoding="base64">
+ _AQAAAAAAAAAAgAAAAAAAABUAAAAAAAAAHQAAAAAAAAA=eF4z0zPRM9A1NDK30E03NDU1NTZISTUGADAoBLk=AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAAVQAAAAAAAAA=eF5jYGCYPHPmyv0MYHBhLwNWcMAWjY+mDqYPJg5TD6Nh8hdg4mj6cZl3Ad08HPrQ6QdQ+gYufWjiMBrd/7jcfQNNHAYw/I/DfegALD4TCOwBtac3VQ==AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAAdAAAAAAAAAA=eF51zcENgzAQRNFJJ9QADVCKO2FrICGBmyVOVIFLcSdEgr+KbIXT08xqtJKNy7LtUkI9LzM5olV6nyt9H06PYyjyz+R/UK9S7xuyYfQ7hqnUUG/2lRH1+W+PAZub3jDfqJnd6UMRE31Gz33bdev+BQ63tnQ=AQAAAAAAAAAAgAAAAAAAAJABAAAAAAAAPgEAAAAAAAA=eF5LSwMCpUKHbx8dz4Po9lpZi3QgLcH5cwGIXjHlCmcGkLZQ2FAEoo+v7r4NosPM010ygfTzQ05rQXSFn5xYFpDmuPWzDkTPSr36HETrftwQkA2kD9b07ATRwRwZSjkgfZOdu0F0jbz8FxAtuPpXTC6QXm527SiItjm0US8PSF/07Z0OojNuZjDkA2mGVJdMED3jg/wlEG1c89uqAEifZ7++GERnTd7EUwik2eT7SkH00lWZ90C0s5mrexGQfnQQ6BEg3ez7R6IYSKvevN4Ioo+nbH4FojM/9AWXAGnemqw9IHoDu5tqKZAOmqzYB6K/yf39BqJnrroRXwak7cy2nADRTw72G5YD6U7f7FkgWv+mG3MFkL6aopQDoms//L0CotVqbtpWgt27dRmIrpw8gb8KSCvK51SAaACU1KHdAQAAAAAAAAAAgAAAAAAAAMQAAAAAAAAADAAAAAAAAAA=eF5jYBj6AAAAxAABAQAAAAAAAAAAgAAAAAAAADEAAAAAAAAADQAAAAAAAAA=eF5jYCANMAIAADIAAg==AQAAAAAAAAAAgAAAAAAAALAEAAAAAAAADQEAAAAAAAA=eF510itswlAUBmB0NdjOFo3ekMyCrh6WWbCrrsFsAkMNpgnJkpWELGOM7JGWsAeUxwap3mzR67lX/v2P/Nqc5y0U8qNT9GtJ+6qa713wwHR+60GP/O+B31xLDEgeH/yvLl+GJP8tuBH3W6fGHak7Ai/b7mNqjUk/9+Belj3uP5A+J+AVFVPS/xP4NLWyjmZkrmfwhmr0hcz7Cp60pfAb2UMIfinjtiKynzm4jgXZG/q5auid7BNdurTdD7LnT3C1thhd7/8L/ETGMpfkLujNC4kVuVcMLlWLPrq+4xr8KGVTdH3fDfiZKrwld0dX4zo78h6+waNQAl2/kx/wkhp4T94PuiprH8D/Ad6FwlM=AQAAAAAAAAAAgAAAAAAAABADAAAAAAAAZQAAAAAAAAA=eF5dxUkCgQAAAMDsCSFE9qVsWf//OQedZi4TBH81rnODm9ziNnc45C5H3OM+DzjmIY94zAlPeMozTnnOC854ySte84a3vOM9H/jIJ8654DNf+Mo3vnPJD37yi9/84W/1D7vfCWI=AQAAAAAAAAAAgAAAAAAAAIgBAAAAAAAAaQAAAAAAAAA=eF4txdEGwgAAAMDJJEmSJEmSZGaSSZKZzCRJMkmSZP//ET1093KN4C900y233XHXPfc98NAjjz3x1DPPvfDSkWMnXnnt1BtvvfPemXMfXLj00SefffHVN1e+++GnX377469r/wAXYwmTAQAAAAAAAAAAgAAAAAAAADEAAAAAAAAADAAAAAAAAAA=eF5jZiYRAAAOjACU
+ </AppendedData>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_cell_properties_cfg2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_cell_properties_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..f277ddc3f3b0c71f875381b15cdf51c81292e2ce
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_cell_properties_cfg2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_cell_properties_val2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_cell_properties_val2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..1b5cd21f22d7492ba501fc9c8d97705ddd79b69f
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_cell_properties_val2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_cfg2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_cfg2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..b9a46a8d7a5cf8b23da627a3e967add852bb98cb
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_cfg2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_ele2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_ele2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..17f8d48e8a6912ef8b775c662fbc4c7cc584b587
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_ele2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_ele_g2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_ele_g2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..2b33c30e84f158903bbceed486a10cca158be188
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_ele_g2.bin differ
diff --git a/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_nod2.bin b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_nod2.bin
new file mode 100644
index 0000000000000000000000000000000000000000..8dbf409e28bd18db43d2339c3501af176b5fbd99
Binary files /dev/null and b/Tests/Data/Parabolic/T/1D_dirichlet_source-term/2/line_1_lines_1e2_partitioned_msh_nod2.bin differ