From 365cdff54367aa4c850a442a5f36eada8d6084b2 Mon Sep 17 00:00:00 2001
From: Martin Beck <martin.beck@iws.uni-stuttgart.de>
Date: Tue, 13 Oct 2015 11:35:03 +0200
Subject: [PATCH] [el2p] Introduced input files, increased simulation time

Tests now use input files instead of command line input.
Simulation time changed to the original setup with 1e6 s of initialisation and 1e6 s of injection. Reference files are adapted accordingly.
Name of the executable for both the parallel and the non-parallel test is now the same.
(reviewed by fetzer)
---
 test/geomechanics/el2p/CMakeLists.txt       |  10 +-
 test/geomechanics/el2p/el2pproblem.hh       |  63 ++-
 test/geomechanics/el2p/test_el2p.cc         | 116 ++---
 test/geomechanics/el2p/test_el2p.input      |  12 +
 test/references/el2p-parallel-reference.vtu | 158 +++----
 test/references/el2p-reference.vtu          | 462 ++++++++++----------
 6 files changed, 412 insertions(+), 409 deletions(-)
 create mode 100644 test/geomechanics/el2p/test_el2p.input

diff --git a/test/geomechanics/el2p/CMakeLists.txt b/test/geomechanics/el2p/CMakeLists.txt
index 644f28b7fb..ac7e639bb5 100644
--- a/test/geomechanics/el2p/CMakeLists.txt
+++ b/test/geomechanics/el2p/CMakeLists.txt
@@ -4,16 +4,16 @@ add_dumux_test(test_el2p test_el2p test_el2p.cc
                python ${CMAKE_SOURCE_DIR}/bin/runtest.py
                  --script fuzzy
                  --files ${CMAKE_SOURCE_DIR}/test/references/el2p-reference.vtu
-                         ${CMAKE_CURRENT_BINARY_DIR}/el2p-00003.vtu
-                 --command "${CMAKE_CURRENT_BINARY_DIR}/test_el2p ${CMAKE_CURRENT_SOURCE_DIR}/grids/test_el2p.dgf 100 200 100")
+                         ${CMAKE_CURRENT_BINARY_DIR}/el2p-00048.vtu
+                 --command "${CMAKE_CURRENT_BINARY_DIR}/test_el2p")
 
 if(MPI_FOUND)
-  add_dumux_test(test_el2p_parallel test_el2p_parallel test_el2p.cc
+  add_dumux_test(test_el2p_parallel test_el2p test_el2p.cc
                  python ${CMAKE_SOURCE_DIR}/bin/runtest.py
                    --script fuzzy
                    --files ${CMAKE_SOURCE_DIR}/test/references/el2p-parallel-reference.vtu
-                           ${CMAKE_CURRENT_BINARY_DIR}/s0004-p0001-el2p-00003.vtu
-                   --command "${MPIEXEC} -np 4 ${CMAKE_CURRENT_BINARY_DIR}/test_el2p_parallel ${CMAKE_CURRENT_SOURCE_DIR}/grids/test_el2p.dgf 100 200 100")
+                           ${CMAKE_CURRENT_BINARY_DIR}/s0004-p0001-el2p-00048.vtu
+                   --command "${MPIEXEC} -np 4 ${CMAKE_CURRENT_BINARY_DIR}/test_el2p")
 endif()
 
 #install sources
diff --git a/test/geomechanics/el2p/el2pproblem.hh b/test/geomechanics/el2p/el2pproblem.hh
index b1b1413caa..097a5bb920 100644
--- a/test/geomechanics/el2p/el2pproblem.hh
+++ b/test/geomechanics/el2p/el2pproblem.hh
@@ -180,6 +180,7 @@ class El2P_TestProblem : public ImplicitPorousMediaProblem<TypeTag>
 
     typedef typename GET_PROP_TYPE(TypeTag, PrimaryVariables) PrimaryVariables;
     typedef typename GET_PROP_TYPE(TypeTag, BoundaryTypes) BoundaryTypes;
+    typedef typename GET_PROP_TYPE(TypeTag, GridCreator) GridCreator;
     typedef typename GET_PROP_TYPE(TypeTag, TimeManager) TimeManager;
 
     typedef typename GridView::template Codim<0>::Entity Element;
@@ -196,6 +197,7 @@ class El2P_TestProblem : public ImplicitPorousMediaProblem<TypeTag>
     typedef Dune::BlockVector<GlobalPosition> InitialStressField;
 
     typedef typename GET_PROP_TYPE(TypeTag, LocalFEMSpace) LocalFEMSpace;
+    typedef typename GET_PROP_TYPE(TypeTag, PTAG(Problem)) Problem;
 
     typedef typename GET_PROP_TYPE(TypeTag, CO2Table) CO2Table;
     typedef Dumux::CO2<Scalar, CO2Table> CO2;
@@ -209,21 +211,19 @@ public:
      * \param tInitEnd End of initialization period
      */
     El2P_TestProblem(TimeManager &timeManager,
-                    const GridView &gridView,
-                    const Scalar tInitEnd)
+                    const GridView &gridView)
         : ParentType(timeManager, gridView),
-        gridView_(gridView),
-        vertexMapper_(gridView)
+        gridView_(gridView)
     {
         std::cout << "El2P_TestProblem: Initializing the fluid system for the el2p model\n";
 
         // initialize the tables of the fluid system
-        FluidSystem::init(/*Tmin=*/273,
-                          /*Tmax=*/400,
-                          /*nT=*/120,
-                          /*pmin=*/1e5,
-                          /*pmax=*/1e8,
-                          /*np=*/200);
+//         FluidSystem::init(/*Tmin=*/273,
+//                           /*Tmax=*/400,
+//                           /*nT=*/120,
+//                           /*pmin=*/1e5,
+//                           /*pmax=*/1e8,
+//                           /*np=*/200);
 
         // resize the pressure field vector with the number of vertices
         pInit_.resize(gridView.size(dim));
@@ -239,10 +239,31 @@ public:
         // (usually the coupling is switched off for the initialization run)
         coupled_ = false;
         // set initial episode length equal to length of initialization period
+        Scalar tInitEnd = GET_RUNTIME_PARAM(TypeTag, Scalar,TimeManager.TInitEnd);
         this->timeManager().startNextEpisode(tInitEnd);
         // transfer the episode index to spatial parameters
         // (during intialization episode hydraulic different parameters might be applied)
         this->spatialParams().setEpisode(this->timeManager().episodeIndex());
+
+        depthBOR_ = GET_RUNTIME_PARAM(TypeTag, Scalar, Injection.DepthBOR);
+        episodeLength_ = GET_RUNTIME_PARAM(TypeTag, Scalar, TimeManager.EpisodeLength);
+
+        dt_ = GET_RUNTIME_PARAM(TypeTag, Scalar, TimeManager.DtInitial);
+    }
+
+    void init()
+    {
+        if (this->timeManager().time() < 1e-8)
+        {
+            // set the initial approximated hydrostatic pressure distribution
+            // based on an averaged brine density
+            // or based on a pressure polynomial
+            this->initializePressure();
+            // output is written
+            this->setOutput(true);
+        }
+
+        ParentType::init();
     }
 
     // note: pInit is < 0 (just due to geomechanics sign convention applied here)
@@ -300,8 +321,7 @@ public:
         VertexIterator vEndIt = gridView_.template end<dim>();
         for(; vIt != vEndIt; ++vIt)
         {
-            int vIdxGlobal = vertexMapper_.index(*vIt);
-            //
+            int vIdxGlobal = this->vertexMapper().index(*vIt);
             pInit_[vIdxGlobal] = -this->model().curSol().base()[vIdxGlobal*2][0];
         }
     }
@@ -675,17 +695,28 @@ public:
     void episodeEnd()
     {
         this->timeManager().startNextEpisode(episodeLength_);
+        // At the end of the initializationRun
+        if (this->timeManager().time() == GET_RUNTIME_PARAM(TypeTag, Scalar,TimeManager.TInitEnd))
+        {
+            this->timeManager().setTimeStepSize(dt_);
+
+            this->setCoupled(true);
+            // pressure field resulting from the initialization period is applied for the initial
+            // and the Dirichlet boundary conditions
+            this->setPressure();
+            // output is written
+            this->setOutput(true);
+        }
     }
 
 private:
     static constexpr Scalar eps_ = 3e-6;
-    static constexpr Scalar depthBOR_ = 2000;
+    Scalar depthBOR_;
     static constexpr Scalar brineDensity_ = 1059;
-    static constexpr Scalar episodeLength_ = 1e5;
+    Scalar episodeLength_;// = GET_RUNTIME_PARAM(TypeTag, Scalar, TimeManager.EpisodeLength);
 
     std::vector<Scalar> pInit_;
     GridView gridView_;
-    VertexMapper vertexMapper_;
     Scalar dt_;
 public:
     bool initializationRun_, coupled_, output_;
@@ -868,7 +899,7 @@ public:
 
 private:
     static constexpr Scalar eps_ = 3e-6;
-    static constexpr Scalar depthBOR_ = 2000;
+    Scalar depthBOR_;
     std::vector<Scalar> pInit_;
     GridView gridView_;
     VertexMapper vertexMapper_;
diff --git a/test/geomechanics/el2p/test_el2p.cc b/test/geomechanics/el2p/test_el2p.cc
index ccfa4e1102..cbf5f31128 100644
--- a/test/geomechanics/el2p/test_el2p.cc
+++ b/test/geomechanics/el2p/test_el2p.cc
@@ -26,90 +26,50 @@
 #if HAVE_DUNE_PDELAB
 
 #include <dune/common/precision.hh>
-#include <dumux/common/start.hh>
+#include <dune/common/version.hh>
 #include "el2pproblem.hh"
+#include <dumux/common/start.hh>
 
