diff --git a/ProcessLib/ComponentTransport/Tests.cmake b/ProcessLib/ComponentTransport/Tests.cmake
index 080b3baa30d60b4e72dd29c8e64f293fccdd9e70..ee4b9fc8bf2f929669b4396228078bf40bd6dda0 100644
--- a/ProcessLib/ComponentTransport/Tests.cmake
+++ b/ProcessLib/ComponentTransport/Tests.cmake
@@ -929,3 +929,4 @@ endif()
 
 OgsTest(PROJECTFILE Parabolic/ComponentTransport/MassFlux/only_grad_c.prj RUNTIME 1)
 OgsTest(PROJECTFILE Parabolic/ComponentTransport/MassFlux/grad_c_and_grad_p.prj RUNTIME 3)
+OgsTest(PROJECTFILE Parabolic/ComponentTransport/MassFlux/only_grad_p.prj RUNTIME 1)
diff --git a/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p.prj b/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p.prj
new file mode 100644
index 0000000000000000000000000000000000000000..c8bda3d5fbdda0583f79b843570ceed83e6a3614
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p.prj
@@ -0,0 +1,326 @@
+<?xml version='1.0' encoding='ISO-8859-1'?>
+<OpenGeoSysProject>
+    <meshes>
+        <mesh>only_grad_p.vtu</mesh>
+        <mesh>only_grad_p_left.vtu</mesh>
+        <mesh>only_grad_p_right.vtu</mesh>
+    </meshes>
+    <processes>
+        <process>
+            <name>HC</name>
+            <type>ComponentTransport</type>
+            <integration_order>4</integration_order>
+            <process_variables>
+                <concentration>C</concentration>
+                <pressure>pressure</pressure>
+            </process_variables>
+            <secondary_variables>
+                <secondary_variable internal_name="darcy_velocity" output_name="darcy_velocity"/>
+                <secondary_variable internal_name="CFlux" output_name="CFlux"/>
+            </secondary_variables>
+            <specific_body_force>0 0 0</specific_body_force>
+        </process>
+    </processes>
+    <media>
+        <medium id="0">
+            <phases>
+                <phase>
+                    <type>AqueousLiquid</type>
+                    <components>
+                        <component>
+                            <name>C</name>
+                            <properties>
+                                <property>
+                                    <name>pore_diffusion</name>
+                                    <type>Constant</type>
+                                    <value>1.e-9</value>
+                                </property>
+                                <property>
+                                    <name>retardation_factor</name>
+                                    <type>Constant</type>
+                                    <value>1.0</value>
+                                </property>
+                                <property>
+                                    <name>decay_rate</name>
+                                    <type>Parameter</type>
+                                    <parameter_name>decay</parameter_name>
+                                </property>
+                            </properties>
+                        </component>
+                    </components>
+                    <properties>
+                        <property>
+                            <name>density</name>
+                            <type>Constant</type>
+                            <value>1</value>
+                        </property>
+                        <property>
+                            <name>viscosity</name>
+                            <type>Constant</type>
+                            <value>1.0</value>
+                        </property>
+                    </properties>
+                </phase>
+            </phases>
+            <properties>
+                <property>
+                    <name>permeability</name>
+                    <type>Constant</type>
+                    <value>1.e-9</value>
+                </property>
+                <property>
+                    <name>porosity</name>
+                    <type>Constant</type>
+                    <value>0.15</value>
+                </property>
+                <property>
+                    <name>longitudinal_dispersivity</name>
+                    <type>Constant</type>
+                    <value>0</value>
+                </property>
+                <property>
+                    <name>transversal_dispersivity</name>
+                    <type>Constant</type>
+                    <value>0.0</value>
+                </property>
+            </properties>
+        </medium>
+        <medium id="1">
+            <phases>
+                <phase>
+                    <type>AqueousLiquid</type>
+                    <components>
+                        <component>
+                            <name>C</name>
+                            <properties>
+                                <property>
+                                    <name>pore_diffusion</name>
+                                    <type>Constant</type>
+                                    <value>1.e-9</value>
+                                </property>
+                                <property>
+                                    <name>retardation_factor</name>
+                                    <type>Constant</type>
+                                    <value>1.0</value>
+                                </property>
+                                <property>
+                                    <name>decay_rate</name>
+                                    <type>Parameter</type>
+                                    <parameter_name>decay</parameter_name>
+                                </property>
+                            </properties>
+                        </component>
+                    </components>
+                    <properties>
+                        <property>
+                            <name>density</name>
+                            <type>Constant</type>
+                            <value>1</value>
+                        </property>
+                        <property>
+                            <name>viscosity</name>
+                            <type>Constant</type>
+                            <value>1.0</value>
+                        </property>
+                    </properties>
+                </phase>
+            </phases>
+            <properties>
+                <property>
+                    <name>permeability</name>
+                    <type>Constant</type>
+                    <value>1.e-9</value>
+                </property>
+                <property>
+                    <name>porosity</name>
+                    <type>Constant</type>
+                    <value>0.65</value>
+                </property>
+                <property>
+                    <name>longitudinal_dispersivity</name>
+                    <type>Constant</type>
+                    <value>0</value>
+                </property>
+                <property>
+                    <name>transversal_dispersivity</name>
+                    <type>Constant</type>
+                    <value>0.0</value>
+                </property>
+            </properties>
+        </medium>
+    </media>
+    <parameters>
+        <parameter>
+            <name>decay</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <parameter>
+            <name>C_right</name>
+            <type>Constant</type>
+            <value>1</value>
+        </parameter>
+        <parameter>
+            <name>C_left</name>
+            <type>Constant</type>
+            <value>1</value>
+        </parameter>
+        <parameter>
+            <name>C_ini</name>
+            <type>MeshNode</type>
+            <field_name>C_ini</field_name>
+        </parameter>
+        <parameter>
+            <name>p_left</name>
+            <type>Constant</type>
+            <value>0.3</value>
+        </parameter>
+        <parameter>
+            <name>p_right</name>
+            <type>Constant</type>
+            <value>0</value>
+        </parameter>
+        <parameter>
+            <name>p_ini</name>
+            <type>MeshNode</type>
+            <field_name>p_ini</field_name>
+        </parameter>
+    </parameters>
+    <process_variables>
+        <process_variable>
+            <name>pressure</name>
+            <components>1</components>
+            <order>1</order>
+            <initial_condition>p_ini</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <mesh>only_grad_p_left</mesh>
+                    <type>Dirichlet</type>
+                    <parameter>p_left</parameter>
+                </boundary_condition>
+                <boundary_condition>
+                    <mesh>only_grad_p_right</mesh>
+                    <type>Dirichlet</type>
+                    <parameter>p_right</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+        </process_variable>
+        <process_variable>
+            <name>C</name>
+            <components>1</components>
+            <order>1</order>
+            <initial_condition>C_ini</initial_condition>
+            <boundary_conditions>
+                <boundary_condition>
+                    <mesh>only_grad_p_right</mesh>
+                    <type>Dirichlet</type>
+                    <parameter>C_right</parameter>
+                </boundary_condition>
+                <boundary_condition>
+                    <mesh>only_grad_p_left</mesh>
+                    <type>Dirichlet</type>
+                    <parameter>C_left</parameter>
+                </boundary_condition>
+            </boundary_conditions>
+        </process_variable>
+    </process_variables>
+    <time_loop>
+        <processes>
+            <process ref="HC">
+                <nonlinear_solver>basic_picard</nonlinear_solver>
+                <convergence_criterion>
+                    <type>PerComponentDeltaX</type>
+                    <norm_type>INFINITY_N</norm_type>
+                    <!--<reltols>1e-12 1e-10</reltols>-->
+                    <abstols>1e-15 1e-15</abstols>
+                </convergence_criterion>
+                <time_discretization>
+                    <type>BackwardEuler</type>
+                </time_discretization>
+                <time_stepping>
+                    <type>FixedTimeStepping</type>
+                    <t_initial> 0.0 </t_initial>
+                    <t_end> 4e8 </t_end>
+                    <timesteps>
+                        <pair>
+                            <repeat>4</repeat>
+                            <delta_t>1e8</delta_t>
+                        </pair>
+                    </timesteps>
+                </time_stepping>
+            </process>
+        </processes>
+        <output>
+            <type>VTK</type>
+            <prefix>only_grad_p</prefix>
+            <timesteps>
+                <pair>
+                    <repeat>1000</repeat>
+                    <each_steps>1</each_steps>
+                </pair>
+            </timesteps>
+            <variables>
+            </variables>
+            <output_extrapolation_residuals>true</output_extrapolation_residuals>
+        </output>
+    </time_loop>
+    <nonlinear_solvers>
+        <nonlinear_solver>
+            <name>basic_picard</name>
+            <type>Picard</type>
+            <max_iter>50</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>
+            <petsc>
+                <prefix>hc</prefix>
+                <parameters>-hc_ksp_type bcgs -hc_pc_type bjacobi -hc_ksp_rtol 1e-8 -hc_ksp_max_it 20000</parameters>
+            </petsc>
+        </linear_solver>
+    </linear_solvers>
+    <test_definition>
+        <vtkdiff>
+            <file>only_grad_p_ts_4_t_400000000.000000.vtu</file>
+            <field>pressure</field>
+            <absolute_tolerance>4.5e-9</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <file>only_grad_p_ts_4_t_400000000.000000.vtu</file>
+            <field>C</field>
+            <absolute_tolerance>6.3e-15</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <file>only_grad_p_ts_4_t_400000000.000000.vtu</file>
+            <field>CFlux</field>
+            <absolute_tolerance>5.6e-17</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <file>only_grad_p_ts_4_t_400000000.000000.vtu</file>
+            <field>CFlux_residual</field>
+            <absolute_tolerance>1.7e-17</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <file>only_grad_p_ts_4_t_400000000.000000.vtu</file>
+            <field>darcy_velocity</field>
+            <absolute_tolerance>5.6e-17</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+        <vtkdiff>
+            <file>only_grad_p_ts_4_t_400000000.000000.vtu</file>
+            <field>darcy_velocity_residual</field>
+            <absolute_tolerance>1.7e-17</absolute_tolerance>
+            <relative_tolerance>0</relative_tolerance>
+        </vtkdiff>
+    </test_definition>
+</OpenGeoSysProject>
diff --git a/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p.vtu b/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..24a5a7f0de2d580f9bf1843f1c5316fb19bd8541
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p.vtu
@@ -0,0 +1,47 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="44" NumberOfCells="10">
+      <PointData Scalars="C_ini">
+        <DataArray type="Float64" Name="C_ini" format="binary" RangeMin="1" RangeMax="1">
+          AQAAAACAAABgAQAAEgAAAA==eJxjYACBD/YMozRNaABeMTQV
+        </DataArray>
+        <DataArray type="Float64" Name="p_ini" format="binary" RangeMin="-4.470348341500951e-9" RangeMax="0.3">
+          AQAAAACAAABgAQAAXAAAAA==eJwzNgaBy/ZAQnSd+0X79LS0d4E7ztmnpaVtcn94CkSX2XIdB9F6xsaH7c+dOfMncMc+EJ1jy7XdfvbMmZJBO9bZmwMNCNwxz/7///+GQOY+41FzweYCAJipnj0=
+        </DataArray>
+        <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="0" RangeMax="43">
+          AQAAAACAAABgAQAAVQAAAA==eJwtxUkCgQAAAMAoFJEtJZGiiP//z8HMZYLgb+KpQ0eeee6FYydeeuXUa2+ceeud9z746NwnFy59duWLa199c+O7W3d++Oneg19+e/THX/8AvN8Dsw==
+        </DataArray>
+      </PointData>
+      <CellData Scalars="MaterialIDs">
+        <DataArray type="Int32" Name="MaterialIDs" format="binary" RangeMin="0" RangeMax="1">
+          AQAAAACAAAAoAAAADgAAAA==eJxjYMAEjFgwAABkAAY=
+        </DataArray>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="0.8079604065183404">
+          AQAAAACAAAAgBAAArAAAAA==eJx90rEJwzAUhGHXmSKLmGSFrOMqjYdwlyFSZADVrtSkERYIjMgSqR48H/+zG+MPccLHDUP43Ov0vrEn8u3xWYPzmfy1LN8gp5Bf518N8hv5c7zswb2d3OUd/tu+wRO59QDnM7n1ADmF3HqA/EZuPcC9nTzq46SXRK57cOczue7B5RRy3YPLb+S6B3dvJ4/+/+SdyHUP2oO67kF7UNc9aA/qugft4eB/jXzcvQ==
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              0
+            </Value>
+            <Value index="1">
+              0.80796040652
+            </Value>
+          </InformationKey>
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="43">
+          AQAAAACAAACAAgAAbgAAAA==eJxdyzcCggAQAEGMmEHBgGJWDPz/fxZeNdtsNUnyrxOfxmfxTXwbv8ZvuC5ujtvh7rgeboGrcA9cH5fh9rgnboDLcQdcgxvilrga98KluBXuiHvjRrgCd8J9cGNciTvjvrgJbo274Nr4D5lpBrk=
+        </DataArray>
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="8" RangeMax="80">
+          AQAAAACAAABQAAAAIwAAAA==eJzjYIAAASgtAaUVoLQGlDaA0hZQ2gFKe0DpACgNADdQAbk=
+        </DataArray>
+        <DataArray type="UInt8" Name="types" format="binary" RangeMin="11" RangeMax="11">
+          AQAAAACAAAAKAAAACwAAAA==eJzj5oYBAAJnAG8=
+        </DataArray>
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p_left.vtu b/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p_left.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..5b9149032f9cd2c67818d39d0aeab4bbe45c1393
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p_left.vtu
@@ -0,0 +1,53 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="4" NumberOfCells="1">
+      <PointData Scalars="C_ini">
+        <DataArray type="Float64" Name="C_ini" format="binary" RangeMin="1" RangeMax="1">
+          AQAAAACAAAAgAAAAEAAAAA==eJxjYACBD/YMOGgAQWwEvQ==
+        </DataArray>
+        <DataArray type="Float64" Name="p_ini" format="binary" RangeMin="0.3" RangeMax="0.3">
+          AQAAAACAAAAgAAAAEAAAAA==eJwzNgaBy/bGOGgAjsAJEQ==
+        </DataArray>
+        <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="0" RangeMax="33">
+          AQAAAACAAAAgAAAAEwAAAA==eJxjYIAAMSjNDaUVoTQAA+gAQw==
+        </DataArray>
+        <DataArray type="Int64" IdType="1" Name="vtkOriginalPointIds" format="binary" RangeMin="0" RangeMax="6">
+          AQAAAACAAAAgAAAAEwAAAA==eJxjYIAAZijNAqXZoDQAANgADg==
+        </DataArray>
+      </PointData>
+      <CellData Scalars="MaterialIDs">
+        <DataArray type="Int32" Name="MaterialIDs" format="binary" RangeMin="0" RangeMax="0">
+          AQAAAACAAAAEAAAADAAAAA==eJxjYGBgAAAABAAB
+        </DataArray>
+        <DataArray type="Int64" IdType="1" Name="vtkOriginalCellIds" format="binary" RangeMin="2" RangeMax="2">
+          AQAAAACAAAAIAAAACwAAAA==eJxjYoAAAAAYAAM=
+        </DataArray>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="0.1131370824610367">
+          AQAAAACAAABgAAAAFgAAAA==eJxjYCAaODys2mKPh48uDqYBAKMKOQ==
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              0
+            </Value>
+            <Value index="1">
+              0.11313708246
+            </Value>
+          </InformationKey>
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="3">
+          AQAAAACAAAAgAAAAEwAAAA==eJxjYIAARijNDKWZoDQAAHgABw==
+        </DataArray>
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="4" RangeMax="4">
+          AQAAAACAAAAIAAAACwAAAA==eJxjYYAAAAAoAAU=
+        </DataArray>
+        <DataArray type="UInt8" Name="types" format="binary" RangeMin="9" RangeMax="9">
+          AQAAAACAAAABAAAACQAAAA==eJzjBAAACgAK
+        </DataArray>
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p_right.vtu b/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p_right.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..f8296dd5aa4d678abd8e2a84509fc8f42ba680b7
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p_right.vtu
@@ -0,0 +1,53 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="4" NumberOfCells="1">
+      <PointData Scalars="C_ini">
+        <DataArray type="Float64" Name="C_ini" format="binary" RangeMin="1" RangeMax="1">
+          AQAAAACAAAAgAAAAEAAAAA==eJxjYACBD/YMOGgAQWwEvQ==
+        </DataArray>
+        <DataArray type="Float64" Name="p_ini" format="binary" RangeMin="-4.470348341500951e-9" RangeMax="-4.470348341500951e-9">
+          AQAAAACAAAAgAAAAEwAAAA==eJz7//+/obGx8b7/OGgAOsMSFQ==
+        </DataArray>
+        <DataArray type="UInt64" Name="bulk_node_ids" format="binary" RangeMin="10" RangeMax="43">
+          AQAAAACAAAAgAAAAEwAAAA==eJzjYoAABSgtCqW1oTQABwgAaw==
+        </DataArray>
+        <DataArray type="Int64" IdType="1" Name="vtkOriginalPointIds" format="binary" RangeMin="31" RangeMax="34">
+          AQAAAACAAAAgAAAAEwAAAA==eJyTZ4AABSitCKWVoDQACiAAgw==
+        </DataArray>
+      </PointData>
+      <CellData Scalars="MaterialIDs">
+        <DataArray type="Int32" Name="MaterialIDs" format="binary" RangeMin="1" RangeMax="1">
+          AQAAAACAAAAEAAAADAAAAA==eJxjZGBgAAAACAAC
+        </DataArray>
+        <DataArray type="Int64" IdType="1" Name="vtkOriginalCellIds" format="binary" RangeMin="21" RangeMax="21">
+          AQAAAACAAAAIAAAACwAAAA==eJwTZYAAAACwABY=
+        </DataArray>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0.800000011920929" RangeMax="0.8079604065183404">
+          AQAAAACAAABgAAAAIwAAAA==eJxjYGBYMHPmS3sGTIAu7vCwaos9kjiMj64eJg6mAZFhFiE=
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              0.80000001192
+            </Value>
+            <Value index="1">
+              0.80796040652
+            </Value>
+          </InformationKey>
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="3">
+          AQAAAACAAAAgAAAAEwAAAA==eJxjYIAAJijNDKUZoTQAAIgABw==
+        </DataArray>
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="4" RangeMax="4">
+          AQAAAACAAAAIAAAACwAAAA==eJxjYYAAAAAoAAU=
+        </DataArray>
+        <DataArray type="UInt8" Name="types" format="binary" RangeMin="9" RangeMax="9">
+          AQAAAACAAAABAAAACQAAAA==eJzjBAAACgAK
+        </DataArray>
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>
diff --git a/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p_ts_4_t_400000000.000000.vtu b/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p_ts_4_t_400000000.000000.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..d5b077f66cbf6ea9f5206faacc000dab2cb27132
--- /dev/null
+++ b/Tests/Data/Parabolic/ComponentTransport/MassFlux/only_grad_p_ts_4_t_400000000.000000.vtu
@@ -0,0 +1,85 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian" header_type="UInt32" compressor="vtkZLibDataCompressor">
+  <UnstructuredGrid>
+    <Piece NumberOfPoints="44" NumberOfCells="10">
+      <PointData Scalars="pressure">
+        <DataArray type="Float64" Name="pressure" format="binary" RangeMin="-4.470348341500951e-9" RangeMax="0.3">
+          AQAAAACAAABgAQAAXAAAAA==eJwzNgaBy/ZAQnSd+0X79LS0d4E7ztmnpaVtcn94CkSX2XIdB9F6xsaH7c+dOfMncMc+EJ1jy7XdfvbMmZJBO9bZmwMNCNwxz/7///+GQOY+41FzweYCAJipnj0=
+        </DataArray>
+        <DataArray type="Float64" Name="C" format="binary" RangeMin="1" RangeMax="1">
+          AQAAAACAAABgAQAAEgAAAA==eJxjYACBD/YMozRNaABeMTQV
+        </DataArray>
+        <DataArray type="Float64" Name="CFlux" NumberOfComponents="3" format="binary" RangeMin="3.75e-10" RangeMax="3.75e-10">
+          AQAAAACAAAAgBAAAHQAAAA==eJy772h1R/DoT1sGNHB/VHxUfFR8xIsDALkEu1k=
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              3.75e-10
+            </Value>
+            <Value index="1">
+              3.75e-10
+            </Value>
+          </InformationKey>
+        </DataArray>
+        <DataArray type="Float64" Name="darcy_velocity" NumberOfComponents="3" format="binary" RangeMin="3.75e-10" RangeMax="3.75e-10">
+          AQAAAACAAAAgBAAAHQAAAA==eJy772h1R/DoT1sGNHB/VHxUfFR8xIsDALkEu1k=
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              3.75e-10
+            </Value>
+            <Value index="1">
+              3.75e-10
+            </Value>
+          </InformationKey>
+        </DataArray>
+      </PointData>
+      <CellData Scalars="CFlux_residual">
+        <DataArray type="Float64" Name="CFlux_residual" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="0">
+          AQAAAACAAADwAAAADAAAAA==eJxjYBhZAAAA8AAB
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              0
+            </Value>
+            <Value index="1">
+              0
+            </Value>
+          </InformationKey>
+        </DataArray>
+        <DataArray type="Float64" Name="darcy_velocity_residual" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="0">
+          AQAAAACAAADwAAAADAAAAA==eJxjYBhZAAAA8AAB
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              0
+            </Value>
+            <Value index="1">
+              0
+            </Value>
+          </InformationKey>
+        </DataArray>
+      </CellData>
+      <Points>
+        <DataArray type="Float64" Name="Points" NumberOfComponents="3" format="binary" RangeMin="0" RangeMax="0.8079604065183404">
+          AQAAAACAAAAgBAAArAAAAA==eJx90rEJwzAUhGHXmSKLmGSFrOMqjYdwlyFSZADVrtSkERYIjMgSqR48H/+zG+MPccLHDUP43Ov0vrEn8u3xWYPzmfy1LN8gp5Bf518N8hv5c7zswb2d3OUd/tu+wRO59QDnM7n1ADmF3HqA/EZuPcC9nTzq46SXRK57cOczue7B5RRy3YPLb+S6B3dvJ4/+/+SdyHUP2oO67kF7UNc9aA/qugft4eB/jXzcvQ==
+          <InformationKey name="L2_NORM_RANGE" location="vtkDataArray" length="2">
+            <Value index="0">
+              0
+            </Value>
+            <Value index="1">
+              0.80796040652
+            </Value>
+          </InformationKey>
+        </DataArray>
+      </Points>
+      <Cells>
+        <DataArray type="Int64" Name="connectivity" format="binary" RangeMin="0" RangeMax="43">
+          AQAAAACAAACAAgAAbgAAAA==eJxdyzcCggAQAEGMmEHBgGJWDPz/fxZeNdtsNUnyrxOfxmfxTXwbv8ZvuC5ujtvh7rgeboGrcA9cH5fh9rgnboDLcQdcgxvilrga98KluBXuiHvjRrgCd8J9cGNciTvjvrgJbo274Nr4D5lpBrk=
+        </DataArray>
+        <DataArray type="Int64" Name="offsets" format="binary" RangeMin="8" RangeMax="80">
+          AQAAAACAAABQAAAAIwAAAA==eJzjYIAAASgtAaUVoLQGlDaA0hZQ2gFKe0DpACgNADdQAbk=
+        </DataArray>
+        <DataArray type="UInt8" Name="types" format="binary" RangeMin="11" RangeMax="11">
+          AQAAAACAAAAKAAAACwAAAA==eJzj5oYBAAJnAG8=
+        </DataArray>
+      </Cells>
+    </Piece>
+  </UnstructuredGrid>
+</VTKFile>