From c8d099f18e130fdc12d6262b92f1e099bbe6933d Mon Sep 17 00:00:00 2001
From: Thomas Fetzer <thomas.fetzer@iws.uni-stuttgart.de>
Date: Fri, 11 Dec 2015 11:16:18 +0100
Subject: [PATCH] [boundarytypes,non-isothermal] Revert boundarytypes to almost
 original state (see below), update the non-isothermal models and reference
 solutions.

Problems occured if the boundary types are herited or more precise the problem was to reset out of the parent class the boundary types from the childs classes.
---
 dumux/common/boundarytypes.hh                 | 267 +++++--
 .../2cnistokes2p2cnilocaloperator.hh          |   8 +-
 .../common/multidomainboundarytypes.hh        | 219 ------
 .../common/subdomainpropertydefaults.hh       |   5 -
 .../2p2cnicouplinglocalresidual.hh            |  12 +-
 .../stokesncnicouplinglocalresidual.hh        |  35 +-
 .../2cnistokes2p2cniproblem.hh                |  13 +-
 .../2cnistokes2p2cni/2p2cnisubproblem.hh      |   2 +-
 .../2cnistokes2p2cni/stokes2cnisubproblem.hh  |   2 +-
 .../2cnizeroeq2p2cni/2p2cnisubproblem.hh      |   2 +-
 .../2cnizeroeq2p2cni/zeroeq2cnisubproblem.hh  |   2 +-
 .../2cstokes2p2c/2cstokes2p2cproblem.hh       |   2 -
 .../2czeroeq2p2c/2p2csubproblem.hh            |   2 +-
 .../2czeroeq2p2c/zeroeq2csubproblem.hh        |   2 +-
 .../2cnistokes2p2cni-ff-reference.vtu         | 156 ++--
 ...stokes2p2cniboundarylayer-ff-reference.vtu | 162 ++--
 ...stokes2p2cniboundarylayer-pm-reference.vtu | 720 +++++++++---------
 17 files changed, 749 insertions(+), 862 deletions(-)
 delete mode 100644 dumux/multidomain/common/multidomainboundarytypes.hh

diff --git a/dumux/common/boundarytypes.hh b/dumux/common/boundarytypes.hh
index 1e8a7faa74..8d4864c695 100644
--- a/dumux/common/boundarytypes.hh
+++ b/dumux/common/boundarytypes.hh
@@ -40,29 +40,41 @@ public:
     { reset(); }
 
     /*!
-     * \brief Reset the boundary types.
+     * \brief Reset the boundary types for all equations.
      *
      * After this method no equations will be disabled and neither
      * Neumann nor Dirichlet conditions will be evaluated. This
-     * corrosponds to a Neumann zero boundary.
+     * corresponds to a Neumann zero boundary.
      */
     void reset()
     {
-        for (int i=0; i < numEq; ++i) {
-            boundaryInfo_[i].visited = 0;
-
-            boundaryInfo_[i].isDirichlet = 0;
-            boundaryInfo_[i].isNeumann = 0;
-            boundaryInfo_[i].isOutflow = 0;
-            boundaryInfo_[i].isCouplingInflow = 0;
-            boundaryInfo_[i].isCouplingOutflow = 0;
-            boundaryInfo_[i].isMortarCoupling = 0;
-
-            eq2pvIdx_[i] = i;
-            pv2eqIdx_[i] = i;
+        for (int eqIdx=0; eqIdx < numEq; ++eqIdx)
+        {
+            resetEq(eqIdx);
         }
     }
 
+    /*!
+     * \brief Reset the boundary types for one equation.
+     */
+    void resetEq(int eqIdx)
+    {
+          boundaryInfo_[eqIdx].visited = 0;
+
+          boundaryInfo_[eqIdx].isDirichlet = 0;
+          boundaryInfo_[eqIdx].isNeumann = 0;
+          boundaryInfo_[eqIdx].isOutflow = 0;
+          boundaryInfo_[eqIdx].isCouplingInflow = 0;
+          boundaryInfo_[eqIdx].isCouplingOutflow = 0;
+          boundaryInfo_[eqIdx].isMortarCoupling = 0;
+          boundaryInfo_[eqIdx].isCouplingDirichlet = 0;
+          boundaryInfo_[eqIdx].isCouplingNeumann = 0;
+          boundaryInfo_[eqIdx].isCouplingMortar = 0;
+
+          eq2pvIdx_[eqIdx] = eqIdx;
+          pv2eqIdx_[eqIdx] = eqIdx;
+    }
+
     /*!
      * \brief Returns true if the boundary types for a given equation
      *        has been specified.
@@ -75,7 +87,7 @@ public:
     /*!
      * \brief Make sure the boundary conditions are well-posed.
      *
-     * If they are not, an assertation fails and the program aborts!
+     * If they are not, an assertion fails and the program aborts!
      * (if the NDEBUG macro is not defined)
      */
     void checkWellPosed() const
@@ -120,22 +132,47 @@ public:
         }
     }
 
+    /*!
+     * \brief Set all boundary conditions to Dirichlet-like coupling
+     */
+    void setAllCouplingDirichlet()
+    {
+        for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
+        {
+            setCouplingDirichlet(eqIdx);
+        }
+    }
+
+    /*!
+     * \brief Set all boundary conditions to Neumann-like coupling.
+     */
+    void setAllCouplingNeumann()
+    {
+        for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
+        {
+            setCouplingNeumann(eqIdx);
+        }
+    }
+    /*!
+     * \brief Set all boundary conditions to mortar coupling.
+     */
+    void setAllCouplingMortar()
+    {
+        for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
+        {
+            setCouplingDirichlet(eqIdx);
+        }
+    }
+
     /*!
      * \brief Set all boundary conditions to coupling inflow.
      */
     DUNE_DEPRECATED_MSG("setAllCouplingInflow() is deprecated")
     void setAllCouplingInflow()
     {
-        for (int eqIdx = 0; eqIdx < numEq; ++eqIdx) {
-            boundaryInfo_[eqIdx].visited = 1;
-            boundaryInfo_[eqIdx].isDirichlet = 0;
-            boundaryInfo_[eqIdx].isNeumann = 0;
-            boundaryInfo_[eqIdx].isOutflow = 0;
-            boundaryInfo_[eqIdx].isCouplingInflow = 1;
-            boundaryInfo_[eqIdx].isCouplingOutflow = 0;
-            boundaryInfo_[eqIdx].isMortarCoupling = 0;
-
-            Valgrind::SetDefined(boundaryInfo_[eqIdx]);
+        for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
+        {
+            setCouplingInflow(eqIdx);
         }
     }
 
@@ -145,16 +182,9 @@ public:
     DUNE_DEPRECATED_MSG("setAllCouplingOutflow() is deprecated")
     void setAllCouplingOutflow()
     {
-        for (int eqIdx = 0; eqIdx < numEq; ++eqIdx) {
-            boundaryInfo_[eqIdx].visited = 1;
-            boundaryInfo_[eqIdx].isDirichlet = 0;
-            boundaryInfo_[eqIdx].isNeumann = 0;
-            boundaryInfo_[eqIdx].isOutflow = 0;
-            boundaryInfo_[eqIdx].isCouplingInflow = 0;
-            boundaryInfo_[eqIdx].isCouplingOutflow = 1;
-            boundaryInfo_[eqIdx].isMortarCoupling = 0;
-
-            Valgrind::SetDefined(boundaryInfo_[eqIdx]);
+        for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
+        {
+            setCouplingOutflow(eqIdx);
         }
     }
 
@@ -164,16 +194,9 @@ public:
     DUNE_DEPRECATED_MSG("setAllMortarCoupling() is deprecated")
     void setAllMortarCoupling()
     {
-        for (int eqIdx = 0; eqIdx < numEq; ++eqIdx) {
-            boundaryInfo_[eqIdx].visited = 1;
-            boundaryInfo_[eqIdx].isDirichlet = 0;
-            boundaryInfo_[eqIdx].isNeumann = 0;
-            boundaryInfo_[eqIdx].isOutflow = 0;
-            boundaryInfo_[eqIdx].isCouplingInflow = 0;
-            boundaryInfo_[eqIdx].isCouplingOutflow = 0;
-            boundaryInfo_[eqIdx].isMortarCoupling = 1;
-
-            Valgrind::SetDefined(boundaryInfo_[eqIdx]);
+        for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
+        {
+            setCouplingMortar(eqIdx);
         }
     }
 
@@ -185,13 +208,9 @@ public:
      */
     void setNeumann(int eqIdx)
     {
+        resetEq(eqIdx);
         boundaryInfo_[eqIdx].visited = 1;
-        boundaryInfo_[eqIdx].isDirichlet = 0;
         boundaryInfo_[eqIdx].isNeumann = 1;
-        boundaryInfo_[eqIdx].isOutflow = 0;
-        boundaryInfo_[eqIdx].isCouplingInflow = 0;
-        boundaryInfo_[eqIdx].isCouplingOutflow = 0;
-        boundaryInfo_[eqIdx].isMortarCoupling = 0;
 
         Valgrind::SetDefined(boundaryInfo_[eqIdx]);
     }
@@ -207,13 +226,9 @@ public:
      */
     void setDirichlet(int pvIdx, int eqIdx)
     {
+        resetEq(eqIdx);
         boundaryInfo_[eqIdx].visited = 1;
         boundaryInfo_[eqIdx].isDirichlet = 1;
-        boundaryInfo_[eqIdx].isNeumann = 0;
-        boundaryInfo_[eqIdx].isOutflow = 0;
-        boundaryInfo_[eqIdx].isCouplingInflow = 0;
-        boundaryInfo_[eqIdx].isCouplingOutflow = 0;
-        boundaryInfo_[eqIdx].isMortarCoupling = 0;
 
         // update the equation <-> primary variable mapping
         eq2pvIdx_[eqIdx] = pvIdx;
@@ -231,13 +246,48 @@ public:
      */
     void setOutflow(int eqIdx)
     {
+        resetEq(eqIdx);
         boundaryInfo_[eqIdx].visited = 1;
-        boundaryInfo_[eqIdx].isDirichlet = 0;
-        boundaryInfo_[eqIdx].isNeumann = 0;
         boundaryInfo_[eqIdx].isOutflow = 1;
-        boundaryInfo_[eqIdx].isCouplingInflow = 0;
-        boundaryInfo_[eqIdx].isCouplingOutflow = 0;
-        boundaryInfo_[eqIdx].isMortarCoupling = 0;
+
+        Valgrind::SetDefined(boundaryInfo_[eqIdx]);
+    }
+
+    /*!
+     * \brief Set a boundary condition for a single equation to
+     *        a Dirichlet-like coupling condition.
+     */
+    void setCouplingDirichlet(int eqIdx)
+    {
+        resetEq(eqIdx);
+        boundaryInfo_[eqIdx].visited = 1;
+        boundaryInfo_[eqIdx].isCouplingDirichlet = 1;
+
+        Valgrind::SetDefined(boundaryInfo_[eqIdx]);
+    }
+
+    /*!
+     * \brief Set a boundary condition for a single equation to
+     *        a Neumann-like coupling condition.
+     */
+    void setCouplingNeumann(int eqIdx)
+    {
+        resetEq(eqIdx);
+        boundaryInfo_[eqIdx].visited = 1;
+        boundaryInfo_[eqIdx].isCouplingNeumann = 1;
+
+        Valgrind::SetDefined(boundaryInfo_[eqIdx]);
+    }
+
+    /*!
+     * \brief Set a boundary condition for a single equation to
+     *        a mortar coupling condition.
+     */
+    void setCouplingMortar(int eqIdx)
+    {
+        resetEq(eqIdx);
+        boundaryInfo_[eqIdx].visited = 1;
+        boundaryInfo_[eqIdx].isCouplingMortar = 1;
 
         Valgrind::SetDefined(boundaryInfo_[eqIdx]);
     }
@@ -248,13 +298,9 @@ public:
     DUNE_DEPRECATED_MSG("setCouplingInflow() is deprecated")
     void setCouplingInflow(int eqIdx)
     {
+        resetEq(eqIdx);
         boundaryInfo_[eqIdx].visited = 1;
-        boundaryInfo_[eqIdx].isDirichlet = 0;
-        boundaryInfo_[eqIdx].isNeumann = 0;
-        boundaryInfo_[eqIdx].isOutflow = 0;
         boundaryInfo_[eqIdx].isCouplingInflow = 1;
-        boundaryInfo_[eqIdx].isCouplingOutflow = 0;
-        boundaryInfo_[eqIdx].isMortarCoupling = 0;
 
         Valgrind::SetDefined(boundaryInfo_[eqIdx]);
     }
@@ -265,13 +311,9 @@ public:
     DUNE_DEPRECATED_MSG("setCouplingOutflow() is deprecated")
     void setCouplingOutflow(int eqIdx)
     {
+        resetEq(eqIdx);
         boundaryInfo_[eqIdx].visited = 1;
-        boundaryInfo_[eqIdx].isDirichlet = 0;
-        boundaryInfo_[eqIdx].isNeumann = 0;
-        boundaryInfo_[eqIdx].isOutflow = 0;
-        boundaryInfo_[eqIdx].isCouplingInflow = 0;
         boundaryInfo_[eqIdx].isCouplingOutflow = 1;
-        boundaryInfo_[eqIdx].isMortarCoupling = 0;
 
         Valgrind::SetDefined(boundaryInfo_[eqIdx]);
     }
@@ -282,12 +324,8 @@ public:
     DUNE_DEPRECATED_MSG("setMortarCoupling() is deprecated")
     void setMortarCoupling(int eqIdx)
     {
+        resetEq(eqIdx);
         boundaryInfo_[eqIdx].visited = 1;
-        boundaryInfo_[eqIdx].isDirichlet = 0;
-        boundaryInfo_[eqIdx].isNeumann = 0;
-        boundaryInfo_[eqIdx].isOutflow = 0;
-        boundaryInfo_[eqIdx].isCouplingInflow = 0;
-        boundaryInfo_[eqIdx].isCouplingOutflow = 0;
         boundaryInfo_[eqIdx].isMortarCoupling = 1;
 
         Valgrind::SetDefined(boundaryInfo_[eqIdx]);
@@ -442,17 +480,82 @@ public:
         return false;
     }
 
+    /*!
+     * \brief Returns true if an equation is used to specify an
+     *        Dirichlet coupling condition.
+     *
+     * \param eqIdx The index of the equation
+     */
+    bool isCouplingDirichlet(unsigned eqIdx) const
+    { return boundaryInfo_[eqIdx].isCouplingDirichlet; }
+
+    /*!
+     * \brief Returns true if some equation is used to specify an
+     *        Dirichlet coupling condition.
+     */
+    bool hasCouplingDirichlet() const
+    {
+        for (int i = 0; i < numEq; ++i)
+            if (boundaryInfo_[i].isCouplingDirichlet)
+                return true;
+        return false;
+    }
+
+    /*!
+     * \brief Returns true if an equation is used to specify an
+     *        Neumann coupling condition.
+     *
+     * \param eqIdx The index of the equation
+     */
+    bool isCouplingNeumann(unsigned eqIdx) const
+    { return boundaryInfo_[eqIdx].isCouplingNeumann; }
+
+    /*!
+     * \brief Returns true if some equation is used to specify an
+     *        Neumann coupling condition.
+     */
+    bool hasCouplingNeumann() const
+    {
+        for (int i = 0; i < numEq; ++i)
+            if (boundaryInfo_[i].isCouplingNeumann)
+                return true;
+        return false;
+    }
+
     /*!
      * \brief Returns true if an equation is used to specify a
-     *        coupling condition.
+     *        Mortar coupling condition.
      *
      * \param eqIdx The index of the equation
+     */
+    bool isCouplingMortar(unsigned eqIdx) const
+    {
+        return boundaryInfo_[eqIdx].isCouplingMortar;
+    }
+
+    /*!
+     * \brief Returns true if some equation is used to specify an
+     *        Mortar coupling condition.
+     */
+    bool hasCouplingMortar() const
+    {
+        for (int i = 0; i < numEq; ++i)
+            if (boundaryInfo_[i].isCouplingMortar)
+                return true;
+        return false;
+    }
+
+    /*!
+     * \brief Returns true if an equation is used to specify a
+     *        coupling condition.
      *
-     * This only used for correcting the pressure in the Stokes equation.
+     * \param eqIdx The index of the equation
      */
     bool isCoupling(unsigned eqIdx) const
     {
-        return false;
+        return boundaryInfo_[eqIdx].isCouplingDirichlet
+               || boundaryInfo_[eqIdx].isCouplingNeumann
+               || boundaryInfo_[eqIdx].isCouplingMortar;
     }
 
     /*!
@@ -461,6 +564,9 @@ public:
      */
     bool hasCoupling() const
     {
+        for (int i = 0; i < numEq; ++i)
+            if (isCoupling(i))
+                return true;
         return false;
     }
 
@@ -486,7 +592,7 @@ public:
     unsigned eqToDirichletIndex(unsigned eqIdx) const
     { return eq2pvIdx_[eqIdx]; }
 
-private:
+protected:
     // this is a bitfield structure!
     struct __attribute__((__packed__)) {
         unsigned char visited : 1;
@@ -496,6 +602,9 @@ private:
         unsigned char isCouplingInflow : 1;
         unsigned char isCouplingOutflow : 1;
         unsigned char isMortarCoupling : 1;
+        unsigned char isCouplingDirichlet : 1;
+        unsigned char isCouplingNeumann : 1;
+        unsigned char isCouplingMortar : 1;
     } boundaryInfo_[numEq];
 
     unsigned char eq2pvIdx_[numEq];
diff --git a/dumux/multidomain/2cnistokes2p2cni/2cnistokes2p2cnilocaloperator.hh b/dumux/multidomain/2cnistokes2p2cni/2cnistokes2p2cnilocaloperator.hh
index a470c6eda2..9fa2046249 100644
--- a/dumux/multidomain/2cnistokes2p2cni/2cnistokes2p2cnilocaloperator.hh
+++ b/dumux/multidomain/2cnistokes2p2cni/2cnistokes2p2cnilocaloperator.hh
@@ -220,7 +220,7 @@ public:
                                    cParams,
                                    couplingRes1, couplingRes2);
 
-        if (cParams.boundaryTypes2.isCouplingInflow(energyEqIdx2))
+        if (cParams.boundaryTypes2.isCouplingNeumann(energyEqIdx2))
         {
             unsigned int blModel = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, int, BoundaryLayer, Model);
             // only enter here, if a boundary layer model  is used for the computation of the diffusive fluxes
@@ -342,7 +342,7 @@ public:
                                         globalProblem.localResidual1().residual(vertInElem1)[energyEqIdx1]);
             }
         }
-        if (cParams.boundaryTypes2.isCouplingOutflow(energyEqIdx2))
+        if (cParams.boundaryTypes2.isCouplingDirichlet(energyEqIdx2))
         {
             // set residualDarcy[energyEqIdx2] = T in 2p2cnilocalresidual.hh
             couplingRes2.accumulate(lfsu_n.child(energyEqIdx2), vertInElem2,
@@ -378,13 +378,13 @@ public:
             normalMassFlux2[phaseIdx] = -boundaryVars2.volumeFlux(phaseIdx)*
                 cParams.elemVolVarsCur2[vertInElem2].density(phaseIdx);
 
-        if (cParams.boundaryTypes1.isCouplingOutflow(energyEqIdx1))
+        if (cParams.boundaryTypes1.isCouplingDirichlet(energyEqIdx1))
         {
             // set residualStokes[energyIdx1] = T in stokes2cnilocalresidual.hh
             couplingRes1.accumulate(lfsu_s.child(energyEqIdx1), vertInElem1,
                                     -cParams.elemVolVarsCur2[vertInElem2].temperature());
         }
-        if (cParams.boundaryTypes1.isCouplingInflow(energyEqIdx1))
+        if (cParams.boundaryTypes1.isCouplingNeumann(energyEqIdx1))
         {
             if (globalProblem.sdProblem2().isCornerPoint(globalPos2))
             {
diff --git a/dumux/multidomain/common/multidomainboundarytypes.hh b/dumux/multidomain/common/multidomainboundarytypes.hh
deleted file mode 100644
index 61059881f7..0000000000
--- a/dumux/multidomain/common/multidomainboundarytypes.hh
+++ /dev/null
@@ -1,219 +0,0 @@
-// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
-// vi: set et ts=4 sw=4 sts=4:
-/*****************************************************************************
- *   See the file COPYING for full copying permissions.                      *
- *                                                                           *
- *   This program is free software: you can redistribute it and/or modify    *
- *   it under the terms of the GNU General Public License as published by    *
- *   the Free Software Foundation, either version 2 of the License, or       *
- *   (at your option) any later version.                                     *
- *                                                                           *
- *   This program is distributed in the hope that it will be useful,         *
- *   but WITHOUT ANY WARRANTY; without even the implied warranty of          *
- *   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the            *
- *   GNU General Public License for more details.                            *
- *                                                                           *
- *   You should have received a copy of the GNU General Public License       *
- *   along with this program.  If not, see <http://www.gnu.org/licenses/>.   *
- *****************************************************************************/
-/*!
- * \file
- * \brief Class to specify the type of a boundary for multidomain problems.
- */
-#ifndef MULTIDOMAIN_BOUNDARY_TYPES_HH
-#define MULTIDOMAIN_BOUNDARY_TYPES_HH
-
-#include <dumux/common/boundarytypes.hh>
-#include <dumux/common/valgrind.hh>
-
-namespace Dumux
-{
-
-/*!
- * \ingroup BC
- * \brief Class to specify the type of a boundary for multidomain problems.
- */
-template <int numEq>
-class MultidomainBoundaryTypes : public BoundaryTypes<numEq>
-{
-    typedef BoundaryTypes<numEq> ParentType;
-
-public:
-    MultidomainBoundaryTypes()
-    { reset(); }
-
-    /*!
-     * \brief Reset the boundary types.
-     *
-     * After this method no equations will be disabled and neither
-     * neumann nor dirichlet conditions will be evaluated. This
-     * corrosponds to a neumann zero boundary.
-     */
-    void reset()
-    {
-        ParentType::reset();
-        for (int i=0; i < numEq; ++i) {
-            boundaryCouplingInfo_[i].visited = 0;
-            boundaryCouplingInfo_[i].isCouplingDirichlet = 0;
-            boundaryCouplingInfo_[i].isCouplingNeumann = 0;
-        }
-    }
-
-    /*!
-     * \brief Returns true if the boundary types for a given equation
-     *        has been specified.
-     *
-     * \param eqIdx The index of the equation
-     */
-    bool isSet(int eqIdx) const
-    {
-        return ParentType::boundaryInfo_[eqIdx].visited
-               || boundaryCouplingInfo_[eqIdx].visited;
-    }
-
-    /*!
-     * \brief Set all boundary conditions to coupling Dirichlet.
-     */
-    void setAllCouplingDirichlet()
-    {
-        for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
-        {
-            setCouplingDirichlet(eqIdx);
-        }
-    }
-
-    /*!
-     * \brief Set all boundary conditions to coupling Neumann.
-     */
-    void setAllCouplingNeumann()
-    {
-        for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
-        {
-            setCouplingNeumann(eqIdx);
-        }
-    }
-
-    /*!
-     * \brief Set a boundary condition for a single equation to coupling inflow.
-     */
-    void setCouplingDirichlet(int eqIdx)
-    {
-        ParentType::reset();
-        boundaryCouplingInfo_[eqIdx].visited = 1;
-        boundaryCouplingInfo_[eqIdx].isCouplingDirichlet = 1;
-        boundaryCouplingInfo_[eqIdx].isCouplingNeumann = 0;
-
-        Valgrind::SetDefined(boundaryCouplingInfo_[eqIdx]);
-    }
-
-    /*!
-     * \brief Set a boundary condition for a single equation to coupling outflow.
-     */
-    void setCouplingNeumann(int eqIdx)
-    {
-        ParentType::reset();
-        boundaryCouplingInfo_[eqIdx].visited = 1;
-        boundaryCouplingInfo_[eqIdx].isCouplingDirichlet = 0;
-        boundaryCouplingInfo_[eqIdx].isCouplingNeumann = 1;
-
-        Valgrind::SetDefined(boundaryCouplingInfo_[eqIdx]);
-    }
-
-    /*!
-     * \brief Returns true if an equation is used to specify an
-     *        Dirichlet coupling condition.
-     *
-     * \param eqIdx The index of the equation
-     */
-    bool isCouplingDirichlet(unsigned eqIdx) const
-    { return boundaryCouplingInfo_[eqIdx].isCouplingDirichlet; }
-
-    /*!
-     * \brief Returns true if some equation is used to specify an
-     *        Dirichlet coupling condition.
-     */
-    bool hasCouplingDirichlet() const
-    {
-        for (int i = 0; i < numEq; ++i)
-            if (boundaryCouplingInfo_[i].isCouplingDirichlet)
-                return true;
-        return false;
-    }
-
-    /*!
-     * \brief Returns true if an equation is used to specify an
-     *        Neumann coupling condition.
-     *
-     * \param eqIdx The index of the equation
-     */
-    bool isCouplingNeumann(unsigned eqIdx) const
-    { return boundaryCouplingInfo_[eqIdx].isCouplingNeumann; }
-
-    /*!
-     * \brief Returns true if some equation is used to specify an
-     *        Neumann coupling condition.
-     */
-    bool hasCouplingNeumann() const
-    {
-        for (int i = 0; i < numEq; ++i)
-            if (boundaryCouplingInfo_[i].isCouplingNeumann)
-                return true;
-        return false;
-    }
-
-    /*!
-     * \brief Returns true if an equation is used to specify a
-     *        Mortar coupling condition.
-     *
-     * \param eqIdx The index of the equation
-     */
-    bool isCouplingMortar(unsigned eqIdx) const
-    {
-        return false;
-    }
-
-    /*!
-     * \brief Returns true if some equation is used to specify an
-     *        Mortar coupling condition.
-     */
-    bool hasCouplingMortar() const
-    {
-        return false;
-    }
-
-    /*!
-     * \brief Returns true if an equation is used to specify a
-     *        coupling condition.
-     *
-     * \param eqIdx The index of the equation
-     */
-    bool isCoupling(unsigned eqIdx) const
-    {
-        return boundaryCouplingInfo_[eqIdx].isCouplingDirichlet
-               || boundaryCouplingInfo_[eqIdx].isCouplingNeumann;
-    }
-
-    /*!
-     * \brief Returns true if some equation is used to specify an
-     *        coupling condition.
-     */
-    bool hasCoupling() const
-    {
-        for (int i = 0; i < numEq; ++i)
-            if (isCoupling(i))
-                return true;
-        return false;
-    }
-
-private:
-    // this is a bitfield structure!
-    struct __attribute__((__packed__)) {
-        unsigned char visited : 1;
-        unsigned char isCouplingDirichlet : 1;
-        unsigned char isCouplingNeumann : 1;
-    } boundaryCouplingInfo_[numEq];
-};
-
-}
-
-#endif
diff --git a/dumux/multidomain/common/subdomainpropertydefaults.hh b/dumux/multidomain/common/subdomainpropertydefaults.hh
index cb10a0822c..6f36b304e3 100644
--- a/dumux/multidomain/common/subdomainpropertydefaults.hh
+++ b/dumux/multidomain/common/subdomainpropertydefaults.hh
@@ -33,7 +33,6 @@
 #include <dune/pdelab/constraints/conforming.hh>
 
 #include "subdomainproperties.hh"
-#include "multidomainboundarytypes.hh"
 #include "multidomainproperties.hh"
 #include "multidomainlocaloperator.hh"
 #include <dumux/multidomain/couplinglocalresiduals/boxcouplinglocalresidual.hh>
@@ -133,10 +132,6 @@ public:
     { return ParameterTree::unusedNewRunTimeParams(); }
 };
 
-//! Boundary types at a single degree of freedom
-SET_TYPE_PROP(SubDomain, BoundaryTypes,
-              MultidomainBoundaryTypes<GET_PROP_VALUE(TypeTag, NumEq)>);
-
 } // namespace Properties
 } // namespace Dumux
 
diff --git a/dumux/multidomain/couplinglocalresiduals/2p2cnicouplinglocalresidual.hh b/dumux/multidomain/couplinglocalresiduals/2p2cnicouplinglocalresidual.hh
index c6441ab018..8d4a43d144 100644
--- a/dumux/multidomain/couplinglocalresiduals/2p2cnicouplinglocalresidual.hh
+++ b/dumux/multidomain/couplinglocalresiduals/2p2cnicouplinglocalresidual.hh
@@ -109,14 +109,13 @@ public:
                 // loop over the single vertices on the current face
                 for (int faceVertexIdx = 0; faceVertexIdx < numFaceVertices; ++faceVertexIdx)
                 {
-                    const int boundaryFaceIdx = this->fvGeometry_().boundaryFaceIndex(fIdx, faceVertexIdx);
                     const int vIdx = refElement.subEntity(fIdx, 1, faceVertexIdx, dim);
                     // only evaluate, if we consider the same face vertex as in the outer
                     // loop over the element vertices
                     if (vIdx != idx)
                         continue;
 
-                    if (boundaryHasCoupling_(this->bcTypes_(idx)))
+                    if (this->bcTypes_(idx).hasCoupling())
                         evalCouplingVertex_(idx);
                 }
             }