-int main(int argc, char** argv) {
-    try {
-        typedef TTAG(El2P_TestProblem) TypeTag;
-        typedef GET_PROP_TYPE(TypeTag, Grid) Grid;
-        typedef GET_PROP_TYPE(TypeTag, Problem) Problem;
-        typedef GET_PROP_TYPE(TypeTag, TimeManager) TimeManager;
-        typedef Dune::GridPtr<Grid> GridPointer;
-
-        // initialize MPI, finalize is done automatically on exit
-        Dune::MPIHelper::instance(argc, argv);
-
-        if (argc < 5) {
-            std::cout<<"usage: "<<argv[0]<<" grid tEnd_Initialization tEnd dt "<<std::endl;
-            return 1;
-        }
-
-        int argPos = 1;
-        double tEnd(0), dt(0), tInitEnd(0);
-
-        // get the grid name
-        const char *dgfFileName = argv[argPos++];
-        // get the end of the initialization period
-        std::istringstream(argv[argPos++]) >> tInitEnd;
-        // get the end of the real simulation period
-        std::istringstream(argv[argPos++]) >> tEnd;
-        // get the initial time step (applied for initialization and for real simulation)
-        std::istringstream(argv[argPos++]) >> dt;
-
-        // load the grid from file
-        GridPointer gridPtr(dgfFileName);
-        (*gridPtr).loadBalance();
-
-        // Instantiate the time manager
-        TimeManager timeManager;
-
-        // instantiate problem
-        Problem problem(timeManager, gridPtr->leafGridView(), tInitEnd);
-
-        // set the initial approximated hydrostatic pressure distribution
-        // based on an averaged brine density
-        // or based on a pressure polynomial
-        problem.initializePressure();
-
-        // start initialization run to initialize the pressure field correctly
-        timeManager.init(problem, 0.0, // initial time
-                dt, // initial time step
-                tInitEnd); // end of initialization period
-        std::cout<<"tInit: "<<tInitEnd<<" tEnd: "<<tEnd<<" dt: "<<dt<<std::endl;
-        timeManager.run();
-
-        // for the real simulation the coupling between mass balances and momentum equation
-        // is turned on
-        problem.setCoupled(true);
-        // pressure field resulting from the initialization period is applied for the initial
-        // and the Dirichlet boundary conditions
-        problem.setPressure();
-        // output is written
-        problem.setOutput(true);
-        // run the real simulation
-        timeManager.init(problem, tInitEnd, // initial time
-                dt, // initial time step
-                tEnd + tInitEnd); // final time
-
-        timeManager.run();
-
-        return 0;
-    } catch (Dune::Exception &e) {
-        std::cerr << "Dune reported error: " << e << std::endl;
-    } catch (...) {
-        std::cerr << "Unknown exception thrown!" << std::endl;
+/*!
+ * \brief Provides an interface for customizing error messages associated with
+ *        reading in parameters.
+ *
+ * \param progName  The name of the program, that was tried to be started.
+ * \param errorMsg  The error message that was issued by the start function.
+ *                  Comprises the thing that went wrong and a general help message.
+ */
+void usage(const char *progName, const std::string &errorMsg)
+{
+    if (errorMsg.size() > 0) {
+        std::string errorMessageOut = "\nUsage: ";
+        errorMessageOut += progName;
+        errorMessageOut += " [options]\n";
+        errorMessageOut += errorMsg;
+        errorMessageOut += "\n\nThe list of mandatory options for this program is:\n"
+                           "\t-TimeManager.TInitEnd          End of the initialization [s] \n"
+                           "\t-TimeManager.TEnd              End of the simulation [s] \n"
+                           "\t-TimeManager.DtInitial         Initial timestep size [s] \n"
+                           "\t-Grid.File                     Name of the file containing the grid \n"
+                           "\t                               definition in DGF format\n"
+                           "\t-Problem.Name                  String for naming of the output files \n"
+                           "\n";
+
+        std::cout << errorMessageOut << std::endl;
     }
 }
 
+////////////////////////
+// the main function
+////////////////////////
+int main(int argc, char** argv)
+{
+    typedef TTAG(El2P_TestProblem) TypeTag;
+    return Dumux::start<TypeTag>(argc, argv, usage);
+}
+
+
 #else // HAVE_DUNE_PDELAB
 
-#include <iostream>
+#warning You need to have dune-pdelab installed and patched to run this test.
 
-int main()
-{
-#warning You need to have dune-pdelab (>= 2.0) installed to run this test.
-    std::cerr << "You need to have dune-pdelab (>= 2.0) installed to run this test.\n";
-    return 77;
-}
 #endif // HAVE_DUNE_PDELAB
diff --git a/test/geomechanics/el2p/test_el2p.input b/test/geomechanics/el2p/test_el2p.input
new file mode 100644
index 0000000000..fa06d1020b
--- /dev/null
+++ b/test/geomechanics/el2p/test_el2p.input
@@ -0,0 +1,12 @@
+[TimeManager]
+TInitEnd = 1e6       # [s]
+TEnd = 2e6        # [s]
+DtInitial = 10       # [s]
+EpisodeLength = 1e6  # [s]
+
+[Grid]
+File = ./grids/test_el2p.dgf
+Refine = 0
+
+[Injection]
+DepthBOR = 2000     # Depth at the bottom of the scenario [m]
diff --git a/test/references/el2p-parallel-reference.vtu b/test/references/el2p-parallel-reference.vtu
index 6db52d6f97..d6908630c5 100644
--- a/test/references/el2p-parallel-reference.vtu
+++ b/test/references/el2p-parallel-reference.vtu
@@ -11,118 +11,118 @@
           313.15 313.15
         </DataArray>
         <DataArray type="Float32" Name="pW" NumberOfComponents="1" format="ascii">
-          2.08787e+07 2.08787e+07 2.08786e+07 2.08789e+07 1.82821e+07 1.82821e+07 1.82818e+07 1.82815e+07 2.08789e+07 2.08807e+07 1.82815e+07 1.82797e+07
-          2.08786e+07 2.08789e+07 1.82818e+07 1.82815e+07 2.08787e+07 2.08787e+07 1.82821e+07 1.82821e+07 1.56848e+07 1.56848e+07 1.56843e+07 1.56827e+07
-          1.56827e+07 1.57146e+07 1.56843e+07 1.56827e+07 1.56848e+07 1.56848e+07 1.30875e+07 1.30875e+07 1.30873e+07 1.30869e+07 1.30869e+07 1.30853e+07
-          1.30873e+07 1.30869e+07 1.30875e+07 1.30875e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07
+          2.09108e+07 2.09108e+07 2.11522e+07 2.12389e+07 1.83157e+07 1.83157e+07 1.87008e+07 1.90482e+07 2.12389e+07 2.13996e+07 1.90482e+07 1.91824e+07
+          2.11522e+07 2.12389e+07 1.87008e+07 1.90482e+07 2.09108e+07 2.09108e+07 1.83157e+07 1.83157e+07 1.57134e+07 1.57134e+07 1.63926e+07 1.65411e+07
+          1.65411e+07 2.45233e+07 1.63926e+07 1.65411e+07 1.57134e+07 1.57134e+07 1.31046e+07 1.31046e+07 1.34523e+07 1.37949e+07 1.37949e+07 1.39894e+07
+          1.34523e+07 1.37949e+07 1.31046e+07 1.31046e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07
           1.04901e+07 1.04901e+07
         </DataArray>
         <DataArray type="Float32" Name="pN" NumberOfComponents="1" format="ascii">
-          2.08962e+07 2.08962e+07 2.08961e+07 2.08965e+07 1.82997e+07 1.82997e+07 1.82994e+07 1.82991e+07 2.08965e+07 2.08982e+07 1.82991e+07 1.82973e+07
-          2.08961e+07 2.08965e+07 1.82994e+07 1.82991e+07 2.08962e+07 2.08962e+07 1.82997e+07 1.82997e+07 1.57024e+07 1.57024e+07 1.57019e+07 1.57003e+07
-          1.57003e+07 1.57322e+07 1.57019e+07 1.57003e+07 1.57024e+07 1.57024e+07 1.31051e+07 1.31051e+07 1.31048e+07 1.31045e+07 1.31045e+07 1.31029e+07
-          1.31048e+07 1.31045e+07 1.31051e+07 1.31051e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07
+          2.09283e+07 2.09283e+07 2.11697e+07 2.12565e+07 1.83333e+07 1.83333e+07 1.87184e+07 1.90658e+07 2.12565e+07 2.14172e+07 1.90658e+07 1.92e+07
+          2.11697e+07 2.12565e+07 1.87184e+07 1.90658e+07 2.09283e+07 2.09283e+07 1.83333e+07 1.83333e+07 1.5731e+07 1.5731e+07 1.64102e+07 1.65587e+07
+          1.65587e+07 2.45438e+07 1.64102e+07 1.65587e+07 1.5731e+07 1.5731e+07 1.31222e+07 1.31222e+07 1.34699e+07 1.38124e+07 1.38124e+07 1.40069e+07
+          1.34699e+07 1.38124e+07 1.31222e+07 1.31222e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07
           1.05077e+07 1.05077e+07
         </DataArray>
         <DataArray type="Float32" Name="pC" NumberOfComponents="1" format="ascii">
+          17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.8
           17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
-          17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
-          17569.9 17570.6 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
+          17569.9 20467.7 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17570.1
           17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
           17569.9 17569.9
         </DataArray>
         <DataArray type="Float32" Name="SW" NumberOfComponents="1" format="ascii">
           1 1 1 1 1 1 1 1 1 1 1 1
-          1 1 1 1 1 1 1 1 1 1 1 1
-          1 0.999985 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 1 1 1 1 0.999999
+          0.999999 0.938099 1 0.999999 1 1 1 1 1 1 1 0.999994
           1 1 1 1 1 1 1 1 1 1 1 1
           1 1
         </DataArray>
         <DataArray type="Float32" Name="SN" NumberOfComponents="1" format="ascii">
-          0 0 0 0 0 0 0 0 0 0 0 0
-          0 0 0 0 0 0 0 0 0 0 5.53686e-20 1.50343e-12
-          1.50343e-12 1.48314e-05 5.53671e-20 1.50343e-12 0 0 0 0 3.13428e-27 1.3241e-19 1.32409e-19 2.19497e-12
-          3.13427e-27 1.3241e-19 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 0 -2.55006e-16 -2.85954e-11 0 0 -2.85954e-11 -1.90771e-06
+          0 0 -2.55006e-16 -2.85954e-11 0 0 0 0 0 0 1.22683e-11 1.15208e-06
+          1.15208e-06 0.0619015 1.22683e-11 1.15208e-06 0 0 0 0 1.64771e-15 9.82589e-11 9.82589e-11 5.80911e-06
+          1.64771e-15 9.82589e-11 0 0 0 0 0 0 0 0 0 0
           0 0
         </DataArray>
         <DataArray type="Float32" Name="rhoW" NumberOfComponents="1" format="ascii">
-          1056.55 1056.55 1056.55 1056.55 1059.66 1059.66 1059.66 1059.66 1056.55 1056.55 1059.66 1059.66
-          1056.55 1056.55 1059.66 1059.66 1056.55 1056.55 1059.66 1059.66 1062.49 1062.49 1062.49 1062.49
-          1062.49 1062.5 1062.49 1062.49 1062.49 1062.49 1064.99 1064.99 1064.99 1064.99 1064.99 1064.98
-          1064.99 1064.99 1064.99 1064.99 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12
+          1056.56 1056.56 1056.65 1056.68 1059.68 1059.68 1059.82 1059.95 1056.68 1056.74 1059.95 1060
+          1056.65 1056.68 1059.82 1059.95 1056.56 1056.56 1059.68 1059.68 1062.5 1062.5 1062.76 1062.81
+          1062.81 1065.82 1062.76 1062.81 1062.5 1062.5 1064.99 1064.99 1065.13 1065.26 1065.26 1065.34
+          1065.13 1065.26 1064.99 1064.99 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12
           1067.12 1067.12
         </DataArray>
         <DataArray type="Float32" Name="rhoN" NumberOfComponents="1" format="ascii">
-          695.737 695.737 695.736 695.742 675.355 675.355 675.349 675.342 695.742 695.772 675.342 675.303
-          695.736 695.742 675.349 675.342 695.737 695.737 675.355 675.355 673.471 673.471 673.457 673.415
-          673.415 674.245 673.457 673.415 673.471 673.471 669.582 669.582 669.573 669.562 669.562 669.506
-          669.573 669.562 669.582 669.582 660.265 660.265 660.265 660.265 660.265 660.265 660.265 660.265
+          696.287 696.287 700.31 701.729 676.066 676.066 684.062 690.898 701.729 704.326 690.898 693.449
+          700.31 701.729 684.062 690.898 696.287 696.287 676.066 676.066 674.214 674.214 690.891 694.264
+          694.264 806.405 690.891 694.264 674.214 674.214 670.167 670.167 681.6 691.793 691.793 697.196
+          681.6 691.793 670.167 670.167 660.265 660.265 660.265 660.265 660.265 660.265 660.265 660.265
           660.265 660.265
         </DataArray>
         <DataArray type="Float32" Name="u" NumberOfComponents="3" format="ascii">
-          -1.97296e-06 0 0 -4.07123e-14 0 0 -3.93317e-06 -3.93317e-06 0 -3.40134e-14 -1.34961e-05 0
-          -1.7081e-06 0 4.2673e-06 -8.21726e-14 0 1.20408e-06 -1.54798e-06 -1.54798e-06 4.59428e-06 -2.67261e-13 -1.03498e-05 -1.74234e-06
-          -1.34961e-05 -8.72173e-14 0 1.32216e-13 -1.0052e-13 0 -1.03498e-05 -2.77762e-13 -1.74234e-06 4.97113e-13 -2.31617e-13 -9.89639e-05
-          -3.93317e-06 3.93317e-06 0 7.1976e-14 1.34961e-05 0 -1.54798e-06 1.54798e-06 4.59428e-06 1.08941e-13 1.03498e-05 -1.74234e-06
-          -1.97296e-06 0 0 6.90176e-14 0 0 -1.7081e-06 0 4.2673e-06 9.1342e-14 0 1.20408e-06
-          -5.26212e-06 0 1.06304e-05 1.28886e-14 0 1.1194e-05 -9.87207e-06 -9.87207e-06 1.1028e-05 1.18437e-13 -0.000113413 1.17761e-05
-          -0.000113413 -1.04943e-13 1.17761e-05 -1.8588e-13 -5.74986e-13 1.21934e-05 -9.87207e-06 9.87207e-06 1.1028e-05 2.01963e-13 0.000113413 1.17761e-05
-          -5.26212e-06 0 1.06304e-05 8.77731e-14 0 1.1194e-05 -1.15664e-06 0 1.59192e-05 -1.66614e-14 0 2.10784e-05
-          -5.44738e-07 -5.44738e-07 1.74412e-05 -2.46879e-14 -8.99169e-06 2.73609e-05 -8.99169e-06 -2.45514e-13 2.73609e-05 -1.43423e-13 -9.52702e-14 0.000132764
-          -5.44738e-07 5.44738e-07 1.74412e-05 7.85512e-14 8.99169e-06 2.73609e-05 -1.15664e-06 0 1.59192e-05 4.31436e-14 0 2.10784e-05
-          -1.64153e-06 0 1.66183e-05 -2.7253e-14 0 1.83931e-05 -2.37623e-06 -2.37623e-06 1.93307e-05 -3.61824e-14 -6.26453e-06 2.64385e-05
-          -6.26453e-06 -2.31283e-13 2.64385e-05 -8.6403e-14 -1.66826e-13 2.75843e-05 -2.37623e-06 2.37623e-06 1.93307e-05 3.90133e-14 6.26453e-06 2.64385e-05
-          -1.64153e-06 0 1.66183e-05 5.16667e-14 0 1.83931e-05
+          -0.00277155 0 0 -1.54695e-14 0 0 -0.00520181 -0.00520181 0 -2.50259e-14 -0.00758968 0
+          -0.00357935 0 0.0138348 -1.39857e-14 0 0.0124751 -0.00762265 -0.00762265 0.0113754 -2.09942e-14 -0.0123083 0.007474
+          -0.00758968 1.40066e-14 0 -3.70754e-14 1.43706e-14 0 -0.0123083 2.53144e-14 0.007474 -3.84442e-14 2.54465e-14 -0.0187025
+          -0.00520181 0.00520181 0 -3.41893e-14 0.00758968 0 -0.00762265 0.00762265 0.0113754 -4.08003e-14 0.0123083 0.007474
+          -0.00277155 0 0 -2.39843e-14 0 0 -0.00357935 0 0.0138348 -2.61755e-14 0 0.0124751
+          -0.00482172 0 0.0321837 -1.29802e-14 0 0.0329652 -0.0114343 -0.0114343 0.0331809 -1.79138e-14 -0.0385637 0.0340232
+          -0.0385637 4.47982e-14 0.0340232 -4.25656e-14 3.68844e-14 0.0346988 -0.0114343 0.0114343 0.0331809 -5.16361e-14 0.0385637 0.0340232
+          -0.00482172 0 0.0321837 -2.96162e-14 0 0.0329652 -0.00261688 0 0.0481444 -1.56999e-14 0 0.0520194
+          -0.00588075 -0.00588075 0.0547624 -2.07578e-14 -0.0102909 0.0621726 -0.0102909 4.76463e-14 0.0621726 -4.5631e-14 3.23186e-14 0.0929006
+          -0.00588075 0.00588075 0.0547624 -5.04998e-14 0.0102909 0.0621726 -0.00261688 0 0.0481444 -3.14825e-14 0 0.0520194
+          -0.00332892 0 0.0544532 -2.43962e-14 0 0.0582693 -0.00503534 -0.00503534 0.0621787 -3.04972e-14 -0.00716444 0.0707033
+          -0.00716444 4.75369e-14 0.0707033 -4.7659e-14 3.73042e-14 0.0800043 -0.00503534 0.00503534 0.0621787 -5.3383e-14 0.00716444 0.0707033
+          -0.00332892 0 0.0544532 -4.66631e-14 0 0.0582693
         </DataArray>
       </PointData>
       <CellData Scalars="principal stress 1" Vectors="effective stress changes X">
         <DataArray type="Float32" Name="effective stress changes X" NumberOfComponents="3" format="ascii">
-          -26.0752 50.4121 16.8745 -81.3016 91.8238 245.067 -81.3016 -91.8238 245.067 -26.0752 -50.4121 16.8745
-          42.1948 291.982 49.9157 -1367.57 561.716 534.451 -1367.57 -561.716 534.451 42.1948 -291.982 49.9157
-          47.181 289.966 -74.5584 -1370.44 559.941 -575.594 -1370.44 -559.941 -575.594 47.181 -289.966 -74.5584
-          -6.13386 31.1451 -54.113 23.7489 61.6763 -311.181 23.7489 -61.6763 -311.181 -6.13387 -31.1451 -54.113
+          -148564 33867.6 21224.4 -285388 35367.4 93043.3 -285388 -35367.4 93043.3 -148564 -33867.6 21224.4
+          -211828 106177 21728.3 -796583 159075 146568 -796583 -159075 146568 -211828 -106177 21728.3
+          -195985 103540 -51668.6 -778217 157698 -183706 -778217 -157698 -183706 -195985 -103540 -51668.6
+          -112637 28773.9 -59398 -255731 32696.3 -149839 -255731 -32696.3 -149839 -112637 -28773.9 -59398
         </DataArray>
         <DataArray type="Float32" Name="effective stress changes Y" NumberOfComponents="3" format="ascii">
-          50.4121 166.371 -7.28027 91.8238 -81.3016 245.067 -91.8238 -81.3016 -245.067 -50.4121 166.371 7.28027
-          291.982 810.062 282.568 561.716 -1367.57 534.451 -561.716 -1367.57 -534.451 -291.982 810.062 -282.568
-          289.966 795.468 -306.333 559.941 -1370.44 -575.594 -559.941 -1370.44 575.594 -289.966 795.468 306.333
-          31.1451 113.348 -48.6447 61.6763 23.749 -311.181 -61.6763 23.7489 311.181 -31.1451 113.348 48.6447
+          33867.6 110925 36499.6 35367.4 -285388 93043.3 -35367.4 -285388 -93043.3 -33867.6 110925 -36499.6
+          106177 275107 88680.7 159075 -796583 146568 -159075 -796583 -146568 -106177 275107 -88680.7
+          103540 258632 -131632 157698 -778217 -183706 -157698 -778217 183706 -103540 258632 131632
+          28773.9 113530 -102257 32696.3 -255731 -149839 -32696.3 -255731 149839 -28773.9 113530 102257
         </DataArray>
         <DataArray type="Float32" Name="effective stress changes Z" NumberOfComponents="3" format="ascii">
-          16.8745 -7.28027 -21.8807 245.067 245.067 554.605 245.067 -245.067 554.605 16.8745 7.28027 -21.8807
-          49.9157 282.568 -47.3803 534.451 534.451 -1414.8 534.451 -534.451 -1414.8 49.9157 -282.568 -47.3803
-          -74.5584 -306.333 -54.4964 -575.594 -575.594 -1497.1 -575.594 575.594 -1497.1 -74.5584 306.333 -54.4964
-          -54.113 -48.6447 27.557 -311.181 -311.181 640.311 -311.181 311.181 640.311 -54.113 48.6447 27.557
+          21224.4 36499.6 -278484 93043.3 93043.3 -159881 93043.3 -93043.3 -159881 21224.4 -36499.6 -278484
+          21728.3 88680.7 -510507 146568 146568 -1.08846e+06 146568 -146568 -1.08846e+06 21728.3 -88680.7 -510507
+          -51668.6 -131632 -495945 -183706 -183706 -1.12778e+06 -183706 183706 -1.12778e+06 -51668.6 131632 -495945
+          -59398 -102257 -170856 -149839 -149839 -171782 -149839 149839 -171782 -59398 102257 -170856
         </DataArray>
         <DataArray type="Float32" Name="principal stress 1" NumberOfComponents="1" format="ascii">
-          4.28906e+07 4.28909e+07 4.28909e+07 4.28906e+07 3.71722e+07 3.71731e+07 3.71731e+07 3.71722e+07 3.14534e+07 3.14544e+07 3.14544e+07 3.14534e+07
-          2.57343e+07 2.57345e+07 2.57345e+07 2.57343e+07
+          4.32124e+07 4.32436e+07 4.32436e+07 4.32124e+07 3.7775e+07 3.8108e+07 3.8108e+07 3.7775e+07 3.20298e+07 3.23947e+07 3.23947e+07 3.20298e+07
+          2.59776e+07 2.58892e+07 2.58892e+07 2.59776e+07
         </DataArray>
         <DataArray type="Float32" Name="principal stress 2" NumberOfComponents="1" format="ascii">
-          4.28904e+07 4.28903e+07 4.28903e+07 4.28904e+07 3.71714e+07 3.71731e+07 3.71731e+07 3.71714e+07 3.14527e+07 3.14544e+07 3.14544e+07 3.14527e+07
-          2.57343e+07 2.57339e+07 2.57339e+07 2.57343e+07
+          4.29529e+07 4.31181e+07 4.31181e+07 4.29529e+07 3.7288e+07 3.8108e+07 3.8108e+07 3.7288e+07 3.15752e+07 3.23947e+07 3.23947e+07 3.15752e+07
+          2.57514e+07 2.58052e+07 2.58052e+07 2.57514e+07
         </DataArray>
         <DataArray type="Float32" Name="principal stress 3" NumberOfComponents="1" format="ascii">
-          4.28904e+07 4.28903e+07 4.28903e+07 4.28904e+07 3.71713e+07 3.7173e+07 3.7173e+07 3.71713e+07 3.14526e+07 3.14542e+07 3.14542e+07 3.14526e+07
-          2.57342e+07 2.57339e+07 2.57339e+07 2.57342e+07
+          4.28229e+07 4.31181e+07 4.31181e+07 4.28229e+07 3.69894e+07 3.78161e+07 3.78161e+07 3.69894e+07 3.12752e+07 3.20451e+07 3.20451e+07 3.12752e+07
+          2.56932e+07 2.58052e+07 2.58052e+07 2.56932e+07
         </DataArray>
         <DataArray type="Float32" Name="total stresses X" NumberOfComponents="3" format="ascii">
-          4.28904e+07 50.4121 16.8745 4.28903e+07 91.8238 245.067 4.28903e+07 -91.8238 245.067 4.28904e+07 -50.4121 16.8745
-          3.71714e+07 291.982 49.9157 3.71731e+07 561.716 534.451 3.71731e+07 -561.716 534.451 3.71714e+07 -291.982 49.9157
-          3.14527e+07 289.966 -74.5584 3.14544e+07 559.941 -575.594 3.14544e+07 -559.941 -575.594 3.14527e+07 -289.966 -74.5584
-          2.57342e+07 31.1451 -54.113 2.57339e+07 61.6763 -311.181 2.57339e+07 -61.6763 -311.181 2.57342e+07 -31.1451 -54.113
+          4.29529e+07 33867.6 21224.4 4.31181e+07 35367.4 93043.3 4.31181e+07 -35367.4 93043.3 4.29529e+07 -33867.6 21224.4
+          3.7288e+07 106177 21728.3 3.8108e+07 159075 146568 3.8108e+07 -159075 146568 3.7288e+07 -106177 21728.3
+          3.15752e+07 103540 -51668.6 3.23947e+07 157698 -183706 3.23947e+07 -157698 -183706 3.15752e+07 -103540 -51668.6
+          2.57514e+07 28773.9 -59398 2.58052e+07 32696.3 -149839 2.58052e+07 -32696.3 -149839 2.57514e+07 -28773.9 -59398
         </DataArray>
         <DataArray type="Float32" Name="total stresses Y" NumberOfComponents="3" format="ascii">
-          50.4121 4.28906e+07 -7.28027 91.8238 4.28903e+07 245.067 -91.8238 4.28903e+07 -245.067 -50.4121 4.28906e+07 7.28027
-          291.982 3.71722e+07 282.568 561.716 3.71731e+07 534.451 -561.716 3.71731e+07 -534.451 -291.982 3.71722e+07 -282.568
-          289.966 3.14534e+07 -306.333 559.941 3.14544e+07 -575.594 -559.941 3.14544e+07 575.594 -289.966 3.14534e+07 306.333
-          31.1451 2.57343e+07 -48.6447 61.6763 2.57339e+07 -311.181 -61.6763 2.57339e+07 311.181 -31.1451 2.57343e+07 48.6447
+          33867.6 4.32124e+07 36499.6 35367.4 4.31181e+07 93043.3 -35367.4 4.31181e+07 -93043.3 -33867.6 4.32124e+07 -36499.6
+          106177 3.7775e+07 88680.7 159075 3.8108e+07 146568 -159075 3.8108e+07 -146568 -106177 3.7775e+07 -88680.7
+          103540 3.20298e+07 -131632 157698 3.23947e+07 -183706 -157698 3.23947e+07 183706 -103540 3.20298e+07 131632
+          28773.9 2.59776e+07 -102257 32696.3 2.58052e+07 -149839 -32696.3 2.58052e+07 149839 -28773.9 2.59776e+07 102257
         </DataArray>
         <DataArray type="Float32" Name="total stresses Z" NumberOfComponents="3" format="ascii">
-          16.8745 -7.28027 4.28904e+07 245.067 245.067 4.28909e+07 245.067 -245.067 4.28909e+07 16.8745 7.28027 4.28904e+07
-          49.9157 282.568 3.71713e+07 534.451 534.451 3.7173e+07 534.451 -534.451 3.7173e+07 49.9157 -282.568 3.71713e+07
-          -74.5584 -306.333 3.14526e+07 -575.594 -575.594 3.14542e+07 -575.594 575.594 3.14542e+07 -74.5584 306.333 3.14526e+07
-          -54.113 -48.6447 2.57343e+07 -311.181 -311.181 2.57345e+07 -311.181 311.181 2.57345e+07 -54.113 48.6447 2.57343e+07
+          21224.4 36499.6 4.28229e+07 93043.3 93043.3 4.32436e+07 93043.3 -93043.3 4.32436e+07 21224.4 -36499.6 4.28229e+07
+          21728.3 88680.7 3.69894e+07 146568 146568 3.78161e+07 146568 -146568 3.78161e+07 21728.3 -88680.7 3.69894e+07
+          -51668.6 -131632 3.12752e+07 -183706 -183706 3.20451e+07 -183706 183706 3.20451e+07 -51668.6 131632 3.12752e+07
+          -59398 -102257 2.56932e+07 -149839 -149839 2.58892e+07 -149839 149839 2.58892e+07 -59398 102257 2.56932e+07
         </DataArray>
         <DataArray type="Float32" Name="initial stresses X" NumberOfComponents="3" format="ascii">
           4.28905e+07 0 0 4.28905e+07 0 0 4.28905e+07 0 0 4.28905e+07 0 0
@@ -143,28 +143,28 @@
           0 0 2.57343e+07 0 0 2.57343e+07 0 0 2.57343e+07 0 0 2.57343e+07
         </DataArray>
         <DataArray type="Float32" Name="delta pEff" NumberOfComponents="1" format="ascii">
-          -93.4588 -185.879 -185.879 -93.4588 -451.051 2620.35 2620.35 -451.051 -437.656 2663.34 2663.34 -437.656
-          -102.047 -450.662 -450.662 -102.047
+          210886 512906 512906 210886 328052 1.7328e+06 1.7328e+06 328052 318111 1.71984e+06 1.71984e+06 318111
+          129751 326633 326633 129751
         </DataArray>
         <DataArray type="Float32" Name="effectivePressure" NumberOfComponents="1" format="ascii">
-          1.95803e+07 1.95802e+07 1.95802e+07 1.95803e+07 1.6983e+07 1.69861e+07 1.69861e+07 1.6983e+07 1.43857e+07 1.43888e+07 1.43888e+07 1.43857e+07
-          1.17887e+07 1.17884e+07 1.17884e+07 1.17887e+07
+          1.98241e+07 2.01262e+07 2.01262e+07 1.98241e+07 1.73426e+07 1.87474e+07 1.87474e+07 1.73426e+07 1.47271e+07 1.61289e+07 1.61289e+07 1.47271e+07
+          1.19271e+07 1.2124e+07 1.2124e+07 1.19271e+07
         </DataArray>
         <DataArray type="Float32" Name="Pcr_tensile" NumberOfComponents="1" format="ascii">
-          -2.33101e+07 -2.33101e+07 -2.33101e+07 -2.33101e+07 -2.01883e+07 -2.01869e+07 -2.01869e+07 -2.01883e+07 -1.70668e+07 -1.70654e+07 -1.70654e+07 -1.70668e+07
-          -1.39455e+07 -1.39455e+07 -1.39455e+07 -1.39455e+07
+          -2.29988e+07 -2.29919e+07 -2.29919e+07 -2.29988e+07 -1.96467e+07 -1.90688e+07 -1.90688e+07 -1.96467e+07 -1.65481e+07 -1.59163e+07 -1.59163e+07 -1.65481e+07
+          -1.37661e+07 -1.36813e+07 -1.36813e+07 -1.37661e+07
         </DataArray>
         <DataArray type="Float32" Name="Pcr_shear" NumberOfComponents="1" format="ascii">
-          -2.331e+07 -2.33098e+07 -2.33098e+07 -2.331e+07 -2.01878e+07 -2.01869e+07 -2.01869e+07 -2.01878e+07 -1.70664e+07 -1.70653e+07 -1.70653e+07 -1.70664e+07
-          -1.39455e+07 -1.39452e+07 -1.39452e+07 -1.39455e+07
+          -2.28041e+07 -2.29291e+07 -2.29291e+07 -2.28041e+07 -1.92539e+07 -1.89228e+07 -1.89228e+07 -1.92539e+07 -1.61708e+07 -1.57415e+07 -1.57415e+07 -1.61708e+07
+          -1.36239e+07 -1.36393e+07 -1.36393e+07 -1.36239e+07
         </DataArray>
         <DataArray type="Float32" Name="effective Kxx" NumberOfComponents="1" format="ascii">
-          9.99999e-15 9.99997e-15 9.99997e-15 9.99999e-15 9.99993e-15 1.00004e-14 1.00004e-14 9.99993e-15 9.99993e-15 1.00004e-14 1.00004e-14 9.99993e-15
-          9.99999e-15 9.99994e-15 9.99994e-15 9.99999e-15
+          1.00281e-14 1.00651e-14 1.00651e-14 1.00281e-14 1.00398e-14 1.0241e-14 1.0241e-14 1.00398e-14 1.00386e-14 1.02413e-14 1.02413e-14 1.00386e-14
+          1.00151e-14 1.00609e-14 1.00609e-14 1.00151e-14
         </DataArray>
         <DataArray type="Float32" Name="effective Porosity" NumberOfComponents="1" format="ascii">
-          0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
-          0.2 0.2 0.2 0.2
+          0.200025 0.200058 0.200058 0.200025 0.200036 0.200214 0.200214 0.200036 0.200035 0.200215 0.200215 0.200035
+          0.200014 0.200055 0.200055 0.200014
         </DataArray>
       </CellData>
       <Points>
diff --git a/test/references/el2p-reference.vtu b/test/references/el2p-reference.vtu
index e2146dc974..d0a8aa0a5e 100644
--- a/test/references/el2p-reference.vtu
+++ b/test/references/el2p-reference.vtu
@@ -17,40 +17,40 @@
           313.15 313.15 313.15 313.15 313.15
         </DataArray>
         <DataArray type="Float32" Name="pW" NumberOfComponents="1" format="ascii">
-          2.08787e+07 2.08787e+07 2.08787e+07 2.08786e+07 1.82821e+07 1.82821e+07 1.82821e+07 1.82818e+07 2.08787e+07 2.08789e+07 1.82821e+07 1.82815e+07
-          2.08787e+07 2.08786e+07 1.82821e+07 1.82818e+07 2.08787e+07 2.08787e+07 1.82821e+07 1.82821e+07 2.08787e+07 2.08789e+07 1.82821e+07 1.82815e+07
-          2.08807e+07 1.82797e+07 2.08789e+07 1.82815e+07 2.08787e+07 1.82821e+07 2.08787e+07 2.08786e+07 1.82821e+07 1.82818e+07 2.08789e+07 1.82815e+07
-          2.08786e+07 1.82818e+07 2.08787e+07 1.82821e+07 2.08787e+07 2.08787e+07 1.82821e+07 1.82821e+07 2.08787e+07 1.82821e+07 2.08787e+07 1.82821e+07
-          2.08787e+07 1.82821e+07 1.56848e+07 1.56848e+07 1.56848e+07 1.56843e+07 1.56848e+07 1.56827e+07 1.56848e+07 1.56843e+07 1.56848e+07 1.56848e+07
-          1.56848e+07 1.56827e+07 1.57146e+07 1.56827e+07 1.56848e+07 1.56848e+07 1.56843e+07 1.56827e+07 1.56843e+07 1.56848e+07 1.56848e+07 1.56848e+07
-          1.56848e+07 1.56848e+07 1.56848e+07 1.30875e+07 1.30875e+07 1.30875e+07 1.30873e+07 1.30875e+07 1.30869e+07 1.30875e+07 1.30873e+07 1.30875e+07
-          1.30875e+07 1.30875e+07 1.30869e+07 1.30853e+07 1.30869e+07 1.30875e+07 1.30875e+07 1.30873e+07 1.30869e+07 1.30873e+07 1.30875e+07 1.30875e+07
-          1.30875e+07 1.30875e+07 1.30875e+07 1.30875e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07
+          2.09108e+07 2.09108e+07 2.09108e+07 2.11522e+07 1.83157e+07 1.83157e+07 1.83157e+07 1.87008e+07 2.09108e+07 2.12389e+07 1.83157e+07 1.90482e+07
+          2.09108e+07 2.11522e+07 1.83157e+07 1.87008e+07 2.09108e+07 2.09108e+07 1.83157e+07 1.83157e+07 2.09108e+07 2.12389e+07 1.83157e+07 1.90482e+07
+          2.13996e+07 1.91824e+07 2.12389e+07 1.90482e+07 2.09108e+07 1.83157e+07 2.09108e+07 2.11522e+07 1.83157e+07 1.87008e+07 2.12389e+07 1.90482e+07
+          2.11522e+07 1.87008e+07 2.09108e+07 1.83157e+07 2.09108e+07 2.09108e+07 1.83157e+07 1.83157e+07 2.09108e+07 1.83157e+07 2.09108e+07 1.83157e+07
+          2.09108e+07 1.83157e+07 1.57134e+07 1.57134e+07 1.57134e+07 1.63926e+07 1.57134e+07 1.65411e+07 1.57134e+07 1.63926e+07 1.57134e+07 1.57134e+07
+          1.57134e+07 1.65411e+07 2.45233e+07 1.65411e+07 1.57134e+07 1.57134e+07 1.63926e+07 1.65411e+07 1.63926e+07 1.57134e+07 1.57134e+07 1.57134e+07
+          1.57134e+07 1.57134e+07 1.57134e+07 1.31046e+07 1.31046e+07 1.31046e+07 1.34523e+07 1.31046e+07 1.37949e+07 1.31046e+07 1.34523e+07 1.31046e+07
+          1.31046e+07 1.31046e+07 1.37949e+07 1.39894e+07 1.37949e+07 1.31046e+07 1.31046e+07 1.34523e+07 1.37949e+07 1.34523e+07 1.31046e+07 1.31046e+07
+          1.31046e+07 1.31046e+07 1.31046e+07 1.31046e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07
           1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07
           1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07 1.04901e+07
         </DataArray>
         <DataArray type="Float32" Name="pN" NumberOfComponents="1" format="ascii">
-          2.08962e+07 2.08962e+07 2.08962e+07 2.08961e+07 1.82997e+07 1.82997e+07 1.82997e+07 1.82994e+07 2.08962e+07 2.08965e+07 1.82997e+07 1.82991e+07
-          2.08962e+07 2.08961e+07 1.82997e+07 1.82994e+07 2.08962e+07 2.08962e+07 1.82997e+07 1.82997e+07 2.08962e+07 2.08965e+07 1.82997e+07 1.82991e+07
-          2.08982e+07 1.82973e+07 2.08965e+07 1.82991e+07 2.08962e+07 1.82997e+07 2.08962e+07 2.08961e+07 1.82997e+07 1.82994e+07 2.08965e+07 1.82991e+07
-          2.08961e+07 1.82994e+07 2.08962e+07 1.82997e+07 2.08962e+07 2.08962e+07 1.82997e+07 1.82997e+07 2.08962e+07 1.82997e+07 2.08962e+07 1.82997e+07
-          2.08962e+07 1.82997e+07 1.57024e+07 1.57024e+07 1.57024e+07 1.57019e+07 1.57024e+07 1.57003e+07 1.57024e+07 1.57019e+07 1.57024e+07 1.57024e+07
-          1.57024e+07 1.57003e+07 1.57322e+07 1.57003e+07 1.57024e+07 1.57024e+07 1.57019e+07 1.57003e+07 1.57019e+07 1.57024e+07 1.57024e+07 1.57024e+07
-          1.57024e+07 1.57024e+07 1.57024e+07 1.31051e+07 1.31051e+07 1.31051e+07 1.31048e+07 1.31051e+07 1.31045e+07 1.31051e+07 1.31048e+07 1.31051e+07
-          1.31051e+07 1.31051e+07 1.31045e+07 1.31029e+07 1.31045e+07 1.31051e+07 1.31051e+07 1.31048e+07 1.31045e+07 1.31048e+07 1.31051e+07 1.31051e+07
-          1.31051e+07 1.31051e+07 1.31051e+07 1.31051e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07
+          2.09283e+07 2.09283e+07 2.09283e+07 2.11697e+07 1.83333e+07 1.83333e+07 1.83333e+07 1.87184e+07 2.09283e+07 2.12565e+07 1.83333e+07 1.90658e+07
+          2.09283e+07 2.11697e+07 1.83333e+07 1.87184e+07 2.09283e+07 2.09283e+07 1.83333e+07 1.83333e+07 2.09283e+07 2.12565e+07 1.83333e+07 1.90658e+07
+          2.14172e+07 1.92e+07 2.12565e+07 1.90658e+07 2.09283e+07 1.83333e+07 2.09283e+07 2.11697e+07 1.83333e+07 1.87184e+07 2.12565e+07 1.90658e+07
+          2.11697e+07 1.87184e+07 2.09283e+07 1.83333e+07 2.09283e+07 2.09283e+07 1.83333e+07 1.83333e+07 2.09283e+07 1.83333e+07 2.09283e+07 1.83333e+07
+          2.09283e+07 1.83333e+07 1.5731e+07 1.5731e+07 1.5731e+07 1.64102e+07 1.5731e+07 1.65587e+07 1.5731e+07 1.64102e+07 1.5731e+07 1.5731e+07
+          1.5731e+07 1.65587e+07 2.45438e+07 1.65587e+07 1.5731e+07 1.5731e+07 1.64102e+07 1.65587e+07 1.64102e+07 1.5731e+07 1.5731e+07 1.5731e+07
+          1.5731e+07 1.5731e+07 1.5731e+07 1.31222e+07 1.31222e+07 1.31222e+07 1.34699e+07 1.31222e+07 1.38124e+07 1.31222e+07 1.34699e+07 1.31222e+07
+          1.31222e+07 1.31222e+07 1.38124e+07 1.40069e+07 1.38124e+07 1.31222e+07 1.31222e+07 1.34699e+07 1.38124e+07 1.34699e+07 1.31222e+07 1.31222e+07
+          1.31222e+07 1.31222e+07 1.31222e+07 1.31222e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07
           1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07
           1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07 1.05077e+07
         </DataArray>
         <DataArray type="Float32" Name="pC" NumberOfComponents="1" format="ascii">
           17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
           17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
+          17569.9 17569.8 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
           17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
           17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
+          17569.9 17569.9 20467.7 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
           17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
-          17569.9 17569.9 17570.6 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
-          17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
-          17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
+          17569.9 17569.9 17569.9 17570.1 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
           17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
           17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9 17569.9
           17569.9 17569.9 17569.9 17569.9 17569.9
@@ -60,220 +60,220 @@
           1 1 1 1 1 1 1 1 1 1 1 1
           1 1 1 1 1 1 1 1 1 1 1 1
           1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 1 1 1 1 0.999999 1 1 1 1
+          1 0.999999 0.938099 0.999999 1 1 1 0.999999 1 1 1 1
           1 1 1 1 1 1 1 1 1 1 1 1
-          1 1 0.999985 1 1 1 1 1 1 1 1 1
-          1 1 1 1 1 1 1 1 1 1 1 1
-          1 1 1 1 1 1 1 1 1 1 1 1
+          1 1 1 0.999994 1 1 1 1 1 1 1 1
           1 1 1 1 1 1 1 1 1 1 1 1
           1 1 1 1 1 1 1 1 1 1 1 1
           1 1 1 1 1
         </DataArray>
         <DataArray type="Float32" Name="SN" NumberOfComponents="1" format="ascii">
-          0 0 0 0 0 0 0 0 0 0 0 0
-          0 0 0 0 0 0 0 0 0 0 0 0
-          0 0 0 0 0 0 0 0 0 0 0 0
-          0 0 0 0 0 0 0 0 0 0 0 0
-          0 0 0 0 0 5.53698e-20 0 1.50343e-12 0 5.53697e-20 0 0
-          0 1.50343e-12 1.48314e-05 1.50343e-12 0 0 5.53699e-20 1.50343e-12 5.53704e-20 0 0 0
-          0 0 0 0 0 0 3.13457e-27 0 1.32414e-19 0 3.13452e-27 0
-          0 0 1.32414e-19 2.19497e-12 1.32413e-19 0 0 3.13455e-27 1.32414e-19 3.13456e-27 0 0
+          0 0 0 0 0 0 0 -2.55006e-16 0 0 0 -2.85954e-11
+          0 0 0 -2.55006e-16 0 0 0 0 0 0 0 -2.85954e-11
+          0 -1.90771e-06 0 -2.85954e-11 0 0 0 0 0 -2.55006e-16 0 -2.85954e-11
+          0 -2.55006e-16 0 0 0 0 0 0 0 0 0 0
+          0 0 0 0 0 1.22683e-11 0 1.15208e-06 0 1.22683e-11 0 0
+          0 1.15208e-06 0.0619015 1.15208e-06 0 0 1.22683e-11 1.15208e-06 1.22683e-11 0 0 0
+          0 0 0 0 0 0 1.64771e-15 0 9.82589e-11 0 1.64771e-15 0
+          0 0 9.82589e-11 5.80911e-06 9.82589e-11 0 0 1.64771e-15 9.82589e-11 1.64771e-15 0 0
           0 0 0 0 0 0 0 0 0 0 0 0
           0 0 0 0 0 0 0 0 0 0 0 0
           0 0 0 0 0
         </DataArray>
         <DataArray type="Float32" Name="rhoW" NumberOfComponents="1" format="ascii">
-          1056.55 1056.55 1056.55 1056.55 1059.66 1059.66 1059.66 1059.66 1056.55 1056.55 1059.66 1059.66
-          1056.55 1056.55 1059.66 1059.66 1056.55 1056.55 1059.66 1059.66 1056.55 1056.55 1059.66 1059.66
-          1056.55 1059.66 1056.55 1059.66 1056.55 1059.66 1056.55 1056.55 1059.66 1059.66 1056.55 1059.66
-          1056.55 1059.66 1056.55 1059.66 1056.55 1056.55 1059.66 1059.66 1056.55 1059.66 1056.55 1059.66
-          1056.55 1059.66 1062.49 1062.49 1062.49 1062.49 1062.49 1062.49 1062.49 1062.49 1062.49 1062.49
-          1062.49 1062.49 1062.5 1062.49 1062.49 1062.49 1062.49 1062.49 1062.49 1062.49 1062.49 1062.49
-          1062.49 1062.49 1062.49 1064.99 1064.99 1064.99 1064.99 1064.99 1064.99 1064.99 1064.99 1064.99
-          1064.99 1064.99 1064.99 1064.98 1064.99 1064.99 1064.99 1064.99 1064.99 1064.99 1064.99 1064.99
+          1056.56 1056.56 1056.56 1056.65 1059.68 1059.68 1059.68 1059.82 1056.56 1056.68 1059.68 1059.95
+          1056.56 1056.65 1059.68 1059.82 1056.56 1056.56 1059.68 1059.68 1056.56 1056.68 1059.68 1059.95
+          1056.74 1060 1056.68 1059.95 1056.56 1059.68 1056.56 1056.65 1059.68 1059.82 1056.68 1059.95
+          1056.65 1059.82 1056.56 1059.68 1056.56 1056.56 1059.68 1059.68 1056.56 1059.68 1056.56 1059.68
+          1056.56 1059.68 1062.5 1062.5 1062.5 1062.76 1062.5 1062.81 1062.5 1062.76 1062.5 1062.5
+          1062.5 1062.81 1065.82 1062.81 1062.5 1062.5 1062.76 1062.81 1062.76 1062.5 1062.5 1062.5
+          1062.5 1062.5 1062.5 1064.99 1064.99 1064.99 1065.13 1064.99 1065.26 1064.99 1065.13 1064.99
+          1064.99 1064.99 1065.26 1065.34 1065.26 1064.99 1064.99 1065.13 1065.26 1065.13 1064.99 1064.99
           1064.99 1064.99 1064.99 1064.99 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12
           1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12 1067.12
           1067.12 1067.12 1067.12 1067.12 1067.12
         </DataArray>
         <DataArray type="Float32" Name="rhoN" NumberOfComponents="1" format="ascii">
-          695.737 695.737 695.737 695.736 675.355 675.355 675.355 675.349 695.737 695.742 675.355 675.342
-          695.737 695.736 675.355 675.349 695.737 695.737 675.355 675.355 695.737 695.742 675.355 675.342
-          695.772 675.303 695.742 675.342 695.737 675.355 695.737 695.736 675.355 675.349 695.742 675.342
-          695.736 675.349 695.737 675.355 695.737 695.737 675.355 675.355 695.737 675.355 695.737 675.355
-          695.737 675.355 673.471 673.471 673.471 673.457 673.471 673.415 673.471 673.457 673.471 673.471
-          673.471 673.415 674.245 673.415 673.471 673.471 673.457 673.415 673.457 673.471 673.471 673.471
-          673.471 673.471 673.471 669.582 669.582 669.582 669.573 669.582 669.562 669.582 669.573 669.582
-          669.582 669.582 669.562 669.506 669.562 669.582 669.582 669.573 669.562 669.573 669.582 669.582
-          669.582 669.582 669.582 669.582 660.265 660.265 660.265 660.265 660.265 660.265 660.265 660.265
+          696.287 696.287 696.287 700.31 676.066 676.066 676.066 684.062 696.287 701.729 676.066 690.898
+          696.287 700.31 676.066 684.062 696.287 696.287 676.066 676.066 696.287 701.729 676.066 690.898
+          704.326 693.449 701.729 690.898 696.287 676.066 696.287 700.31 676.066 684.062 701.729 690.898
+          700.31 684.062 696.287 676.066 696.287 696.287 676.066 676.066 696.287 676.066 696.287 676.066
+          696.287 676.066 674.214 674.214 674.214 690.891 674.214 694.264 674.214 690.891 674.214 674.214
+          674.214 694.264 806.405 694.264 674.214 674.214 690.891 694.264 690.891 674.214 674.214 674.214
+          674.214 674.214 674.214 670.167 670.167 670.167 681.6 670.167 691.793 670.167 681.6 670.167
+          670.167 670.167 691.793 697.196 691.793 670.167 670.167 681.6 691.793 681.6 670.167 670.167
+          670.167 670.167 670.167 670.167 660.265 660.265 660.265 660.265 660.265 660.265 660.265 660.265
           660.265 660.265 660.265 660.265 660.265 660.265 660.265 660.265 660.265 660.265 660.265 660.265
           660.265 660.265 660.265 660.265 660.265
         </DataArray>
         <DataArray type="Float32" Name="u" NumberOfComponents="3" format="ascii">
-          0 0 0 -1.97296e-06 0 0 0 -1.97296e-06 0 -3.93317e-06 -3.93317e-06 0
-          0 0 4.50451e-06 -1.7081e-06 0 4.2673e-06 0 -1.7081e-06 4.2673e-06 -1.54798e-06 -1.54798e-06 4.59428e-06
-          -3.78082e-14 0 0 -4.76066e-14 -1.34961e-05 0 -4.02437e-14 0 1.20408e-06 -1.46697e-13 -1.03498e-05 -1.74234e-06
-          1.97296e-06 0 0 3.93317e-06 -3.93317e-06 0 1.7081e-06 0 4.2673e-06 1.54798e-06 -1.54798e-06 4.59428e-06
-          0 0 0 0 -1.97296e-06 0 0 0 4.50451e-06 0 -1.7081e-06 4.2673e-06
-          0 1.57442e-14 0 -1.34961e-05 1.97486e-14 0 0 3.52791e-14 1.20408e-06 -1.03498e-05 3.76843e-14 -1.74234e-06
-          -1.08073e-13 -2.99072e-13 0 -1.09318e-14 4.11799e-15 -9.89639e-05 1.34961e-05 4.64771e-14 0 1.03498e-05 1.75278e-13 -1.74234e-06
-          0 4.83087e-14 0 0 9.33911e-14 1.20408e-06 0 1.97296e-06 0 -3.93317e-06 3.93317e-06 0
-          0 1.7081e-06 4.2673e-06 -1.54798e-06 1.54798e-06 4.59428e-06 -2.08581e-14 1.34961e-05 0 -9.11819e-15 1.03498e-05 -1.74234e-06
-          3.93317e-06 3.93317e-06 0 1.54798e-06 1.54798e-06 4.59428e-06 0 1.97296e-06 0 0 1.7081e-06 4.2673e-06
-          0 0 0 -1.97296e-06 0 0 0 0 4.50451e-06 -1.7081e-06 0 4.2673e-06
-          -5.18215e-15 0 0 1.78245e-14 0 1.20408e-06 1.97296e-06 0 0 1.7081e-06 0 4.2673e-06
-          0 0 0 0 0 4.50451e-06 0 0 1.01828e-05 -5.26212e-06 0 1.06304e-05
-          0 -5.26212e-06 1.06304e-05 -9.87207e-06 -9.87207e-06 1.1028e-05 9.51674e-14 0 1.1194e-05 2.62086e-13 -0.000113413 1.17761e-05
-          5.26212e-06 0 1.06304e-05 9.87207e-06 -9.87207e-06 1.1028e-05 0 0 1.01828e-05 0 -5.26212e-06 1.06304e-05
-          0 9.25283e-14 1.1194e-05 -0.000113413 3.09797e-13 1.17761e-05 -5.42592e-14 -8.91101e-13 1.21934e-05 0.000113413 2.83736e-13 1.17761e-05
-          0 1.12155e-13 1.1194e-05 0 5.26212e-06 1.06304e-05 -9.87207e-06 9.87207e-06 1.1028e-05 2.36023e-13 0.000113413 1.17761e-05
-          9.87207e-06 9.87207e-06 1.1028e-05 0 5.26212e-06 1.06304e-05 0 0 1.01828e-05 -5.26212e-06 0 1.06304e-05
-          1.14802e-13 0 1.1194e-05 5.26212e-06 0 1.06304e-05 0 0 1.01828e-05 0 0 1.44722e-05
-          -1.15664e-06 0 1.59192e-05 0 -1.15664e-06 1.59192e-05 -5.44738e-07 -5.44738e-07 1.74412e-05 1.04919e-13 0 2.10784e-05
-          1.88325e-13 -8.99169e-06 2.73609e-05 1.15664e-06 0 1.59192e-05 5.44738e-07 -5.44738e-07 1.74412e-05 0 0 1.44722e-05
-          0 -1.15664e-06 1.59192e-05 0 3.16928e-14 2.10784e-05 -8.99169e-06 2.09631e-14 2.73609e-05 1.76345e-14 1.71377e-14 0.000132764
-          8.99169e-06 3.56923e-13 2.73609e-05 0 1.74448e-13 2.10784e-05 0 1.15664e-06 1.59192e-05 -5.44738e-07 5.44738e-07 1.74412e-05
-          5.24308e-13 8.99169e-06 2.73609e-05 5.44738e-07 5.44738e-07 1.74412e-05 0 1.15664e-06 1.59192e-05 0 0 1.44722e-05
-          -1.15664e-06 0 1.59192e-05 2.47647e-13 0 2.10784e-05 1.15664e-06 0 1.59192e-05 0 0 1.44722e-05
-          0 0 1.55884e-05 -1.64153e-06 0 1.66183e-05 0 -1.64153e-06 1.66183e-05 -2.37623e-06 -2.37623e-06 1.93307e-05
-          1.57007e-13 0 1.83931e-05 1.89476e-13 -6.26453e-06 2.64385e-05 1.64153e-06 0 1.66183e-05 2.37623e-06 -2.37623e-06 1.93307e-05
-          0 0 1.55884e-05 0 -1.64153e-06 1.66183e-05 0 3.97844e-14 1.83931e-05 -6.26453e-06 4.03098e-14 2.64385e-05
-          1.73529e-13 7.35824e-14 2.75843e-05 6.26453e-06 2.77773e-13 2.64385e-05 0 2.52172e-13 1.83931e-05 0 1.64153e-06 1.66183e-05
-          -2.37623e-06 2.37623e-06 1.93307e-05 4.26958e-13 6.26453e-06 2.64385e-05 2.37623e-06 2.37623e-06 1.93307e-05 0 1.64153e-06 1.66183e-05
-          0 0 1.55884e-05 -1.64153e-06 0 1.66183e-05 3.6937e-13 0 1.83931e-05 1.64153e-06 0 1.66183e-05
-          0 0 1.55884e-05
+          0 0 0 -0.00277155 0 0 0 -0.00277155 0 -0.00520181 -0.00520181 0
+          0 0 0.0145612 -0.00357935 0 0.0138348 0 -0.00357935 0.0138348 -0.00762265 -0.00762265 0.0113754
+          -1.63979e-14 0 0 -2.58893e-14 -0.00758968 0 -1.53258e-14 0 0.0124751 -2.46952e-14 -0.0123083 0.007474
+          0.00277155 0 0 0.00520181 -0.00520181 0 0.00357935 0 0.0138348 0.00762265 -0.00762265 0.0113754
+          0 0 0 0 -0.00277155 0 0 0 0.0145612 0 -0.00357935 0.0138348
+          0 -2.14586e-14 0 -0.00758968 -3.10759e-14 0 0 -2.34212e-14 0.0124751 -0.0123083 -3.68178e-14 0.007474
+          -2.97767e-14 -3.52575e-14 0 -2.89874e-14 -4.48231e-14 -0.0187025 0.00758968 -2.58354e-14 0 0.0123083 -3.11301e-14 0.007474
+          0 -1.82393e-14 0 0 -2.01583e-14 0.0124751 0 0.00277155 0 -0.00520181 0.00520181 0
+          0 0.00357935 0.0138348 -0.00762265 0.00762265 0.0113754 -2.06788e-14 0.00758968 0 -1.89846e-14 0.0123083 0.007474
+          0.00520181 0.00520181 0 0.00762265 0.00762265 0.0113754 0 0.00277155 0 0 0.00357935 0.0138348
+          0 0 0 -0.00277155 0 0 0 0 0.0145612 -0.00357935 0 0.0138348
+          -1.31558e-14 0 0 -1.2093e-14 0 0.0124751 0.00277155 0 0 0.00357935 0 0.0138348
+          0 0 0 0 0 0.0145612 0 0 0.0313194 -0.00482172 0 0.0321837
+          0 -0.00482172 0.0321837 -0.0114343 -0.0114343 0.0331809 -1.18752e-14 0 0.0329652 -1.77758e-14 -0.0385637 0.0340232
+          0.00482172 0 0.0321837 0.0114343 -0.0114343 0.0331809 0 0 0.0313194 0 -0.00482172 0.0321837
+          0 -2.57592e-14 0.0329652 -0.0385637 -4.02443e-14 0.0340232 -2.43873e-14 -6.51967e-14 0.0346988 0.0385637 -3.56469e-14 0.0340232
+          0 -2.29542e-14 0.0329652 0 0.00482172 0.0321837 -0.0114343 0.0114343 0.0331809 -1.31924e-14 0.0385637 0.0340232
+          0.0114343 0.0114343 0.0331809 0 0.00482172 0.0321837 0 0 0.0313194 -0.00482172 0 0.0321837
+          -9.08114e-15 0 0.0329652 0.00482172 0 0.0321837 0 0 0.0313194 0 0 0.0448354
+          -0.00261688 0 0.0481444 0 -0.00261688 0.0481444 -0.00588075 -0.00588075 0.0547624 -1.02362e-14 0 0.0520194
+          -1.48656e-14 -0.0102909 0.0621726 0.00261688 0 0.0481444 0.00588075 -0.00588075 0.0547624 0 0 0.0448354
+          0 -0.00261688 0.0481444 0 -3.12155e-14 0.0520194 -0.0102909 -4.69702e-14 0.0621726 -1.37621e-14 -7.12923e-14 0.0929006
+          0.0102909 -3.94686e-14 0.0621726 0 -2.74933e-14 0.0520194 0 0.00261688 0.0481444 -0.00588075 0.00588075 0.0547624
+          -7.37208e-15 0.0102909 0.0621726 0.00588075 0.00588075 0.0547624 0 0.00261688 0.0481444 0 0 0.0448354
+          -0.00261688 0 0.0481444 -6.50998e-15 0 0.0520194 0.00261688 0 0.0481444 0 0 0.0448354
+          0 0 0.0506485 -0.00332892 0 0.0544532 0 -0.00332892 0.0544532 -0.00503534 -0.00503534 0.0621787
+          -1.4797e-14 0 0.0582693 -1.64628e-14 -0.00716444 0.0707033 0.00332892 0 0.0544532 0.00503534 -0.00503534 0.0621787
+          0 0 0.0506485 0 -0.00332892 0.0544532 0 -4.76022e-14 0.0582693 -0.00716444 -5.69408e-14 0.0707033
+          -1.52038e-14 -6.82184e-14 0.0800043 0.00716444 -5.13696e-14 0.0707033 0 -4.23153e-14 0.0582693 0 0.00332892 0.0544532
+          -0.00503534 0.00503534 0.0621787 -1.08813e-14 0.00716444 0.0707033 0.00503534 0.00503534 0.0621787 0 0.00332892 0.0544532
+          0 0 0.0506485 -0.00332892 0 0.0544532 -9.51683e-15 0 0.0582693 0.00332892 0 0.0544532
+          0 0 0.0506485
         </DataArray>
       </PointData>
       <CellData Scalars="principal stress 1" Vectors="effective stress changes X">
         <DataArray type="Float32" Name="effective stress changes X" NumberOfComponents="3" format="ascii">
-          46.9628 9.00043 -6.84956 -26.0752 50.4121 16.8744 -26.0752 -50.4121 -16.8745 46.9628 -9.00043 6.84956
-          166.371 50.4121 -7.28028 -81.3016 91.8238 245.067 -81.3016 -91.8238 -245.067 166.371 -50.4121 7.28028
-          166.371 -50.4121 -7.28028 -81.3016 -91.8238 245.067 -81.3016 91.8238 -245.067 166.371 50.4121 7.28028
-          46.9628 -9.00043 -6.84956 -26.0752 -50.4121 16.8745 -26.0752 50.4121 -16.8745 46.9628 9.00043 6.84956
-          111.855 22.2491 27.3578 42.1948 291.982 49.9157 42.1948 -291.982 -49.9157 111.855 -22.2491 -27.3578
-          810.062 291.982 282.568 -1367.57 561.716 534.451 -1367.57 -561.716 -534.451 810.062 -291.982 -282.568
-          810.062 -291.982 282.568 -1367.57 -561.716 534.451 -1367.57 561.716 -534.451 810.062 291.982 -282.568
-          111.855 -22.2491 27.3578 42.1948 -291.982 49.9157 42.1948 291.982 -49.9157 111.855 22.2491 -27.3578
-          104.829 19.9902 -43.1174 47.181 289.966 -74.5584 47.181 -289.966 74.5584 104.829 -19.9902 43.1174
-          795.468 289.966 -306.333 -1370.44 559.941 -575.594 -1370.44 -559.941 575.594 795.468 -289.966 306.333
-          795.468 -289.966 -306.333 -1370.44 -559.941 -575.594 -1370.44 559.941 575.594 795.468 289.966 306.333
-          104.829 -19.9902 -43.1174 47.181 -289.966 -74.5584 47.181 289.966 74.5584 104.829 19.9902 43.1174
-          40.3198 0.613978 -10.987 -6.13386 31.1451 -54.113 -6.13386 -31.1451 54.113 40.3198 -0.613978 10.987
-          113.348 31.1451 -48.6447 23.749 61.6763 -311.181 23.749 -61.6763 311.181 113.348 -31.1451 48.6447
-          113.348 -31.1451 -48.6447 23.7489 -61.6763 -311.181 23.749 61.6763 311.181 113.348 31.1451 48.6447
-          40.3198 -0.613978 -10.987 -6.13387 -31.1451 -54.113 -6.13385 31.1451 54.113 40.3198 0.613978 10.987
+          70451.1 32367.8 16035.9 -148564 33867.6 21224.4 -148564 -33867.6 -21224.4 70451.1 -32367.8 -16035.9
+          110925 33867.6 36499.6 -285388 35367.4 93043.3 -285388 -35367.4 -93043.3 110925 -33867.6 -36499.6
+          110925 -33867.6 36499.6 -285388 -35367.4 93043.3 -285388 35367.4 -93043.3 110925 33867.6 -36499.6
+          70451.1 -32367.8 16035.9 -148564 -33867.6 21224.4 -148564 33867.6 -21224.4 70451.1 32367.8 -16035.9
+          103381 53279.6 15945.6 -211828 106177 21728.3 -211828 -106177 -21728.3 103381 -53279.6 -15945.6
+          275107 106177 88680.7 -796583 159075 146568 -796583 -159075 -146568 275107 -106177 -88680.7
+          275107 -106177 88680.7 -796583 -159075 146568 -796583 159075 -146568 275107 106177 -88680.7
+          103381 -53279.6 15945.6 -211828 -106177 21728.3 -211828 106177 -21728.3 103381 53279.6 -15945.6
+          94582.4 49382.4 -48867.3 -195985 103540 -51668.6 -195985 -103540 51668.6 94582.4 -49382.4 48867.3
+          258632 103540 -131632 -778217 157698 -183706 -778217 -157698 183706 258632 -103540 131632
+          258632 -103540 -131632 -778217 -157698 -183706 -778217 157698 183706 258632 103540 131632
+          94582.4 -49382.4 -48867.3 -195985 -103540 -51668.6 -195985 103540 51668.6 94582.4 49382.4 48867.3
+          97437.4 24851.4 -53976.3 -112637 28773.9 -59398 -112637 -28773.9 59398 97437.4 -24851.4 53976.3
+          113530 28773.9 -102257 -255731 32696.3 -149839 -255731 -32696.3 149839 113530 -28773.9 102257
+          113530 -28773.9 -102257 -255731 -32696.3 -149839 -255731 32696.3 149839 113530 28773.9 102257
+          97437.4 -24851.4 -53976.3 -112637 -28773.9 -59398 -112637 28773.9 59398 97437.4 24851.4 53976.3
         </DataArray>
         <DataArray type="Float32" Name="effective stress changes Y" NumberOfComponents="3" format="ascii">
-          9.00043 46.9628 -6.84956 50.4121 166.371 -7.28028 -50.4121 166.371 -7.28028 -9.00043 46.9628 -6.84956
-          50.4121 -26.0752 16.8744 91.8238 -81.3016 245.067 -91.8238 -81.3016 245.067 -50.4121 -26.0752 16.8745
-          -50.4121 -26.0752 -16.8745 -91.8238 -81.3016 -245.067 91.8238 -81.3016 -245.067 50.4121 -26.0752 -16.8745
-          -9.00043 46.9628 6.84956 -50.4121 166.371 7.28028 50.4121 166.371 7.28028 9.00043 46.9628 6.84956
-          22.2491 111.855 27.3578 291.982 810.062 282.568 -291.982 810.062 282.568 -22.2491 111.855 27.3578
-          291.982 42.1948 49.9157 561.716 -1367.57 534.451 -561.716 -1367.57 534.451 -291.982 42.1948 49.9157
-          -291.982 42.1948 -49.9157 -561.716 -1367.57 -534.451 561.716 -1367.57 -534.451 291.982 42.1948 -49.9157
-          -22.2491 111.855 -27.3578 -291.982 810.062 -282.568 291.982 810.062 -282.568 22.2491 111.855 -27.3578
-          19.9902 104.829 -43.1174 289.966 795.468 -306.333 -289.966 795.468 -306.333 -19.9902 104.829 -43.1174
-          289.966 47.181 -74.5584 559.941 -1370.44 -575.594 -559.941 -1370.44 -575.594 -289.966 47.181 -74.5584
-          -289.966 47.181 74.5584 -559.941 -1370.44 575.594 559.941 -1370.44 575.594 289.966 47.181 74.5584
-          -19.9902 104.829 43.1174 -289.966 795.468 306.333 289.966 795.468 306.333 19.9902 104.829 43.1174
-          0.613978 40.3198 -10.987 31.1451 113.348 -48.6447 -31.1451 113.348 -48.6447 -0.613978 40.3198 -10.987
-          31.1451 -6.13386 -54.113 61.6763 23.749 -311.181 -61.6763 23.749 -311.181 -31.1451 -6.13386 -54.113
-          -31.1451 -6.13386 54.113 -61.6763 23.749 311.181 61.6763 23.749 311.181 31.1451 -6.13385 54.113
-          -0.613978 40.3198 10.987 -31.1451 113.348 48.6447 31.1451 113.348 48.6447 0.613978 40.3198 10.987
+          32367.8 70451.1 16035.9 33867.6 110925 36499.6 -33867.6 110925 36499.6 -32367.8 70451.1 16035.9
+          33867.6 -148564 21224.4 35367.4 -285388 93043.3 -35367.4 -285388 93043.3 -33867.6 -148564 21224.4
+          -33867.6 -148564 -21224.4 -35367.4 -285388 -93043.3 35367.4 -285388 -93043.3 33867.6 -148564 -21224.4
+          -32367.8 70451.1 -16035.9 -33867.6 110925 -36499.6 33867.6 110925 -36499.6 32367.8 70451.1 -16035.9
+          53279.6 103381 15945.6 106177 275107 88680.7 -106177 275107 88680.7 -53279.6 103381 15945.6
+          106177 -211828 21728.3 159075 -796583 146568 -159075 -796583 146568 -106177 -211828 21728.3
+          -106177 -211828 -21728.3 -159075 -796583 -146568 159075 -796583 -146568 106177 -211828 -21728.3
+          -53279.6 103381 -15945.6 -106177 275107 -88680.7 106177 275107 -88680.7 53279.6 103381 -15945.6
+          49382.4 94582.4 -48867.3 103540 258632 -131632 -103540 258632 -131632 -49382.4 94582.4 -48867.3
+          103540 -195985 -51668.6 157698 -778217 -183706 -157698 -778217 -183706 -103540 -195985 -51668.6
+          -103540 -195985 51668.6 -157698 -778217 183706 157698 -778217 183706 103540 -195985 51668.6
+          -49382.4 94582.4 48867.3 -103540 258632 131632 103540 258632 131632 49382.4 94582.4 48867.3
+          24851.4 97437.4 -53976.3 28773.9 113530 -102257 -28773.9 113530 -102257 -24851.4 97437.4 -53976.3
+          28773.9 -112637 -59398 32696.3 -255731 -149839 -32696.3 -255731 -149839 -28773.9 -112637 -59398
+          -28773.9 -112637 59398 -32696.3 -255731 149839 32696.3 -255731 149839 28773.9 -112637 59398
+          -24851.4 97437.4 53976.3 -28773.9 113530 102257 28773.9 113530 102257 24851.4 97437.4 53976.3
         </DataArray>
         <DataArray type="Float32" Name="effective stress changes Z" NumberOfComponents="3" format="ascii">
-          -6.84956 -6.84956 -87.0152 16.8744 -7.28028 -21.8807 -16.8745 -7.28028 -21.8807 6.84956 -6.84956 -87.0152
-          -7.28028 16.8744 -21.8807 245.067 245.067 554.605 -245.067 245.067 554.605 7.28028 16.8745 -21.8807
-          -7.28028 -16.8745 -21.8807 245.067 -245.067 554.605 -245.067 -245.067 554.605 7.28028 -16.8745 -21.8807
-          -6.84956 6.84956 -87.0152 16.8745 7.28028 -21.8807 -16.8745 7.28028 -21.8807 6.84956 6.84956 -87.0152
-          27.3578 27.3578 -104.288 49.9157 282.568 -47.3803 -49.9157 282.568 -47.3803 -27.3578 27.3578 -104.288
-          282.568 49.9157 -47.3803 534.451 534.451 -1414.8 -534.451 534.451 -1414.8 -282.568 49.9157 -47.3803
-          282.568 -49.9157 -47.3803 534.451 -534.451 -1414.8 -534.451 -534.451 -1414.8 -282.568 -49.9157 -47.3803
-          27.3578 -27.3578 -104.288 49.9157 -282.568 -47.3803 -49.9157 -282.568 -47.3803 -27.3578 -27.3578 -104.288
-          -43.1174 -43.1174 -85.7491 -74.5584 -306.333 -54.4964 74.5584 -306.333 -54.4964 43.1174 -43.1174 -85.7491
-          -306.333 -74.5584 -54.4964 -575.594 -575.594 -1497.1 575.594 -575.594 -1497.1 306.333 -74.5584 -54.4964
-          -306.333 74.5584 -54.4964 -575.594 575.594 -1497.1 575.594 575.594 -1497.1 306.333 74.5584 -54.4964
-          -43.1174 43.1174 -85.7491 -74.5584 306.333 -54.4964 74.5584 306.333 -54.4964 43.1174 43.1174 -85.7491
-          -10.987 -10.987 -10.2948 -54.113 -48.6447 27.557 54.113 -48.6447 27.557 10.987 -10.987 -10.2948
-          -48.6447 -54.113 27.557 -311.181 -311.181 640.311 311.181 -311.181 640.311 48.6447 -54.113 27.557
-          -48.6447 54.113 27.557 -311.181 311.181 640.311 311.181 311.181 640.311 48.6447 54.113 27.557
-          -10.987 10.987 -10.2948 -54.113 48.6447 27.557 54.113 48.6447 27.557 10.987 10.987 -10.2948
+          16035.9 16035.9 -293456 21224.4 36499.6 -278484 -21224.4 36499.6 -278484 -16035.9 16035.9 -293456
+          36499.6 21224.4 -278484 93043.3 93043.3 -159881 -93043.3 93043.3 -159881 -36499.6 21224.4 -278484
+          36499.6 -21224.4 -278484 93043.3 -93043.3 -159881 -93043.3 -93043.3 -159881 -36499.6 -21224.4 -278484
+          16035.9 -16035.9 -293456 21224.4 -36499.6 -278484 -21224.4 -36499.6 -278484 -16035.9 -16035.9 -293456
+          15945.6 15945.6 -410217 21728.3 88680.7 -510507 -21728.3 88680.7 -510507 -15945.6 15945.6 -410217
+          88680.7 21728.3 -510507 146568 146568 -1.08846e+06 -146568 146568 -1.08846e+06 -88680.7 21728.3 -510507
+          88680.7 -21728.3 -510507 146568 -146568 -1.08846e+06 -146568 -146568 -1.08846e+06 -88680.7 -21728.3 -510507
+          15945.6 -15945.6 -410217 21728.3 -88680.7 -510507 -21728.3 -88680.7 -510507 -15945.6 -15945.6 -410217
+          -48867.3 -48867.3 -364281 -51668.6 -131632 -495945 51668.6 -131632 -495945 48867.3 -48867.3 -364281
+          -131632 -51668.6 -495945 -183706 -183706 -1.12778e+06 183706 -183706 -1.12778e+06 131632 -51668.6 -495945
+          -131632 51668.6 -495945 -183706 183706 -1.12778e+06 183706 183706 -1.12778e+06 131632 51668.6 -495945
+          -48867.3 48867.3 -364281 -51668.6 131632 -495945 51668.6 131632 -495945 48867.3 48867.3 -364281
+          -53976.3 -53976.3 -116107 -59398 -102257 -170856 59398 -102257 -170856 53976.3 -53976.3 -116107
+          -102257 -59398 -170856 -149839 -149839 -171782 149839 -149839 -171782 102257 -59398 -170856
+          -102257 59398 -170856 -149839 149839 -171782 149839 149839 -171782 102257 59398 -170856
+          -53976.3 53976.3 -116107 -59398 102257 -170856 59398 102257 -170856 53976.3 53976.3 -116107
         </DataArray>
         <DataArray type="Float32" Name="principal stress 1" NumberOfComponents="1" format="ascii">
-          4.28906e+07 4.28906e+07 4.28906e+07 4.28906e+07 4.28906e+07 4.28911e+07 4.28911e+07 4.28906e+07 4.28906e+07 4.28911e+07 4.28911e+07 4.28906e+07
-          4.28906e+07 4.28906e+07 4.28906e+07 4.28906e+07 3.71718e+07 3.71723e+07 3.71723e+07 3.71718e+07 3.71723e+07 3.71741e+07 3.71741e+07 3.71723e+07
-          3.71723e+07 3.71741e+07 3.71741e+07 3.71723e+07 3.71718e+07 3.71723e+07 3.71723e+07 3.71718e+07 3.14531e+07 3.14536e+07 3.14536e+07 3.14531e+07
-          3.14536e+07 3.14555e+07 3.14555e+07 3.14536e+07 3.14536e+07 3.14555e+07 3.14555e+07 3.14536e+07 3.14531e+07 3.14536e+07 3.14536e+07 3.14531e+07
-          2.57343e+07 2.57344e+07 2.57344e+07 2.57343e+07 2.57344e+07 2.57348e+07 2.57348e+07 2.57344e+07 2.57344e+07 2.57348e+07 2.57348e+07 2.57344e+07
-          2.57343e+07 2.57344e+07 2.57344e+07 2.57343e+07
+          4.30393e+07 4.32124e+07 4.32124e+07 4.30393e+07 4.32124e+07 4.32436e+07 4.32436e+07 4.32124e+07 4.32124e+07 4.32436e+07 4.32436e+07 4.32124e+07
+          4.30393e+07 4.32124e+07 4.32124e+07 4.30393e+07 3.74082e+07 3.7775e+07 3.7775e+07 3.74082e+07 3.7775e+07 3.8108e+07 3.8108e+07 3.7775e+07
+          3.7775e+07 3.8108e+07 3.8108e+07 3.7775e+07 3.74082e+07 3.7775e+07 3.7775e+07 3.74082e+07 3.1676e+07 3.20298e+07 3.20298e+07 3.1676e+07
+          3.20298e+07 3.23947e+07 3.23947e+07 3.20298e+07 3.20298e+07 3.23947e+07 3.23947e+07 3.20298e+07 3.1676e+07 3.20298e+07 3.20298e+07 3.1676e+07
+          2.58752e+07 2.59776e+07 2.59776e+07 2.58752e+07 2.59776e+07 2.58892e+07 2.58892e+07 2.59776e+07 2.59776e+07 2.58892e+07 2.58892e+07 2.59776e+07
+          2.58752e+07 2.59776e+07 2.59776e+07 2.58752e+07
         </DataArray>
         <DataArray type="Float32" Name="principal stress 2" NumberOfComponents="1" format="ascii">
-          4.28905e+07 4.28904e+07 4.28904e+07 4.28905e+07 4.28904e+07 4.28902e+07 4.28902e+07 4.28904e+07 4.28904e+07 4.28902e+07 4.28902e+07 4.28904e+07
-          4.28905e+07 4.28904e+07 4.28904e+07 4.28905e+07 3.71718e+07 3.71713e+07 3.71713e+07 3.71718e+07 3.71713e+07 3.71725e+07 3.71725e+07 3.71713e+07
-          3.71713e+07 3.71725e+07 3.71725e+07 3.71713e+07 3.71718e+07 3.71713e+07 3.71713e+07 3.71718e+07 3.14531e+07 3.14526e+07 3.14526e+07 3.14531e+07
-          3.14526e+07 3.14538e+07 3.14538e+07 3.14526e+07 3.14526e+07 3.14538e+07 3.14538e+07 3.14526e+07 3.14531e+07 3.14526e+07 3.14526e+07 3.14531e+07
-          2.57343e+07 2.57343e+07 2.57343e+07 2.57343e+07 2.57343e+07 2.57338e+07 2.57338e+07 2.57343e+07 2.57343e+07 2.57338e+07 2.57338e+07 2.57343e+07
-          2.57343e+07 2.57343e+07 2.57343e+07 2.57343e+07
+          4.30393e+07 4.29529e+07 4.29529e+07 4.30393e+07 4.29529e+07 4.31181e+07 4.31181e+07 4.29529e+07 4.29529e+07 4.31181e+07 4.31181e+07 4.29529e+07
+          4.30393e+07 4.29529e+07 4.29529e+07 4.30393e+07 3.74082e+07 3.7288e+07 3.7288e+07 3.74082e+07 3.7288e+07 3.8108e+07 3.8108e+07 3.7288e+07
+          3.7288e+07 3.8108e+07 3.8108e+07 3.7288e+07 3.74082e+07 3.7288e+07 3.7288e+07 3.74082e+07 3.1676e+07 3.15752e+07 3.15752e+07 3.1676e+07
+          3.15752e+07 3.23947e+07 3.23947e+07 3.15752e+07 3.15752e+07 3.23947e+07 3.23947e+07 3.15752e+07 3.1676e+07 3.15752e+07 3.15752e+07 3.1676e+07
+          2.58752e+07 2.57514e+07 2.57514e+07 2.58752e+07 2.57514e+07 2.58052e+07 2.58052e+07 2.57514e+07 2.57514e+07 2.58052e+07 2.58052e+07 2.57514e+07
+          2.58752e+07 2.57514e+07 2.57514e+07 2.58752e+07
         </DataArray>
         <DataArray type="Float32" Name="principal stress 3" NumberOfComponents="1" format="ascii">
-          4.28904e+07 4.28904e+07 4.28904e+07 4.28904e+07 4.28904e+07 4.28902e+07 4.28902e+07 4.28904e+07 4.28904e+07 4.28902e+07 4.28902e+07 4.28904e+07
-          4.28904e+07 4.28904e+07 4.28904e+07 4.28904e+07 3.71716e+07 3.71713e+07 3.71713e+07 3.71716e+07 3.71712e+07 3.71725e+07 3.71725e+07 3.71712e+07
-          3.71712e+07 3.71725e+07 3.71725e+07 3.71712e+07 3.71716e+07 3.71713e+07 3.71713e+07 3.71716e+07 3.14529e+07 3.14525e+07 3.14525e+07 3.14529e+07
-          3.14525e+07 3.14537e+07 3.14537e+07 3.14525e+07 3.14525e+07 3.14537e+07 3.14537e+07 3.14525e+07 3.14529e+07 3.14525e+07 3.14525e+07 3.14529e+07
-          2.57343e+07 2.57342e+07 2.57342e+07 2.57343e+07 2.57342e+07 2.57337e+07 2.57337e+07 2.57342e+07 2.57342e+07 2.57337e+07 2.57337e+07 2.57342e+07
-          2.57343e+07 2.57342e+07 2.57342e+07 2.57343e+07
+          4.26754e+07 4.28229e+07 4.28229e+07 4.26754e+07 4.28229e+07 4.31181e+07 4.31181e+07 4.28229e+07 4.28229e+07 4.31181e+07 4.31181e+07 4.28229e+07
+          4.26754e+07 4.28229e+07 4.28229e+07 4.26754e+07 3.68946e+07 3.69894e+07 3.69894e+07 3.68946e+07 3.69894e+07 3.78161e+07 3.78161e+07 3.69894e+07
+          3.69894e+07 3.78161e+07 3.78161e+07 3.69894e+07 3.68946e+07 3.69894e+07 3.69894e+07 3.68946e+07 3.12172e+07 3.12752e+07 3.12752e+07 3.12172e+07
+          3.12752e+07 3.20451e+07 3.20451e+07 3.12752e+07 3.12752e+07 3.20451e+07 3.20451e+07 3.12752e+07 3.12172e+07 3.12752e+07 3.12752e+07 3.12172e+07
+          2.56617e+07 2.56932e+07 2.56932e+07 2.56617e+07 2.56932e+07 2.58052e+07 2.58052e+07 2.56932e+07 2.56932e+07 2.58052e+07 2.58052e+07 2.56932e+07
+          2.56617e+07 2.56932e+07 2.56932e+07 2.56617e+07
         </DataArray>
         <DataArray type="Float32" Name="total stresses X" NumberOfComponents="3" format="ascii">
-          4.28905e+07 9.00043 -6.84956 4.28904e+07 50.4121 16.8744 4.28904e+07 -50.4121 -16.8745 4.28905e+07 -9.00043 6.84956
-          4.28906e+07 50.4121 -7.28028 4.28903e+07 91.8238 245.067 4.28903e+07 -91.8238 -245.067 4.28906e+07 -50.4121 7.28028
-          4.28906e+07 -50.4121 -7.28028 4.28903e+07 -91.8238 245.067 4.28903e+07 91.8238 -245.067 4.28906e+07 50.4121 7.28028
-          4.28905e+07 -9.00043 -6.84956 4.28904e+07 -50.4121 16.8745 4.28904e+07 50.4121 -16.8745 4.28905e+07 9.00043 6.84956
-          3.71718e+07 22.2491 27.3578 3.71714e+07 291.982 49.9157 3.71714e+07 -291.982 -49.9157 3.71718e+07 -22.2491 -27.3578
-          3.71722e+07 291.982 282.568 3.71731e+07 561.716 534.451 3.71731e+07 -561.716 -534.451 3.71722e+07 -291.982 -282.568
-          3.71722e+07 -291.982 282.568 3.71731e+07 -561.716 534.451 3.71731e+07 561.716 -534.451 3.71722e+07 291.982 -282.568
-          3.71718e+07 -22.2491 27.3578 3.71714e+07 -291.982 49.9157 3.71714e+07 291.982 -49.9157 3.71718e+07 22.2491 -27.3578
-          3.14531e+07 19.9902 -43.1174 3.14527e+07 289.966 -74.5584 3.14527e+07 -289.966 74.5584 3.14531e+07 -19.9902 43.1174
-          3.14534e+07 289.966 -306.333 3.14544e+07 559.941 -575.594 3.14544e+07 -559.941 575.594 3.14534e+07 -289.966 306.333
-          3.14534e+07 -289.966 -306.333 3.14544e+07 -559.941 -575.594 3.14544e+07 559.941 575.594 3.14534e+07 289.966 306.333
-          3.14531e+07 -19.9902 -43.1174 3.14527e+07 -289.966 -74.5584 3.14527e+07 289.966 74.5584 3.14531e+07 19.9902 43.1174
-          2.57343e+07 0.613978 -10.987 2.57342e+07 31.1451 -54.113 2.57342e+07 -31.1451 54.113 2.57343e+07 -0.613978 10.987
-          2.57343e+07 31.1451 -48.6447 2.57339e+07 61.6763 -311.181 2.57339e+07 -61.6763 311.181 2.57343e+07 -31.1451 48.6447
-          2.57343e+07 -31.1451 -48.6447 2.57339e+07 -61.6763 -311.181 2.57339e+07 61.6763 311.181 2.57343e+07 31.1451 48.6447
-          2.57343e+07 -0.613978 -10.987 2.57342e+07 -31.1451 -54.113 2.57342e+07 31.1451 54.113 2.57343e+07 0.613978 10.987
+          4.30393e+07 32367.8 16035.9 4.29529e+07 33867.6 21224.4 4.29529e+07 -33867.6 -21224.4 4.30393e+07 -32367.8 -16035.9
+          4.32124e+07 33867.6 36499.6 4.31181e+07 35367.4 93043.3 4.31181e+07 -35367.4 -93043.3 4.32124e+07 -33867.6 -36499.6
+          4.32124e+07 -33867.6 36499.6 4.31181e+07 -35367.4 93043.3 4.31181e+07 35367.4 -93043.3 4.32124e+07 33867.6 -36499.6
+          4.30393e+07 -32367.8 16035.9 4.29529e+07 -33867.6 21224.4 4.29529e+07 33867.6 -21224.4 4.30393e+07 32367.8 -16035.9
+          3.74082e+07 53279.6 15945.6 3.7288e+07 106177 21728.3 3.7288e+07 -106177 -21728.3 3.74082e+07 -53279.6 -15945.6
+          3.7775e+07 106177 88680.7 3.8108e+07 159075 146568 3.8108e+07 -159075 -146568 3.7775e+07 -106177 -88680.7
+          3.7775e+07 -106177 88680.7 3.8108e+07 -159075 146568 3.8108e+07 159075 -146568 3.7775e+07 106177 -88680.7
+          3.74082e+07 -53279.6 15945.6 3.7288e+07 -106177 21728.3 3.7288e+07 106177 -21728.3 3.74082e+07 53279.6 -15945.6
+          3.1676e+07 49382.4 -48867.3 3.15752e+07 103540 -51668.6 3.15752e+07 -103540 51668.6 3.1676e+07 -49382.4 48867.3
+          3.20298e+07 103540 -131632 3.23947e+07 157698 -183706 3.23947e+07 -157698 183706 3.20298e+07 -103540 131632
+          3.20298e+07 -103540 -131632 3.23947e+07 -157698 -183706 3.23947e+07 157698 183706 3.20298e+07 103540 131632
+          3.1676e+07 -49382.4 -48867.3 3.15752e+07 -103540 -51668.6 3.15752e+07 103540 51668.6 3.1676e+07 49382.4 48867.3
+          2.58752e+07 24851.4 -53976.3 2.57514e+07 28773.9 -59398 2.57514e+07 -28773.9 59398 2.58752e+07 -24851.4 53976.3
+          2.59776e+07 28773.9 -102257 2.58052e+07 32696.3 -149839 2.58052e+07 -32696.3 149839 2.59776e+07 -28773.9 102257
+          2.59776e+07 -28773.9 -102257 2.58052e+07 -32696.3 -149839 2.58052e+07 32696.3 149839 2.59776e+07 28773.9 102257
+          2.58752e+07 -24851.4 -53976.3 2.57514e+07 -28773.9 -59398 2.57514e+07 28773.9 59398 2.58752e+07 24851.4 53976.3
         </DataArray>
         <DataArray type="Float32" Name="total stresses Y" NumberOfComponents="3" format="ascii">
-          9.00043 4.28905e+07 -6.84956 50.4121 4.28906e+07 -7.28028 -50.4121 4.28906e+07 -7.28028 -9.00043 4.28905e+07 -6.84956
-          50.4121 4.28904e+07 16.8744 91.8238 4.28903e+07 245.067 -91.8238 4.28903e+07 245.067 -50.4121 4.28904e+07 16.8745
-          -50.4121 4.28904e+07 -16.8745 -91.8238 4.28903e+07 -245.067 91.8238 4.28903e+07 -245.067 50.4121 4.28904e+07 -16.8745
-          -9.00043 4.28905e+07 6.84956 -50.4121 4.28906e+07 7.28028 50.4121 4.28906e+07 7.28028 9.00043 4.28905e+07 6.84956
-          22.2491 3.71718e+07 27.3578 291.982 3.71722e+07 282.568 -291.982 3.71722e+07 282.568 -22.2491 3.71718e+07 27.3578
-          291.982 3.71714e+07 49.9157 561.716 3.71731e+07 534.451 -561.716 3.71731e+07 534.451 -291.982 3.71714e+07 49.9157
-          -291.982 3.71714e+07 -49.9157 -561.716 3.71731e+07 -534.451 561.716 3.71731e+07 -534.451 291.982 3.71714e+07 -49.9157
-          -22.2491 3.71718e+07 -27.3578 -291.982 3.71722e+07 -282.568 291.982 3.71722e+07 -282.568 22.2491 3.71718e+07 -27.3578
-          19.9902 3.14531e+07 -43.1174 289.966 3.14534e+07 -306.333 -289.966 3.14534e+07 -306.333 -19.9902 3.14531e+07 -43.1174
-          289.966 3.14527e+07 -74.5584 559.941 3.14544e+07 -575.594 -559.941 3.14544e+07 -575.594 -289.966 3.14527e+07 -74.5584
-          -289.966 3.14527e+07 74.5584 -559.941 3.14544e+07 575.594 559.941 3.14544e+07 575.594 289.966 3.14527e+07 74.5584
-          -19.9902 3.14531e+07 43.1174 -289.966 3.14534e+07 306.333 289.966 3.14534e+07 306.333 19.9902 3.14531e+07 43.1174
-          0.613978 2.57343e+07 -10.987 31.1451 2.57343e+07 -48.6447 -31.1451 2.57343e+07 -48.6447 -0.613978 2.57343e+07 -10.987
-          31.1451 2.57342e+07 -54.113 61.6763 2.57339e+07 -311.181 -61.6763 2.57339e+07 -311.181 -31.1451 2.57342e+07 -54.113
-          -31.1451 2.57342e+07 54.113 -61.6763 2.57339e+07 311.181 61.6763 2.57339e+07 311.181 31.1451 2.57342e+07 54.113
-          -0.613978 2.57343e+07 10.987 -31.1451 2.57343e+07 48.6447 31.1451 2.57343e+07 48.6447 0.613978 2.57343e+07 10.987
+          32367.8 4.30393e+07 16035.9 33867.6 4.32124e+07 36499.6 -33867.6 4.32124e+07 36499.6 -32367.8 4.30393e+07 16035.9
+          33867.6 4.29529e+07 21224.4 35367.4 4.31181e+07 93043.3 -35367.4 4.31181e+07 93043.3 -33867.6 4.29529e+07 21224.4
+          -33867.6 4.29529e+07 -21224.4 -35367.4 4.31181e+07 -93043.3 35367.4 4.31181e+07 -93043.3 33867.6 4.29529e+07 -21224.4
+          -32367.8 4.30393e+07 -16035.9 -33867.6 4.32124e+07 -36499.6 33867.6 4.32124e+07 -36499.6 32367.8 4.30393e+07 -16035.9
+          53279.6 3.74082e+07 15945.6 106177 3.7775e+07 88680.7 -106177 3.7775e+07 88680.7 -53279.6 3.74082e+07 15945.6
+          106177 3.7288e+07 21728.3 159075 3.8108e+07 146568 -159075 3.8108e+07 146568 -106177 3.7288e+07 21728.3
+          -106177 3.7288e+07 -21728.3 -159075 3.8108e+07 -146568 159075 3.8108e+07 -146568 106177 3.7288e+07 -21728.3
+          -53279.6 3.74082e+07 -15945.6 -106177 3.7775e+07 -88680.7 106177 3.7775e+07 -88680.7 53279.6 3.74082e+07 -15945.6
+          49382.4 3.1676e+07 -48867.3 103540 3.20298e+07 -131632 -103540 3.20298e+07 -131632 -49382.4 3.1676e+07 -48867.3
+          103540 3.15752e+07 -51668.6 157698 3.23947e+07 -183706 -157698 3.23947e+07 -183706 -103540 3.15752e+07 -51668.6
+          -103540 3.15752e+07 51668.6 -157698 3.23947e+07 183706 157698 3.23947e+07 183706 103540 3.15752e+07 51668.6
+          -49382.4 3.1676e+07 48867.3 -103540 3.20298e+07 131632 103540 3.20298e+07 131632 49382.4 3.1676e+07 48867.3
+          24851.4 2.58752e+07 -53976.3 28773.9 2.59776e+07 -102257 -28773.9 2.59776e+07 -102257 -24851.4 2.58752e+07 -53976.3
+          28773.9 2.57514e+07 -59398 32696.3 2.58052e+07 -149839 -32696.3 2.58052e+07 -149839 -28773.9 2.57514e+07 -59398
+          -28773.9 2.57514e+07 59398 -32696.3 2.58052e+07 149839 32696.3 2.58052e+07 149839 28773.9 2.57514e+07 59398
+          -24851.4 2.58752e+07 53976.3 -28773.9 2.59776e+07 102257 28773.9 2.59776e+07 102257 24851.4 2.58752e+07 53976.3
         </DataArray>
         <DataArray type="Float32" Name="total stresses Z" NumberOfComponents="3" format="ascii">
-          -6.84956 -6.84956 4.28904e+07 16.8744 -7.28028 4.28904e+07 -16.8745 -7.28028 4.28904e+07 6.84956 -6.84956 4.28904e+07
-          -7.28028 16.8744 4.28904e+07 245.067 245.067 4.28909e+07 -245.067 245.067 4.28909e+07 7.28028 16.8745 4.28904e+07
-          -7.28028 -16.8745 4.28904e+07 245.067 -245.067 4.28909e+07 -245.067 -245.067 4.28909e+07 7.28028 -16.8745 4.28904e+07
-          -6.84956 6.84956 4.28904e+07 16.8745 7.28028 4.28904e+07 -16.8745 7.28028 4.28904e+07 6.84956 6.84956 4.28904e+07
-          27.3578 27.3578 3.71716e+07 49.9157 282.568 3.71713e+07 -49.9157 282.568 3.71713e+07 -27.3578 27.3578 3.71716e+07
-          282.568 49.9157 3.71713e+07 534.451 534.451 3.7173e+07 -534.451 534.451 3.7173e+07 -282.568 49.9157 3.71713e+07
-          282.568 -49.9157 3.71713e+07 534.451 -534.451 3.7173e+07 -534.451 -534.451 3.7173e+07 -282.568 -49.9157 3.71713e+07
-          27.3578 -27.3578 3.71716e+07 49.9157 -282.568 3.71713e+07 -49.9157 -282.568 3.71713e+07 -27.3578 -27.3578 3.71716e+07
-          -43.1174 -43.1174 3.14529e+07 -74.5584 -306.333 3.14526e+07 74.5584 -306.333 3.14526e+07 43.1174 -43.1174 3.14529e+07
-          -306.333 -74.5584 3.14526e+07 -575.594 -575.594 3.14542e+07 575.594 -575.594 3.14542e+07 306.333 -74.5584 3.14526e+07
-          -306.333 74.5584 3.14526e+07 -575.594 575.594 3.14542e+07 575.594 575.594 3.14542e+07 306.333 74.5584 3.14526e+07
-          -43.1174 43.1174 3.14529e+07 -74.5584 306.333 3.14526e+07 74.5584 306.333 3.14526e+07 43.1174 43.1174 3.14529e+07
-          -10.987 -10.987 2.57343e+07 -54.113 -48.6447 2.57343e+07 54.113 -48.6447 2.57343e+07 10.987 -10.987 2.57343e+07
-          -48.6447 -54.113 2.57343e+07 -311.181 -311.181 2.57345e+07 311.181 -311.181 2.57345e+07 48.6447 -54.113 2.57343e+07
-          -48.6447 54.113 2.57343e+07 -311.181 311.181 2.57345e+07 311.181 311.181 2.57345e+07 48.6447 54.113 2.57343e+07
-          -10.987 10.987 2.57343e+07 -54.113 48.6447 2.57343e+07 54.113 48.6447 2.57343e+07 10.987 10.987 2.57343e+07
+          16035.9 16035.9 4.26754e+07 21224.4 36499.6 4.28229e+07 -21224.4 36499.6 4.28229e+07 -16035.9 16035.9 4.26754e+07
+          36499.6 21224.4 4.28229e+07 93043.3 93043.3 4.32436e+07 -93043.3 93043.3 4.32436e+07 -36499.6 21224.4 4.28229e+07
+          36499.6 -21224.4 4.28229e+07 93043.3 -93043.3 4.32436e+07 -93043.3 -93043.3 4.32436e+07 -36499.6 -21224.4 4.28229e+07
+          16035.9 -16035.9 4.26754e+07 21224.4 -36499.6 4.28229e+07 -21224.4 -36499.6 4.28229e+07 -16035.9 -16035.9 4.26754e+07
+          15945.6 15945.6 3.68946e+07 21728.3 88680.7 3.69894e+07 -21728.3 88680.7 3.69894e+07 -15945.6 15945.6 3.68946e+07
+          88680.7 21728.3 3.69894e+07 146568 146568 3.78161e+07 -146568 146568 3.78161e+07 -88680.7 21728.3 3.69894e+07
+          88680.7 -21728.3 3.69894e+07 146568 -146568 3.78161e+07 -146568 -146568 3.78161e+07 -88680.7 -21728.3 3.69894e+07
+          15945.6 -15945.6 3.68946e+07 21728.3 -88680.7 3.69894e+07 -21728.3 -88680.7 3.69894e+07 -15945.6 -15945.6 3.68946e+07
+          -48867.3 -48867.3 3.12172e+07 -51668.6 -131632 3.12752e+07 51668.6 -131632 3.12752e+07 48867.3 -48867.3 3.12172e+07
+          -131632 -51668.6 3.12752e+07 -183706 -183706 3.20451e+07 183706 -183706 3.20451e+07 131632 -51668.6 3.12752e+07
+          -131632 51668.6 3.12752e+07 -183706 183706 3.20451e+07 183706 183706 3.20451e+07 131632 51668.6 3.12752e+07
+          -48867.3 48867.3 3.12172e+07 -51668.6 131632 3.12752e+07 51668.6 131632 3.12752e+07 48867.3 48867.3 3.12172e+07
+          -53976.3 -53976.3 2.56617e+07 -59398 -102257 2.56932e+07 59398 -102257 2.56932e+07 53976.3 -53976.3 2.56617e+07
+          -102257 -59398 2.56932e+07 -149839 -149839 2.58892e+07 149839 -149839 2.58892e+07 102257 -59398 2.56932e+07
+          -102257 59398 2.56932e+07 -149839 149839 2.58892e+07 149839 149839 2.58892e+07 102257 59398 2.56932e+07
+          -53976.3 53976.3 2.56617e+07 -59398 102257 2.56932e+07 59398 102257 2.56932e+07 53976.3 53976.3 2.56617e+07
         </DataArray>
         <DataArray type="Float32" Name="initial stresses X" NumberOfComponents="3" format="ascii">
           4.28905e+07 0 0 4.28905e+07 0 0 4.28905e+07 0 0 4.28905e+07 0 0
@@ -330,52 +330,52 @@
           0 0 2.57343e+07 0 0 2.57343e+07 0 0 2.57343e+07 0 0 2.57343e+07
         </DataArray>
         <DataArray type="Float32" Name="delta pEff" NumberOfComponents="1" format="ascii">
-          -50.2657 -93.4588 -93.4588 -50.2657 -93.4588 -185.879 -185.879 -93.4588 -93.4588 -185.879 -185.879 -93.4588
-          -50.2657 -93.4588 -93.4588 -50.2657 -103.762 -451.051 -451.051 -103.762 -451.051 2620.35 2620.35 -451.051
-          -451.051 2620.35 2620.35 -451.051 -103.762 -451.051 -451.051 -103.762 -96.0599 -437.656 -437.656 -96.0599
-          -437.656 2663.34 2663.34 -437.656 -437.656 2663.34 2663.34 -437.656 -96.0599 -437.656 -437.656 -96.06
-          -30.375 -102.047 -102.047 -30.375 -102.047 -450.662 -450.662 -102.047 -102.047 -450.662 -450.662 -102.047
-          -30.375 -102.047 -102.047 -30.375
+          78304.7 210886 210886 78304.7 210886 512906 512906 210886 210886 512906 512906 210886
+          78304.7 210886 210886 78304.7 133030 328052 328052 133030 328052 1.7328e+06 1.7328e+06 328052
+          328052 1.7328e+06 1.7328e+06 328052 133030 328052 328052 133030 128366 318111 318111 128366
+          318111 1.71984e+06 1.71984e+06 318111 318111 1.71984e+06 1.71984e+06 318111 128366 318111 318111 128366
+          43466.4 129751 129751 43466.4 129751 326633 326633 129751 129751 326633 326633 129751
+          43466.4 129751 129751 43466.4
         </DataArray>
         <DataArray type="Float32" Name="effectivePressure" NumberOfComponents="1" format="ascii">
-          1.95803e+07 1.95803e+07 1.95803e+07 1.95803e+07 1.95803e+07 1.95802e+07 1.95802e+07 1.95803e+07 1.95803e+07 1.95802e+07 1.95802e+07 1.95803e+07
-          1.95803e+07 1.95803e+07 1.95803e+07 1.95803e+07 1.69834e+07 1.6983e+07 1.6983e+07 1.69834e+07 1.6983e+07 1.69861e+07 1.69861e+07 1.6983e+07
-          1.6983e+07 1.69861e+07 1.69861e+07 1.6983e+07 1.69834e+07 1.6983e+07 1.6983e+07 1.69834e+07 1.43861e+07 1.43857e+07 1.43857e+07 1.43861e+07
-          1.43857e+07 1.43888e+07 1.43888e+07 1.43857e+07 1.43857e+07 1.43888e+07 1.43888e+07 1.43857e+07 1.43861e+07 1.43857e+07 1.43857e+07 1.43861e+07
-          1.17888e+07 1.17887e+07 1.17887e+07 1.17888e+07 1.17887e+07 1.17884e+07 1.17884e+07 1.17887e+07 1.17887e+07 1.17884e+07 1.17884e+07 1.17887e+07
-          1.17888e+07 1.17887e+07 1.17887e+07 1.17888e+07
+          1.96916e+07 1.98241e+07 1.98241e+07 1.96916e+07 1.98241e+07 2.01262e+07 2.01262e+07 1.98241e+07 1.98241e+07 2.01262e+07 2.01262e+07 1.98241e+07
+          1.96916e+07 1.98241e+07 1.98241e+07 1.96916e+07 1.71476e+07 1.73426e+07 1.73426e+07 1.71476e+07 1.73426e+07 1.87474e+07 1.87474e+07 1.73426e+07
+          1.73426e+07 1.87474e+07 1.87474e+07 1.73426e+07 1.71476e+07 1.73426e+07 1.73426e+07 1.71476e+07 1.45374e+07 1.47271e+07 1.47271e+07 1.45374e+07
+          1.47271e+07 1.61289e+07 1.61289e+07 1.47271e+07 1.47271e+07 1.61289e+07 1.61289e+07 1.47271e+07 1.45374e+07 1.47271e+07 1.47271e+07 1.45374e+07
+          1.18408e+07 1.19271e+07 1.19271e+07 1.18408e+07 1.19271e+07 1.2124e+07 1.2124e+07 1.19271e+07 1.19271e+07 1.2124e+07 1.2124e+07 1.19271e+07
+          1.18408e+07 1.19271e+07 1.19271e+07 1.18408e+07
         </DataArray>
         <DataArray type="Float32" Name="Pcr_tensile" NumberOfComponents="1" format="ascii">
-          -2.33101e+07 -2.33101e+07 -2.33101e+07 -2.33101e+07 -2.33101e+07 -2.331e+07 -2.331e+07 -2.33101e+07 -2.33101e+07 -2.331e+07 -2.331e+07 -2.33101e+07
-          -2.33101e+07 -2.33101e+07 -2.33101e+07 -2.33101e+07 -2.01882e+07 -2.01882e+07 -2.01882e+07 -2.01882e+07 -2.01882e+07 -2.01864e+07 -2.01864e+07 -2.01882e+07
-          -2.01882e+07 -2.01864e+07 -2.01864e+07 -2.01882e+07 -2.01882e+07 -2.01882e+07 -2.01882e+07 -2.01882e+07 -1.70668e+07 -1.70668e+07 -1.70668e+07 -1.70668e+07
-          -1.70667e+07 -1.70648e+07 -1.70648e+07 -1.70667e+07 -1.70667e+07 -1.70648e+07 -1.70648e+07 -1.70667e+07 -1.70668e+07 -1.70668e+07 -1.70668e+07 -1.70668e+07
-          -1.39455e+07 -1.39455e+07 -1.39455e+07 -1.39455e+07 -1.39455e+07 -1.39454e+07 -1.39454e+07 -1.39455e+07 -1.39455e+07 -1.39454e+07 -1.39454e+07 -1.39455e+07
-          -1.39455e+07 -1.39455e+07 -1.39455e+07 -1.39455e+07
+          -2.29838e+07 -2.29988e+07 -2.29988e+07 -2.29838e+07 -2.29988e+07 -2.29919e+07 -2.29919e+07 -2.29988e+07 -2.29988e+07 -2.29919e+07 -2.29919e+07 -2.29988e+07
+          -2.29838e+07 -2.29988e+07 -2.29988e+07 -2.29838e+07 -1.9747e+07 -1.96467e+07 -1.96467e+07 -1.9747e+07 -1.96467e+07 -1.90688e+07 -1.90688e+07 -1.96467e+07
+          -1.96467e+07 -1.90688e+07 -1.90688e+07 -1.96467e+07 -1.9747e+07 -1.96467e+07 -1.96467e+07 -1.9747e+07 -1.66798e+07 -1.65481e+07 -1.65481e+07 -1.66798e+07
+          -1.65481e+07 -1.59163e+07 -1.59163e+07 -1.65481e+07 -1.65481e+07 -1.59163e+07 -1.59163e+07 -1.65481e+07 -1.66798e+07 -1.65481e+07 -1.65481e+07 -1.66798e+07
+          -1.38209e+07 -1.37661e+07 -1.37661e+07 -1.38209e+07 -1.37661e+07 -1.36813e+07 -1.36813e+07 -1.37661e+07 -1.37661e+07 -1.36813e+07 -1.36813e+07 -1.37661e+07
+          -1.38209e+07 -1.37661e+07 -1.37661e+07 -1.38209e+07
         </DataArray>
         <DataArray type="Float32" Name="Pcr_shear" NumberOfComponents="1" format="ascii">
-          -2.331e+07 -2.331e+07 -2.331e+07 -2.331e+07 -2.331e+07 -2.33095e+07 -2.33095e+07 -2.331e+07 -2.331e+07 -2.33095e+07 -2.33095e+07 -2.331e+07
-          -2.331e+07 -2.331e+07 -2.331e+07 -2.331e+07 -2.01881e+07 -2.01877e+07 -2.01877e+07 -2.01881e+07 -2.01876e+07 -2.01856e+07 -2.01856e+07 -2.01876e+07
-          -2.01876e+07 -2.01856e+07 -2.01856e+07 -2.01876e+07 -2.01881e+07 -2.01877e+07 -2.01877e+07 -2.01881e+07 -1.70667e+07 -1.70662e+07 -1.70662e+07 -1.70667e+07
-          -1.70662e+07 -1.7064e+07 -1.7064e+07 -1.70662e+07 -1.70662e+07 -1.7064e+07 -1.7064e+07 -1.70662e+07 -1.70667e+07 -1.70662e+07 -1.70662e+07 -1.70667e+07
-          -1.39455e+07 -1.39454e+07 -1.39454e+07 -1.39455e+07 -1.39454e+07 -1.39448e+07 -1.39448e+07 -1.39454e+07 -1.39454e+07 -1.39448e+07 -1.39448e+07 -1.39454e+07
-          -1.39455e+07 -1.39454e+07 -1.39454e+07 -1.39455e+07
+          -2.28019e+07 -2.28041e+07 -2.28041e+07 -2.28019e+07 -2.28041e+07 -2.29291e+07 -2.29291e+07 -2.28041e+07 -2.28041e+07 -2.29291e+07 -2.29291e+07 -2.28041e+07
+          -2.28019e+07 -2.28041e+07 -2.28041e+07 -2.28019e+07 -1.94902e+07 -1.92539e+07 -1.92539e+07 -1.94902e+07 -1.92539e+07 -1.89228e+07 -1.89228e+07 -1.92539e+07
+          -1.92539e+07 -1.89228e+07 -1.89228e+07 -1.92539e+07 -1.94902e+07 -1.92539e+07 -1.92539e+07 -1.94902e+07 -1.64503e+07 -1.61708e+07 -1.61708e+07 -1.64503e+07
+          -1.61708e+07 -1.57415e+07 -1.57415e+07 -1.61708e+07 -1.61708e+07 -1.57415e+07 -1.57415e+07 -1.61708e+07 -1.64503e+07 -1.61708e+07 -1.61708e+07 -1.64503e+07
+          -1.37141e+07 -1.36239e+07 -1.36239e+07 -1.37141e+07 -1.36239e+07 -1.36393e+07 -1.36393e+07 -1.36239e+07 -1.36239e+07 -1.36393e+07 -1.36393e+07 -1.36239e+07
+          -1.37141e+07 -1.36239e+07 -1.36239e+07 -1.37141e+07
         </DataArray>
         <DataArray type="Float32" Name="effective Kxx" NumberOfComponents="1" format="ascii">
-          1e-14 9.99999e-15 9.99999e-15 1e-14 9.99999e-15 9.99997e-15 9.99997e-15 9.99999e-15 9.99999e-15 9.99997e-15 9.99997e-15 9.99999e-15
-          1e-14 9.99999e-15 9.99999e-15 1e-14 9.99999e-15 9.99993e-15 9.99993e-15 9.99999e-15 9.99993e-15 1.00004e-14 1.00004e-14 9.99993e-15
-          9.99993e-15 1.00004e-14 1.00004e-14 9.99993e-15 9.99999e-15 9.99993e-15 9.99993e-15 9.99999e-15 9.99999e-15 9.99993e-15 9.99993e-15 9.99999e-15
-          9.99993e-15 1.00004e-14 1.00004e-14 9.99993e-15 9.99993e-15 1.00004e-14 1.00004e-14 9.99993e-15 9.99999e-15 9.99993e-15 9.99993e-15 9.99999e-15
-          9.99999e-15 9.99999e-15 9.99999e-15 9.99999e-15 9.99999e-15 9.99994e-15 9.99994e-15 9.99999e-15 9.99999e-15 9.99994e-15 9.99994e-15 9.99999e-15
-          9.99999e-15 9.99999e-15 9.99999e-15 9.99999e-15
+          1.00136e-14 1.00281e-14 1.00281e-14 1.00136e-14 1.00281e-14 1.00651e-14 1.00651e-14 1.00281e-14 1.00281e-14 1.00651e-14 1.00651e-14 1.00281e-14
+          1.00136e-14 1.00281e-14 1.00281e-14 1.00136e-14 1.00181e-14 1.00398e-14 1.00398e-14 1.00181e-14 1.00398e-14 1.0241e-14 1.0241e-14 1.00398e-14
+          1.00398e-14 1.0241e-14 1.0241e-14 1.00398e-14 1.00181e-14 1.00398e-14 1.00398e-14 1.00181e-14 1.00156e-14 1.00386e-14 1.00386e-14 1.00156e-14
+          1.00386e-14 1.02413e-14 1.02413e-14 1.00386e-14 1.00386e-14 1.02413e-14 1.02413e-14 1.00386e-14 1.00156e-14 1.00386e-14 1.00386e-14 1.00156e-14
+          9.99301e-15 1.00151e-14 1.00151e-14 9.99301e-15 1.00151e-14 1.00609e-14 1.00609e-14 1.00151e-14 1.00151e-14 1.00609e-14 1.00609e-14 1.00151e-14
+          9.99301e-15 1.00151e-14 1.00151e-14 9.99301e-15
         </DataArray>
         <DataArray type="Float32" Name="effective Porosity" NumberOfComponents="1" format="ascii">
-          0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
-          0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
-          0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
-          0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
-          0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
-          0.2 0.2 0.2 0.2
+          0.200012 0.200025 0.200025 0.200012 0.200025 0.200058 0.200058 0.200025 0.200025 0.200058 0.200058 0.200025
+          0.200012 0.200025 0.200025 0.200012 0.200016 0.200036 0.200036 0.200016 0.200036 0.200214 0.200214 0.200036
+          0.200036 0.200214 0.200214 0.200036 0.200016 0.200036 0.200036 0.200016 0.200014 0.200035 0.200035 0.200014
+          0.200035 0.200215 0.200215 0.200035 0.200035 0.200215 0.200215 0.200035 0.200014 0.200035 0.200035 0.200014
+          0.199994 0.200014 0.200014 0.199994 0.200014 0.200055 0.200055 0.200014 0.200014 0.200055 0.200055 0.200014
+          0.199994 0.200014 0.200014 0.199994
         </DataArray>
       </CellData>
       <Points>
-- 
GitLab