@@ -248,6 +247,8 @@ public:
      * \brief Set the Dirichlet-like conditions for the coupling
      *        and replace the existing residual
      *
+     * TODO integrate this function into the evalBoundary_?
+     *
      * \param scvIdx Sub control vertex index for the coupling condition
      */
     void evalCouplingVertex_(const int scvIdx)
@@ -255,21 +256,22 @@ public:
         const VolumeVariables &volVars = this->curVolVars_()[scvIdx];
 
         // set pressure as part of the momentum coupling
-        if (this->bcTypes_(scvIdx).isCouplingOutflow(massBalanceIdx))
+        if (this->bcTypes_(scvIdx).isCouplingDirichlet(massBalanceIdx))
             this->residual_[scvIdx][massBalanceIdx] = volVars.pressure(nPhaseIdx);
 
         // set mass fraction;
-        if (this->bcTypes_(scvIdx).isCouplingOutflow(contiWEqIdx))
+        if (this->bcTypes_(scvIdx).isCouplingDirichlet(contiWEqIdx))
             this->residual_[scvIdx][contiWEqIdx] = volVars.massFraction(nPhaseIdx, wCompIdx);
 
         // set temperature as part of the coupling
-        if (this->bcTypes_(scvIdx).isCouplingOutflow(energyEqIdx))
+        if (this->bcTypes_(scvIdx).isCouplingDirichlet(energyEqIdx))
             this->residual_[scvIdx][energyEqIdx] = volVars.temperature();
     }
 
     /*!
      * \brief Check if one of the boundary conditions is coupling.
      */
+    DUNE_DEPRECATED_MSG("boundaryHasCoupling_ is deprecated")
     bool boundaryHasCoupling_(const BoundaryTypes& bcTypes) const
     {
         for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
diff --git a/dumux/multidomain/couplinglocalresiduals/stokesncnicouplinglocalresidual.hh b/dumux/multidomain/couplinglocalresiduals/stokesncnicouplinglocalresidual.hh
index 2a68740f67..7042411d7d 100644
--- a/dumux/multidomain/couplinglocalresiduals/stokesncnicouplinglocalresidual.hh
+++ b/dumux/multidomain/couplinglocalresiduals/stokesncnicouplinglocalresidual.hh
@@ -175,11 +175,12 @@ namespace Dumux
                         asImp_()->evalNeumannSegment_(&intersection, idx, boundaryFaceIdx, boundaryVars);
                         asImp_()->evalOutflowSegment_(&intersection, idx, boundaryFaceIdx, boundaryVars);
 
+                        // TODO: check if this code is used
                         //for the corner points, the boundary flux across the vertical non-coupling boundary faces
                         //has to be calculated to fulfill the mass balance
                         //convert suddomain intersection into multidomain intersection and check whether it is an outer boundary
                         if(!GridView::Grid::multiDomainIntersection(intersection).neighbor()
-                           && (this->boundaryHasMortarCoupling_(this->bcTypes_(idx)) || this->momentumBalanceHasNeumann_(this->bcTypes_(idx))))
+                           && (bcTypes.hasCouplingMortar() || this->momentumBalanceHasNeumann_(this->bcTypes_(idx))))
                         {
                             const GlobalPosition& globalPos = this->fvGeometry_().subContVol[idx].global;
                             //problem specific function, in problem orientation of interface is known
@@ -211,11 +212,11 @@ namespace Dumux
                             }
                         }
                         // Beavers-Joseph condition at the coupling boundary/interface
-                        if(boundaryHasCoupling_(bcTypes))
+                        if(bcTypes.hasCoupling())
                         {
                             evalBeaversJoseph_(&intersection, idx, boundaryFaceIdx, boundaryVars);
                         }
-                        if(boundaryHasCoupling_(bcTypes) || boundaryHasMortarCoupling_(bcTypes))
+                        if(bcTypes.hasCoupling() || bcTypes.hasCouplingMortar())
                         {
                             asImp_()->evalCouplingVertex_(&intersection, idx, boundaryFaceIdx, boundaryVars);
                         }
@@ -271,17 +272,19 @@ namespace Dumux
             const VolumeVariables &volVars = this->curVolVars_()[scvIdx];
             const BoundaryTypes &bcTypes = this->bcTypes_(scvIdx);
 
+            // TODO revise comment
             // set velocity normal to the interface
-            if (bcTypes.isCouplingInflow(momentumYIdx))
-                this->residual_[scvIdx][momentumYIdx] =
-                volVars.velocity() *
-                boundaryVars.face().normal /
-                boundaryVars.face().normal.two_norm();
-            Valgrind::CheckDefined(this->residual_[scvIdx][momentumYIdx]);
+            if (bcTypes.isCouplingNeumann(momentumYIdx))
+            {
+                this->residual_[scvIdx][momentumYIdx] = volVars.velocity()
+                                                        * boundaryVars.face().normal
+                                                        / boundaryVars.face().normal.two_norm();
+                Valgrind::CheckDefined(this->residual_[scvIdx][momentumYIdx]);
+            }
 
             // add pressure correction - required for pressure coupling,
             // if p.n comes from the pm
-            if (bcTypes.isCouplingOutflow(momentumYIdx) || bcTypes.isMortarCoupling(momentumYIdx))
+            if (bcTypes.isCouplingDirichlet(momentumYIdx) || bcTypes.isCouplingMortar(momentumYIdx))
             {
                 DimVector pressureCorrection(isIt->centerUnitOuterNormal());
                 pressureCorrection *= volVars.pressure();
@@ -294,7 +297,7 @@ namespace Dumux
             {
                 int eqIdx =  dim + compIdx;
                 if (eqIdx != massBalanceIdx) {
-                    if (bcTypes.isCouplingOutflow(eqIdx))
+                    if (bcTypes.isCouplingDirichlet(eqIdx))
                     {
                         if(useMoles)
                             this->residual_[scvIdx][eqIdx] = volVars.moleFraction(compIdx);
@@ -305,7 +308,7 @@ namespace Dumux
                 }
             }
             // set temperature
-            if (bcTypes.isCouplingOutflow(energyEqIdx))
+            if (bcTypes.isCouplingDirichlet(energyEqIdx))
             {
                 this->residual_[scvIdx][energyEqIdx] = volVars.temperature();
                 Valgrind::CheckDefined(this->residual_[scvIdx][energyEqIdx]);
@@ -330,9 +333,10 @@ namespace Dumux
             assert(beaversJosephCoeff > 0.0);
             const DimVector& elementUnitNormal = isIt->centerUnitOuterNormal();
 
+            // TODO revise comment
             // implementation as NEUMANN condition /////////////////////////////////////////////
             // (v.n)n
-            if (bcTypes.isCouplingOutflow(momentumXIdx))
+            if (bcTypes.isCouplingDirichlet(momentumXIdx))
             {
                 const Scalar normalComp = boundaryVars.velocity()*elementUnitNormal;
                 DimVector normalV = elementUnitNormal;
@@ -358,9 +362,10 @@ namespace Dumux
 //                     this->residual_[scvIdx][dimIdx] += priVars[dimIdx]*
 //                                                     boundaryFaceArea;
             }
-            if (bcTypes.isCouplingInflow(momentumXIdx))
+            if (bcTypes.isCouplingNeumann(momentumXIdx))
             {
                 assert(beaversJosephCoeff > 0);
+                // TODO revise comment and move upwards
                 ///////////////////////////////////////////////////////////////////////////////////////////
                 //IMPLEMENTATION AS DIRICHLET CONDITION
                 //tangential componment: vx = sqrt K /alpha * (grad v n(unity))t
@@ -379,6 +384,7 @@ namespace Dumux
         }
 
         // return true, if at least one equation on the boundary has a  coupling condition
+        DUNE_DEPRECATED_MSG("boundaryHasCoupling_() is deprecated")
         bool boundaryHasCoupling_(const BoundaryTypes& bcTypes) const
         {
             for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
@@ -388,6 +394,7 @@ namespace Dumux
         }
 
         // return true, if at least one equation on the boundary has a mortar coupling condition
+        DUNE_DEPRECATED_MSG("boundaryHasMortarCoupling_() is deprecated")
         bool boundaryHasMortarCoupling_(const BoundaryTypes& bcTypes) const
         {
             for (int eqIdx = 0; eqIdx < numEq; ++eqIdx)
diff --git a/test/multidomain/2cnistokes2p2cni/2cnistokes2p2cniproblem.hh b/test/multidomain/2cnistokes2p2cni/2cnistokes2p2cniproblem.hh
index 4658d35a33..c68cfd1a2a 100644
--- a/test/multidomain/2cnistokes2p2cni/2cnistokes2p2cniproblem.hh
+++ b/test/multidomain/2cnistokes2p2cni/2cnistokes2p2cniproblem.hh
@@ -159,7 +159,6 @@ class TwoCNIStokesTwoPTwoCNIProblem : public MultiDomainProblem<TypeTag>
     typedef typename MDGrid::SubDomainGrid SDGrid;
 
     typedef typename MDGrid::Traits::template Codim<0>::Entity MDElement;
-    typedef typename MDGrid::Traits::template Codim<0>::EntityPointer MDElementPointer;
     typedef typename Stokes2cniGridView::template Codim<0>::Entity SDElement1;
     typedef typename TwoPTwoCNIGridView::template Codim<0>::Entity SDElement2;
     typedef typename SDGrid::Traits::template Codim<0>::EntityPointer SDElementPointer;
@@ -384,10 +383,8 @@ public:
                  mdGrid.leafSubDomainInterfaceBegin(stokes2cni_, twoPtwoCNI_); ifIt != endIfIt; ++ifIt)
         {
             const int firstFaceIdx = ifIt->indexInFirstCell();
-            const MDElementPointer mdElementPointer1 = ifIt->firstCell(); // ATTENTION!!!
-            const MDElement& mdElement1 = *mdElementPointer1;     // Entity pointer has to be copied before.
-            const SDElementPointer sdElementPointer1 = this->sdElementPointer1(mdElement1);
-            const SDElement1& sdElement1 = *sdElementPointer1;
+            const MDElement& mdElement1 = ifIt->firstCell();
+            const SDElement1& sdElement1 = this->sdElementPointer1(mdElement1);
             fvGeometry1.update(this->sdGridView1(), sdElement1);
 
             const Dune::ReferenceElement<typename MDGrid::ctype,dim>& referenceElement1 =
@@ -520,10 +517,8 @@ public:
                  mdGrid.leafSubDomainInterfaceBegin(stokes2cni_, twoPtwoCNI_); ifIt != endIfIt; ++ifIt)
         {
             const int secondFaceIdx = ifIt->indexInSecondCell();
-            const MDElementPointer mdElementPointer2 = ifIt->secondCell(); // ATTENTION!!!
-            const MDElement& mdElement2 = *mdElementPointer2;     // Entity pointer has to be copied before.
-            const SDElementPointer sdElementPointer2 = this->sdElementPointer2(mdElement2);
-            const SDElement2& sdElement2 = *sdElementPointer2;
+            const MDElement& mdElement2 = ifIt->secondCell();
+            const SDElement2& sdElement2 = this->sdElementPointer2(mdElement2);
             fvGeometry2.update(this->sdGridView2(), sdElement2);
 
             const Dune::ReferenceElement<typename MDGrid::ctype,dim>& referenceElement2 =
diff --git a/test/multidomain/2cnistokes2p2cni/2p2cnisubproblem.hh b/test/multidomain/2cnistokes2p2cni/2p2cnisubproblem.hh
index b1bcaadbc3..67c8894ab4 100644
--- a/test/multidomain/2cnistokes2p2cni/2p2cnisubproblem.hh
+++ b/test/multidomain/2cnistokes2p2cni/2p2cnisubproblem.hh
@@ -263,7 +263,7 @@ public:
             {
                 if (globalPos[0] > runUpDistanceX_ - eps_)
                 {
-                    values.setAllCouplingInflow();
+                    values.setAllCouplingNeumann();
                 }
                 else
                     values.setAllNeumann();
diff --git a/test/multidomain/2cnistokes2p2cni/stokes2cnisubproblem.hh b/test/multidomain/2cnistokes2p2cni/stokes2cnisubproblem.hh
index c90bd7f466..499f42d69b 100644
--- a/test/multidomain/2cnistokes2p2cni/stokes2cnisubproblem.hh
+++ b/test/multidomain/2cnistokes2p2cni/stokes2cnisubproblem.hh
@@ -272,7 +272,7 @@ public:
 
             if (globalPos[0] > runUpDistanceX_-eps_ && time > initializationTime_)
             {
-                values.setAllCouplingOutflow();
+                values.setAllCouplingDirichlet();
             }
         }
 
diff --git a/test/multidomain/2cnizeroeq2p2cni/2p2cnisubproblem.hh b/test/multidomain/2cnizeroeq2p2cni/2p2cnisubproblem.hh
index 7a05762161..60cdd580dd 100644
--- a/test/multidomain/2cnizeroeq2p2cni/2p2cnisubproblem.hh
+++ b/test/multidomain/2cnizeroeq2p2cni/2p2cnisubproblem.hh
@@ -233,7 +233,7 @@ public:
             if (globalPos[0] > runUpDistanceX1_ - eps_
                 && globalPos[0] < runUpDistanceX2_ + eps_)
             {
-                values.setAllCouplingInflow();
+                values.setAllCouplingNeumann();
             }
             else
                 values.setAllNeumann();
diff --git a/test/multidomain/2cnizeroeq2p2cni/zeroeq2cnisubproblem.hh b/test/multidomain/2cnizeroeq2p2cni/zeroeq2cnisubproblem.hh
index 507c48f1b0..0496944231 100644
--- a/test/multidomain/2cnizeroeq2p2cni/zeroeq2cnisubproblem.hh
+++ b/test/multidomain/2cnizeroeq2p2cni/zeroeq2cnisubproblem.hh
@@ -223,7 +223,7 @@ public:
             if (globalPos[0] > runUpDistanceX1_ - eps_
                 && globalPos[0] < runUpDistanceX2_ + eps_)
             {
-                values.setAllCouplingOutflow();
+                values.setAllCouplingDirichlet();
             }
         }
 
diff --git a/test/multidomain/2cstokes2p2c/2cstokes2p2cproblem.hh b/test/multidomain/2cstokes2p2c/2cstokes2p2cproblem.hh
index 58e19ff016..e399ed5fad 100644
--- a/test/multidomain/2cstokes2p2c/2cstokes2p2cproblem.hh
+++ b/test/multidomain/2cstokes2p2c/2cstokes2p2cproblem.hh
@@ -152,10 +152,8 @@ class TwoCStokesTwoPTwoCProblem : public MultiDomainProblem<TypeTag>
     typedef typename MDGrid::LeafGridView MDGridView;
     enum { dim = MDGridView::dimension };
 
-    typedef typename MDGrid::Traits::template Codim<0>::EntityPointer MDElementPointer;
     typedef typename MDGrid::Traits::template Codim<0>::Entity MDElement;
     typedef typename MDGrid::SubDomainGrid SDGrid;
-    typedef typename SDGrid::template Codim<0>::EntityPointer SDElementPointer;
 
     typedef typename GET_PROP_TYPE(Stokes2cTypeTag, GridView) Stokes2cGridView;
     typedef typename GET_PROP_TYPE(TwoPTwoCTypeTag, GridView) TwoPTwoCGridView;
diff --git a/test/multidomain/2czeroeq2p2c/2p2csubproblem.hh b/test/multidomain/2czeroeq2p2c/2p2csubproblem.hh
index 7f7f929565..61e187008a 100644
--- a/test/multidomain/2czeroeq2p2c/2p2csubproblem.hh
+++ b/test/multidomain/2czeroeq2p2c/2p2csubproblem.hh
@@ -226,7 +226,7 @@ public:
         if (onUpperBoundary_(globalPos))
         {
             if (globalPos[0] > runUpDistanceX_ - eps_)
-                values.setAllCouplingInflow();
+                values.setAllCouplingNeumann();
             else
                 values.setAllNeumann();
         }
diff --git a/test/multidomain/2czeroeq2p2c/zeroeq2csubproblem.hh b/test/multidomain/2czeroeq2p2c/zeroeq2csubproblem.hh
index 185c23c345..bd86c454ba 100644
--- a/test/multidomain/2czeroeq2p2c/zeroeq2csubproblem.hh
+++ b/test/multidomain/2czeroeq2p2c/zeroeq2csubproblem.hh
@@ -233,7 +233,7 @@ public:
             values.setNeumann(transportEqIdx);
 
             if (globalPos[0] > runUpDistanceX_-eps_)
-                values.setAllCouplingOutflow();
+                values.setAllCouplingDirichlet();
         }
         if (onLeftBoundary_(globalPos))
         {
diff --git a/test/references/2cnistokes2p2cni-ff-reference.vtu b/test/references/2cnistokes2p2cni-ff-reference.vtu
index 15db143745..22c0138ac2 100644
--- a/test/references/2cnistokes2p2cni-ff-reference.vtu
+++ b/test/references/2cnistokes2p2cni-ff-reference.vtu
@@ -4,9 +4,9 @@
     <Piece NumberOfCells="120" NumberOfPoints="144">
       <PointData Scalars="temperature" Vectors="v">
         <DataArray type="Float32" Name="temperature" NumberOfComponents="1" format="ascii">
-          297.361 297.358 298.15 298.149 297.356 298.149 297.354 298.148 297.352 298.147 297.352 298.147
-          297.352 298.146 297.352 298.146 297.353 298.145 297.354 298.144 297.355 298.144 297.356 298.143
-          297.357 298.142 297.358 298.142 297.359 298.141 297.365 298.141 298.15 298.15 298.15 298.15
+          297.454 297.451 298.15 298.149 297.448 298.149 297.446 298.148 297.445 298.148 297.445 298.147
+          297.445 298.147 297.445 298.146 297.446 298.145 297.447 298.145 297.448 298.144 297.449 298.144
+          297.45 298.143 297.451 298.143 297.452 298.142 297.457 298.142 298.15 298.15 298.15 298.15
           298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15
           298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15
           298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15
@@ -32,23 +32,23 @@
           100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
         </DataArray>
         <DataArray type="Float32" Name="delP" NumberOfComponents="1" format="ascii">
-          0.00202823 0.00189852 0.0020099 0.00188019 0.00176355 0.00174687 0.00162778 0.001615 0.00149164 0.00148087 0.00135569 0.00134753
-          0.00121974 0.00121294 0.00108398 0.00107885 0.000948261 0.000944022 0.00081267 0.000809516 0.00067711 0.000674576 0.000541637 0.000539837
-          0.000406179 0.000404845 0.000270781 0.000269967 0.000135386 0.000134962 0 0 0.00203203 0.00189634 0.00176095 0.00162532
-          0.0014898 0.0013542 0.00121867 0.00108314 0.000947662 0.000812207 0.000676782 0.000541391 0.000406009 0.000270667 0.000135308 0
-          0.00202281 0.00188698 0.00175254 0.0016178 0.00148329 0.00134872 0.00121403 0.00107937 0.000944531 0.000809756 0.000674815 0.000539941
-          0.000404943 0.000270002 0.000134975 0 0.00203097 0.00189542 0.00175994 0.00162434 0.00148893 0.00135338 0.00121801 0.00108254
-          0.00094721 0.000811804 0.000676502 0.000541153 0.000405867 0.000270561 0.000135278 0 0.00202645 0.00189072 0.00175571 0.00162037
-          0.00148544 0.00135029 0.00121536 0.00108029 0.000945324 0.000810271 0.000675267 0.000540211 0.000405178 0.000270115 0.000135067 0
-          0.00202979 0.00189428 0.00175893 0.00162342 0.0014881 0.00135264 0.00121738 0.00108197 0.000946745 0.0008114 0.000676195 0.000540897
-          0.000405703 0.000270437 0.000135246 0 0.00202797 0.00189245 0.00175726 0.00162185 0.00148674 0.00135141 0.00121632 0.00108106
-          0.000945976 0.000810762 0.000675677 0.000540495 0.000405405 0.000270245 0.000135147 0 0.00202823 0.00189271 0.00175772 0.00162221
-          0.00148703 0.00135163 0.00121649 0.00108119 0.000946082 0.000810848 0.000675747 0.000540551 0.000405445 0.000270272 0.000135149 0
+          0.00202825 0.00189845 0.00201022 0.00188042 0.00176347 0.00174708 0.0016277 0.00161515 0.00149157 0.001481 0.00135564 0.00134763
+          0.0012197 0.00121303 0.00108395 0.00107891 0.000948233 0.000944082 0.000812647 0.000809563 0.000677092 0.000674617 0.000541623 0.000539868
+          0.000406168 0.00040487 0.000270773 0.000269983 0.000135381 0.000134971 0 0 0.00203199 0.0018963 0.00176091 0.00162528
+          0.00148976 0.00135418 0.00121865 0.00108312 0.000947646 0.000812194 0.000676771 0.000541382 0.000406002 0.000270662 0.000135305 0
+          0.00202292 0.0018871 0.00175265 0.0016179 0.00148338 0.0013488 0.0012141 0.00107943 0.000944579 0.000809795 0.000674848 0.000539967
+          0.000404962 0.000270015 0.000134982 0 0.00203094 0.0018954 0.00175992 0.00162432 0.00148891 0.00135337 0.00121801 0.00108253
+          0.000947205 0.0008118 0.000676499 0.000541152 0.000405865 0.00027056 0.000135277 0 0.00202651 0.00189078 0.00175577 0.00162042
+          0.00148549 0.00135034 0.0012154 0.00108033 0.000945357 0.000810299 0.00067529 0.00054023 0.000405192 0.000270125 0.000135071 0
+          0.00202979 0.00189428 0.00175893 0.00162342 0.00148811 0.00135265 0.00121738 0.00108198 0.00094675 0.000811406 0.000676199 0.000540901
+          0.000405705 0.000270439 0.000135247 0 0.002028 0.00189248 0.00175729 0.00162188 0.00148676 0.00135144 0.00121635 0.00108108
+          0.000945996 0.00081078 0.000675692 0.000540508 0.000405414 0.000270251 0.00013515 0 0.00202825 0.00189274 0.00175774 0.00162223
+          0.00148705 0.00135165 0.00121651 0.00108121 0.0009461 0.000810865 0.00067576 0.000540562 0.000405453 0.000270278 0.000135152 0
         </DataArray>
         <DataArray type="Float32" Name="x_gas^H2O" NumberOfComponents="1" format="ascii">
-          0.0302368 0.0302319 0.0127981 0.0128109 0.0302266 0.0128237 0.0302229 0.0128365 0.0302207 0.0128493 0.0302199 0.012862
-          0.03022 0.0128748 0.0302208 0.0128875 0.0302222 0.0129002 0.0302239 0.0129128 0.0302258 0.0129255 0.0302277 0.0129382
-          0.0302296 0.0129509 0.0302315 0.0129636 0.0302335 0.0129763 0.0302431 0.0129827 0.0127981 0.0127981 0.0127981 0.0127981
+          0.0304051 0.0304002 0.0127981 0.012811 0.0303949 0.012824 0.0303913 0.0128369 0.0303892 0.0128498 0.0303884 0.0128627
+          0.0303885 0.0128755 0.0303894 0.0128884 0.0303908 0.0129012 0.0303925 0.012914 0.0303944 0.0129268 0.0303963 0.0129396
+          0.0303982 0.0129524 0.0304 0.0129653 0.030402 0.0129781 0.0304114 0.0129845 0.0127981 0.0127981 0.0127981 0.0127981
           0.0127981 0.0127981 0.0127982 0.0127982 0.0127982 0.0127983 0.0127983 0.0127983 0.0127984 0.0127984 0.0127985 0.0127985
           0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981
           0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981
@@ -56,27 +56,27 @@
           0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981
           0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981
           0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981
-          0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127998 0.0128018 0.0128035
-          0.0128049 0.0128061 0.012807 0.0128077 0.0128084 0.0128089 0.0128093 0.0128096 0.0128098 0.01281 0.01281 0.0127981
+          0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127998 0.0128017 0.0128034
+          0.0128048 0.0128059 0.0128068 0.0128075 0.0128082 0.0128086 0.012809 0.0128094 0.0128096 0.0128097 0.0128097 0.0127981
         </DataArray>
         <DataArray type="Float32" Name="X_gas^H2O" NumberOfComponents="1" format="ascii">
-          0.0190269 0.0190237 0.008 0.00800808 0.0190204 0.00801613 0.019018 0.00802417 0.0190166 0.00803219 0.0190161 0.00804019
-          0.0190162 0.00804818 0.0190167 0.00805617 0.0190176 0.00806414 0.0190187 0.0080721 0.0190199 0.00808006 0.0190211 0.00808803
-          0.0190223 0.00809601 0.0190235 0.008104 0.0190248 0.00811199 0.0190309 0.00811598 0.008 0.008 0.00800001 0.00800002
-          0.00800003 0.00800004 0.00800006 0.00800007 0.0080001 0.00800012 0.00800014 0.00800017 0.0080002 0.00800024 0.00800027 0.00800029
+          0.019134 0.0191309 0.008 0.00800816 0.0191275 0.00801629 0.0191252 0.0080244 0.0191239 0.0080325 0.0191234 0.00804059
+          0.0191234 0.00804866 0.019124 0.00805672 0.0191249 0.00806477 0.019126 0.00807281 0.0191272 0.00808085 0.0191284 0.0080889
+          0.0191296 0.00809696 0.0191308 0.00810503 0.0191321 0.00811309 0.019138 0.00811713 0.008 0.008 0.00800001 0.00800002
+          0.00800003 0.00800004 0.00800006 0.00800007 0.0080001 0.00800012 0.00800015 0.00800018 0.00800021 0.00800024 0.00800028 0.00800029
           0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
           0.008 0.008 0.008 0.00799999 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
           0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
           0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
           0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
           0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
-          0.008 0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.008 0.0080011 0.00800233 0.00800341
-          0.0080043 0.00800502 0.0080056 0.00800608 0.00800647 0.00800678 0.00800704 0.00800724 0.00800738 0.00800747 0.00800749 0.008
+          0.008 0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.008 0.00800108 0.00800229 0.00800335
+          0.00800422 0.00800492 0.00800549 0.00800596 0.00800634 0.00800665 0.0080069 0.00800709 0.00800723 0.00800732 0.00800734 0.008
         </DataArray>
         <DataArray type="Float32" Name="x_gas^Air" NumberOfComponents="1" format="ascii">
-          0.969763 0.969768 0.987202 0.987189 0.969773 0.987176 0.969777 0.987163 0.969779 0.987151 0.96978 0.987138
-          0.96978 0.987125 0.969779 0.987113 0.969778 0.9871 0.969776 0.987087 0.969774 0.987074 0.969772 0.987062
-          0.96977 0.987049 0.969769 0.987036 0.969766 0.987024 0.969757 0.987017 0.987202 0.987202 0.987202 0.987202
+          0.969595 0.9696 0.987202 0.987189 0.969605 0.987176 0.969609 0.987163 0.969611 0.98715 0.969612 0.987137
+          0.969611 0.987125 0.969611 0.987112 0.969609 0.987099 0.969607 0.987086 0.969606 0.987073 0.969604 0.98706
+          0.969602 0.987048 0.9696 0.987035 0.969598 0.987022 0.969589 0.987015 0.987202 0.987202 0.987202 0.987202
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987201
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202
@@ -85,12 +85,12 @@
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.9872 0.987198 0.987197
-          0.987195 0.987194 0.987193 0.987192 0.987192 0.987191 0.987191 0.98719 0.98719 0.98719 0.98719 0.987202
+          0.987195 0.987194 0.987193 0.987192 0.987192 0.987191 0.987191 0.987191 0.98719 0.98719 0.98719 0.987202
         </DataArray>
         <DataArray type="Float32" Name="X_gas^Air" NumberOfComponents="1" format="ascii">
-          0.980973 0.980976 0.992 0.991992 0.98098 0.991984 0.980982 0.991976 0.980983 0.991968 0.980984 0.99196
-          0.980984 0.991952 0.980983 0.991944 0.980982 0.991936 0.980981 0.991928 0.98098 0.99192 0.980979 0.991912
-          0.980978 0.991904 0.980977 0.991896 0.980975 0.991888 0.980969 0.991884 0.992 0.992 0.992 0.992
+          0.980866 0.980869 0.992 0.991992 0.980873 0.991984 0.980875 0.991976 0.980876 0.991967 0.980877 0.991959
+          0.980877 0.991951 0.980876 0.991943 0.980875 0.991935 0.980874 0.991927 0.980873 0.991919 0.980872 0.991911
+          0.98087 0.991903 0.980869 0.991895 0.980868 0.991887 0.980862 0.991883 0.992 0.992 0.992 0.992
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992
@@ -99,26 +99,26 @@
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.991999 0.991998 0.991997
-          0.991996 0.991995 0.991994 0.991994 0.991994 0.991993 0.991993 0.991993 0.991993 0.991993 0.991993 0.992
+          0.991996 0.991995 0.991995 0.991994 0.991994 0.991993 0.991993 0.991993 0.991993 0.991993 0.991993 0.992
         </DataArray>
         <DataArray type="Float32" Name="h" NumberOfComponents="1" format="ascii">
-          71297.5 71286.9 44886.1 44905.4 71275.5 44924.7 71267.7 44943.9 71263 44963 71261.2 44982.1
-          71261.4 45001.2 71263.2 45020.3 71266.1 45039.4 71269.8 45058.4 71273.9 45077.4 71278 45096.5
-          71282 45115.5 71286 45134.6 71290.4 45153.7 71311.1 45163.3 44886.1 44886.1 44886.1 44886.1
-          44886.2 44886.2 44886.2 44886.3 44886.3 44886.4 44886.4 44886.5 44886.6 44886.7 44886.8 44886.8
+          71659 71648.4 44886.1 44905.7 71637.1 44925.2 71629.4 44944.6 71624.8 44964 71623.1 44983.4
+          71623.4 45002.8 71625.3 45022.1 71628.3 45041.4 71632 45060.7 71636.2 45080 71640.2 45099.3
+          71644.2 45118.6 71648.1 45137.9 71652.5 45157.3 71672.6 45167 44886.1 44886.1 44886.1 44886.1
+          44886.2 44886.2 44886.2 44886.3 44886.3 44886.4 44886.5 44886.5 44886.6 44886.7 44886.8 44886.8
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
-          44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44888.8 44891.9 44894.5
-          44896.7 44898.5 44899.9 44901.1 44902.1 44902.9 44903.5 44904 44904.3 44904.6 44904.6 44886.1
+          44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44888.8 44891.8 44894.4
+          44896.5 44898.3 44899.7 44900.8 44901.8 44902.5 44903.1 44903.6 44904 44904.2 44904.2 44886.1
         </DataArray>
         <DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
-          1.15798 1.15799 1.16259 1.16258 1.15801 1.16258 1.15802 1.16258 1.15802 1.16257 1.15802 1.16257
-          1.15802 1.16257 1.15802 1.16257 1.15802 1.16256 1.15801 1.16256 1.15801 1.16256 1.158 1.16255
-          1.158 1.16255 1.15799 1.16255 1.15799 1.16254 1.15796 1.16254 1.16259 1.16259 1.16259 1.16259
+          1.15755 1.15756 1.16259 1.16258 1.15757 1.16258 1.15758 1.16258 1.15759 1.16257 1.15759 1.16257
+          1.15759 1.16257 1.15759 1.16256 1.15758 1.16256 1.15758 1.16256 1.15757 1.16255 1.15757 1.16255
+          1.15756 1.16255 1.15756 1.16254 1.15755 1.16254 1.15753 1.16254 1.16259 1.16259 1.16259 1.16259
           1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259
           1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259
           1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259
@@ -130,9 +130,9 @@
           1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16259
         </DataArray>
         <DataArray type="Float32" Name="mu" NumberOfComponents="1" format="ascii">
-          1.79268e-05 1.79267e-05 1.81222e-05 1.81221e-05 1.79266e-05 1.81219e-05 1.79265e-05 1.81218e-05 1.79265e-05 1.81216e-05 1.79265e-05 1.81215e-05
-          1.79265e-05 1.81214e-05 1.79265e-05 1.81212e-05 1.79265e-05 1.81211e-05 1.79265e-05 1.81209e-05 1.79266e-05 1.81208e-05 1.79266e-05 1.81206e-05
-          1.79266e-05 1.81205e-05 1.79267e-05 1.81203e-05 1.79267e-05 1.81202e-05 1.79269e-05 1.81201e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05
+          1.79296e-05 1.79295e-05 1.81222e-05 1.81221e-05 1.79294e-05 1.81219e-05 1.79294e-05 1.81218e-05 1.79293e-05 1.81217e-05 1.79293e-05 1.81215e-05
+          1.79293e-05 1.81214e-05 1.79293e-05 1.81212e-05 1.79294e-05 1.81211e-05 1.79294e-05 1.81209e-05 1.79294e-05 1.81208e-05 1.79294e-05 1.81207e-05
+          1.79295e-05 1.81205e-05 1.79295e-05 1.81204e-05 1.79295e-05 1.81202e-05 1.79297e-05 1.81202e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05
           1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05
           1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05
           1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05
@@ -144,38 +144,38 @@
           1.81222e-05 1.81222e-05 1.81222e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81222e-05
         </DataArray>
         <DataArray type="Float32" Name="v" NumberOfComponents="3" format="ascii">
-          0.000913725 1.24857e-05 0 0.000915844 1.31614e-05 0 1.11851 0 0 1.11893 -0.000212801 0
-          0.000917975 1.31999e-05 0 1.1197 -0.000226546 0 0.000918633 1.31479e-05 0 1.12031 -0.000177815 0
-          0.00091916 1.30915e-05 0 1.12085 -0.000139228 0 0.000919439 1.30219e-05 0 1.12123 -9.917e-05 0
-          0.000919677 1.29714e-05 0 1.12154 -7.26501e-05 0 0.000919808 1.29215e-05 0 1.12175 -4.90274e-05 0
-          0.000919923 1.28844e-05 0 1.12193 -3.35985e-05 0 0.000919986 1.28482e-05 0 1.12205 -2.03253e-05 0
-          0.000920041 1.28199e-05 0 1.12214 -1.15699e-05 0 0.00092007 1.28089e-05 0 1.12221 -4.239e-06 0
-          0.000920096 1.27999e-05 0 1.12226 8.88714e-07 0 0.000920103 1.27808e-05 0 1.12229 5.22126e-06 0
-          0.000920121 1.28486e-05 0 1.12231 9.30377e-06 0 0.000920109 1.24843e-05 0 1.12232 1.02407e-05 0
-          2.1001 0 0 2.10015 -0.0014171 0 2.10037 -0.00155935 0 2.10071 -0.0012866 0
-          2.10104 -0.00104255 0 2.10134 -0.000781599 0 2.1016 -0.000600314 0 2.1018 -0.000438247 0
-          2.10196 -0.000328567 0 2.10208 -0.000234721 0 2.10217 -0.000170839 0 2.10224 -0.000116801 0
-          2.10229 -7.85714e-05 0 2.10233 -4.58596e-05 0 2.10235 -1.96306e-05 0 2.10236 3.55784e-06 0
-          2.88064 0 0 2.88051 -0.000678473 0 2.88031 -0.00109956 0 2.88029 -0.00117452 0
-          2.88035 -0.00107508 0 2.88048 -0.000911039 0 2.88062 -0.00074125 0 2.88076 -0.000583201 0
-          2.88089 -0.000451297 0 2.881 -0.000340476 0 2.8811 -0.000253855 0 2.88117 -0.000183064 0
-          2.88122 -0.000128311 0 2.88126 -8.19835e-05 0 2.88129 -4.33622e-05 0 2.8813 -5.6877e-06 0
-          3.37718 0 0 3.37714 -0.00093452 0 3.37702 -0.001318 0 3.37688 -0.00142967 0
-          3.37677 -0.00137673 0 3.3767 -0.00122482 0 3.37667 -0.00104628 0 3.37666 -0.000856242 0
-          3.37666 -0.000688328 0 3.37668 -0.000535048 0 3.37669 -0.000410736 0 3.3767 -0.000302732 0
-          3.37671 -0.00021677 0 3.37672 -0.0001409 0 3.37673 -7.66445e-05 0 3.37673 -1.33419e-05 0
-          3.48529 0 0 3.4852 -0.00027326 0 3.485 -0.000566851 0 3.48479 -0.000733801 0
-          3.48458 -0.000789047 0 3.48441 -0.000764817 0 3.48427 -0.000692255 0 3.48417 -0.00059635 0
-          3.48409 -0.000494753 0 3.48404 -0.000397058 0 3.484 -0.000309984 0 3.48397 -0.000233629 0
-          3.48395 -0.000169464 0 3.48394 -0.000113007 0 3.48394 -6.39343e-05 0 3.48393 -1.52791e-05 0
-          3.07414 0 0 3.07409 -0.000364916 0 3.07397 -0.000530221 0 3.0738 -0.000622876 0
-          3.07361 -0.000651474 0 3.07341 -0.000632812 0 3.07323 -0.000581394 0 3.07308 -0.000509791 0
-          3.07294 -0.000432188 0 3.07284 -0.000353048 0 3.07275 -0.000280871 0 3.07268 -0.000214126 0
-          3.07263 -0.000156489 0 3.0726 -0.000103807 0 3.07257 -5.67097e-05 0 3.07256 -1.04414e-05 0
-          1.98047 0 0 1.98042 0.000143717 0 1.98031 0.000120605 0 1.98018 8.16373e-05 0
-          1.98003 4.56076e-05 0 1.97989 2.23357e-05 0 1.97975 7.75314e-06 0 1.97962 1.09233e-06 0
-          1.97951 -1.11118e-06 0 1.97942 -1.036e-06 0 1.97935 4.54722e-08 0 1.9793 5.70455e-07 0
-          1.97925 8.07197e-07 0 1.97922 -6.03331e-08 0 1.9792 -9.97972e-07 0 1.9792 -3.51497e-06 0
+          0.000913702 1.26213e-05 0 0.000915789 1.32841e-05 0 1.11851 0 0 1.11892 -0.000209087 0
+          0.000917887 1.33234e-05 0 1.11968 -0.000222597 0 0.000918533 1.32731e-05 0 1.12029 -0.000174603 0
+          0.000919051 1.32183e-05 0 1.12081 -0.000136623 0 0.000919326 1.31501e-05 0 1.12119 -9.71949e-05 0
+          0.00091956 1.31005e-05 0 1.12149 -7.11025e-05 0 0.000919689 1.3051e-05 0 1.1217 -4.78628e-05 0
+          0.000919803 1.30138e-05 0 1.12187 -3.26939e-05 0 0.000919864 1.29769e-05 0 1.12199 -1.96483e-05 0
+          0.000919919 1.29486e-05 0 1.12208 -1.10618e-05 0 0.000919947 1.29372e-05 0 1.12215 -3.87895e-06 0
+          0.000919972 1.2927e-05 0 1.1222 1.13979e-06 0 0.00091998 1.29071e-05 0 1.12223 5.36724e-06 0
+          0.000919997 1.29718e-05 0 1.12225 9.3302e-06 0 0.000919985 1.26166e-05 0 1.12226 1.02498e-05 0
+          2.1001 0 0 2.10015 -0.00139383 0 2.10037 -0.00153358 0 2.1007 -0.00126497 0
+          2.10103 -0.00102476 0 2.10132 -0.000767907 0 2.10157 -0.000589532 0 2.10177 -0.000430066 0
+          2.10193 -0.000322187 0 2.10205 -0.000229893 0 2.10214 -0.000167106 0 2.10221 -0.000114032 0
+          2.10226 -7.65407e-05 0 2.10229 -4.45356e-05 0 2.10231 -1.89405e-05 0 2.10233 3.54856e-06 0
+          2.88064 0 0 2.88051 -0.000667733 0 2.88032 -0.00108192 0 2.88029 -0.00115551 0
+          2.88035 -0.00105756 0 2.88048 -0.000896081 0 2.88062 -0.000729 0 2.88076 -0.00057349 0
+          2.88089 -0.000443734 0 2.881 -0.000334727 0 2.88109 -0.000249546 0 2.88116 -0.000179942 0
+          2.88121 -0.000126128 0 2.88125 -8.05891e-05 0 2.88128 -4.26262e-05 0 2.88129 -5.55978e-06 0
+          3.37718 0 0 3.37714 -0.000919092 0 3.37702 -0.00129624 0 3.37688 -0.00140597 0
+          3.37678 -0.00135378 0 3.37671 -0.00120422 0 3.37668 -0.00102853 0 3.37667 -0.000841551 0
+          3.37667 -0.000676384 0 3.37668 -0.00052563 0 3.3767 -0.000403408 0 3.37671 -0.000297247 0
+          3.37672 -0.000212794 0 3.37673 -0.000138284 0 3.37674 -7.52133e-05 0 3.37674 -1.30958e-05 0
+          3.48529 0 0 3.4852 -0.000269103 0 3.48501 -0.000557919 0 3.4848 -0.000722098 0
+          3.48459 -0.000776375 0 3.48442 -0.000752475 0 3.48429 -0.000681041 0 3.48419 -0.000586658 0
+          3.48411 -0.0004867 0 3.48406 -0.00039059 0 3.48402 -0.000304938 0 3.48399 -0.000229828 0
+          3.48398 -0.000166705 0 3.48397 -0.00011116 0 3.48396 -6.28707e-05 0 3.48396 -1.49938e-05 0
+          3.07414 0 0 3.0741 -0.00035881 0 3.07397 -0.000521416 0 3.0738 -0.000612538 0
+          3.07361 -0.000640632 0 3.07343 -0.000622229 0 3.07325 -0.000571618 0 3.0731 -0.000501164 0
+          3.07296 -0.000424832 0 3.07286 -0.000346995 0 3.07277 -0.000276021 0 3.07271 -0.000210399 0
+          3.07266 -0.000153745 0 3.07262 -0.000101971 0 3.0726 -5.56977e-05 0 3.07259 -1.0245e-05 0
+          1.98047 0 0 1.98042 0.000140974 0 1.98031 0.000118159 0 1.98018 7.98149e-05 0
+          1.98004 4.43909e-05 0 1.9799 2.15179e-05 0 1.97976 7.20266e-06 0 1.97964 6.77249e-07 0
+          1.97953 -1.45598e-06 0 1.97944 -1.34568e-06 0 1.97937 -2.37399e-07 0 1.97932 3.25492e-07 0
+          1.97927 6.11835e-07 0 1.97924 -1.8378e-07 0 1.97923 -1.04098e-06 0 1.97922 -3.44668e-06 0
           0.00134 0 0 0.00134 0 0 0.00134 0 0 0.00134 0 0
           0.00134 0 0 0.00134 0 0 0.00134 0 0 0.00134 0 0
           0.00134 0 0 0.00134 0 0 0.00134 0 0 0.00134 0 0
diff --git a/test/references/2cnistokes2p2cniboundarylayer-ff-reference.vtu b/test/references/2cnistokes2p2cniboundarylayer-ff-reference.vtu
index a3a34cac8c..2dc94521b7 100644
--- a/test/references/2cnistokes2p2cniboundarylayer-ff-reference.vtu
+++ b/test/references/2cnistokes2p2cniboundarylayer-ff-reference.vtu
@@ -4,9 +4,9 @@
     <Piece NumberOfCells="120" NumberOfPoints="144">
       <PointData Scalars="temperature" Vectors="v">
         <DataArray type="Float32" Name="temperature" NumberOfComponents="1" format="ascii">
-          296.449 296.445 298.15 298.149 296.441 298.148 296.437 298.146 296.432 298.145 296.428 298.144
-          296.424 298.143 296.42 298.141 296.416 298.14 296.412 298.139 296.409 298.138 296.406 298.136
-          296.403 298.135 296.4 298.134 296.398 298.133 296.396 298.132 298.15 298.15 298.15 298.15
+          296.027 296.025 298.15 298.149 296.023 298.147 296.021 298.146 296.018 298.144 296.016 298.143
+          296.013 298.141 296.011 298.14 296.009 298.138 296.007 298.137 296.005 298.135 296.004 298.134
+          296.002 298.132 296.001 298.13 296 298.129 295.999 298.128 298.15 298.15 298.15 298.15
           298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15
           298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15
           298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15 298.15
@@ -32,51 +32,51 @@
           100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
         </DataArray>
         <DataArray type="Float32" Name="delP" NumberOfComponents="1" format="ascii">
-          0.00202834 0.00189957 0.00200667 0.0018779 0.00176431 0.00174487 0.00162856 0.00161358 0.00149221 0.00147959 0.00135619 0.0013466
-          0.0012201 0.0012121 0.00108429 0.00107822 0.000948479 0.000943453 0.000812857 0.000809096 0.000677236 0.000674203 0.000541745 0.000539572
-          0.000406247 0.000404626 0.000270832 0.000269834 0.000135405 0.000134884 0 0 0.00203252 0.00189681 0.00176133 0.00162569
-          0.00149008 0.00135447 0.00121887 0.00108332 0.000947789 0.000812322 0.000676859 0.00054146 0.000406052 0.000270702 0.000135323 0
-          0.00202165 0.00188575 0.00175146 0.00161682 0.00148243 0.00134799 0.00121339 0.00107885 0.000944072 0.000809395 0.000674507 0.000539713
-          0.000404763 0.000269891 0.000134915 0 0.00203123 0.00189566 0.00176013 0.0016245 0.00148905 0.00135348 0.00121809 0.00108259
-          0.000947247 0.000811832 0.00067652 0.000541167 0.000405875 0.000270567 0.00013528 0 0.00202591 0.00189012 0.00175516 0.00161981
-          0.00148494 0.00134983 0.00121495 0.00107993 0.000945012 0.000810005 0.000675048 0.000540038 0.000405049 0.00027003 0.000135025 0
-          0.00202984 0.00189429 0.00175893 0.0016234 0.00148808 0.00135259 0.00121733 0.00108192 0.000946696 0.000811348 0.000676154 0.000540857
-          0.000405677 0.000270417 0.00013524 0 0.00202769 0.00189214 0.00175697 0.00162154 0.00148646 0.00135114 0.00121608 0.00108083
-          0.000945783 0.000810586 0.000675537 0.000540376 0.000405321 0.000270184 0.00013512 0 0.002028 0.00189245 0.00175751 0.00162198
-          0.00148682 0.00135141 0.0012163 0.001081 0.000945921 0.0008107 0.000675631 0.000540452 0.000405378 0.000270226 0.000135135 0
+          0.00202833 0.00189996 0.00200513 0.00187677 0.00176475 0.00174389 0.00162895 0.00161284 0.00149254 0.00147896 0.00135647 0.00134611
+          0.00122033 0.00121166 0.00108448 0.00107787 0.000948635 0.000943152 0.000812985 0.000808856 0.000677342 0.000673999 0.000541831 0.000539415
+          0.000406315 0.000404502 0.00027088 0.000269753 0.000135441 0.000134836 0 0 0.00203277 0.00189703 0.00176155 0.00162587
+          0.00149025 0.00135462 0.001219 0.00108343 0.000947882 0.000812399 0.000676924 0.000541512 0.000406092 0.00027073 0.000135338 0
+          0.00202109 0.00188515 0.00175093 0.00161633 0.00148201 0.00134762 0.00121306 0.00107857 0.000943836 0.0008092 0.000674345 0.000539587
+          0.000404667 0.000269828 0.000134881 0 0.00203137 0.00189578 0.00176024 0.00162458 0.00148913 0.00135353 0.00121814 0.00108263
+          0.000947282 0.000811857 0.000676542 0.000541182 0.000405887 0.000270575 0.000135284 0 0.00202564 0.00188982 0.00175488 0.00161953
+          0.0014847 0.00134959 0.00121475 0.00107974 0.000944853 0.000809867 0.000674936 0.000539946 0.000404982 0.000269984 0.000135003 0
+          0.00202986 0.00189431 0.00175894 0.00162339 0.00148807 0.00135258 0.00121732 0.00108189 0.000946679 0.000811326 0.000676139 0.000540841
+          0.000405668 0.000270408 0.000135238 0 0.00202755 0.00189198 0.00175682 0.00162139 0.00148633 0.001351 0.00121596 0.00108072
+          0.000945687 0.000810497 0.000675466 0.000540316 0.000405278 0.000270154 0.000135107 0 0.00202789 0.00189232 0.00175742 0.00162186
+          0.00148671 0.0013513 0.0012162 0.00108091 0.000945843 0.000810626 0.000675573 0.000540402 0.000405343 0.000270201 0.000135123 0
         </DataArray>
         <DataArray type="Float32" Name="x_gas^H2O" NumberOfComponents="1" format="ascii">
-          0.0286211 0.0286154 0.0127981 0.0128097 0.028608 0.0128213 0.0286001 0.0128329 0.0285923 0.0128445 0.0285848 0.012856
-          0.0285777 0.0128676 0.0285709 0.0128793 0.0285645 0.0128911 0.0285584 0.012903 0.0285525 0.0129149 0.028547 0.0129268
-          0.0285419 0.0129387 0.0285371 0.0129507 0.0285331 0.0129626 0.0285303 0.0129686 0.0127981 0.0127981 0.0127981 0.0127981
-          0.0127981 0.0127981 0.0127981 0.0127982 0.0127982 0.0127982 0.0127983 0.0127983 0.0127984 0.0127984 0.0127985 0.0127985
+          0.0279006 0.0278977 0.0127981 0.0128092 0.0278937 0.0128203 0.0278894 0.0128313 0.0278851 0.0128423 0.027881 0.0128533
+          0.0278771 0.0128644 0.0278735 0.0128756 0.0278701 0.0128868 0.0278669 0.0128981 0.0278639 0.0129095 0.027861 0.0129208
+          0.0278583 0.0129322 0.0278559 0.0129436 0.0278539 0.012955 0.0278525 0.0129607 0.0127981 0.0127981 0.0127981 0.0127981
+          0.0127981 0.0127981 0.0127981 0.0127982 0.0127982 0.0127982 0.0127983 0.0127983 0.0127983 0.0127984 0.0127984 0.0127985
           0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981
           0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981
           0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981
           0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981
           0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981
           0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981
-          0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0128001 0.0128025 0.0128045
-          0.0128062 0.0128076 0.0128087 0.0128097 0.0128104 0.0128111 0.0128116 0.012812 0.0128123 0.0128125 0.0128126 0.0127981
+          0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0127981 0.0128003 0.0128028 0.012805
+          0.0128069 0.0128084 0.0128097 0.0128107 0.0128115 0.0128122 0.0128128 0.0128133 0.0128136 0.0128138 0.0128139 0.0127981
         </DataArray>
         <DataArray type="Float32" Name="X_gas^H2O" NumberOfComponents="1" format="ascii">
-          0.0179991 0.0179954 0.008 0.00800732 0.0179907 0.00801461 0.0179857 0.00802189 0.0179807 0.00802915 0.017976 0.00803639
-          0.0179714 0.00804369 0.0179671 0.00805105 0.017963 0.00805847 0.0179592 0.00806592 0.0179555 0.00807339 0.017952 0.00808088
-          0.0179487 0.00808837 0.0179457 0.00809587 0.0179431 0.00810337 0.0179413 0.00810712 0.008 0.008 0.00800001 0.00800002
-          0.00800002 0.00800004 0.00800005 0.00800007 0.00800009 0.00800011 0.00800013 0.00800016 0.00800019 0.00800022 0.00800025 0.00800026
+          0.0175411 0.0175393 0.008 0.00800698 0.0175367 0.00801394 0.017534 0.00802089 0.0175313 0.00802781 0.0175287 0.00803473
+          0.0175262 0.00804167 0.0175239 0.00804868 0.0175218 0.00805575 0.0175197 0.00806285 0.0175178 0.00806998 0.017516 0.00807712
+          0.0175143 0.00808428 0.0175127 0.00809144 0.0175114 0.0080986 0.0175106 0.00810218 0.008 0.008 0.00800001 0.00800001
+          0.00800002 0.00800003 0.00800005 0.00800006 0.00800008 0.0080001 0.00800013 0.00800015 0.00800018 0.00800021 0.00800024 0.00800025
           0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
           0.008 0.008 0.008 0.00799999 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
           0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
           0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
           0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
-          0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
-          0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.00800002 0.008 0.0080013 0.00800276 0.00800405
-          0.00800512 0.00800599 0.0080067 0.00800729 0.00800777 0.00800817 0.00800849 0.00800875 0.00800894 0.00800906 0.00800911 0.008
+          0.008 0.008 0.008 0.00800001 0.008 0.008 0.008 0.008 0.008 0.008 0.008 0.008
+          0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.00800001 0.00800002 0.008 0.0080014 0.00800299 0.00800438
+          0.00800554 0.0080065 0.00800728 0.00800793 0.00800846 0.00800891 0.00800927 0.00800955 0.00800976 0.0080099 0.00800995 0.008
         </DataArray>
         <DataArray type="Float32" Name="x_gas^Air" NumberOfComponents="1" format="ascii">
-          0.971379 0.971385 0.987202 0.98719 0.971392 0.987179 0.9714 0.987167 0.971408 0.987156 0.971415 0.987144
-          0.971422 0.987132 0.971429 0.987121 0.971436 0.987109 0.971442 0.987097 0.971447 0.987085 0.971453 0.987073
-          0.971458 0.987061 0.971463 0.987049 0.971467 0.987037 0.97147 0.987031 0.987202 0.987202 0.987202 0.987202
+          0.972099 0.972102 0.987202 0.987191 0.972106 0.98718 0.972111 0.987169 0.972115 0.987158 0.972119 0.987147
+          0.972123 0.987136 0.972126 0.987124 0.97213 0.987113 0.972133 0.987102 0.972136 0.987091 0.972139 0.987079
+          0.972142 0.987068 0.972144 0.987056 0.972146 0.987045 0.972147 0.987039 0.987202 0.987202 0.987202 0.987202
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202
@@ -84,13 +84,13 @@
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202
           0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202
-          0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.9872 0.987198 0.987195
-          0.987194 0.987192 0.987191 0.98719 0.98719 0.987189 0.987188 0.987188 0.987188 0.987188 0.987187 0.987202
+          0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.987202 0.9872 0.987197 0.987195
+          0.987193 0.987192 0.98719 0.987189 0.987188 0.987188 0.987187 0.987187 0.987186 0.987186 0.987186 0.987202
         </DataArray>
         <DataArray type="Float32" Name="X_gas^Air" NumberOfComponents="1" format="ascii">
-          0.982001 0.982005 0.992 0.991993 0.982009 0.991985 0.982014 0.991978 0.982019 0.991971 0.982024 0.991964
-          0.982029 0.991956 0.982033 0.991949 0.982037 0.991942 0.982041 0.991934 0.982045 0.991927 0.982048 0.991919
-          0.982051 0.991912 0.982054 0.991904 0.982057 0.991897 0.982059 0.991893 0.992 0.992 0.992 0.992
+          0.982459 0.982461 0.992 0.991993 0.982463 0.991986 0.982466 0.991979 0.982469 0.991972 0.982471 0.991965
+          0.982474 0.991958 0.982476 0.991951 0.982478 0.991944 0.98248 0.991937 0.982482 0.99193 0.982484 0.991923
+          0.982486 0.991916 0.982487 0.991909 0.982489 0.991901 0.982489 0.991898 0.992 0.992 0.992 0.992
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992
@@ -99,26 +99,26 @@
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992
           0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.992 0.991999 0.991997 0.991996
-          0.991995 0.991994 0.991993 0.991993 0.991992 0.991992 0.991992 0.991991 0.991991 0.991991 0.991991 0.992
+          0.991994 0.991993 0.991993 0.991992 0.991992 0.991991 0.991991 0.99199 0.99199 0.99199 0.99199 0.992
         </DataArray>
         <DataArray type="Float32" Name="h" NumberOfComponents="1" format="ascii">
-          67805.9 67793.4 44886.1 44902.9 67777.3 44919.7 67760.1 44936.5 67743.2 44953.1 67726.9 44969.8
-          67711.4 44986.6 67696.7 45003.5 67682.7 45020.5 67669.4 45037.6 67656.8 45054.8 67644.8 45072
-          67633.6 45089.2 67623.3 45106.4 67614.4 45123.6 67608.5 45132.2 44886.1 44886.1 44886.1 44886.1
-          44886.2 44886.2 44886.2 44886.3 44886.3 44886.4 44886.4 44886.5 44886.5 44886.6 44886.7 44886.7
+          66235.7 66229.3 44886.1 44901.8 66220.5 44917.5 66211.1 44933.2 66201.7 44948.8 66192.8 44964.3
+          66184.4 44980 66176.5 44995.8 66169 45011.7 66162 45027.7 66155.4 45043.7 66149.1 45059.8
+          66143.3 45075.9 66138 45092 66133.5 45108.1 66130.6 45116.2 44886.1 44886.1 44886.1 44886.1
+          44886.2 44886.2 44886.2 44886.2 44886.3 44886.3 44886.4 44886.4 44886.5 44886.6 44886.7 44886.6
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
           44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1
-          44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.2 44886.1 44889.3 44892.9 44896.1
-          44898.7 44900.9 44902.7 44904.1 44905.3 44906.3 44907.1 44907.7 44908.2 44908.5 44908.6 44886.1
+          44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.1 44886.2 44886.1 44889.6 44893.5 44896.9
+          44899.8 44902.1 44904.1 44905.7 44907 44908.1 44909 44909.7 44910.2 44910.6 44910.7 44886.1
         </DataArray>
         <DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
-          1.16226 1.16227 1.16259 1.16259 1.16229 1.16259 1.16232 1.16259 1.16234 1.16259 1.16236 1.16259
-          1.16238 1.16259 1.1624 1.16258 1.16241 1.16258 1.16243 1.16258 1.16245 1.16258 1.16246 1.16258
-          1.16247 1.16258 1.16249 1.16258 1.1625 1.16258 1.16251 1.16258 1.16259 1.16259 1.16259 1.16259
+          1.16423 1.16424 1.16259 1.16259 1.16425 1.16259 1.16426 1.16259 1.16428 1.16259 1.16429 1.16259
+          1.1643 1.16259 1.16431 1.16259 1.16432 1.16259 1.16433 1.1626 1.16433 1.1626 1.16434 1.1626
+          1.16435 1.1626 1.16436 1.1626 1.16436 1.1626 1.16437 1.1626 1.16259 1.16259 1.16259 1.16259
           1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259
           1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259
           1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259 1.16259
@@ -130,9 +130,9 @@
           1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16258 1.16259
         </DataArray>
         <DataArray type="Float32" Name="mu" NumberOfComponents="1" format="ascii">
-          1.78983e-05 1.78982e-05 1.81222e-05 1.81221e-05 1.78981e-05 1.81219e-05 1.78979e-05 1.81217e-05 1.78978e-05 1.81216e-05 1.78977e-05 1.81214e-05
-          1.78975e-05 1.81213e-05 1.78974e-05 1.81211e-05 1.78973e-05 1.81209e-05 1.78972e-05 1.81208e-05 1.78971e-05 1.81206e-05 1.7897e-05 1.81204e-05
-          1.78969e-05 1.81203e-05 1.78968e-05 1.81201e-05 1.78967e-05 1.81199e-05 1.78967e-05 1.81198e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05
+          1.78849e-05 1.78849e-05 1.81222e-05 1.81221e-05 1.78848e-05 1.81219e-05 1.78847e-05 1.81217e-05 1.78846e-05 1.81216e-05 1.78846e-05 1.81214e-05
+          1.78845e-05 1.81212e-05 1.78844e-05 1.8121e-05 1.78844e-05 1.81209e-05 1.78843e-05 1.81207e-05 1.78842e-05 1.81205e-05 1.78842e-05 1.81203e-05
+          1.78841e-05 1.81202e-05 1.78841e-05 1.812e-05 1.78841e-05 1.81198e-05 1.7884e-05 1.81197e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05
           1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05
           1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05
           1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05 1.81222e-05
@@ -144,38 +144,38 @@
           1.81222e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81221e-05 1.81222e-05
         </DataArray>
         <DataArray type="Float32" Name="v" NumberOfComponents="3" format="ascii">
-          0.000913908 6.5757e-05 0 0.000916501 6.54864e-05 0 1.11851 0 0 1.11901 -0.000191569 0
-          0.000918941 6.54e-05 0 1.1199 -0.000216258 0 0.000919633 6.49947e-05 0 1.12061 -0.000160287 0
-          0.000920281 6.47231e-05 0 1.12124 -0.000117459 0 0.000920591 6.44368e-05 0 1.12168 -7.21678e-05 0
-          0.000920879 6.43523e-05 0 1.12204 -4.21656e-05 0 0.000921032 6.42917e-05 0 1.12229 -1.51763e-05 0
-          0.000921176 6.42294e-05 0 1.12249 2.45076e-06 0 0.000921251 6.41697e-05 0 1.12263 1.77885e-05 0
-          0.000921324 6.41124e-05 0 1.12275 2.80133e-05 0 0.000921361 6.40592e-05 0 1.12282 3.69028e-05 0
-          0.000921399 6.40088e-05 0 1.12289 4.32361e-05 0 0.000921412 6.39898e-05 0 1.12292 4.8873e-05 0
-          0.000921439 6.38708e-05 0 1.12295 5.41611e-05 0 0.000921427 6.44729e-05 0 1.12296 5.71181e-05 0
-          2.1001 0 0 2.10017 -0.00161615 0 2.10044 -0.00177051 0 2.10082 -0.00145144 0
-          2.10122 -0.00117062 0 2.10157 -0.00086881 0 2.10186 -0.000660704 0 2.1021 -0.000473429 0
-          2.10229 -0.000347347 0 2.10243 -0.000238597 0 2.10254 -0.000164715 0 2.10262 -0.00010143 0
-          2.10268 -5.64272e-05 0 2.10273 -1.69298e-05 0 2.10276 1.53434e-05 0 2.10277 4.52532e-05 0
-          2.88064 0 0 2.88049 -0.000762778 0 2.88027 -0.00124236 0 2.88025 -0.00132361 0
-          2.88032 -0.00120854 0 2.88047 -0.00101843 0 2.88064 -0.000823411 0 2.88081 -0.000641011 0
-          2.88097 -0.000489413 0 2.8811 -0.000361333 0 2.88121 -0.000261367 0 2.88129 -0.000179054 0
-          2.88136 -0.000115242 0 2.8814 -6.07998e-05 0 2.88144 -1.51889e-05 0 2.88145 2.91898e-05 0
-          3.37718 0 0 3.37713 -0.00107484 0 3.377 -0.00151097 0 3.37684 -0.00163489 0
-          3.37672 -0.00157115 0 3.37665 -0.00139383 0 3.37661 -0.00118724 0 3.37661 -0.000967113 0
-          3.37662 -0.000773212 0 3.37663 -0.000595612 0 3.37665 -0.000451692 0 3.37667 -0.000325974 0
-          3.37668 -0.000225687 0 3.3767 -0.000136556 0 3.37671 -6.07348e-05 0 3.37672 1.41093e-05 0
-          3.48529 0 0 3.48519 -0.000307139 0 3.48496 -0.000642364 0 3.48472 -0.000831702 0
-          3.48448 -0.000893871 0 3.48428 -0.000864381 0 3.48412 -0.000780326 0 3.48401 -0.000669104 0
-          3.48392 -0.000551963 0 3.48387 -0.000438932 0 3.48382 -0.000338417 0 3.4838 -0.000249861 0
-          3.48378 -0.000175408 0 3.48377 -0.000109571 0 3.48377 -5.22398e-05 0 3.48377 4.67529e-06 0
-          3.07414 0 0 3.07409 -0.000422106 0 3.07395 -0.00061078 0 3.07375 -0.000715606 0
-          3.07353 -0.000746987 0 3.07331 -0.000724133 0 3.07311 -0.000663989 0 3.07293 -0.000580614 0
-          3.07278 -0.000490634 0 3.07265 -0.000398726 0 3.07256 -0.000314977 0 3.07248 -0.000237222 0
-          3.07243 -0.00016993 0 3.07239 -0.000108037 0 3.07237 -5.24447e-05 0 3.07236 2.22957e-06 0
-          1.98047 0 0 1.98041 0.000171143 0 1.98028 0.000145652 0 1.98013 0.000101384 0
-          1.97996 5.98385e-05 0 1.9798 3.31558e-05 0 1.97964 1.60532e-05 0 1.9795 8.26939e-06 0
-          1.97937 5.32146e-06 0 1.97927 5.17517e-06 0 1.97919 5.99339e-06 0 1.97913 6.32331e-06 0
-          1.97908 6.22937e-06 0 1.97905 5.03645e-06 0 1.97903 3.91182e-06 0 1.97902 1.77727e-06 0
+          0.000914023 6.2468e-05 0 0.000916759 6.27426e-05 0 1.11851 0 0 1.11904 -0.000208196 0
+          0.000919355 6.26682e-05 0 1.11998 -0.000233126 0 0.000920098 6.22975e-05 0 1.12074 -0.000173526 0
+          0.00092078 6.20333e-05 0 1.1214 -0.00012783 0 0.000921109 6.17618e-05 0 1.12187 -7.96635e-05 0
+          0.000921411 6.16272e-05 0 1.12225 -4.78658e-05 0 0.000921571 6.1582e-05 0 1.12252 -1.92743e-05 0
+          0.000921721 6.15216e-05 0 1.12273 -6.39622e-07 0 0.000921799 6.14656e-05 0 1.12288 1.5553e-05 0
+          0.000921874 6.1411e-05 0 1.123 2.63232e-05 0 0.00092191 6.13592e-05 0 1.12308 3.56722e-05 0
+          0.000921949 6.13096e-05 0 1.12314 4.23078e-05 0 0.000921959 6.12775e-05 0 1.12318 4.81534e-05 0
+          0.000921987 6.11883e-05 0 1.12321 5.37155e-05 0 0.00092198 6.14791e-05 0 1.12322 5.54225e-05 0
+          2.1001 0 0 2.10017 -0.00172429 0 2.10045 -0.00188896 0 2.10086 -0.00154979 0
+          2.10128 -0.00125064 0 2.10165 -0.000929676 0 2.10196 -0.000708148 0 2.10221 -0.000509061 0
+          2.10241 -0.000374919 0 2.10256 -0.000259316 0 2.10268 -0.000180695 0 2.10277 -0.000113331 0
+          2.10283 -6.53273e-05 0 2.10288 -2.30341e-05 0 2.10291 1.15648e-05 0 2.10292 4.40748e-05 0
+          2.88064 0 0 2.88048 -0.000809572 0 2.88024 -0.00131946 0 2.88022 -0.00140602 0
+          2.88029 -0.00128344 0 2.88045 -0.00108136 0 2.88063 -0.000874004 0 2.88081 -0.000680297 0
+          2.88098 -0.000519323 0 2.88111 -0.000383508 0 2.88123 -0.000277574 0 2.88132 -0.000190534 0
+          2.88139 -0.000123202 0 2.88144 -6.6033e-05 0 2.88147 -1.83534e-05 0 2.88149 2.77276e-05 0
+          3.37718 0 0 3.37713 -0.00114483 0 3.37698 -0.00160838 0 3.37682 -0.00174 0
+          3.37669 -0.00167196 0 3.37661 -0.00148338 0 3.37658 -0.00126363 0 3.37657 -0.0010297 0
+          3.37657 -0.00082359 0 3.37659 -0.000634948 0 3.37661 -0.000482057 0 3.37663 -0.000348576 0
+          3.37665 -0.000242098 0 3.37666 -0.000147494 0 3.37667 -6.70048e-05 0 3.37668 1.2469e-05 0
+          3.48529 0 0 3.48518 -0.000323428 0 3.48494 -0.000678991 0 3.48468 -0.000879959 0
+          3.48442 -0.000945896 0 3.48422 -0.000914613 0 3.48405 -0.000825414 0 3.48393 -0.000707505 0
+          3.48383 -0.000583365 0 3.48377 -0.000463725 0 3.48373 -0.000357438 0 3.4837 -0.000263977 0
+          3.48368 -0.000185575 0 3.48367 -0.000116445 0 3.48367 -5.64092e-05 0 3.48367 3.11334e-06 0
+          3.07414 0 0 3.07409 -0.000449736 0 3.07394 -0.000649757 0 3.07373 -0.000760754 0
+          3.07349 -0.000793851 0 3.07326 -0.000769477 0 3.07304 -0.00070552 0 3.07285 -0.000616976 0
+          3.07269 -0.000521419 0 3.07256 -0.00042388 0 3.07246 -0.000334998 0 3.07238 -0.000252532 0
+          3.07232 -0.000181182 0 3.07228 -0.0001156 0 3.07226 -5.6714e-05 0 3.07225 1.19443e-06 0
+          1.98047 0 0 1.9804 0.000185419 0 1.98027 0.000159102 0 1.98011 0.000112266 0
+          1.97993 6.80647e-05 0 1.97976 3.95752e-05 0 1.97958 2.11606e-05 0 1.97944 1.26034e-05 0
+          1.9793 9.11205e-06 0 1.9792 8.55727e-06 0 1.97911 8.96824e-06 0 1.97904 8.81083e-06 0
+          1.97899 8.14678e-06 0 1.97896 6.27105e-06 0 1.97894 4.41143e-06 0 1.97893 1.36795e-06 0
           0.00134 0 0 0.00134 0 0 0.00134 0 0 0.00134 0 0
           0.00134 0 0 0.00134 0 0 0.00134 0 0 0.00134 0 0
           0.00134 0 0 0.00134 0 0 0.00134 0 0 0.00134 0 0
diff --git a/test/references/2cnistokes2p2cniboundarylayer-pm-reference.vtu b/test/references/2cnistokes2p2cniboundarylayer-pm-reference.vtu
index e41c8fd218..df324448ea 100644
--- a/test/references/2cnistokes2p2cniboundarylayer-pm-reference.vtu
+++ b/test/references/2cnistokes2p2cniboundarylayer-pm-reference.vtu
@@ -4,24 +4,24 @@
     <Piece NumberOfCells="180" NumberOfPoints="208">
       <PointData Scalars="Sn" Vectors="velocityW">
         <DataArray type="Float32" Name="Sn" NumberOfComponents="1" format="ascii">
-          0.0100516 0.0100516 0.0101095 0.0101095 0.0100516 0.0101095 0.0100516 0.0101095 0.0100516 0.0101095 0.0100516 0.0101095
-          0.0100516 0.0101095 0.0100516 0.0101095 0.0100516 0.0101095 0.0100516 0.0101095 0.0100516 0.0101095 0.0100516 0.0101095
-          0.0100516 0.0101095 0.0100516 0.0101095 0.0100516 0.0101095 0.0100516 0.0101095 0.0101528 0.0101528 0.0101528 0.0101528
-          0.0101528 0.0101528 0.0101528 0.0101528 0.0101528 0.0101528 0.0101528 0.0101528 0.0101528 0.0101527 0.0101528 0.0101527
-          0.0102012 0.0102012 0.0102012 0.0102012 0.0102012 0.0102012 0.0102012 0.0102012 0.0102012 0.0102012 0.0102011 0.0102011
-          0.0102011 0.0102011 0.0102011 0.0102011 0.0103938 0.0103938 0.0103938 0.0103938 0.0103938 0.0103938 0.0103938 0.0103938
-          0.0103938 0.0103938 0.0103938 0.0103938 0.0103938 0.0103938 0.0103938 0.0103937 0.0113488 0.0113488 0.0113488 0.0113488
-          0.0113488 0.0113488 0.0113488 0.0113488 0.0113488 0.0113488 0.0113488 0.0113488 0.0113487 0.0113487 0.0113487 0.0113487
-          0.0127517 0.0127517 0.0127517 0.0127517 0.0127517 0.0127517 0.0127517 0.0127517 0.0127517 0.0127517 0.0127517 0.0127517
-          0.0127517 0.0127516 0.0127516 0.0127516 0.0143782 0.0143782 0.0143782 0.0143781 0.0143781 0.0143781 0.0143781 0.0143781
-          0.0143781 0.0143781 0.0143781 0.0143781 0.0143781 0.0143781 0.0143781 0.0143781 0.0164691 0.0164691 0.0164691 0.0164691
-          0.0164691 0.0164691 0.0164691 0.016469 0.016469 0.016469 0.016469 0.016469 0.016469 0.016469 0.016469 0.016469
-          0.022698 0.022698 0.0226979 0.0226979 0.0226979 0.0226979 0.0226979 0.0226978 0.0226978 0.0226978 0.0226978 0.0226978
-          0.0226978 0.0226978 0.0226977 0.0226977 0.0475076 0.0475075 0.0475074 0.0475074 0.0475073 0.0475073 0.0475073 0.0475073
-          0.0475073 0.0475073 0.0475072 0.0475072 0.0475072 0.0475072 0.0475071 0.0475069 0.096405 0.0964046 0.0964045 0.0964044
-          0.0964044 0.0964044 0.0964044 0.0964044 0.0964044 0.0964044 0.0964044 0.0964044 0.0964044 0.0964044 0.0964042 0.0964038
-          0.173248 0.173247 0.173247 0.173247 0.173247 0.173247 0.173247 0.173247 0.173248 0.173248 0.173248 0.173248
-          0.173248 0.173248 0.173248 0.173247
+          0.0100515 0.0100515 0.0101094 0.0101094 0.0100515 0.0101094 0.0100515 0.0101094 0.0100515 0.0101094 0.0100515 0.0101094
+          0.0100515 0.0101094 0.0100515 0.0101094 0.0100515 0.0101094 0.0100515 0.0101094 0.0100515 0.0101094 0.0100515 0.0101094
+          0.0100515 0.0101094 0.0100515 0.0101094 0.0100515 0.0101094 0.0100515 0.0101094 0.0101526 0.0101526 0.0101526 0.0101526
+          0.0101526 0.0101526 0.0101526 0.0101526 0.0101526 0.0101526 0.0101526 0.0101526 0.0101526 0.0101526 0.0101526 0.0101526
+          0.0102007 0.0102007 0.0102007 0.0102007 0.0102007 0.0102007 0.0102007 0.0102007 0.0102007 0.0102007 0.0102007 0.0102007
+          0.0102007 0.0102007 0.0102007 0.0102007 0.0103922 0.0103922 0.0103922 0.0103922 0.0103922 0.0103922 0.0103922 0.0103922
+          0.0103922 0.0103922 0.0103922 0.0103922 0.0103922 0.0103922 0.0103922 0.0103922 0.0113453 0.0113453 0.0113453 0.0113453
+          0.0113453 0.0113452 0.0113452 0.0113452 0.0113452 0.0113452 0.0113452 0.0113452 0.0113452 0.0113452 0.0113452 0.0113452
+          0.0127476 0.0127476 0.0127476 0.0127476 0.0127475 0.0127475 0.0127475 0.0127475 0.0127475 0.0127475 0.0127475 0.0127475
+          0.0127475 0.0127475 0.0127475 0.0127475 0.014373 0.014373 0.014373 0.014373 0.014373 0.014373 0.014373 0.0143729
+          0.0143729 0.0143729 0.0143729 0.0143729 0.0143729 0.0143729 0.0143729 0.0143729 0.0164612 0.0164612 0.0164612 0.0164612
+          0.0164612 0.0164612 0.0164612 0.0164612 0.0164612 0.0164612 0.0164612 0.0164612 0.0164612 0.0164612 0.0164612 0.0164611
+          0.0226457 0.0226457 0.0226456 0.0226456 0.0226456 0.0226456 0.0226456 0.0226456 0.0226455 0.0226455 0.0226455 0.0226455
+          0.0226455 0.0226455 0.0226454 0.0226454 0.0473786 0.0473785 0.0473784 0.0473784 0.0473783 0.0473783 0.0473782 0.0473782
+          0.0473782 0.0473781 0.0473781 0.0473781 0.0473781 0.047378 0.0473779 0.0473778 0.0961579 0.0961576 0.0961574 0.0961573
+          0.0961573 0.0961572 0.0961572 0.0961571 0.0961571 0.096157 0.096157 0.0961569 0.0961569 0.0961568 0.0961567 0.0961563
+          0.172869 0.172868 0.172868 0.172868 0.172868 0.172868 0.172868 0.172868 0.172868 0.172868 0.172867 0.172867
+          0.172867 0.172867 0.172867 0.172866
         </DataArray>
         <DataArray type="Float32" Name="Sw" NumberOfComponents="1" format="ascii">
           0.989948 0.989948 0.989891 0.989891 0.989948 0.989891 0.989948 0.989891 0.989948 0.989891 0.989948 0.989891
@@ -29,19 +29,19 @@
           0.989948 0.989891 0.989948 0.989891 0.989948 0.989891 0.989948 0.989891 0.989847 0.989847 0.989847 0.989847
           0.989847 0.989847 0.989847 0.989847 0.989847 0.989847 0.989847 0.989847 0.989847 0.989847 0.989847 0.989847
           0.989799 0.989799 0.989799 0.989799 0.989799 0.989799 0.989799 0.989799 0.989799 0.989799 0.989799 0.989799
-          0.989799 0.989799 0.989799 0.989799 0.989606 0.989606 0.989606 0.989606 0.989606 0.989606 0.989606 0.989606
-          0.989606 0.989606 0.989606 0.989606 0.989606 0.989606 0.989606 0.989606 0.988651 0.988651 0.988651 0.988651
-          0.988651 0.988651 0.988651 0.988651 0.988651 0.988651 0.988651 0.988651 0.988651 0.988651 0.988651 0.988651
-          0.987248 0.987248 0.987248 0.987248 0.987248 0.987248 0.987248 0.987248 0.987248 0.987248 0.987248 0.987248
-          0.987248 0.987248 0.987248 0.987248 0.985622 0.985622 0.985622 0.985622 0.985622 0.985622 0.985622 0.985622
-          0.985622 0.985622 0.985622 0.985622 0.985622 0.985622 0.985622 0.985622 0.983531 0.983531 0.983531 0.983531
-          0.983531 0.983531 0.983531 0.983531 0.983531 0.983531 0.983531 0.983531 0.983531 0.983531 0.983531 0.983531
-          0.977302 0.977302 0.977302 0.977302 0.977302 0.977302 0.977302 0.977302 0.977302 0.977302 0.977302 0.977302
-          0.977302 0.977302 0.977302 0.977302 0.952492 0.952493 0.952493 0.952493 0.952493 0.952493 0.952493 0.952493
-          0.952493 0.952493 0.952493 0.952493 0.952493 0.952493 0.952493 0.952493 0.903595 0.903595 0.903596 0.903596
-          0.903596 0.903596 0.903596 0.903596 0.903596 0.903596 0.903596 0.903596 0.903596 0.903596 0.903596 0.903596
-          0.826752 0.826753 0.826753 0.826753 0.826753 0.826753 0.826753 0.826753 0.826752 0.826752 0.826752 0.826752
-          0.826752 0.826752 0.826752 0.826753
+          0.989799 0.989799 0.989799 0.989799 0.989608 0.989608 0.989608 0.989608 0.989608 0.989608 0.989608 0.989608
+          0.989608 0.989608 0.989608 0.989608 0.989608 0.989608 0.989608 0.989608 0.988655 0.988655 0.988655 0.988655
+          0.988655 0.988655 0.988655 0.988655 0.988655 0.988655 0.988655 0.988655 0.988655 0.988655 0.988655 0.988655
+          0.987252 0.987252 0.987252 0.987252 0.987252 0.987252 0.987252 0.987252 0.987252 0.987252 0.987252 0.987252
+          0.987252 0.987252 0.987252 0.987252 0.985627 0.985627 0.985627 0.985627 0.985627 0.985627 0.985627 0.985627
+          0.985627 0.985627 0.985627 0.985627 0.985627 0.985627 0.985627 0.985627 0.983539 0.983539 0.983539 0.983539
+          0.983539 0.983539 0.983539 0.983539 0.983539 0.983539 0.983539 0.983539 0.983539 0.983539 0.983539 0.983539
+          0.977354 0.977354 0.977354 0.977354 0.977354 0.977354 0.977354 0.977354 0.977354 0.977354 0.977354 0.977355
+          0.977355 0.977355 0.977355 0.977355 0.952621 0.952622 0.952622 0.952622 0.952622 0.952622 0.952622 0.952622
+          0.952622 0.952622 0.952622 0.952622 0.952622 0.952622 0.952622 0.952622 0.903842 0.903842 0.903843 0.903843
+          0.903843 0.903843 0.903843 0.903843 0.903843 0.903843 0.903843 0.903843 0.903843 0.903843 0.903843 0.903844
+          0.827131 0.827132 0.827132 0.827132 0.827132 0.827132 0.827132 0.827132 0.827132 0.827132 0.827133 0.827133
+          0.827133 0.827133 0.827133 0.827134
         </DataArray>
         <DataArray type="Float32" Name="pn" NumberOfComponents="1" format="ascii">
           101149 101149 100833 100833 101149 100833 101149 100833 101149 100833 101149 100833
@@ -64,124 +64,124 @@
           100000 100000 100000 100000
         </DataArray>
         <DataArray type="Float32" Name="pw" NumberOfComponents="1" format="ascii">
-          101139 101139 100813 100813 101139 100813 101139 100813 101139 100813 101139 100813
-          101139 100813 101139 100813 101139 100813 101139 100813 101139 100813 101139 100813
-          101139 100813 101139 100813 101139 100813 101139 100813 100516 100516 100516 100516
-          100516 100516 100516 100516 100516 100516 100516 100516 100516 100516 100516 100516
+          101140 101140 100813 100813 101140 100813 101140 100813 101140 100813 101140 100813
+          101140 100813 101140 100813 101140 100813 101140 100813 101140 100813 101140 100813
+          101140 100813 101140 100813 101140 100813 101140 100813 100517 100517 100517 100517
+          100517 100517 100517 100517 100517 100517 100517 100517 100517 100517 100517 100517
           100247 100247 100247 100247 100247 100247 100247 100247 100247 100247 100247 100247
           100247 100247 100247 100247 100002 100002 100002 100002 100002 100002 100002 100002
-          100002 100002 100002 100002 100002 100002 100002 100002 99778.8 99778.8 99778.8 99778.8
-          99778.8 99778.8 99778.8 99778.8 99778.8 99778.8 99778.8 99778.8 99778.8 99778.8 99778.8 99778.8
-          99576.2 99576.2 99576.2 99576.2 99576.2 99576.2 99576.2 99576.2 99576.2 99576.2 99576.2 99576.2
-          99576.2 99576.2 99576.2 99576.2 99392 99392 99392 99392 99392 99392 99392 99392
-          99392 99392 99392 99392 99392 99392 99392 99392 99224.6 99224.6 99224.6 99224.6
-          99224.6 99224.6 99224.6 99224.6 99224.6 99224.6 99224.6 99224.6 99224.6 99224.6 99224.6 99224.6
-          99072.4 99072.4 99072.4 99072.4 99072.4 99072.4 99072.4 99072.4 99072.4 99072.4 99072.4 99072.4
-          99072.4 99072.4 99072.4 99072.4 98934 98934 98934 98934 98934 98934 98934 98934
-          98934 98934 98934 98934 98934 98934 98934 98934 98808.1 98808.1 98808.1 98808.1
-          98808.1 98808.1 98808.1 98808.1 98808.1 98808.1 98808.1 98808.1 98808.1 98808.1 98808.1 98808.1
-          98693.7 98693.7 98693.7 98693.7 98693.7 98693.7 98693.7 98693.7 98693.7 98693.7 98693.7 98693.7
-          98693.7 98693.7 98693.7 98693.7
+          100002 100002 100002 100002 100002 100002 100002 100002 99779.3 99779.3 99779.3 99779.3
+          99779.3 99779.3 99779.3 99779.3 99779.3 99779.3 99779.3 99779.3 99779.3 99779.3 99779.3 99779.3
+          99576.7 99576.7 99576.7 99576.7 99576.7 99576.7 99576.7 99576.7 99576.7 99576.7 99576.7 99576.7
+          99576.7 99576.7 99576.7 99576.7 99392.5 99392.5 99392.5 99392.5 99392.5 99392.5 99392.5 99392.5
+          99392.5 99392.5 99392.5 99392.5 99392.5 99392.5 99392.5 99392.5 99225.1 99225.1 99225.1 99225.1
+          99225.1 99225.1 99225.1 99225.1 99225.1 99225.1 99225.1 99225.1 99225.1 99225.1 99225.1 99225.1
+          99072.8 99072.8 99072.8 99072.8 99072.8 99072.8 99072.8 99072.8 99072.8 99072.8 99072.8 99072.8
+          99072.8 99072.8 99072.8 99072.8 98934.4 98934.4 98934.4 98934.4 98934.4 98934.4 98934.4 98934.4
+          98934.4 98934.4 98934.4 98934.4 98934.4 98934.4 98934.4 98934.4 98808.6 98808.6 98808.6 98808.6
+          98808.6 98808.6 98808.6 98808.6 98808.6 98808.6 98808.6 98808.6 98808.6 98808.6 98808.6 98808.6
+          98694.2 98694.2 98694.2 98694.2 98694.2 98694.2 98694.2 98694.2 98694.2 98694.2 98694.2 98694.2
+          98694.2 98694.2 98694.2 98694.2
         </DataArray>
         <DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii">
-          9.38355 9.38348 19.8537 19.8539 9.3835 19.854 9.38352 19.8538 9.3835 19.8537 9.38348 19.8537
-          9.38349 19.8539 9.38362 19.8533 9.38348 19.8526 9.38362 19.853 9.38371 19.8533 9.38347 19.8532
-          9.38364 19.853 9.3835 19.8529 9.38353 19.8532 9.38339 19.853 27.6332 27.633 27.6331 27.6332
-          27.633 27.6322 27.6324 27.6319 27.6319 27.6317 27.6317 27.6313 27.6315 27.6309 27.6312 27.631
-          36.289 36.2889 36.2886 36.2883 36.2878 36.2867 36.2863 36.2859 36.2855 36.2851 36.2848 36.2845
-          36.284 36.2831 36.2836 36.283 70.2572 70.2567 70.2563 70.2544 70.2531 70.2525 70.2511 70.2494
-          70.2476 70.2467 70.2459 70.2443 70.2437 70.2426 70.2421 70.2417 227.444 227.444 227.443 227.443
-          227.442 227.441 227.44 227.438 227.437 227.436 227.435 227.434 227.433 227.433 227.432 227.432
-          425.61 425.61 425.61 425.609 425.609 425.608 425.608 425.607 425.607 425.606 425.605 425.605
-          425.605 425.604 425.604 425.604 608.915 608.915 608.915 608.914 608.914 608.914 608.913 608.913
-          608.913 608.912 608.912 608.912 608.911 608.911 608.911 608.911 776.029 776.028 776.028 776.028
-          776.028 776.027 776.027 776.027 776.027 776.027 776.026 776.026 776.026 776.026 776.026 776.025
-          928.076 928.076 928.075 928.075 928.075 928.075 928.075 928.075 928.075 928.074 928.074 928.074
-          928.074 928.074 928.073 928.073 1066.32 1066.32 1066.32 1066.32 1066.32 1066.31 1066.31 1066.31
-          1066.31 1066.31 1066.31 1066.31 1066.31 1066.31 1066.31 1066.31 1192 1192 1192 1192
-          1192 1192 1192 1192 1192 1192 1192 1192 1192 1192 1192 1191.99
-          1306.26 1306.26 1306.26 1306.26 1306.26 1306.26 1306.26 1306.26 1306.26 1306.26 1306.26 1306.26
-          1306.26 1306.26 1306.26 1306.26
+          9.37424 9.37434 19.8369 19.8372 9.37419 19.837 9.37421 19.8367 9.37439 19.8368 9.37424 19.8369
+          9.37401 19.8364 9.37421 19.8366 9.37415 19.837 9.37428 19.8368 9.37417 19.8362 9.37424 19.8365
+          9.37402 19.8368 9.37407 19.837 9.3741 19.8369 9.37416 19.8364 27.6016 27.601 27.6011 27.6013
+          27.6012 27.601 27.6013 27.6006 27.6008 27.6005 27.6009 27.6007 27.6005 27.6006 27.6013 27.6008
+          36.2131 36.2135 36.2128 36.2127 36.2131 36.2121 36.2125 36.212 36.2119 36.2113 36.2119 36.2114
+          36.211 36.2111 36.2116 36.2112 69.9762 69.9757 69.9749 69.975 69.9747 69.9739 69.9728 69.9725
+          69.9721 69.9717 69.9712 69.9709 69.9703 69.9703 69.9696 69.9703 226.893 226.893 226.893 226.893
+          226.892 226.892 226.891 226.891 226.89 226.889 226.889 226.888 226.888 226.888 226.888 226.887
+          425.081 425.081 425.081 425.08 425.08 425.08 425.079 425.079 425.079 425.078 425.078 425.078
+          425.077 425.077 425.077 425.077 608.406 608.406 608.406 608.406 608.406 608.405 608.405 608.405
+          608.405 608.404 608.404 608.404 608.404 608.403 608.403 608.403 775.537 775.536 775.536 775.536
+          775.536 775.535 775.535 775.535 775.535 775.535 775.535 775.534 775.534 775.534 775.534 775.533
+          927.591 927.591 927.59 927.59 927.59 927.59 927.59 927.589 927.589 927.589 927.589 927.589
+          927.589 927.589 927.588 927.588 1065.84 1065.84 1065.84 1065.84 1065.84 1065.84 1065.84 1065.84
+          1065.84 1065.84 1065.84 1065.84 1065.84 1065.84 1065.83 1065.83 1191.53 1191.53 1191.53 1191.53
+          1191.53 1191.52 1191.52 1191.52 1191.52 1191.52 1191.52 1191.52 1191.52 1191.52 1191.52 1191.52
+          1305.8 1305.8 1305.8 1305.8 1305.8 1305.8 1305.8 1305.8 1305.8 1305.8 1305.8 1305.8
+          1305.8 1305.8 1305.8 1305.8
         </DataArray>
         <DataArray type="Float32" Name="rhoW" NumberOfComponents="1" format="ascii">
           997.056 997.056 997.057 997.057 997.056 997.057 997.056 997.057 997.056 997.057 997.056 997.057
           997.056 997.057 997.056 997.057 997.056 997.057 997.056 997.057 997.056 997.057 997.056 997.057
-          997.056 997.057 997.056 997.057 997.056 997.057 997.056 997.057 997.057 997.057 997.057 997.057
-          997.057 997.057 997.057 997.057 997.057 997.057 997.057 997.057 997.057 997.057 997.057 997.057
-          997.06 997.06 997.06 997.06 997.06 997.06 997.06 997.06 997.06 997.06 997.06 997.06
-          997.06 997.06 997.06 997.06 997.066 997.066 997.066 997.066 997.066 997.066 997.066 997.066
-          997.066 997.066 997.066 997.066 997.066 997.066 997.066 997.066 997.079 997.079 997.079 997.079
-          997.079 997.079 997.079 997.079 997.079 997.079 997.079 997.079 997.079 997.079 997.079 997.079
-          997.102 997.102 997.102 997.102 997.102 997.102 997.102 997.102 997.102 997.102 997.103 997.103
-          997.103 997.103 997.103 997.103 997.138 997.138 997.138 997.138 997.138 997.138 997.138 997.138
-          997.139 997.139 997.139 997.139 997.139 997.139 997.139 997.139 997.188 997.188 997.188 997.188
-          997.188 997.188 997.189 997.189 997.189 997.189 997.19 997.19 997.19 997.19 997.19 997.19
-          997.25 997.25 997.25 997.251 997.251 997.252 997.252 997.252 997.253 997.253 997.253 997.254
-          997.254 997.254 997.254 997.254 997.322 997.322 997.323 997.323 997.324 997.324 997.325 997.326
-          997.326 997.326 997.327 997.327 997.328 997.328 997.328 997.328 997.399 997.4 997.4 997.401
-          997.402 997.403 997.403 997.404 997.405 997.406 997.406 997.407 997.407 997.408 997.408 997.408
-          997.477 997.478 997.479 997.48 997.481 997.482 997.483 997.484 997.485 997.486 997.487 997.488
-          997.488 997.489 997.49 997.49
+          997.056 997.057 997.056 997.057 997.056 997.057 997.056 997.057 997.058 997.058 997.058 997.058
+          997.058 997.058 997.058 997.058 997.058 997.058 997.058 997.058 997.058 997.058 997.058 997.058
+          997.061 997.061 997.061 997.061 997.061 997.061 997.061 997.061 997.061 997.061 997.061 997.061
+          997.061 997.061 997.061 997.061 997.068 997.068 997.068 997.068 997.068 997.068 997.068 997.068
+          997.068 997.068 997.068 997.068 997.068 997.068 997.068 997.068 997.083 997.083 997.083 997.083
+          997.083 997.083 997.083 997.083 997.083 997.083 997.083 997.083 997.083 997.083 997.083 997.083
+          997.109 997.109 997.109 997.109 997.109 997.109 997.109 997.109 997.109 997.109 997.109 997.109
+          997.109 997.109 997.109 997.109 997.152 997.152 997.152 997.152 997.152 997.152 997.152 997.152
+          997.152 997.152 997.152 997.152 997.152 997.152 997.152 997.152 997.211 997.211 997.212 997.212
+          997.212 997.212 997.212 997.212 997.212 997.212 997.212 997.212 997.212 997.212 997.213 997.213
+          997.288 997.288 997.288 997.288 997.288 997.289 997.289 997.289 997.289 997.289 997.289 997.29
+          997.29 997.29 997.29 997.29 997.378 997.378 997.378 997.378 997.379 997.379 997.379 997.379
+          997.38 997.38 997.38 997.38 997.38 997.381 997.381 997.381 997.476 997.476 997.476 997.477
+          997.477 997.478 997.478 997.478 997.479 997.479 997.479 997.48 997.48 997.48 997.48 997.48
+          997.578 997.578 997.579 997.579 997.58 997.58 997.581 997.581 997.582 997.582 997.583 997.583
+          997.583 997.584 997.584 997.584
         </DataArray>
         <DataArray type="Float32" Name="rhoN" NumberOfComponents="1" format="ascii">
-          1.16769 1.16769 1.16401 1.16401 1.16769 1.16401 1.16769 1.16401 1.16769 1.16401 1.16769 1.16401
-          1.16769 1.16401 1.16769 1.16401 1.16769 1.16401 1.16769 1.16401 1.16769 1.16401 1.16769 1.16401
-          1.16769 1.16401 1.16769 1.16401 1.16769 1.16401 1.16769 1.16401 1.16066 1.16066 1.16066 1.16066
+          1.1677 1.1677 1.16402 1.16402 1.1677 1.16402 1.1677 1.16402 1.1677 1.16402 1.1677 1.16402
+          1.1677 1.16402 1.1677 1.16402 1.1677 1.16402 1.1677 1.16402 1.1677 1.16402 1.1677 1.16402
+          1.1677 1.16402 1.1677 1.16402 1.1677 1.16402 1.1677 1.16402 1.16066 1.16066 1.16066 1.16066
           1.16066 1.16066 1.16066 1.16066 1.16066 1.16066 1.16066 1.16066 1.16066 1.16066 1.16066 1.16066
-          1.15766 1.15766 1.15766 1.15766 1.15766 1.15766 1.15766 1.15766 1.15766 1.15766 1.15766 1.15766
-          1.15766 1.15766 1.15766 1.15766 1.15531 1.15531 1.15531 1.15531 1.15531 1.15531 1.15531 1.15531
-          1.15531 1.15531 1.15531 1.15531 1.15531 1.15531 1.15531 1.15531 1.15478 1.15478 1.15478 1.15478
-          1.15478 1.15478 1.15478 1.15478 1.15478 1.15478 1.15478 1.15478 1.15478 1.15478 1.15478 1.15478
-          1.15515 1.15515 1.15515 1.15515 1.15515 1.15515 1.15515 1.15516 1.15516 1.15516 1.15516 1.15516
-          1.15516 1.15516 1.15516 1.15516 1.1558 1.1558 1.1558 1.15581 1.15581 1.15581 1.15581 1.15582
-          1.15582 1.15582 1.15582 1.15583 1.15583 1.15583 1.15583 1.15583 1.15673 1.15673 1.15673 1.15673
-          1.15674 1.15674 1.15675 1.15675 1.15676 1.15676 1.15676 1.15677 1.15677 1.15677 1.15678 1.15678
-          1.1579 1.1579 1.15791 1.15791 1.15792 1.15793 1.15794 1.15794 1.15795 1.15796 1.15796 1.15797
-          1.15797 1.15798 1.15798 1.15798 1.15927 1.15927 1.15928 1.15929 1.1593 1.15931 1.15932 1.15933
-          1.15934 1.15935 1.15936 1.15937 1.15937 1.15938 1.15938 1.15939 1.16074 1.16075 1.16076 1.16078
-          1.16079 1.16081 1.16082 1.16084 1.16085 1.16086 1.16088 1.16089 1.1609 1.16091 1.16092 1.16092
-          1.16226 1.16227 1.16229 1.16232 1.16234 1.16236 1.16238 1.1624 1.16241 1.16243 1.16245 1.16246
-          1.16247 1.16249 1.1625 1.16251
+          1.15767 1.15767 1.15767 1.15767 1.15767 1.15767 1.15767 1.15767 1.15767 1.15767 1.15767 1.15767
+          1.15767 1.15767 1.15767 1.15767 1.15534 1.15534 1.15534 1.15534 1.15534 1.15534 1.15534 1.15534
+          1.15534 1.15534 1.15534 1.15534 1.15534 1.15534 1.15534 1.15534 1.15484 1.15484 1.15484 1.15484
+          1.15484 1.15484 1.15484 1.15484 1.15484 1.15484 1.15484 1.15484 1.15484 1.15484 1.15484 1.15484
+          1.15528 1.15528 1.15528 1.15528 1.15528 1.15528 1.15528 1.15529 1.15529 1.15529 1.15529 1.15529
+          1.15529 1.15529 1.15529 1.15529 1.15606 1.15606 1.15606 1.15606 1.15606 1.15606 1.15606 1.15606
+          1.15607 1.15607 1.15607 1.15607 1.15607 1.15607 1.15607 1.15607 1.15717 1.15718 1.15718 1.15718
+          1.15718 1.15718 1.15718 1.15719 1.15719 1.15719 1.15719 1.15719 1.15719 1.15719 1.1572 1.1572
+          1.15862 1.15862 1.15862 1.15862 1.15863 1.15863 1.15863 1.15864 1.15864 1.15864 1.15864 1.15865
+          1.15865 1.15865 1.15865 1.15865 1.16033 1.16033 1.16033 1.16034 1.16035 1.16035 1.16036 1.16036
+          1.16037 1.16037 1.16037 1.16038 1.16038 1.16038 1.16039 1.16039 1.16223 1.16223 1.16224 1.16225
+          1.16226 1.16226 1.16227 1.16228 1.16229 1.16229 1.1623 1.1623 1.16231 1.16231 1.16232 1.16232
+          1.16423 1.16424 1.16425 1.16426 1.16428 1.16429 1.1643 1.16431 1.16432 1.16433 1.16433 1.16434
+          1.16435 1.16436 1.16436 1.16437
         </DataArray>
         <DataArray type="Float32" Name="mobW" NumberOfComponents="1" format="ascii">
-          1123.57 1123.57 1123.53 1123.53 1123.57 1123.53 1123.57 1123.53 1123.57 1123.53 1123.57 1123.53
-          1123.57 1123.53 1123.57 1123.53 1123.57 1123.53 1123.57 1123.53 1123.57 1123.53 1123.57 1123.53
-          1123.57 1123.53 1123.57 1123.53 1123.57 1123.53 1123.57 1123.53 1123.43 1123.43 1123.43 1123.43
-          1123.43 1123.43 1123.43 1123.42 1123.42 1123.42 1123.42 1123.42 1123.42 1123.42 1123.42 1123.42
-          1123.15 1123.15 1123.15 1123.15 1123.15 1123.15 1123.15 1123.15 1123.15 1123.15 1123.15 1123.15
-          1123.15 1123.15 1123.15 1123.15 1122.51 1122.51 1122.51 1122.51 1122.5 1122.5 1122.5 1122.5
-          1122.5 1122.5 1122.5 1122.5 1122.49 1122.49 1122.49 1122.49 1121.08 1121.08 1121.08 1121.08
-          1121.07 1121.07 1121.07 1121.06 1121.06 1121.06 1121.05 1121.05 1121.05 1121.05 1121.05 1121.05
-          1118.33 1118.33 1118.33 1118.32 1118.31 1118.31 1118.3 1118.29 1118.28 1118.28 1118.27 1118.27
-          1118.26 1118.26 1118.25 1118.25 1113.83 1113.83 1113.82 1113.81 1113.8 1113.78 1113.77 1113.76
-          1113.74 1113.73 1113.72 1113.71 1113.7 1113.69 1113.68 1113.68 1107.08 1107.08 1107.06 1107.04
-          1107.02 1106.99 1106.97 1106.94 1106.92 1106.9 1106.88 1106.86 1106.84 1106.83 1106.82 1106.82
-          1092.11 1092.1 1092.08 1092.04 1092 1091.96 1091.92 1091.88 1091.85 1091.81 1091.78 1091.75
-          1091.72 1091.7 1091.69 1091.68 1009.49 1009.47 1009.43 1009.38 1009.33 1009.27 1009.22 1009.17
-          1009.12 1009.07 1009.03 1008.99 1008.95 1008.92 1008.9 1008.89 787.543 787.52 787.469 787.41
-          787.348 787.288 787.23 787.175 787.122 787.072 787.025 786.98 786.939 786.903 786.874 786.863
-          584.051 584.008 583.949 583.887 583.825 583.766 583.71 583.656 583.605 583.557 583.51 583.467
-          583.426 583.388 583.356 583.337
+          1123.57 1123.57 1123.52 1123.52 1123.57 1123.52 1123.57 1123.52 1123.57 1123.52 1123.57 1123.52
+          1123.57 1123.52 1123.57 1123.52 1123.57 1123.52 1123.57 1123.52 1123.57 1123.52 1123.57 1123.52
+          1123.57 1123.52 1123.57 1123.52 1123.57 1123.52 1123.57 1123.52 1123.41 1123.41 1123.41 1123.41
+          1123.41 1123.41 1123.41 1123.41 1123.41 1123.41 1123.41 1123.41 1123.41 1123.41 1123.41 1123.41
+          1123.1 1123.1 1123.1 1123.1 1123.1 1123.1 1123.1 1123.1 1123.1 1123.1 1123.1 1123.1
+          1123.1 1123.1 1123.1 1123.1 1122.36 1122.36 1122.36 1122.36 1122.36 1122.36 1122.36 1122.36
+          1122.36 1122.36 1122.36 1122.36 1122.36 1122.36 1122.36 1122.36 1120.74 1120.74 1120.74 1120.74
+          1120.74 1120.74 1120.73 1120.73 1120.73 1120.73 1120.73 1120.73 1120.73 1120.73 1120.73 1120.73
+          1117.6 1117.6 1117.6 1117.6 1117.6 1117.59 1117.59 1117.59 1117.59 1117.58 1117.58 1117.58
+          1117.58 1117.58 1117.58 1117.58 1112.44 1112.44 1112.44 1112.43 1112.43 1112.42 1112.42 1112.41
+          1112.41 1112.4 1112.4 1112.39 1112.39 1112.39 1112.38 1112.38 1104.67 1104.67 1104.66 1104.65
+          1104.64 1104.63 1104.62 1104.61 1104.6 1104.59 1104.58 1104.57 1104.57 1104.56 1104.56 1104.56
+          1088.37 1088.36 1088.35 1088.34 1088.32 1088.3 1088.28 1088.26 1088.25 1088.23 1088.22 1088.21
+          1088.19 1088.19 1088.18 1088.18 1004.71 1004.71 1004.69 1004.66 1004.64 1004.61 1004.58 1004.56
+          1004.53 1004.51 1004.49 1004.47 1004.46 1004.44 1004.43 1004.43 782.669 782.657 782.631 782.6
+          782.567 782.536 782.506 782.477 782.451 782.426 782.402 782.38 782.36 782.342 782.329 782.324
+          579.145 579.124 579.091 579.057 579.023 578.99 578.959 578.93 578.903 578.877 578.852 578.829
+          578.808 578.789 578.772 578.764
         </DataArray>
         <DataArray type="Float32" Name="mobN" NumberOfComponents="1" format="ascii">
-          8.48318e-05 8.48304e-05 0.000407054 0.000407064 8.48309e-05 0.000407067 8.48313e-05 0.000407061 8.48308e-05 0.000407055 8.48305e-05 0.000407054
-          8.48306e-05 0.000407062 8.48331e-05 0.000407037 8.48304e-05 0.000407006 8.4833e-05 0.000407026 8.48347e-05 0.00040704 8.48302e-05 0.000407034
-          8.48334e-05 0.000407027 8.48308e-05 0.000407022 8.48314e-05 0.000407032 8.48288e-05 0.000407025 0.000814835 0.000814818 0.000814824 0.00081483
-          0.000814822 0.000814773 0.000814783 0.000814753 0.000814748 0.000814737 0.000814737 0.000814717 0.000814725 0.000814691 0.000814706 0.000814696
-          0.00144597 0.00144596 0.00144594 0.00144591 0.00144587 0.00144578 0.00144574 0.00144571 0.00144568 0.00144565 0.00144562 0.00144559
-          0.00144555 0.00144548 0.00144552 0.00144547 0.00586129 0.00586121 0.00586114 0.00586081 0.00586056 0.00586046 0.00586021 0.0058599
-          0.00585959 0.00585944 0.00585929 0.005859 0.0058589 0.00585871 0.00585861 0.00585854 0.0761735 0.0761734 0.0761729 0.0761725
-          0.0761718 0.076171 0.0761702 0.0761693 0.0761685 0.0761676 0.0761669 0.0761663 0.0761656 0.0761652 0.0761649 0.0761648
-          0.33643 0.336429 0.336429 0.336428 0.336427 0.336426 0.336425 0.336424 0.336423 0.336422 0.336421 0.33642
-          0.33642 0.336419 0.336419 0.336418 0.885312 0.885312 0.885311 0.88531 0.88531 0.885309 0.885308 0.885308
-          0.885307 0.885306 0.885306 0.885305 0.885304 0.885303 0.885303 0.885302 1.9969 1.9969 1.9969 1.9969
-          1.9969 1.9969 1.9969 1.99691 1.99691 1.99691 1.99691 1.99691 1.99691 1.99691 1.99691 1.99691
-          8.13562 8.13557 8.13554 8.13553 8.13553 8.13554 8.13554 8.13554 8.13555 8.13555 8.13555 8.13555
-          8.13554 8.13553 8.13549 8.13543 77.4793 77.4788 77.4787 77.4788 77.479 77.4792 77.4795 77.4797
-          77.4799 77.4801 77.4803 77.4804 77.4805 77.4806 77.4803 77.4797 438.182 438.179 438.18 438.182
-          438.184 438.187 438.189 438.192 438.194 438.196 438.198 438.2 438.201 438.202 438.202 438.198
-          1632.79 1632.78 1632.79 1632.8 1632.82 1632.83 1632.84 1632.86 1632.87 1632.88 1632.89 1632.9
-          1632.91 1632.92 1632.93 1632.91
+          8.46559e-05 8.46578e-05 0.000406335 0.000406348 8.46551e-05 0.000406339 8.46554e-05 0.000406324 8.46589e-05 0.00040633 8.4656e-05 0.000406332
+          8.46517e-05 0.000406312 8.46553e-05 0.000406322 8.46543e-05 0.00040634 8.46567e-05 0.000406328 8.46547e-05 0.000406303 8.46559e-05 0.000406318
+          8.46519e-05 0.00040633 8.46528e-05 0.000406338 8.46533e-05 0.000406333 8.46545e-05 0.000406313 0.000812873 0.000812841 0.000812844 0.000812856
+          0.000812849 0.000812837 0.000812859 0.000812815 0.000812826 0.000812807 0.000812833 0.00081282 0.000812809 0.000812812 0.000812857 0.000812823
+          0.00143961 0.00143965 0.00143958 0.00143958 0.00143961 0.00143953 0.00143956 0.00143952 0.00143951 0.00143946 0.00143951 0.00143947
+          0.00143943 0.00143944 0.00143948 0.00143945 0.00581149 0.0058114 0.00581126 0.00581128 0.00581122 0.00581108 0.00581088 0.00581084
+          0.00581077 0.00581068 0.0058106 0.00581054 0.00581045 0.00581044 0.00581032 0.00581043 0.0757578 0.0757578 0.0757577 0.0757574
+          0.075757 0.0757566 0.0757562 0.0757558 0.0757554 0.0757549 0.0757546 0.0757543 0.075754 0.0757537 0.0757536 0.0757535
+          0.335392 0.335392 0.335391 0.335391 0.33539 0.33539 0.335389 0.335388 0.335388 0.335387 0.335386 0.335386
+          0.335385 0.335385 0.335384 0.335384 0.883209 0.883209 0.883208 0.883207 0.883206 0.883205 0.883205 0.883204
+          0.883203 0.883202 0.883202 0.883201 0.8832 0.883199 0.883198 0.883198 1.99217 1.99216 1.99216 1.99216
+          1.99216 1.99216 1.99216 1.99216 1.99216 1.99216 1.99216 1.99216 1.99216 1.99216 1.99215 1.99215
+          8.06818 8.06812 8.06808 8.06806 8.06805 8.06804 8.06803 8.06801 8.068 8.06799 8.06798 8.06796
+          8.06795 8.06793 8.06788 8.06782 76.9567 76.9562 76.9559 76.9559 76.9559 76.9559 76.9559 76.9559
+          76.9559 76.9559 76.9559 76.9558 76.9558 76.9557 76.9554 76.9548 435.833 435.829 435.829 435.83
+          435.831 435.831 435.832 435.833 435.833 435.834 435.835 435.835 435.835 435.835 435.834 435.83
+          1626.2 1626.18 1626.19 1626.2 1626.2 1626.21 1626.21 1626.22 1626.22 1626.23 1626.23 1626.24
+          1626.24 1626.24 1626.25 1626.22
         </DataArray>
         <DataArray type="Float32" Name="X_liquid^H2O" NumberOfComponents="1" format="ascii">
           0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978 0.999978
@@ -204,64 +204,64 @@
           0.999978 0.999978 0.999978 0.999978
         </DataArray>
         <DataArray type="Float32" Name="X_liquid^Air" NumberOfComponents="1" format="ascii">
-          2.18024e-05 2.18024e-05 2.17327e-05 2.17327e-05 2.18024e-05 2.17327e-05 2.18024e-05 2.17327e-05 2.18024e-05 2.17327e-05 2.18024e-05 2.17327e-05
-          2.18024e-05 2.17327e-05 2.18024e-05 2.17327e-05 2.18024e-05 2.17327e-05 2.18024e-05 2.17327e-05 2.18024e-05 2.17327e-05 2.18024e-05 2.17327e-05
-          2.18024e-05 2.17327e-05 2.18024e-05 2.17327e-05 2.18024e-05 2.17327e-05 2.18024e-05 2.17327e-05 2.167e-05 2.167e-05 2.167e-05 2.167e-05
-          2.167e-05 2.167e-05 2.167e-05 2.167e-05 2.167e-05 2.167e-05 2.167e-05 2.167e-05 2.167e-05 2.167e-05 2.167e-05 2.167e-05
-          2.16159e-05 2.16159e-05 2.16159e-05 2.16159e-05 2.16159e-05 2.16159e-05 2.16159e-05 2.16159e-05 2.16159e-05 2.1616e-05 2.1616e-05 2.1616e-05
-          2.1616e-05 2.1616e-05 2.1616e-05 2.1616e-05 2.15783e-05 2.15783e-05 2.15783e-05 2.15783e-05 2.15783e-05 2.15784e-05 2.15784e-05 2.15784e-05
-          2.15784e-05 2.15784e-05 2.15785e-05 2.15785e-05 2.15785e-05 2.15785e-05 2.15785e-05 2.15785e-05 2.15828e-05 2.15828e-05 2.15828e-05 2.15829e-05
-          2.15829e-05 2.1583e-05 2.1583e-05 2.15831e-05 2.15831e-05 2.15832e-05 2.15832e-05 2.15833e-05 2.15833e-05 2.15833e-05 2.15833e-05 2.15833e-05
-          2.16159e-05 2.16159e-05 2.1616e-05 2.16161e-05 2.16162e-05 2.16163e-05 2.16164e-05 2.16165e-05 2.16166e-05 2.16167e-05 2.16168e-05 2.16169e-05
-          2.1617e-05 2.1617e-05 2.16171e-05 2.16171e-05 2.16695e-05 2.16695e-05 2.16696e-05 2.16698e-05 2.167e-05 2.16702e-05 2.16704e-05 2.16706e-05
-          2.16708e-05 2.1671e-05 2.16712e-05 2.16714e-05 2.16715e-05 2.16716e-05 2.16717e-05 2.16717e-05 2.17452e-05 2.17453e-05 2.17455e-05 2.17458e-05
-          2.17461e-05 2.17465e-05 2.17469e-05 2.17473e-05 2.17476e-05 2.1748e-05 2.17483e-05 2.17486e-05 2.17488e-05 2.1749e-05 2.17492e-05 2.17492e-05
-          2.18419e-05 2.1842e-05 2.18424e-05 2.18429e-05 2.18435e-05 2.18442e-05 2.18448e-05 2.18454e-05 2.18459e-05 2.18465e-05 2.1847e-05 2.18474e-05
-          2.18478e-05 2.18482e-05 2.18484e-05 2.18485e-05 2.19555e-05 2.19557e-05 2.19564e-05 2.19573e-05 2.19582e-05 2.19591e-05 2.196e-05 2.19609e-05
-          2.19617e-05 2.19625e-05 2.19632e-05 2.19639e-05 2.19645e-05 2.19651e-05 2.19655e-05 2.19656e-05 2.208e-05 2.20805e-05 2.20816e-05 2.20829e-05
-          2.20843e-05 2.20856e-05 2.20868e-05 2.2088e-05 2.20892e-05 2.20903e-05 2.20913e-05 2.20923e-05 2.20932e-05 2.2094e-05 2.20946e-05 2.20949e-05
-          2.22093e-05 2.22106e-05 2.22123e-05 2.22141e-05 2.2216e-05 2.22177e-05 2.22194e-05 2.22209e-05 2.22224e-05 2.22239e-05 2.22252e-05 2.22265e-05
-          2.22277e-05 2.22288e-05 2.22298e-05 2.22304e-05
+          2.18025e-05 2.18025e-05 2.17329e-05 2.17329e-05 2.18025e-05 2.17329e-05 2.18025e-05 2.17329e-05 2.18025e-05 2.17329e-05 2.18025e-05 2.17329e-05
+          2.18025e-05 2.17329e-05 2.18025e-05 2.17329e-05 2.18025e-05 2.17329e-05 2.18025e-05 2.17329e-05 2.18025e-05 2.17329e-05 2.18025e-05 2.17329e-05
+          2.18025e-05 2.17329e-05 2.18025e-05 2.17329e-05 2.18025e-05 2.17329e-05 2.18025e-05 2.17329e-05 2.16703e-05 2.16703e-05 2.16703e-05 2.16703e-05
+          2.16703e-05 2.16703e-05 2.16703e-05 2.16703e-05 2.16703e-05 2.16703e-05 2.16703e-05 2.16703e-05 2.16703e-05 2.16703e-05 2.16703e-05 2.16703e-05
+          2.16168e-05 2.16168e-05 2.16168e-05 2.16168e-05 2.16168e-05 2.16168e-05 2.16168e-05 2.16168e-05 2.16168e-05 2.16168e-05 2.16168e-05 2.16168e-05
+          2.16168e-05 2.16168e-05 2.16168e-05 2.16168e-05 2.15805e-05 2.15805e-05 2.15805e-05 2.15805e-05 2.15805e-05 2.15805e-05 2.15805e-05 2.15805e-05
+          2.15805e-05 2.15805e-05 2.15805e-05 2.15805e-05 2.15805e-05 2.15805e-05 2.15805e-05 2.15805e-05 2.15879e-05 2.15879e-05 2.15879e-05 2.15879e-05
+          2.15879e-05 2.15879e-05 2.15879e-05 2.1588e-05 2.1588e-05 2.1588e-05 2.1588e-05 2.1588e-05 2.1588e-05 2.1588e-05 2.1588e-05 2.1588e-05
+          2.16268e-05 2.16268e-05 2.16268e-05 2.16268e-05 2.16269e-05 2.16269e-05 2.16269e-05 2.1627e-05 2.1627e-05 2.1627e-05 2.16271e-05 2.16271e-05
+          2.16271e-05 2.16272e-05 2.16272e-05 2.16272e-05 2.16904e-05 2.16904e-05 2.16904e-05 2.16905e-05 2.16906e-05 2.16907e-05 2.16908e-05 2.16909e-05
+          2.16909e-05 2.1691e-05 2.16911e-05 2.16911e-05 2.16912e-05 2.16912e-05 2.16913e-05 2.16913e-05 2.17819e-05 2.17819e-05 2.1782e-05 2.17822e-05
+          2.17823e-05 2.17825e-05 2.17826e-05 2.17828e-05 2.17829e-05 2.17831e-05 2.17832e-05 2.17833e-05 2.17834e-05 2.17835e-05 2.17836e-05 2.17836e-05
+          2.19012e-05 2.19013e-05 2.19015e-05 2.19017e-05 2.1902e-05 2.19023e-05 2.19026e-05 2.19029e-05 2.19031e-05 2.19033e-05 2.19036e-05 2.19038e-05
+          2.1904e-05 2.19041e-05 2.19042e-05 2.19042e-05 2.20448e-05 2.2045e-05 2.20453e-05 2.20457e-05 2.20462e-05 2.20467e-05 2.20471e-05 2.20475e-05
+          2.20479e-05 2.20483e-05 2.20487e-05 2.2049e-05 2.20493e-05 2.20495e-05 2.20497e-05 2.20498e-05 2.22067e-05 2.22069e-05 2.22075e-05 2.22082e-05
+          2.2209e-05 2.22097e-05 2.22103e-05 2.22109e-05 2.22115e-05 2.22121e-05 2.22126e-05 2.22131e-05 2.22135e-05 2.22139e-05 2.22143e-05 2.22144e-05
+          2.238e-05 2.23807e-05 2.23817e-05 2.23827e-05 2.23837e-05 2.23847e-05 2.23857e-05 2.23865e-05 2.23874e-05 2.23881e-05 2.23889e-05 2.23896e-05
+          2.23902e-05 2.23908e-05 2.23913e-05 2.23916e-05
         </DataArray>
         <DataArray type="Float32" Name="X_gas^H2O" NumberOfComponents="1" format="ascii">
-          0.0197275 0.0197275 0.0197883 0.0197883 0.0197275 0.0197883 0.0197275 0.0197883 0.0197275 0.0197883 0.0197275 0.0197883
-          0.0197275 0.0197883 0.0197275 0.0197883 0.0197275 0.0197883 0.0197275 0.0197883 0.0197275 0.0197882 0.0197275 0.0197882
-          0.0197275 0.0197882 0.0197275 0.0197882 0.0197275 0.0197882 0.0197275 0.0197882 0.0198411 0.0198411 0.0198411 0.0198411
-          0.0198411 0.0198411 0.019841 0.019841 0.019841 0.019841 0.019841 0.019841 0.019841 0.019841 0.019841 0.019841
-          0.0198806 0.0198806 0.0198806 0.0198806 0.0198806 0.0198805 0.0198805 0.0198805 0.0198805 0.0198804 0.0198804 0.0198804
-          0.0198804 0.0198804 0.0198804 0.0198804 0.0198932 0.0198932 0.0198932 0.0198931 0.0198931 0.019893 0.0198929 0.0198929
-          0.0198928 0.0198927 0.0198927 0.0198926 0.0198926 0.0198926 0.0198925 0.0198925 0.0198458 0.0198457 0.0198457 0.0198455
-          0.0198454 0.0198452 0.0198451 0.0198449 0.0198447 0.0198446 0.0198444 0.0198443 0.0198442 0.0198441 0.0198441 0.0198441
-          0.0197391 0.0197391 0.0197389 0.0197387 0.0197383 0.019738 0.0197376 0.0197373 0.0197369 0.0197366 0.0197363 0.0197361
-          0.0197358 0.0197356 0.0197355 0.0197355 0.0195716 0.0195715 0.0195711 0.0195706 0.01957 0.0195693 0.0195687 0.019568
-          0.0195673 0.0195667 0.0195662 0.0195657 0.0195652 0.0195649 0.0195647 0.0195646 0.0193388 0.0193386 0.0193379 0.019337
-          0.0193359 0.0193347 0.0193336 0.0193325 0.0193314 0.0193303 0.0193294 0.0193285 0.0193277 0.0193271 0.0193267 0.0193266
-          0.0190475 0.0190471 0.019046 0.0190444 0.0190426 0.0190408 0.019039 0.0190372 0.0190355 0.019034 0.0190325 0.0190311
-          0.0190299 0.0190288 0.0190282 0.0190279 0.0187138 0.018713 0.0187111 0.0187086 0.0187059 0.0187032 0.0187006 0.0186981
-          0.0186957 0.0186935 0.0186913 0.0186893 0.0186875 0.0186859 0.0186848 0.0186844 0.0183578 0.0183563 0.0183532 0.0183496
-          0.0183458 0.0183421 0.0183386 0.0183352 0.018332 0.018329 0.0183261 0.0183233 0.0183208 0.0183186 0.0183168 0.0183161
-          0.0179991 0.0179954 0.0179907 0.0179857 0.0179807 0.017976 0.0179714 0.0179671 0.017963 0.0179592 0.0179555 0.017952
-          0.0179487 0.0179457 0.0179431 0.0179413
+          0.0197274 0.0197274 0.019788 0.019788 0.0197274 0.019788 0.0197274 0.019788 0.0197274 0.019788 0.0197274 0.019788
+          0.0197274 0.019788 0.0197274 0.019788 0.0197274 0.019788 0.0197274 0.019788 0.0197274 0.019788 0.0197274 0.019788
+          0.0197274 0.019788 0.0197274 0.019788 0.0197274 0.019788 0.0197274 0.019788 0.0198402 0.0198402 0.0198402 0.0198402
+          0.0198402 0.0198402 0.0198402 0.0198402 0.0198402 0.0198402 0.0198402 0.0198402 0.0198402 0.0198402 0.0198402 0.0198402
+          0.0198781 0.0198781 0.0198781 0.0198781 0.0198781 0.019878 0.019878 0.019878 0.019878 0.019878 0.019878 0.019878
+          0.019878 0.019878 0.019878 0.019878 0.0198865 0.0198865 0.0198864 0.0198864 0.0198864 0.0198864 0.0198864 0.0198864
+          0.0198863 0.0198863 0.0198863 0.0198863 0.0198863 0.0198863 0.0198863 0.0198863 0.0198297 0.0198297 0.0198297 0.0198296
+          0.0198296 0.0198295 0.0198295 0.0198294 0.0198294 0.0198293 0.0198293 0.0198292 0.0198292 0.0198292 0.0198292 0.0198292
+          0.0197051 0.0197051 0.019705 0.0197049 0.0197048 0.0197047 0.0197046 0.0197044 0.0197043 0.0197042 0.0197041 0.019704
+          0.0197039 0.0197039 0.0197038 0.0197038 0.0195069 0.0195069 0.0195067 0.0195065 0.0195063 0.019506 0.0195057 0.0195055
+          0.0195052 0.019505 0.0195048 0.0195046 0.0195044 0.0195043 0.0195042 0.0195042 0.0192275 0.0192274 0.0192271 0.0192267
+          0.0192262 0.0192257 0.0192252 0.0192247 0.0192243 0.0192239 0.0192235 0.0192232 0.0192228 0.0192226 0.0192224 0.0192224
+          0.0188721 0.0188719 0.0188713 0.0188706 0.0188698 0.0188689 0.0188681 0.0188673 0.0188666 0.0188659 0.0188652 0.0188646
+          0.0188641 0.0188637 0.0188634 0.0188633 0.0184574 0.018457 0.018456 0.0184548 0.0184535 0.0184522 0.0184509 0.0184497
+          0.0184486 0.0184475 0.0184465 0.0184456 0.0184448 0.018444 0.0184435 0.0184433 0.0180062 0.0180054 0.0180038 0.0180019
+          0.0179999 0.017998 0.0179962 0.0179945 0.0179929 0.0179914 0.01799 0.0179886 0.0179874 0.0179863 0.0179855 0.0179851
+          0.0175411 0.0175393 0.0175367 0.017534 0.0175313 0.0175287 0.0175262 0.0175239 0.0175218 0.0175197 0.0175178 0.017516
+          0.0175143 0.0175127 0.0175114 0.0175106
         </DataArray>
         <DataArray type="Float32" Name="X_gas^Air" NumberOfComponents="1" format="ascii">
-          0.980272 0.980272 0.980212 0.980212 0.980272 0.980212 0.980272 0.980212 0.980272 0.980212 0.980272 0.980212
-          0.980272 0.980212 0.980272 0.980212 0.980272 0.980212 0.980272 0.980212 0.980272 0.980212 0.980272 0.980212
-          0.980272 0.980212 0.980272 0.980212 0.980272 0.980212 0.980272 0.980212 0.980159 0.980159 0.980159 0.980159
-          0.980159 0.980159 0.980159 0.980159 0.980159 0.980159 0.980159 0.980159 0.980159 0.980159 0.980159 0.980159
-          0.980119 0.980119 0.980119 0.980119 0.980119 0.980119 0.980119 0.98012 0.98012 0.98012 0.98012 0.98012
-          0.98012 0.98012 0.98012 0.98012 0.980107 0.980107 0.980107 0.980107 0.980107 0.980107 0.980107 0.980107
-          0.980107 0.980107 0.980107 0.980107 0.980107 0.980107 0.980107 0.980107 0.980154 0.980154 0.980154 0.980154
-          0.980155 0.980155 0.980155 0.980155 0.980155 0.980155 0.980156 0.980156 0.980156 0.980156 0.980156 0.980156
-          0.980261 0.980261 0.980261 0.980261 0.980262 0.980262 0.980262 0.980263 0.980263 0.980263 0.980264 0.980264
-          0.980264 0.980264 0.980264 0.980264 0.980428 0.980429 0.980429 0.980429 0.98043 0.980431 0.980431 0.980432
-          0.980433 0.980433 0.980434 0.980434 0.980435 0.980435 0.980435 0.980435 0.980661 0.980661 0.980662 0.980663
-          0.980664 0.980665 0.980666 0.980668 0.980669 0.98067 0.980671 0.980671 0.980672 0.980673 0.980673 0.980673
-          0.980953 0.980953 0.980954 0.980956 0.980957 0.980959 0.980961 0.980963 0.980964 0.980966 0.980968 0.980969
-          0.98097 0.980971 0.980972 0.980972 0.981286 0.981287 0.981289 0.981291 0.981294 0.981297 0.981299 0.981302
-          0.981304 0.981307 0.981309 0.981311 0.981313 0.981314 0.981315 0.981316 0.981642 0.981644 0.981647 0.98165
-          0.981654 0.981658 0.981661 0.981665 0.981668 0.981671 0.981674 0.981677 0.981679 0.981681 0.981683 0.981684
-          0.982001 0.982005 0.982009 0.982014 0.982019 0.982024 0.982029 0.982033 0.982037 0.982041 0.982045 0.982048
-          0.982051 0.982054 0.982057 0.982059
+          0.980273 0.980273 0.980212 0.980212 0.980273 0.980212 0.980273 0.980212 0.980273 0.980212 0.980273 0.980212
+          0.980273 0.980212 0.980273 0.980212 0.980273 0.980212 0.980273 0.980212 0.980273 0.980212 0.980273 0.980212
+          0.980273 0.980212 0.980273 0.980212 0.980273 0.980212 0.980273 0.980212 0.98016 0.98016 0.98016 0.98016
+          0.98016 0.98016 0.98016 0.98016 0.98016 0.98016 0.98016 0.98016 0.98016 0.98016 0.98016 0.98016
+          0.980122 0.980122 0.980122 0.980122 0.980122 0.980122 0.980122 0.980122 0.980122 0.980122 0.980122 0.980122
+          0.980122 0.980122 0.980122 0.980122 0.980114 0.980114 0.980114 0.980114 0.980114 0.980114 0.980114 0.980114
+          0.980114 0.980114 0.980114 0.980114 0.980114 0.980114 0.980114 0.980114 0.98017 0.98017 0.98017 0.98017
+          0.98017 0.98017 0.980171 0.980171 0.980171 0.980171 0.980171 0.980171 0.980171 0.980171 0.980171 0.980171
+          0.980295 0.980295 0.980295 0.980295 0.980295 0.980295 0.980295 0.980296 0.980296 0.980296 0.980296 0.980296
+          0.980296 0.980296 0.980296 0.980296 0.980493 0.980493 0.980493 0.980493 0.980494 0.980494 0.980494 0.980494
+          0.980495 0.980495 0.980495 0.980495 0.980496 0.980496 0.980496 0.980496 0.980773 0.980773 0.980773 0.980773
+          0.980774 0.980774 0.980775 0.980775 0.980776 0.980776 0.980776 0.980777 0.980777 0.980777 0.980778 0.980778
+          0.981128 0.981128 0.981129 0.981129 0.98113 0.981131 0.981132 0.981133 0.981133 0.981134 0.981135 0.981135
+          0.981136 0.981136 0.981137 0.981137 0.981543 0.981543 0.981544 0.981545 0.981547 0.981548 0.981549 0.98155
+          0.981551 0.981552 0.981553 0.981554 0.981555 0.981556 0.981556 0.981557 0.981994 0.981995 0.981996 0.981998
+          0.982 0.982002 0.982004 0.982005 0.982007 0.982009 0.98201 0.982011 0.982013 0.982014 0.982015 0.982015
+          0.982459 0.982461 0.982463 0.982466 0.982469 0.982471 0.982474 0.982476 0.982478 0.98248 0.982482 0.982484
+          0.982486 0.982487 0.982489 0.982489
         </DataArray>
         <DataArray type="Float32" Name="x_liquid^H2O" NumberOfComponents="1" format="ascii">
           0.999986 0.999986 0.999986 0.999986 0.999986 0.999986 0.999986 0.999986 0.999986 0.999986 0.999986 0.999986
@@ -284,64 +284,64 @@
           0.999986 0.999986 0.999986 0.999986
         </DataArray>
         <DataArray type="Float32" Name="x_liquid^Air" NumberOfComponents="1" format="ascii">
-          1.35628e-05 1.35628e-05 1.35194e-05 1.35194e-05 1.35628e-05 1.35194e-05 1.35628e-05 1.35194e-05 1.35628e-05 1.35194e-05 1.35628e-05 1.35194e-05
-          1.35628e-05 1.35194e-05 1.35628e-05 1.35194e-05 1.35628e-05 1.35194e-05 1.35628e-05 1.35194e-05 1.35628e-05 1.35194e-05 1.35628e-05 1.35194e-05
-          1.35628e-05 1.35194e-05 1.35628e-05 1.35194e-05 1.35628e-05 1.35194e-05 1.35628e-05 1.35194e-05 1.34804e-05 1.34804e-05 1.34804e-05 1.34804e-05
-          1.34804e-05 1.34804e-05 1.34804e-05 1.34804e-05 1.34804e-05 1.34804e-05 1.34804e-05 1.34804e-05 1.34804e-05 1.34804e-05 1.34804e-05 1.34804e-05
-          1.34467e-05 1.34467e-05 1.34467e-05 1.34467e-05 1.34467e-05 1.34468e-05 1.34468e-05 1.34468e-05 1.34468e-05 1.34468e-05 1.34468e-05 1.34468e-05
-          1.34468e-05 1.34468e-05 1.34468e-05 1.34468e-05 1.34233e-05 1.34233e-05 1.34233e-05 1.34234e-05 1.34234e-05 1.34234e-05 1.34234e-05 1.34234e-05
-          1.34234e-05 1.34234e-05 1.34234e-05 1.34235e-05 1.34235e-05 1.34235e-05 1.34235e-05 1.34235e-05 1.34261e-05 1.34262e-05 1.34262e-05 1.34262e-05
-          1.34262e-05 1.34262e-05 1.34263e-05 1.34263e-05 1.34263e-05 1.34264e-05 1.34264e-05 1.34264e-05 1.34264e-05 1.34265e-05 1.34265e-05 1.34265e-05
-          1.34467e-05 1.34468e-05 1.34468e-05 1.34468e-05 1.34469e-05 1.3447e-05 1.3447e-05 1.34471e-05 1.34472e-05 1.34472e-05 1.34473e-05 1.34474e-05
-          1.34474e-05 1.34474e-05 1.34475e-05 1.34475e-05 1.34801e-05 1.34801e-05 1.34802e-05 1.34803e-05 1.34804e-05 1.34805e-05 1.34807e-05 1.34808e-05
-          1.34809e-05 1.3481e-05 1.34811e-05 1.34813e-05 1.34813e-05 1.34814e-05 1.34815e-05 1.34815e-05 1.35272e-05 1.35272e-05 1.35273e-05 1.35275e-05
-          1.35278e-05 1.3528e-05 1.35282e-05 1.35285e-05 1.35287e-05 1.35289e-05 1.35291e-05 1.35293e-05 1.35294e-05 1.35295e-05 1.35296e-05 1.35297e-05
-          1.35873e-05 1.35874e-05 1.35876e-05 1.3588e-05 1.35884e-05 1.35887e-05 1.35891e-05 1.35895e-05 1.35898e-05 1.35902e-05 1.35905e-05 1.35908e-05
-          1.3591e-05 1.35912e-05 1.35914e-05 1.35914e-05 1.3658e-05 1.36581e-05 1.36586e-05 1.36591e-05 1.36597e-05 1.36602e-05 1.36608e-05 1.36613e-05
-          1.36619e-05 1.36623e-05 1.36628e-05 1.36632e-05 1.36636e-05 1.3664e-05 1.36642e-05 1.36643e-05 1.37354e-05 1.37358e-05 1.37365e-05 1.37373e-05
-          1.37381e-05 1.37389e-05 1.37397e-05 1.37404e-05 1.37412e-05 1.37418e-05 1.37425e-05 1.37431e-05 1.37436e-05 1.37441e-05 1.37445e-05 1.37447e-05
-          1.38159e-05 1.38167e-05 1.38178e-05 1.38189e-05 1.382e-05 1.38211e-05 1.38221e-05 1.38231e-05 1.38241e-05 1.38249e-05 1.38258e-05 1.38266e-05
-          1.38273e-05 1.3828e-05 1.38286e-05 1.3829e-05
+          1.35628e-05 1.35628e-05 1.35195e-05 1.35195e-05 1.35628e-05 1.35195e-05 1.35628e-05 1.35195e-05 1.35628e-05 1.35195e-05 1.35628e-05 1.35195e-05
+          1.35628e-05 1.35195e-05 1.35628e-05 1.35195e-05 1.35628e-05 1.35195e-05 1.35628e-05 1.35195e-05 1.35628e-05 1.35195e-05 1.35628e-05 1.35195e-05
+          1.35628e-05 1.35195e-05 1.35628e-05 1.35195e-05 1.35628e-05 1.35195e-05 1.35628e-05 1.35195e-05 1.34806e-05 1.34806e-05 1.34806e-05 1.34806e-05
+          1.34806e-05 1.34806e-05 1.34806e-05 1.34806e-05 1.34806e-05 1.34806e-05 1.34806e-05 1.34806e-05 1.34806e-05 1.34806e-05 1.34806e-05 1.34806e-05
+          1.34473e-05 1.34473e-05 1.34473e-05 1.34473e-05 1.34473e-05 1.34473e-05 1.34473e-05 1.34473e-05 1.34473e-05 1.34473e-05 1.34473e-05 1.34473e-05
+          1.34473e-05 1.34473e-05 1.34473e-05 1.34473e-05 1.34247e-05 1.34247e-05 1.34247e-05 1.34247e-05 1.34247e-05 1.34247e-05 1.34247e-05 1.34247e-05
+          1.34247e-05 1.34247e-05 1.34247e-05 1.34247e-05 1.34247e-05 1.34247e-05 1.34247e-05 1.34247e-05 1.34293e-05 1.34293e-05 1.34293e-05 1.34293e-05
+          1.34293e-05 1.34293e-05 1.34293e-05 1.34294e-05 1.34294e-05 1.34294e-05 1.34294e-05 1.34294e-05 1.34294e-05 1.34294e-05 1.34294e-05 1.34294e-05
+          1.34535e-05 1.34535e-05 1.34535e-05 1.34535e-05 1.34536e-05 1.34536e-05 1.34536e-05 1.34536e-05 1.34537e-05 1.34537e-05 1.34537e-05 1.34537e-05
+          1.34537e-05 1.34537e-05 1.34537e-05 1.34537e-05 1.34931e-05 1.34931e-05 1.34931e-05 1.34932e-05 1.34932e-05 1.34933e-05 1.34933e-05 1.34934e-05
+          1.34934e-05 1.34935e-05 1.34935e-05 1.34935e-05 1.34936e-05 1.34936e-05 1.34936e-05 1.34936e-05 1.355e-05 1.355e-05 1.35501e-05 1.35502e-05
+          1.35503e-05 1.35504e-05 1.35505e-05 1.35506e-05 1.35507e-05 1.35507e-05 1.35508e-05 1.35509e-05 1.3551e-05 1.3551e-05 1.3551e-05 1.3551e-05
+          1.36242e-05 1.36243e-05 1.36244e-05 1.36246e-05 1.36247e-05 1.36249e-05 1.36251e-05 1.36253e-05 1.36254e-05 1.36256e-05 1.36257e-05 1.36258e-05
+          1.36259e-05 1.3626e-05 1.36261e-05 1.36261e-05 1.37136e-05 1.37137e-05 1.37139e-05 1.37141e-05 1.37144e-05 1.37147e-05 1.3715e-05 1.37153e-05
+          1.37155e-05 1.37157e-05 1.3716e-05 1.37162e-05 1.37163e-05 1.37165e-05 1.37166e-05 1.37166e-05 1.38142e-05 1.38144e-05 1.38148e-05 1.38152e-05
+          1.38157e-05 1.38161e-05 1.38165e-05 1.38169e-05 1.38173e-05 1.38176e-05 1.38179e-05 1.38182e-05 1.38185e-05 1.38188e-05 1.3819e-05 1.3819e-05
+          1.39221e-05 1.39225e-05 1.39231e-05 1.39238e-05 1.39244e-05 1.3925e-05 1.39256e-05 1.39261e-05 1.39266e-05 1.39271e-05 1.39276e-05 1.3928e-05
+          1.39284e-05 1.39288e-05 1.39291e-05 1.39293e-05
         </DataArray>
         <DataArray type="Float32" Name="x_gas^H2O" NumberOfComponents="1" format="ascii">
-          0.0313371 0.0313371 0.0314324 0.0314324 0.0313371 0.0314324 0.0313371 0.0314324 0.0313371 0.0314324 0.0313371 0.0314324
-          0.0313371 0.0314324 0.0313371 0.0314324 0.0313371 0.0314324 0.0313371 0.0314324 0.0313371 0.0314324 0.0313371 0.0314324
-          0.0313371 0.0314324 0.0313371 0.0314324 0.0313371 0.0314324 0.0313371 0.0314324 0.0315153 0.0315153 0.0315153 0.0315153
-          0.0315153 0.0315153 0.0315153 0.0315152 0.0315152 0.0315152 0.0315152 0.0315152 0.0315152 0.0315152 0.0315152 0.0315152
-          0.0315774 0.0315774 0.0315773 0.0315773 0.0315773 0.0315773 0.0315772 0.0315772 0.0315771 0.0315771 0.0315771 0.031577
-          0.031577 0.031577 0.031577 0.031577 0.0315971 0.0315971 0.0315971 0.031597 0.0315969 0.0315968 0.0315967 0.0315966
-          0.0315965 0.0315964 0.0315963 0.0315962 0.0315961 0.0315961 0.0315961 0.031596 0.0315226 0.0315226 0.0315225 0.0315223
-          0.0315221 0.0315218 0.0315216 0.0315213 0.0315211 0.0315208 0.0315206 0.0315204 0.0315202 0.0315201 0.03152 0.03152
-          0.0313553 0.0313552 0.0313549 0.0313545 0.0313541 0.0313535 0.031353 0.0313524 0.0313519 0.0313514 0.0313509 0.0313505
-          0.0313501 0.0313498 0.0313496 0.0313496 0.0310923 0.0310921 0.0310915 0.0310907 0.0310898 0.0310887 0.0310877 0.0310866
-          0.0310856 0.0310847 0.0310838 0.031083 0.0310823 0.0310818 0.0310814 0.0310813 0.0307267 0.0307264 0.0307254 0.0307239
-          0.0307222 0.0307204 0.0307186 0.0307168 0.0307151 0.0307135 0.030712 0.0307106 0.0307094 0.0307084 0.0307078 0.0307076
-          0.0302692 0.0302686 0.0302668 0.0302643 0.0302615 0.0302586 0.0302558 0.0302531 0.0302504 0.0302479 0.0302456 0.0302434
-          0.0302415 0.0302399 0.0302388 0.0302384 0.0297448 0.0297436 0.0297406 0.0297367 0.0297325 0.0297283 0.0297242 0.0297202
-          0.0297165 0.0297129 0.0297095 0.0297064 0.0297035 0.0297011 0.0296993 0.0296987 0.0291853 0.029183 0.0291781 0.0291723
-          0.0291664 0.0291606 0.0291551 0.0291498 0.0291447 0.0291399 0.0291354 0.0291311 0.0291271 0.0291236 0.0291209 0.0291196
-          0.0286211 0.0286154 0.028608 0.0286001 0.0285923 0.0285848 0.0285777 0.0285709 0.0285645 0.0285584 0.0285525 0.028547
-          0.0285419 0.0285371 0.0285331 0.0285303
+          0.031337 0.031337 0.0314319 0.0314319 0.031337 0.0314319 0.031337 0.0314319 0.031337 0.0314319 0.031337 0.0314319
+          0.031337 0.0314319 0.031337 0.0314319 0.031337 0.0314319 0.031337 0.0314319 0.031337 0.0314319 0.031337 0.0314319
+          0.031337 0.0314319 0.031337 0.0314319 0.031337 0.0314319 0.031337 0.0314319 0.0315139 0.0315139 0.0315139 0.0315139
+          0.0315139 0.0315139 0.0315139 0.0315139 0.0315139 0.0315139 0.0315139 0.0315139 0.0315139 0.0315139 0.0315139 0.0315139
+          0.0315734 0.0315734 0.0315733 0.0315733 0.0315733 0.0315733 0.0315733 0.0315733 0.0315733 0.0315733 0.0315733 0.0315733
+          0.0315733 0.0315733 0.0315732 0.0315732 0.0315865 0.0315865 0.0315865 0.0315865 0.0315865 0.0315864 0.0315864 0.0315864
+          0.0315863 0.0315863 0.0315863 0.0315862 0.0315862 0.0315862 0.0315862 0.0315862 0.0314975 0.0314974 0.0314974 0.0314973
+          0.0314973 0.0314972 0.0314971 0.031497 0.0314969 0.0314969 0.0314968 0.0314967 0.0314967 0.0314966 0.0314966 0.0314966
+          0.0313019 0.0313018 0.0313017 0.0313016 0.0313014 0.0313012 0.031301 0.0313008 0.0313006 0.0313005 0.0313003 0.0313001
+          0.0313 0.0312999 0.0312999 0.0312999 0.0309907 0.0309906 0.0309904 0.0309901 0.0309897 0.0309893 0.0309889 0.0309885
+          0.0309881 0.0309877 0.0309874 0.0309871 0.0309869 0.0309867 0.0309865 0.0309865 0.0305519 0.0305518 0.0305513 0.0305507
+          0.0305499 0.0305491 0.0305484 0.0305476 0.0305469 0.0305463 0.0305457 0.0305451 0.0305446 0.0305443 0.030544 0.0305439
+          0.0299936 0.0299933 0.0299924 0.0299913 0.02999 0.0299886 0.0299873 0.0299861 0.029985 0.0299839 0.0299829 0.0299819
+          0.0299811 0.0299804 0.0299799 0.0299798 0.0293418 0.0293412 0.0293397 0.0293378 0.0293357 0.0293336 0.0293316 0.0293297
+          0.029328 0.0293263 0.0293247 0.0293233 0.029322 0.0293208 0.02932 0.0293197 0.0286323 0.0286311 0.0286286 0.0286256
+          0.0286225 0.0286195 0.0286166 0.028614 0.0286114 0.028609 0.0286068 0.0286047 0.0286028 0.0286011 0.0285997 0.0285991
+          0.0279006 0.0278977 0.0278937 0.0278894 0.0278851 0.027881 0.0278771 0.0278735 0.0278701 0.0278669 0.0278639 0.027861
+          0.0278583 0.0278559 0.0278539 0.0278525
         </DataArray>
         <DataArray type="Float32" Name="x_gas^Air" NumberOfComponents="1" format="ascii">
           0.968663 0.968663 0.968568 0.968568 0.968663 0.968568 0.968663 0.968568 0.968663 0.968568 0.968663 0.968568
           0.968663 0.968568 0.968663 0.968568 0.968663 0.968568 0.968663 0.968568 0.968663 0.968568 0.968663 0.968568
-          0.968663 0.968568 0.968663 0.968568 0.968663 0.968568 0.968663 0.968568 0.968485 0.968485 0.968485 0.968485
-          0.968485 0.968485 0.968485 0.968485 0.968485 0.968485 0.968485 0.968485 0.968485 0.968485 0.968485 0.968485
-          0.968423 0.968423 0.968423 0.968423 0.968423 0.968423 0.968423 0.968423 0.968423 0.968423 0.968423 0.968423
-          0.968423 0.968423 0.968423 0.968423 0.968403 0.968403 0.968403 0.968403 0.968403 0.968403 0.968403 0.968403
-          0.968404 0.968404 0.968404 0.968404 0.968404 0.968404 0.968404 0.968404 0.968477 0.968477 0.968477 0.968478
-          0.968478 0.968478 0.968478 0.968479 0.968479 0.968479 0.968479 0.96848 0.96848 0.96848 0.96848 0.96848
-          0.968645 0.968645 0.968645 0.968645 0.968646 0.968646 0.968647 0.968648 0.968648 0.968649 0.968649 0.96865
-          0.96865 0.96865 0.96865 0.96865 0.968908 0.968908 0.968908 0.968909 0.96891 0.968911 0.968912 0.968913
-          0.968914 0.968915 0.968916 0.968917 0.968918 0.968918 0.968919 0.968919 0.969273 0.969274 0.969275 0.969276
-          0.969278 0.96928 0.969281 0.969283 0.969285 0.969287 0.969288 0.969289 0.969291 0.969292 0.969292 0.969292
-          0.969731 0.969731 0.969733 0.969736 0.969738 0.969741 0.969744 0.969747 0.96975 0.969752 0.969754 0.969757
-          0.969759 0.96976 0.969761 0.969762 0.970255 0.970256 0.970259 0.970263 0.970267 0.970272 0.970276 0.97028
-          0.970284 0.970287 0.97029 0.970294 0.970296 0.970299 0.970301 0.970301 0.970815 0.970817 0.970822 0.970828
-          0.970834 0.970839 0.970845 0.97085 0.970855 0.97086 0.970865 0.970869 0.970873 0.970876 0.970879 0.97088
-          0.971379 0.971385 0.971392 0.9714 0.971408 0.971415 0.971422 0.971429 0.971436 0.971442 0.971447 0.971453
-          0.971458 0.971463 0.971467 0.97147
+          0.968663 0.968568 0.968663 0.968568 0.968663 0.968568 0.968663 0.968568 0.968486 0.968486 0.968486 0.968486
+          0.968486 0.968486 0.968486 0.968486 0.968486 0.968486 0.968486 0.968486 0.968486 0.968486 0.968486 0.968486
+          0.968427 0.968427 0.968427 0.968427 0.968427 0.968427 0.968427 0.968427 0.968427 0.968427 0.968427 0.968427
+          0.968427 0.968427 0.968427 0.968427 0.968413 0.968413 0.968413 0.968414 0.968414 0.968414 0.968414 0.968414
+          0.968414 0.968414 0.968414 0.968414 0.968414 0.968414 0.968414 0.968414 0.968503 0.968503 0.968503 0.968503
+          0.968503 0.968503 0.968503 0.968503 0.968503 0.968503 0.968503 0.968503 0.968503 0.968503 0.968503 0.968503
+          0.968698 0.968698 0.968698 0.968698 0.968699 0.968699 0.968699 0.968699 0.968699 0.9687 0.9687 0.9687
+          0.9687 0.9687 0.9687 0.9687 0.969009 0.969009 0.96901 0.96901 0.96901 0.969011 0.969011 0.969012
+          0.969012 0.969012 0.969013 0.969013 0.969013 0.969013 0.969013 0.969014 0.969448 0.969448 0.969449 0.969449
+          0.96945 0.969451 0.969452 0.969452 0.969453 0.969454 0.969454 0.969455 0.969455 0.969456 0.969456 0.969456
+          0.970006 0.970007 0.970008 0.970009 0.97001 0.970011 0.970013 0.970014 0.970015 0.970016 0.970017 0.970018
+          0.970019 0.97002 0.97002 0.97002 0.970658 0.970659 0.97066 0.970662 0.970664 0.970666 0.970668 0.97067
+          0.970672 0.970674 0.970675 0.970677 0.970678 0.970679 0.97068 0.97068 0.971368 0.971369 0.971371 0.971374
+          0.971377 0.971381 0.971383 0.971386 0.971389 0.971391 0.971393 0.971395 0.971397 0.971399 0.9714 0.971401
+          0.972099 0.972102 0.972106 0.972111 0.972115 0.972119 0.972123 0.972126 0.97213 0.972133 0.972136 0.972139
+          0.972142 0.972144 0.972146 0.972147
         </DataArray>
         <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
           0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41 0.41
@@ -366,22 +366,22 @@
         <DataArray type="Float32" Name="temperature" NumberOfComponents="1" format="ascii">
           298.15 298.15 298.148 298.148 298.15 298.148 298.15 298.148 298.15 298.148 298.15 298.148
           298.15 298.148 298.15 298.148 298.15 298.148 298.15 298.148 298.15 298.148 298.15 298.148
-          298.15 298.148 298.15 298.148 298.15 298.148 298.15 298.148 298.145 298.145 298.145 298.145
-          298.145 298.145 298.145 298.145 298.145 298.145 298.145 298.145 298.145 298.144 298.144 298.144
-          298.134 298.134 298.134 298.134 298.134 298.134 298.134 298.134 298.134 298.134 298.134 298.134
-          298.134 298.134 298.134 298.134 298.109 298.109 298.109 298.109 298.109 298.109 298.109 298.109
-          298.109 298.109 298.109 298.109 298.109 298.109 298.109 298.109 298.059 298.059 298.058 298.058
-          298.058 298.058 298.058 298.058 298.058 298.058 298.057 298.057 298.057 298.057 298.057 298.057
-          297.969 297.969 297.968 297.968 297.968 297.968 297.967 297.967 297.967 297.966 297.966 297.966
-          297.966 297.966 297.966 297.966 297.827 297.827 297.827 297.827 297.826 297.826 297.825 297.824
-          297.824 297.823 297.823 297.822 297.822 297.822 297.822 297.821 297.63 297.629 297.629 297.628
-          297.627 297.626 297.625 297.624 297.623 297.622 297.622 297.621 297.62 297.62 297.619 297.619
-          297.379 297.379 297.378 297.376 297.375 297.373 297.372 297.37 297.369 297.367 297.366 297.365
-          297.364 297.363 297.362 297.362 297.088 297.087 297.086 297.083 297.081 297.079 297.076 297.074
-          297.072 297.07 297.068 297.066 297.065 297.064 297.063 297.062 296.772 296.771 296.768 296.765
-          296.762 296.758 296.755 296.752 296.749 296.747 296.744 296.742 296.739 296.737 296.736 296.735
-          296.449 296.445 296.441 296.437 296.432 296.428 296.424 296.42 296.416 296.412 296.409 296.406
-          296.403 296.4 296.398 296.396
+          298.15 298.148 298.15 298.148 298.15 298.148 298.15 298.148 298.144 298.144 298.144 298.144
+          298.144 298.144 298.144 298.144 298.144 298.144 298.144 298.144 298.144 298.144 298.144 298.144
+          298.132 298.132 298.132 298.132 298.132 298.132 298.132 298.132 298.132 298.132 298.132 298.132
+          298.132 298.132 298.132 298.132 298.104 298.103 298.103 298.103 298.103 298.103 298.103 298.103
+          298.103 298.103 298.103 298.103 298.103 298.103 298.103 298.103 298.045 298.045 298.045 298.045
+          298.045 298.045 298.045 298.045 298.045 298.045 298.045 298.045 298.045 298.045 298.045 298.045
+          297.94 297.94 297.94 297.94 297.94 297.94 297.94 297.939 297.939 297.939 297.939 297.939
+          297.939 297.939 297.939 297.939 297.773 297.773 297.772 297.772 297.772 297.772 297.772 297.771
+          297.771 297.771 297.771 297.771 297.771 297.77 297.77 297.77 297.534 297.534 297.534 297.534
+          297.533 297.533 297.532 297.532 297.532 297.531 297.531 297.531 297.53 297.53 297.53 297.53
+          297.227 297.227 297.226 297.225 297.225 297.224 297.223 297.223 297.222 297.221 297.221 297.22
+          297.22 297.219 297.219 297.219 296.861 296.861 296.86 296.859 296.858 296.857 296.855 296.854
+          296.853 296.852 296.852 296.851 296.85 296.849 296.849 296.849 296.455 296.455 296.453 296.451
+          296.45 296.448 296.446 296.445 296.443 296.442 296.44 296.439 296.438 296.437 296.436 296.436
+          296.027 296.025 296.023 296.021 296.018 296.016 296.013 296.011 296.009 296.007 296.005 296.004
+          296.002 296.001 296 295.999
         </DataArray>
         <DataArray type="Float32" Name="phase presence" NumberOfComponents="1" format="ascii">
           3 3 3 3 3 3 3 3 3 3 3 3
@@ -404,112 +404,112 @@
           3 3 3 3
         </DataArray>
         <DataArray type="Float32" Name="velocityW" NumberOfComponents="3" format="ascii">
-          -3.56641e-11 -1.24106e-10 0 -7.76879e-11 -1.24612e-10 0 -3.38168e-11 -4.20353e-10 0 -8.9698e-11 -4.24211e-10 0
-          -1.48634e-10 -1.37998e-10 0 -1.68339e-10 -4.41232e-10 0 -2.00821e-10 -1.4449e-10 0 -2.20091e-10 -4.64923e-10 0
-          -2.5825e-10 -1.57096e-10 0 -2.75364e-10 -4.86643e-10 0 -2.98301e-10 -1.62105e-10 0 -3.20197e-10 -5.09969e-10 0
-          -3.23112e-10 -1.81783e-10 0 -3.47714e-10 -5.32533e-10 0 -3.44617e-10 -1.91275e-10 0 -3.7458e-10 -5.6069e-10 0
-          -3.3566e-10 -2.13891e-10 0 -3.70096e-10 -5.98372e-10 0 -3.14953e-10 -2.28106e-10 0 -3.42589e-10 -6.29672e-10 0
-          -2.87728e-10 -2.43479e-10 0 -3.04214e-10 -6.57082e-10 0 -2.48556e-10 -2.41699e-10 0 -2.5345e-10 -6.76886e-10 0
-          -1.96387e-10 -2.46269e-10 0 -2.06916e-10 -6.86587e-10 0 -1.35373e-10 -2.55403e-10 0 -1.5923e-10 -7.00332e-10 0
-          -6.51031e-11 -2.77986e-10 0 -7.99129e-11 -7.10421e-10 0 -2.47953e-11 -2.76644e-10 0 -1.92573e-11 -7.15827e-10 0
-          -4.30258e-11 -1.13008e-09 0 -1.11868e-10 -1.13618e-09 0 -2.21442e-10 -1.16443e-09 0 -2.77349e-10 -1.20842e-09 0
-          -3.35457e-10 -1.25053e-09 0 -3.80208e-10 -1.31068e-09 0 -4.10146e-10 -1.36968e-09 0 -4.62049e-10 -1.43825e-09 0
-          -4.53732e-10 -1.5093e-09 0 -4.10532e-10 -1.56156e-09 0 -3.68656e-10 -1.61053e-09 0 -2.96671e-10 -1.6542e-09 0
-          -2.37168e-10 -1.68818e-09 0 -1.80805e-10 -1.71974e-09 0 -9.06422e-11 -1.73972e-09 0 -3.85482e-11 -1.7379e-09 0
-          -4.28831e-11 -2.17716e-09 0 -1.36324e-10 -2.18184e-09 0 -2.87379e-10 -2.21472e-09 0 -3.51507e-10 -2.29365e-09 0
-          -4.24066e-10 -2.36226e-09 0 -5.00339e-10 -2.45091e-09 0 -5.3949e-10 -2.56413e-09 0 -5.93849e-10 -2.6804e-09 0
-          -5.56182e-10 -2.78285e-09 0 -4.8949e-10 -2.8609e-09 0 -4.45075e-10 -2.93534e-09 0 -3.71239e-10 -3.00983e-09 0
-          -3.02877e-10 -3.08797e-09 0 -2.37534e-10 -3.12213e-09 0 -1.05124e-10 -3.15381e-09 0 -2.70378e-12 -3.1578e-09 0
-          -4.9946e-11 -3.01506e-09 0 -1.51014e-10 -3.04146e-09 0 -3.66972e-10 -3.0957e-09 0 -4.63625e-10 -3.1899e-09 0
-          -5.22116e-10 -3.32126e-09 0 -6.31349e-10 -3.4567e-09 0 -6.93592e-10 -3.59784e-09 0 -7.2863e-10 -3.72382e-09 0
-          -6.65968e-10 -3.8325e-09 0 -5.86983e-10 -3.95691e-09 0 -5.42759e-10 -4.06425e-09 0 -4.86121e-10 -4.16517e-09 0
-          -3.70568e-10 -4.24949e-09 0 -2.54678e-10 -4.30284e-09 0 -1.3297e-10 -4.32824e-09 0 -3.50355e-11 -4.35346e-09 0
-          -1.0247e-10 -3.19095e-09 0 -2.09473e-10 -3.23763e-09 0 -4.17117e-10 -3.34023e-09 0 -5.63497e-10 -3.46902e-09 0
-          -6.65762e-10 -3.66851e-09 0 -7.4651e-10 -3.8614e-09 0 -7.76597e-10 -4.03783e-09 0 -7.76333e-10 -4.18944e-09 0
-          -7.39952e-10 -4.32666e-09 0 -6.79352e-10 -4.48185e-09 0 -6.11473e-10 -4.62032e-09 0 -5.30977e-10 -4.73728e-09 0
-          -4.1669e-10 -4.8262e-09 0 -2.86787e-10 -4.91778e-09 0 -1.50691e-10 -4.96354e-09 0 -7.31763e-11 -4.98643e-09 0
-          -1.3012e-10 -2.89044e-09 0 -2.55214e-10 -2.94727e-09 0 -4.79986e-10 -3.10112e-09 0 -6.4503e-10 -3.31725e-09 0
-          -7.47179e-10 -3.56688e-09 0 -7.9758e-10 -3.81077e-09 0 -8.07323e-10 -4.03756e-09 0 -7.86296e-10 -4.24425e-09 0
-          -7.44591e-10 -4.42847e-09 0 -6.88844e-10 -4.59762e-09 0 -6.18815e-10 -4.75558e-09 0 -5.31516e-10 -4.90007e-09 0
-          -4.27166e-10 -5.02799e-09 0 -3.02083e-10 -5.13689e-09 0 -1.55695e-10 -5.20905e-09 0 -7.85087e-11 -5.23399e-09 0
-          -1.45041e-10 -2.20405e-09 0 -2.74381e-10 -2.28792e-09 0 -4.96596e-10 -2.50887e-09 0 -6.41174e-10 -2.80567e-09 0
-          -7.09897e-10 -3.12199e-09 0 -7.21991e-10 -3.42076e-09 0 -6.99472e-10 -3.68484e-09 0 -6.57629e-10 -3.91159e-09 0
-          -6.07706e-10 -4.10997e-09 0 -5.55621e-10 -4.28822e-09 0 -4.99788e-10 -4.45451e-09 0 -4.35631e-10 -4.61426e-09 0
-          -3.58392e-10 -4.76624e-09 0 -2.60165e-10 -4.89833e-09 0 -1.37895e-10 -4.99315e-09 0 -7.13467e-11 -5.02817e-09 0
-          -1.34334e-10 -1.24087e-09 0 -2.43511e-10 -1.36339e-09 0 -4.11009e-10 -1.67391e-09 0 -4.8413e-10 -2.06069e-09 0
-          -4.82162e-10 -2.44368e-09 0 -4.36612e-10 -2.78124e-09 0 -3.72438e-10 -3.06529e-09 0 -3.09632e-10 -3.29292e-09 0
-          -2.59577e-10 -3.48681e-09 0 -2.22957e-10 -3.66064e-09 0 -1.99247e-10 -3.82615e-09 0 -1.80548e-10 -3.99108e-09 0
-          -1.62054e-10 -4.155e-09 0 -1.30493e-10 -4.31216e-09 0 -7.42121e-11 -4.43107e-09 0 -3.99645e-11 -4.47647e-09 0
-          -8.25133e-11 3.06505e-09 0 -1.23024e-10 2.87949e-09 0 -1.47934e-10 2.47388e-09 0 -8.76472e-11 2.01363e-09 0
-          1.75436e-11 1.58837e-09 0 1.31696e-10 1.23599e-09 0 2.37092e-10 9.64705e-10 0 3.13858e-10 7.50312e-10 0
-          3.497e-10 5.77496e-10 0 3.54383e-10 4.25329e-10 0 3.29006e-10 2.78514e-10 0 2.8177e-10 1.25244e-10 0
-          2.10815e-10 -3.67346e-11 0 1.31202e-10 -2.00502e-10 0 6.04518e-11 -3.35835e-10 0 2.636e-11 -3.99015e-10 0
-          3.26395e-11 1.53182e-08 0 1.45696e-10 1.50419e-08 0 3.74555e-10 1.45739e-08 0 5.95307e-10 1.40919e-08 0
-          7.94117e-10 1.36682e-08 0 9.65717e-10 1.3332e-08 0 1.08745e-09 1.30946e-08 0 1.15027e-09 1.29159e-08 0
-          1.16146e-09 1.27742e-08 0 1.12327e-09 1.26527e-08 0 1.0414e-09 1.25331e-08 0 9.20377e-10 1.24086e-08 0
-          7.54159e-10 1.22717e-08 0 5.44452e-10 1.21301e-08 0 2.84836e-10 1.20007e-08 0 1.38421e-10 1.19349e-08 0
-          2.55931e-10 3.89235e-08 0 5.41134e-10 3.85663e-08 0 9.89996e-10 3.817e-08 0 1.27958e-09 3.7724e-08 0
-          1.49668e-09 3.73426e-08 0 1.66269e-09 3.70345e-08 0 1.77e-09 3.68389e-08 0 1.79859e-09 3.67191e-08 0
-          1.77102e-09 3.66144e-08 0 1.6977e-09 3.65227e-08 0 1.58208e-09 3.64284e-08 0 1.42736e-09 3.63347e-08 0
-          1.21624e-09 3.62455e-08 0 9.34451e-10 3.61477e-08 0 5.57205e-10 3.60747e-08 0 3.39342e-10 3.61001e-08 0
-          5.05585e-10 5.40328e-08 0 8.35071e-10 5.36458e-08 0 1.29006e-09 5.33383e-08 0 1.53014e-09 5.29221e-08 0
-          1.71001e-09 5.25684e-08 0 1.8416e-09 5.22718e-08 0 1.93118e-09 5.2102e-08 0 1.93447e-09 5.20124e-08 0
-          1.88444e-09 5.19276e-08 0 1.80113e-09 5.18508e-08 0 1.68291e-09 5.17696e-08 0 1.5389e-09 5.16887e-08 0
-          1.34172e-09 5.16277e-08 0 1.06935e-09 5.15557e-08 0 7.3647e-10 5.15286e-08 0 5.57852e-10 5.1636e-08 0
+          -3.43867e-11 -1.68618e-10 0 -3.81932e-11 -1.84462e-10 0 -6.30506e-11 -5.05723e-10 0 -4.38349e-11 -5.04795e-10 0
+          -6.9889e-11 -1.73591e-10 0 -6.1387e-11 -5.08891e-10 0 -1.05558e-10 -1.67396e-10 0 -1.07626e-10 -5.20193e-10 0
+          -1.34029e-10 -1.7107e-10 0 -1.46785e-10 -5.24848e-10 0 -1.60846e-10 -1.83728e-10 0 -1.77657e-10 -5.37268e-10 0
+          -1.67747e-10 -1.88318e-10 0 -1.74661e-10 -5.53815e-10 0 -1.72646e-10 -1.88831e-10 0 -1.86475e-10 -5.6703e-10 0
+          -1.79521e-10 -2.03556e-10 0 -2.05851e-10 -5.85014e-10 0 -1.5921e-10 -2.20046e-10 0 -1.68819e-10 -5.98054e-10 0
+          -1.50105e-10 -2.13484e-10 0 -1.6393e-10 -6.17866e-10 0 -1.37497e-10 -2.35648e-10 0 -1.55103e-10 -6.24774e-10 0
+          -9.51318e-11 -2.35419e-10 0 -8.98215e-11 -6.32334e-10 0 -6.77782e-11 -2.28342e-10 0 -7.50663e-11 -6.38227e-10 0
+          -3.39253e-11 -2.40928e-10 0 -4.48883e-11 -6.47664e-10 0 -1.16279e-11 -2.40864e-10 0 -8.53608e-12 -6.50509e-10 0
+          -3.35342e-11 -1.29601e-09 0 -4.06684e-11 -1.29895e-09 0 -9.73514e-11 -1.30166e-09 0 -1.38976e-10 -1.32693e-09 0
+          -1.64301e-10 -1.35862e-09 0 -2.18161e-10 -1.38118e-09 0 -2.22239e-10 -1.4283e-09 0 -2.29864e-10 -1.45712e-09 0
+          -2.32862e-10 -1.49761e-09 0 -2.14102e-10 -1.51123e-09 0 -1.96247e-10 -1.5427e-09 0 -1.57059e-10 -1.54799e-09 0
+          -1.20273e-10 -1.56029e-09 0 -9.99866e-11 -1.58758e-09 0 -5.78897e-11 -1.60614e-09 0 -1.82255e-11 -1.61002e-09 0
+          -6.61361e-11 -2.45327e-09 0 -4.34786e-11 -2.44528e-09 0 -1.05496e-10 -2.44928e-09 0 -2.06266e-10 -2.48792e-09 0
+          -2.27191e-10 -2.54241e-09 0 -2.89322e-10 -2.59114e-09 0 -2.96918e-10 -2.64428e-09 0 -2.96276e-10 -2.70295e-09 0
+          -2.92156e-10 -2.74164e-09 0 -2.51528e-10 -2.79984e-09 0 -2.29059e-10 -2.84657e-09 0 -1.72823e-10 -2.86521e-09 0
+          -1.56603e-10 -2.89104e-09 0 -1.50341e-10 -2.92228e-09 0 -7.83446e-11 -2.93845e-09 0 -1.76041e-11 -2.94344e-09 0
+          -5.10719e-12 -3.37118e-09 0 -3.30067e-11 -3.3772e-09 0 -1.60414e-10 -3.4347e-09 0 -2.75071e-10 -3.48771e-09 0
+          -3.05516e-10 -3.55204e-09 0 -3.53431e-10 -3.62626e-09 0 -3.81156e-10 -3.68178e-09 0 -3.61524e-10 -3.76319e-09 0
+          -3.59535e-10 -3.81779e-09 0 -3.62694e-10 -3.89437e-09 0 -2.81695e-10 -3.95418e-09 0 -2.10417e-10 -4.01193e-09 0
+          -2.04882e-10 -4.06202e-09 0 -1.6683e-10 -4.08706e-09 0 -8.33102e-11 -4.08527e-09 0 -3.21784e-11 -4.09054e-09 0
+          -5.28702e-11 -3.69318e-09 0 -1.12153e-10 -3.73132e-09 0 -2.34325e-10 -3.80998e-09 0 -3.33695e-10 -3.88244e-09 0
+          -3.89815e-10 -3.97743e-09 0 -4.25912e-10 -4.08155e-09 0 -4.48547e-10 -4.17625e-09 0 -4.42388e-10 -4.26585e-09 0
+          -4.2847e-10 -4.3603e-09 0 -4.0449e-10 -4.42585e-09 0 -3.48951e-10 -4.48553e-09 0 -2.8969e-10 -4.57128e-09 0
+          -2.40997e-10 -4.62998e-09 0 -1.67984e-10 -4.66143e-09 0 -7.86297e-11 -4.67751e-09 0 -3.97125e-11 -4.68468e-09 0
+          -7.97303e-11 -3.61868e-09 0 -1.53591e-10 -3.64858e-09 0 -2.91811e-10 -3.73071e-09 0 -3.99602e-10 -3.8469e-09 0
+          -4.66698e-10 -3.98128e-09 0 -5.0062e-10 -4.12111e-09 0 -5.07336e-10 -4.25212e-09 0 -4.95281e-10 -4.36729e-09 0
+          -4.69022e-10 -4.4713e-09 0 -4.31179e-10 -4.56015e-09 0 -3.87048e-10 -4.64271e-09 0 -3.28614e-10 -4.71866e-09 0
+          -2.61156e-10 -4.78486e-09 0 -1.83005e-10 -4.83342e-09 0 -9.16099e-11 -4.87234e-09 0 -4.53115e-11 -4.88632e-09 0
+          -9.35897e-11 -3.12763e-09 0 -1.80579e-10 -3.17209e-09 0 -3.34526e-10 -3.29344e-09 0 -4.43803e-10 -3.45909e-09 0
+          -5.04985e-10 -3.63885e-09 0 -5.25955e-10 -3.81303e-09 0 -5.17418e-10 -3.96876e-09 0 -4.92055e-10 -4.10395e-09 0
+          -4.57583e-10 -4.21903e-09 0 -4.16736e-10 -4.32183e-09 0 -3.7149e-10 -4.41359e-09 0 -3.18171e-10 -4.49678e-09 0
+          -2.55382e-10 -4.57511e-09 0 -1.81496e-10 -4.64089e-09 0 -9.49612e-11 -4.68761e-09 0 -4.86899e-11 -4.70539e-09 0
+          -9.8924e-11 -2.3625e-09 0 -1.86966e-10 -2.42746e-09 0 -3.35956e-10 -2.59743e-09 0 -4.27137e-10 -2.81885e-09 0
+          -4.64739e-10 -3.04838e-09 0 -4.6226e-10 -3.25666e-09 0 -4.32868e-10 -3.4338e-09 0 -3.93566e-10 -3.57795e-09 0
+          -3.53244e-10 -3.69893e-09 0 -3.15177e-10 -3.80011e-09 0 -2.78075e-10 -3.8969e-09 0 -2.39426e-10 -3.98505e-09 0
+          -1.96583e-10 -4.07274e-09 0 -1.43859e-10 -4.15155e-09 0 -7.80795e-11 -4.21354e-09 0 -4.12293e-11 -4.23844e-09 0
+          -8.03366e-11 1.67509e-09 0 -1.4719e-10 1.5812e-09 0 -2.5018e-10 1.35861e-09 0 -2.9877e-10 1.08636e-09 0
+          -2.97963e-10 8.13275e-10 0 -2.62949e-10 5.79534e-10 0 -2.10885e-10 3.89537e-10 0 -1.61511e-10 2.44837e-10 0
+          -1.20007e-10 1.27085e-10 0 -9.28261e-11 2.95757e-11 0 -7.43822e-11 -6.04444e-11 0 -6.06659e-11 -1.48476e-10 0
+          -5.40526e-11 -2.40362e-10 0 -4.51741e-11 -3.2924e-10 0 -3.10367e-11 -4.02706e-10 0 -1.91759e-11 -4.32327e-10 0
+          -1.52374e-11 1.34737e-08 0 -1.49734e-11 1.33506e-08 0 -1.32398e-11 1.31045e-08 0 -3.69057e-12 1.27988e-08 0
+          3.57742e-11 1.24959e-08 0 9.80083e-11 1.22396e-08 0 1.64955e-10 1.205e-08 0 2.16361e-10 1.19168e-08 0
+          2.43987e-10 1.18168e-08 0 2.54209e-10 1.17254e-08 0 2.47099e-10 1.16506e-08 0 2.22958e-10 1.15677e-08 0
+          1.8466e-10 1.14801e-08 0 1.35165e-10 1.13931e-08 0 7.08694e-11 1.13213e-08 0 3.68992e-11 1.12848e-08 0
+          1.4732e-10 3.63805e-08 0 2.29421e-10 3.62907e-08 0 3.24015e-10 3.60959e-08 0 3.55747e-10 3.57702e-08 0
+          4.06372e-10 3.54478e-08 0 4.77494e-10 3.51692e-08 0 5.49299e-10 3.49928e-08 0 5.89264e-10 3.48863e-08 0
+          5.92674e-10 3.48066e-08 0 5.7909e-10 3.47365e-08 0 5.48443e-10 3.46652e-08 0 4.95868e-10 3.45882e-08 0
+          4.28389e-10 3.45166e-08 0 3.38623e-10 3.44463e-08 0 2.0524e-10 3.43965e-08 0 1.2517e-10 3.43911e-08 0
+          3.50004e-10 5.10695e-08 0 4.36754e-10 5.10246e-08 0 5.1328e-10 5.08795e-08 0 5.24483e-10 5.05445e-08 0
+          5.67322e-10 5.0215e-08 0 6.33413e-10 4.99264e-08 0 7.05339e-10 4.97589e-08 0 7.30714e-10 4.96692e-08 0
+          7.20909e-10 4.95964e-08 0 6.9352e-10 4.9542e-08 0 6.4881e-10 4.94676e-08 0 5.94802e-10 4.93908e-08 0
+          5.25702e-10 4.93306e-08 0 4.29484e-10 4.92699e-08 0 3.01794e-10 4.9238e-08 0 2.29973e-10 4.92716e-08 0
         </DataArray>
         <DataArray type="Float32" Name="velocityN" NumberOfComponents="3" format="ascii">
-          -1.51179e-18 2.1257e-10 0 -3.16051e-17 2.12566e-10 0 -8.46915e-16 6.19375e-10 0 -6.75687e-16 6.19386e-10 0
-          -2.18706e-17 2.12568e-10 0 3.52814e-17 6.1939e-10 0 3.43615e-17 2.12569e-10 0 6.07991e-16 6.19382e-10 0
-          2.73942e-17 2.12568e-10 0 6.82001e-16 6.19375e-10 0 -4.41587e-17 2.12567e-10 0 -2.39115e-16 6.19374e-10 0
-          -2.98717e-17 2.12567e-10 0 8.54739e-16 6.19385e-10 0 1.9796e-16 2.12574e-10 0 3.22481e-15 6.19357e-10 0
-          1.97285e-17 2.12567e-10 0 6.05024e-16 6.19314e-10 0 -2.72329e-16 2.12574e-10 0 -2.02747e-15 6.19343e-10 0
-          2.22132e-17 2.12578e-10 0 -3.69963e-16 6.19363e-10 0 6.88462e-17 2.12567e-10 0 7.43571e-16 6.1935e-10 0
-          1.80252e-17 2.12575e-10 0 7.66652e-16 6.19345e-10 0 2.01388e-17 2.12568e-10 0 -1.9924e-16 6.19336e-10 0
-          3.91003e-17 2.1257e-10 0 -1.52365e-16 6.1935e-10 0 2.01498e-16 2.12563e-10 0 1.1261e-15 6.19336e-10 0
-          2.40364e-15 1.53399e-09 0 8.40378e-16 1.53398e-09 0 -5.01863e-16 1.53399e-09 0 8.97185e-16 1.53399e-09 0
-          5.60206e-15 1.53398e-09 0 3.85979e-15 1.53392e-09 0 2.24986e-15 1.53394e-09 0 4.01336e-15 1.53387e-09 0
-          1.73122e-15 1.53383e-09 0 5.48803e-16 1.53384e-09 0 1.6671e-15 1.53386e-09 0 1.58242e-15 1.53383e-09 0
-          3.2516e-15 1.53383e-09 0 1.60321e-15 1.53378e-09 0 -4.35128e-16 1.53382e-09 0 2.84304e-15 1.53379e-09 0
-          4.80027e-15 2.63302e-09 0 4.76466e-15 2.63299e-09 0 7.69537e-15 2.63298e-09 0 1.17886e-14 2.63296e-09 0
-          1.69176e-14 2.63291e-09 0 1.60009e-14 2.63276e-09 0 1.17751e-14 2.63273e-09 0 1.14097e-14 2.63267e-09 0
-          9.60396e-15 2.63264e-09 0 8.13929e-15 2.6326e-09 0 8.99627e-15 2.63257e-09 0 9.93589e-15 2.63252e-09 0
-          1.47016e-14 2.63249e-09 0 5.94759e-15 2.63237e-09 0 1.59192e-15 2.63244e-09 0 1.13907e-14 2.63237e-09 0
-          3.815e-14 3.84132e-09 0 3.56013e-14 3.84126e-09 0 8.71446e-14 3.84122e-09 0 1.25951e-13 3.84103e-09 0
-          8.61818e-14 3.84088e-09 0 8.8184e-14 3.84075e-09 0 1.2435e-13 3.84061e-09 0 1.34458e-13 3.84045e-09 0
-          1.0478e-13 3.84029e-09 0 7.39544e-14 3.84021e-09 0 9.78119e-14 3.84013e-09 0 8.8118e-14 3.83998e-09 0
-          7.19154e-14 3.8399e-09 0 6.38236e-14 3.83973e-09 0 3.59744e-14 3.83974e-09 0 3.42877e-14 3.83965e-09 0
-          2.53967e-13 4.26753e-09 0 4.23734e-13 4.26745e-09 0 7.33257e-13 4.26729e-09 0 9.65092e-13 4.267e-09 0
-          1.07742e-12 4.26669e-09 0 1.18602e-12 4.26639e-09 0 1.3166e-12 4.26601e-09 0 1.33347e-12 4.26558e-09 0
-          1.25801e-12 4.26518e-09 0 1.16418e-12 4.26483e-09 0 1.05954e-12 4.26452e-09 0 9.41621e-13 4.26422e-09 0
-          7.90982e-13 4.26399e-09 0 5.94497e-13 4.26379e-09 0 3.50087e-13 4.26366e-09 0 2.08533e-13 4.26359e-09 0
-          5.77387e-13 3.58722e-09 0 1.06497e-12 3.58716e-09 0 1.85377e-12 3.58695e-09 0 2.39882e-12 3.58656e-09 0
-          2.80692e-12 3.58595e-09 0 3.06383e-12 3.58518e-09 0 3.1961e-12 3.58433e-09 0 3.2105e-12 3.58345e-09 0
-          3.14403e-12 3.58259e-09 0 3.00837e-12 3.58175e-09 0 2.7692e-12 3.58098e-09 0 2.4944e-12 3.58035e-09 0
-          2.15017e-12 3.57983e-09 0 1.65019e-12 3.57952e-09 0 9.32022e-13 3.57941e-09 0 4.99691e-13 3.5794e-09 0
-          1.60201e-12 2.4345e-09 0 2.52038e-12 2.43511e-09 0 3.71335e-12 2.43541e-09 0 4.10209e-12 2.4347e-09 0
-          4.31822e-12 2.43344e-09 0 4.49801e-12 2.43195e-09 0 4.60705e-12 2.43033e-09 0 4.61966e-12 2.42866e-09 0
-          4.53985e-12 2.42701e-09 0 4.38761e-12 2.42544e-09 0 4.20015e-12 2.42398e-09 0 3.99756e-12 2.42266e-09 0
-          3.84143e-12 2.4216e-09 0 3.60203e-12 2.4211e-09 0 2.52435e-12 2.42163e-09 0 1.63613e-12 2.4224e-09 0
-          7.31522e-12 8.6467e-10 0 8.31891e-12 8.71907e-10 0 7.97511e-12 8.74642e-10 0 6.14299e-12 8.72427e-10 0
-          5.62211e-12 8.69667e-10 0 5.60617e-12 8.6695e-10 0 5.63054e-12 8.6416e-10 0 5.63321e-12 8.61345e-10 0
-          5.59186e-12 8.58535e-10 0 5.51328e-12 8.55878e-10 0 5.43354e-12 8.53377e-10 0 5.44928e-12 8.51015e-10 0
-          6.11423e-12 8.4855e-10 0 8.37706e-12 8.46545e-10 0 9.08454e-12 8.49989e-10 0 8.13729e-12 8.58368e-10 0
-          4.81267e-11 -4.17388e-09 0 4.48756e-11 -4.11534e-09 0 3.19791e-11 -4.1053e-09 0 2.03628e-11 -4.1163e-09 0
-          1.824e-11 -4.1225e-09 0 1.80697e-11 -4.12694e-09 0 1.80434e-11 -4.13128e-09 0 1.8026e-11 -4.13562e-09 0
-          1.79952e-11 -4.13988e-09 0 1.79533e-11 -4.14398e-09 0 1.79298e-11 -4.14785e-09 0 1.81105e-11 -4.1517e-09 0
-          2.06171e-11 -4.15764e-09 0 3.42514e-11 -4.16942e-09 0 4.99186e-11 -4.15767e-09 0 5.42744e-11 -4.09062e-09 0
-          6.24121e-10 -1.77648e-08 0 5.08873e-10 -1.70664e-08 0 2.9261e-10 -1.70619e-08 0 1.79873e-10 -1.71534e-08 0
-          1.67376e-10 -1.71718e-08 0 1.66386e-10 -1.71788e-08 0 1.6607e-10 -1.71852e-08 0 1.65872e-10 -1.71915e-08 0
-          1.65718e-10 -1.71975e-08 0 1.65612e-10 -1.72034e-08 0 1.65647e-10 -1.72091e-08 0 1.6658e-10 -1.72153e-08 0
-          1.81197e-10 -1.7235e-08 0 3.12682e-10 -1.73397e-08 0 5.68408e-10 -1.73374e-08 0 7.06416e-10 -1.65395e-08 0
-          5.33379e-09 -4.67087e-08 0 3.61929e-09 -4.02722e-08 0 1.44466e-09 -4.11842e-08 0 9.68491e-10 -4.16782e-08 0
-          9.49084e-10 -4.17171e-08 0 9.44584e-10 -4.17284e-08 0 9.42691e-10 -4.17375e-08 0 9.41219e-10 -4.17461e-08 0
-          9.40055e-10 -4.17546e-08 0 9.39235e-10 -4.17627e-08 0 9.38825e-10 -4.17707e-08 0 9.418e-10 -4.17815e-08 0
-          9.64079e-10 -4.18255e-08 0 1.51931e-09 -4.23978e-08 0 4.05476e-09 -4.34676e-08 0 6.05422e-09 -3.61141e-08 0
-          2.92915e-08 -6.61935e-08 0 1.69076e-08 -5.46135e-08 0 4.05537e-09 -5.64311e-08 0 3.57769e-09 -5.72547e-08 0
-          3.56192e-09 -5.73042e-08 0 3.55053e-09 -5.73184e-08 0 3.54207e-09 -5.7329e-08 0 3.53507e-09 -5.73388e-08 0
-          3.52924e-09 -5.73487e-08 0 3.52481e-09 -5.73579e-08 0 3.52137e-09 -5.73672e-08 0 3.52223e-09 -5.73813e-08 0
-          3.53505e-09 -5.74384e-08 0 4.10517e-09 -5.83937e-08 0 1.89853e-08 -6.05201e-08 0 3.3306e-08 -4.72901e-08 0
+          -2.07276e-16 2.12134e-10 0 5.89067e-18 2.12139e-10 0 -1.15941e-15 6.18277e-10 0 -5.66668e-18 6.18297e-10 0
+          1.54999e-16 2.12132e-10 0 1.28539e-15 6.18281e-10 0 -7.7848e-17 2.12133e-10 0 3.11727e-16 6.18263e-10 0
+          -5.81158e-17 2.12142e-10 0 -3.89045e-16 6.18275e-10 0 2.5257e-16 2.12135e-10 0 1.06923e-15 6.18274e-10 0
+          4.50829e-17 2.12124e-10 0 6.76003e-16 6.18242e-10 0 -1.92666e-16 2.12133e-10 0 -1.48858e-15 6.18261e-10 0
+          -7.49084e-17 2.1213e-10 0 -4.11729e-16 6.18282e-10 0 1.25018e-16 2.12137e-10 0 1.9922e-15 6.18271e-10 0
+          4.93203e-17 2.12132e-10 0 4.28939e-16 6.18235e-10 0 -4.36255e-17 2.12135e-10 0 -1.39335e-15 6.18257e-10 0
+          -9.56466e-19 2.12124e-10 0 -1.0068e-15 6.18267e-10 0 -5.69176e-17 2.12127e-10 0 -5.57083e-16 6.18278e-10 0
+          4.75746e-17 2.12128e-10 0 1.25815e-15 6.18272e-10 0 9.38725e-17 2.12131e-10 0 2.66545e-15 6.18248e-10 0
+          3.82399e-15 1.53082e-09 0 2.49411e-15 1.5308e-09 0 -9.61742e-17 1.53079e-09 0 -6.17872e-16 1.53079e-09 0
+          1.68722e-15 1.53079e-09 0 -1.70259e-16 1.53078e-09 0 1.9185e-15 1.53078e-09 0 2.40318e-15 1.53074e-09 0
+          8.36323e-16 1.53077e-09 0 6.75163e-17 1.53074e-09 0 -1.06676e-15 1.53073e-09 0 1.91334e-15 1.53074e-09 0
+          4.75669e-16 1.53074e-09 0 -4.31465e-15 1.53076e-09 0 -5.97406e-16 1.53081e-09 0 6.58798e-15 1.53074e-09 0
+          -4.24251e-15 2.62521e-09 0 5.71804e-15 2.62521e-09 0 7.93298e-15 2.62514e-09 0 -2.77664e-15 2.62515e-09 0
+          6.87377e-15 2.62518e-09 0 6.68338e-15 2.62508e-09 0 3.24028e-15 2.62516e-09 0 6.85777e-15 2.62505e-09 0
+          7.00592e-15 2.62506e-09 0 1.2885e-15 2.62498e-09 0 -3.54778e-16 2.62507e-09 0 8.93194e-15 2.62501e-09 0
+          3.26962e-15 2.62496e-09 0 -5.27395e-15 2.62497e-09 0 -1.30328e-15 2.62508e-09 0 6.22841e-15 2.625e-09 0
+          2.79334e-14 3.82598e-09 0 4.71647e-14 3.82598e-09 0 3.24541e-14 3.82583e-09 0 1.02515e-14 3.82584e-09 0
+          4.67626e-14 3.82587e-09 0 7.41374e-14 3.82571e-09 0 5.23828e-14 3.82567e-09 0 2.99341e-14 3.82561e-09 0
+          3.88375e-14 3.82558e-09 0 3.78497e-14 3.8255e-09 0 3.0972e-14 3.82552e-09 0 3.87083e-14 3.82546e-09 0
+          2.49158e-14 3.82539e-09 0 2.22507e-14 3.8254e-09 0 1.23744e-15 3.8254e-09 0 -3.75349e-14 3.82542e-09 0
+          1.2752e-13 4.23673e-09 0 2.0294e-13 4.23669e-09 0 3.40137e-13 4.23662e-09 0 4.50288e-13 4.23656e-09 0
+          5.68448e-13 4.23646e-09 0 6.39321e-13 4.23629e-09 0 6.11959e-13 4.23612e-09 0 5.75682e-13 4.236e-09 0
+          6.05552e-13 4.23587e-09 0 5.84014e-13 4.23571e-09 0 4.98625e-13 4.2356e-09 0 4.6372e-13 4.2355e-09 0
+          3.92779e-13 4.23541e-09 0 2.91069e-13 4.23536e-09 0 1.356e-13 4.23531e-09 0 1.96485e-14 4.23535e-09 0
+          3.56638e-13 3.5173e-09 0 6.57808e-13 3.51734e-09 0 1.19726e-12 3.51737e-09 0 1.54146e-12 3.51723e-09 0
+          1.71206e-12 3.51698e-09 0 1.7982e-12 3.51666e-09 0 1.8258e-12 3.51632e-09 0 1.8276e-12 3.51596e-09 0
+          1.81126e-12 3.5156e-09 0 1.75203e-12 3.51524e-09 0 1.63957e-12 3.51492e-09 0 1.51551e-12 3.51467e-09 0
+          1.3713e-12 3.51449e-09 0 1.10318e-12 3.51442e-09 0 6.37617e-13 3.51448e-09 0 3.54809e-13 3.51452e-09 0
+          1.42214e-12 2.2896e-09 0 2.18598e-12 2.29032e-09 0 3.07207e-12 2.29093e-09 0 3.18743e-12 2.29074e-09 0
+          3.19709e-12 2.29014e-09 0 3.2319e-12 2.28941e-09 0 3.24399e-12 2.28863e-09 0 3.22474e-12 2.28786e-09 0
+          3.1789e-12 2.28709e-09 0 3.11439e-12 2.28637e-09 0 3.04177e-12 2.28568e-09 0 2.9948e-12 2.28506e-09 0
+          3.04207e-12 2.28457e-09 0 3.04789e-12 2.2845e-09 0 2.24104e-12 2.28529e-09 0 1.48175e-12 2.28614e-09 0
+          7.24796e-12 6.09686e-10 0 8.18113e-12 6.17096e-10 0 7.70695e-12 6.20341e-10 0 5.75481e-12 6.18921e-10 0
+          5.13102e-12 6.17198e-10 0 5.04188e-12 6.15686e-10 0 5.02973e-12 6.14208e-10 0 5.01266e-12 6.12736e-10 0
+          4.98062e-12 6.11305e-10 0 4.94279e-12 6.09936e-10 0 4.91457e-12 6.08635e-10 0 5.00347e-12 6.0734e-10 0
+          5.75284e-12 6.05776e-10 0 8.1096e-12 6.04465e-10 0 8.92472e-12 6.0831e-10 0 8.03618e-12 6.16789e-10 0
+          4.79535e-11 -4.47641e-09 0 4.46959e-11 -4.41784e-09 0 3.18306e-11 -4.40712e-09 0 2.02475e-11 -4.41694e-09 0
+          1.80924e-11 -4.4217e-09 0 1.78753e-11 -4.42446e-09 0 1.78245e-11 -4.42699e-09 0 1.77881e-11 -4.42948e-09 0
+          1.77483e-11 -4.43188e-09 0 1.77132e-11 -4.43418e-09 0 1.77032e-11 -4.43637e-09 0 1.79116e-11 -4.43871e-09 0
+          2.0418e-11 -4.44339e-09 0 3.39247e-11 -4.45412e-09 0 4.94514e-11 -4.44191e-09 0 5.37776e-11 -4.37526e-09 0
+          6.22789e-10 -1.79237e-08 0 5.07962e-10 -1.72272e-08 0 2.9265e-10 -1.7222e-08 0 1.80532e-10 -1.73116e-08 0
+          1.6807e-10 -1.73281e-08 0 1.66966e-10 -1.7333e-08 0 1.66546e-10 -1.73371e-08 0 1.66252e-10 -1.73411e-08 0
+          1.66002e-10 -1.73449e-08 0 1.65813e-10 -1.73486e-08 0 1.65786e-10 -1.73521e-08 0 1.66668e-10 -1.73563e-08 0
+          1.8104e-10 -1.73743e-08 0 3.11072e-10 -1.74768e-08 0 5.64286e-10 -1.7474e-08 0 7.01019e-10 -1.66827e-08 0
+          5.32308e-09 -4.64399e-08 0 3.61265e-09 -4.0021e-08 0 1.44419e-09 -4.09307e-08 0 9.70297e-10 -4.14209e-08 0
+          9.50943e-10 -4.14575e-08 0 9.46112e-10 -4.14665e-08 0 9.43848e-10 -4.14732e-08 0 9.4203e-10 -4.14793e-08 0
+          9.40533e-10 -4.14853e-08 0 9.39401e-10 -4.14909e-08 0 9.38717e-10 -4.14965e-08 0 9.41405e-10 -4.15052e-08 0
+          9.63033e-10 -4.1547e-08 0 1.51272e-09 -4.21125e-08 0 4.02693e-09 -4.31733e-08 0 6.01022e-09 -3.58793e-08 0
+          2.92396e-08 -6.56252e-08 0 1.68764e-08 -5.40768e-08 0 4.04785e-09 -5.58904e-08 0 3.57307e-09 -5.67087e-08 0
+          3.55706e-09 -5.67559e-08 0 3.54499e-09 -5.67677e-08 0 3.53577e-09 -5.67759e-08 0 3.52811e-09 -5.67832e-08 0
+          3.52169e-09 -5.67906e-08 0 3.51677e-09 -5.67974e-08 0 3.5129e-09 -5.68042e-08 0 3.51339e-09 -5.68161e-08 0
+          3.52561e-09 -5.68706e-08 0 4.0907e-09 -5.78157e-08 0 1.88623e-08 -5.99247e-08 0 3.3079e-08 -4.68013e-08 0
         </DataArray>
       </PointData>
       <CellData Scalars="process rank">
-- 
GitLab