From 969732f0ce0d8b72637bf69e884e9b1a1402346f Mon Sep 17 00:00:00 2001
From: Katharina Heck <katharina.heck@iws.uni-stuttgart.de>
Date: Thu, 8 Nov 2018 16:07:04 +0100
Subject: [PATCH] [test][multidomain] add 1p2c_1p2c test to compare diffusion
laws. MaxwellStefan and Fick's Law produce the same solution for 2 components
---
.../stokesdarcy/1p2c_1p2c/CMakeLists.txt | 1 +
.../diffusionlawcomparison/CMakeLists.txt | 27 +
.../1p2c_1p2c/diffusionlawcomparison/main.cc | 299 +++++++
.../diffusionlawcomparison/params.input | 42 +
.../diffusionlawcomparison/problem_darcy.hh | 281 +++++++
.../diffusionlawcomparison/problem_stokes.hh | 322 ++++++++
..._stokes1p2c_comparison_darcy-reference.vtu | 658 +++++++++++++++
...stokes1p2c_comparison_stokes-reference.vtu | 760 ++++++++++++++++++
8 files changed, 2390 insertions(+)
create mode 100644 test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/CMakeLists.txt
create mode 100644 test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/main.cc
create mode 100644 test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/params.input
create mode 100644 test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/problem_darcy.hh
create mode 100644 test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/problem_stokes.hh
create mode 100644 test/references/test_md_darcy1p2c_stokes1p2c_comparison_darcy-reference.vtu
create mode 100644 test/references/test_md_darcy1p2c_stokes1p2c_comparison_stokes-reference.vtu
diff --git a/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/CMakeLists.txt b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/CMakeLists.txt
index 44146a2819..05a41dcbc9 100644
--- a/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/CMakeLists.txt
+++ b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/CMakeLists.txt
@@ -1,2 +1,3 @@
add_subdirectory("verticalflow")
add_subdirectory("horizontalflow")
+add_subdirectory("diffusionlawcomparison")
diff --git a/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/CMakeLists.txt b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/CMakeLists.txt
new file mode 100644
index 0000000000..de3e3fa46a
--- /dev/null
+++ b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/CMakeLists.txt
@@ -0,0 +1,27 @@
+add_input_file_links()
+
+dune_add_test(NAME test_md_darcy1p2c_stokes1p2c_maxwellstefan
+ SOURCES main.cc
+ CMAKE_GUARD HAVE_UMFPACK
+ COMPILE_DEFINITIONS DIFFUSIONTYPE=MaxwellStefansLaw<TypeTag>
+ COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+ CMD_ARGS --script fuzzy
+ --files ${CMAKE_SOURCE_DIR}/test/references/test_md_darcy1p2c_stokes1p2c_comparison_stokes-reference.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/test_md_darcy1p2c_stokes1p2c_comparison_stokes-00020.vtu
+ ${CMAKE_SOURCE_DIR}/test/references/test_md_darcy1p2c_stokes1p2c_comparison_darcy-reference.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/test_md_darcy1p2c_stokes1p2c_comparison_darcy-00020.vtu
+
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_md_darcy1p2c_stokes1p2c_maxwellstefan params.input")
+
+dune_add_test(NAME test_md_darcy1p2c_stokes1p2c_fickslaw
+ SOURCES main.cc
+ CMAKE_GUARD HAVE_UMFPACK
+ COMPILE_DEFINITIONS DIFFUSIONTYPE=FicksLaw<TypeTag>
+ COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py
+ CMD_ARGS --script fuzzy
+ --files ${CMAKE_SOURCE_DIR}/test/references/test_md_darcy1p2c_stokes1p2c_comparison_stokes-reference.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/test_md_darcy1p2c_stokes1p2c_comparison_stokes-00020.vtu
+ ${CMAKE_SOURCE_DIR}/test/references/test_md_darcy1p2c_stokes1p2c_comparison_darcy-reference.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/test_md_darcy1p2c_stokes1p2c_comparison_darcy-00020.vtu
+
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_md_darcy1p2c_stokes1p2c_fickslaw params.input")
diff --git a/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/main.cc b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/main.cc
new file mode 100644
index 0000000000..4f7a786e24
--- /dev/null
+++ b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/main.cc
@@ -0,0 +1,299 @@
+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
+// vi: set et ts=4 sw=4 sts=4:
+/*****************************************************************************
+ * See the file COPYING for full copying permissions. *
+ * *
+ * This program is free software: you can redistribute it and/or modify *
+ * it under the terms of the GNU General Public License as published by *
+ * the Free Software Foundation, either version 2 of the License, or *
+ * (at your option) any later version. *
+ * *
+ * This program is distributed in the hope that it will be useful, *
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of *
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
+ * GNU General Public License for more details. *
+ * *
+ * You should have received a copy of the GNU General Public License *
+ * along with this program. If not, see <http://www.gnu.org/licenses/>. *
+ *****************************************************************************/
+/*!
+ * \file
+ *
+ * \brief A test problem for the coupled Stokes/Darcy problem (1p2c) comparing different diffusion laws
+ */
+#include <config.h>
+
+#include <ctime>
+#include <iostream>
+#include <fstream>
+
+#include <dune/common/parallel/mpihelper.hh>
+#include <dune/common/timer.hh>
+#include <dune/istl/io.hh>
+
+#include <dumux/common/properties.hh>
+#include <dumux/common/parameters.hh>
+#include <dumux/common/dumuxmessage.hh>
+#include <dumux/common/geometry/diameter.hh>
+#include <dumux/linear/seqsolverbackend.hh>
+#include <dumux/assembly/fvassembler.hh>
+#include <dumux/assembly/diffmethod.hh>
+#include <dumux/discretization/method.hh>
+#include <dumux/io/vtkoutputmodule.hh>
+#include <dumux/io/staggeredvtkoutputmodule.hh>
+#include <dumux/io/grid/gridmanager.hh>
+
+#include <dumux/multidomain/staggeredtraits.hh>
+#include <dumux/multidomain/fvassembler.hh>
+#include <dumux/multidomain/newtonsolver.hh>
+
+#include <dumux/multidomain/boundary/stokesdarcy/couplingmanager.hh>
+
+#include "problem_darcy.hh"
+#include "problem_stokes.hh"
+
+namespace Dumux {
+namespace Properties {
+
+template<class TypeTag>
+struct CouplingManager<TypeTag, TTag::StokesOnePTwoCTypeTag>
+{
+ using Traits = StaggeredMultiDomainTraits<TypeTag, TypeTag, Properties::TTag::DarcyOnePTwoCTypeTag>;
+ using type = Dumux::StokesDarcyCouplingManager<Traits>;
+};
+
+template<class TypeTag>
+struct CouplingManager<TypeTag, TTag::DarcyOnePTwoCTypeTag>
+{
+ using Traits = StaggeredMultiDomainTraits<Properties::TTag::StokesOnePTwoCTypeTag, Properties::TTag::StokesOnePTwoCTypeTag, TypeTag>;
+ using type = Dumux::StokesDarcyCouplingManager<Traits>;
+};
+
+} // end namespace Properties
+} // end namespace Dumux
+
+int main(int argc, char** argv) try
+{
+ using namespace Dumux;
+
+ // initialize MPI, finalize is done automatically on exit
+ const auto& mpiHelper = Dune::MPIHelper::instance(argc, argv);
+
+ // print dumux start message
+ if (mpiHelper.rank() == 0)
+ DumuxMessage::print(/*firstCall=*/true);
+
+ // parse command line arguments and input file
+ Parameters::init(argc, argv);
+
+ // Define the sub problem type tags
+ using StokesTypeTag = Properties::TTag::StokesOnePTwoCTypeTag;
+ using DarcyTypeTag = Properties::TTag::DarcyOnePTwoCTypeTag;
+
+ // try to create a grid (from the given grid file or the input file)
+ // for both sub-domains
+ using DarcyGridManager = Dumux::GridManager<GetPropType<DarcyTypeTag, Properties::Grid>>;
+ DarcyGridManager darcyGridManager;
+ darcyGridManager.init("Darcy"); // pass parameter group
+
+ using StokesGridManager = Dumux::GridManager<GetPropType<StokesTypeTag, Properties::Grid>>;
+ StokesGridManager stokesGridManager;
+ stokesGridManager.init("Stokes"); // pass parameter group
+
+ // we compute on the leaf grid view
+ const auto& darcyGridView = darcyGridManager.grid().leafGridView();
+ const auto& stokesGridView = stokesGridManager.grid().leafGridView();
+
+ // create the finite volume grid geometry
+ using StokesFVGridGeometry = GetPropType<StokesTypeTag, Properties::FVGridGeometry>;
+ auto stokesFvGridGeometry = std::make_shared<StokesFVGridGeometry>(stokesGridView);
+ stokesFvGridGeometry->update();
+ using DarcyFVGridGeometry = GetPropType<DarcyTypeTag, Properties::FVGridGeometry>;
+ auto darcyFvGridGeometry = std::make_shared<DarcyFVGridGeometry>(darcyGridView);
+ darcyFvGridGeometry->update();
+
+ using Traits = StaggeredMultiDomainTraits<StokesTypeTag, StokesTypeTag, DarcyTypeTag>;
+
+ // the coupling manager
+ using CouplingManager = StokesDarcyCouplingManager<Traits>;
+ auto couplingManager = std::make_shared<CouplingManager>(stokesFvGridGeometry, darcyFvGridGeometry);
+
+ // the indices
+ constexpr auto stokesCellCenterIdx = CouplingManager::stokesCellCenterIdx;
+ constexpr auto stokesFaceIdx = CouplingManager::stokesFaceIdx;
+ constexpr auto darcyIdx = CouplingManager::darcyIdx;
+
+ // the problems (initial and boundary conditions)
+ using StokesProblem = GetPropType<StokesTypeTag, Properties::Problem>;
+ auto stokesProblem = std::make_shared<StokesProblem>(stokesFvGridGeometry, couplingManager);
+ using DarcyProblem = GetPropType<DarcyTypeTag, Properties::Problem>;
+ auto darcyProblem = std::make_shared<DarcyProblem>(darcyFvGridGeometry, couplingManager);
+
+ // initialize the fluidsystem (tabulation)
+ GetPropType<StokesTypeTag, Properties::FluidSystem>::init();
+
+ // get some time loop parameters
+ using Scalar = GetPropType<StokesTypeTag, Properties::Scalar>;
+ const auto tEnd = getParam<Scalar>("TimeLoop.TEnd");
+ const auto maxDt = getParam<Scalar>("TimeLoop.MaxTimeStepSize");
+ auto dt = getParam<Scalar>("TimeLoop.DtInitial");
+
+ // check if we are about to restart a previously interrupted simulation
+ Scalar restartTime = 0;
+ if (Parameters::getTree().hasKey("Restart") || Parameters::getTree().hasKey("TimeLoop.Restart"))
+ restartTime = getParam<Scalar>("TimeLoop.Restart");
+
+ // instantiate time loop
+ auto timeLoop = std::make_shared<CheckPointTimeLoop<Scalar>>(restartTime, dt, tEnd);
+ timeLoop->setMaxTimeStepSize(maxDt);
+
+ // control the injection period
+ const Scalar injectionBegin = getParam<Scalar>("Stokes.Problem.InjectionBegin");
+ const Scalar injectionEnd = getParam<Scalar>("Stokes.Problem.InjectionEnd");
+
+ if(injectionBegin > 0.0)
+ timeLoop->setCheckPoint({injectionBegin, injectionEnd});
+ else
+ timeLoop->setCheckPoint({injectionEnd});
+
+ // the solution vector
+ Traits::SolutionVector sol;
+ sol[stokesCellCenterIdx].resize(stokesFvGridGeometry->numCellCenterDofs());
+ sol[stokesFaceIdx].resize(stokesFvGridGeometry->numFaceDofs());
+ sol[darcyIdx].resize(darcyFvGridGeometry->numDofs());
+
+ const auto& cellCenterSol = sol[stokesCellCenterIdx];
+ const auto& faceSol = sol[stokesFaceIdx];
+
+ // apply initial solution for instationary problems
+ GetPropType<StokesTypeTag, Properties::SolutionVector> stokesSol;
+ std::get<0>(stokesSol) = cellCenterSol;
+ std::get<1>(stokesSol) = faceSol;
+ stokesProblem->applyInitialSolution(stokesSol);
+ sol[stokesCellCenterIdx] = stokesSol[stokesCellCenterIdx];
+ sol[stokesFaceIdx] = stokesSol[stokesFaceIdx];
+
+ darcyProblem->applyInitialSolution(sol[darcyIdx]);
+ auto solDarcyOld = sol[darcyIdx];
+
+ auto solOld = sol;
+
+ couplingManager->init(stokesProblem, darcyProblem, sol);
+
+ // the grid variables
+ using StokesGridVariables = GetPropType<StokesTypeTag, Properties::GridVariables>;
+ auto stokesGridVariables = std::make_shared<StokesGridVariables>(stokesProblem, stokesFvGridGeometry);
+ stokesGridVariables->init(stokesSol);
+ using DarcyGridVariables = GetPropType<DarcyTypeTag, Properties::GridVariables>;
+ auto darcyGridVariables = std::make_shared<DarcyGridVariables>(darcyProblem, darcyFvGridGeometry);
+ darcyGridVariables->init(sol[darcyIdx]);
+
+ // intialize the vtk output module
+ const auto stokesName = getParam<std::string>("Problem.Name") + "_" + stokesProblem->name();
+ const auto darcyName = getParam<std::string>("Problem.Name") + "_" + darcyProblem->name();
+
+ StaggeredVtkOutputModule<StokesGridVariables, GetPropType<StokesTypeTag, Properties::SolutionVector>> stokesVtkWriter(*stokesGridVariables, stokesSol, stokesName);
+ GetPropType<StokesTypeTag, Properties::IOFields>::initOutputModule(stokesVtkWriter);
+ stokesVtkWriter.write(0.0);
+
+ VtkOutputModule<DarcyGridVariables, GetPropType<DarcyTypeTag, Properties::SolutionVector>> darcyVtkWriter(*darcyGridVariables, sol[darcyIdx], darcyName);
+ GetPropType<DarcyTypeTag, Properties::IOFields>::initOutputModule(darcyVtkWriter);
+ darcyVtkWriter.write(0.0);
+
+ // the assembler with time loop for instationary problem
+ using Assembler = MultiDomainFVAssembler<Traits, CouplingManager, DiffMethod::numeric>;
+ auto assembler = std::make_shared<Assembler>(std::make_tuple(stokesProblem, stokesProblem, darcyProblem),
+ std::make_tuple(stokesFvGridGeometry->cellCenterFVGridGeometryPtr(),
+ stokesFvGridGeometry->faceFVGridGeometryPtr(),
+ darcyFvGridGeometry),
+ std::make_tuple(stokesGridVariables->cellCenterGridVariablesPtr(),
+ stokesGridVariables->faceGridVariablesPtr(),
+ darcyGridVariables),
+ couplingManager,
+ timeLoop);
+
+ // the linear solver
+ using LinearSolver = UMFPackBackend;
+ auto linearSolver = std::make_shared<LinearSolver>();
+
+ // the non-linear solver
+ using NewtonSolver = MultiDomainNewtonSolver<Assembler, LinearSolver, CouplingManager>;
+ NewtonSolver nonLinearSolver(assembler, linearSolver, couplingManager);
+
+ constexpr auto eps = 1e-6;
+
+ // time loop
+ timeLoop->start(); do
+ {
+ // set previous solution for storage evaluations
+ assembler->setPreviousSolution(solOld);
+
+
+ if(timeLoop->time() > injectionBegin - eps && timeLoop->time() < injectionEnd + eps)
+ stokesProblem->setInjectionState(true);
+ else
+ stokesProblem->setInjectionState(false);
+
+ // solve the non-linear system with time step control
+ nonLinearSolver.solve(sol, *timeLoop);
+
+ // make the new solution the old solution
+ solOld = sol;
+ stokesGridVariables->advanceTimeStep();
+ darcyGridVariables->advanceTimeStep();
+
+ // advance to the time loop to the next step
+ timeLoop->advanceTimeStep();
+
+ // write vtk output
+ stokesVtkWriter.write(timeLoop->time());
+ darcyVtkWriter.write(timeLoop->time());
+
+ // report statistics of this time step
+ timeLoop->reportTimeStep();
+
+ // set new dt as suggested by newton solver
+ timeLoop->setTimeStepSize(nonLinearSolver.suggestTimeStepSize(timeLoop->timeStepSize()));
+
+ } while (!timeLoop->finished());
+
+ timeLoop->finalize(stokesGridView.comm());
+ timeLoop->finalize(darcyGridView.comm());
+
+ ////////////////////////////////////////////////////////////
+ // finalize, print dumux message to say goodbye
+ ////////////////////////////////////////////////////////////
+
+ // print dumux end message
+ if (mpiHelper.rank() == 0)
+ {
+ Parameters::print();
+ DumuxMessage::print(/*firstCall=*/false);
+ }
+
+ return 0;
+} // end main
+catch (Dumux::ParameterException &e)
+{
+ std::cerr << std::endl << e << " ---> Abort!" << std::endl;
+ return 1;
+}
+catch (Dune::DGFException & e)
+{
+ std::cerr << "DGF exception thrown (" << e <<
+ "). Most likely, the DGF file name is wrong "
+ "or the DGF file is corrupted, "
+ "e.g. missing hash at end of file or wrong number (dimensions) of entries."
+ << " ---> Abort!" << std::endl;
+ return 2;
+}
+catch (Dune::Exception &e)
+{
+ std::cerr << "Dune reported error: " << e << " ---> Abort!" << std::endl;
+ return 3;
+}
+catch (...)
+{
+ std::cerr << "Unknown exception thrown! ---> Abort!" << std::endl;
+ return 4;
+}
diff --git a/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/params.input b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/params.input
new file mode 100644
index 0000000000..fec9099878
--- /dev/null
+++ b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/params.input
@@ -0,0 +1,42 @@
+[TimeLoop]
+TEnd = 3e6 # s
+DtInitial = 100 # s
+
+[Darcy.Grid]
+UpperRight = 1.0 1.0
+Cells = 20 20
+
+[Stokes.Grid]
+LowerLeft = 0.0 1.0
+UpperRight = 1.0 2.0
+Cells = 20 20
+
+[Stokes.Problem]
+Name = stokes
+Velocity = 1e-6
+Pressure = 1.0e5
+InletMoleFraction = 1e-3
+InjectionBegin = 0
+InjectionEnd = 3e6
+
+[Darcy.Problem]
+Name = darcy
+Pressure = 1.0e5
+InitialMoleFraction = 0.0
+
+[SpatialParams]
+AlphaBeaversJoseph = 1.0
+Permeability = 1e-10 # m^2
+Porosity = 0.3
+Tortuosity = 0.5
+
+[Problem]
+Name = test_md_darcy1p2c_stokes1p2c_comparison
+EnableGravity = false
+EnableInertiaTerms = true
+
+[Vtk]
+AddVelocity = 1
+
+[Assembly.NumericDifference]
+BaseEpsilon = 1e-5
diff --git a/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/problem_darcy.hh b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/problem_darcy.hh
new file mode 100644
index 0000000000..78ae3f3698
--- /dev/null
+++ b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/problem_darcy.hh
@@ -0,0 +1,281 @@
+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
+// vi: set et ts=4 sw=4 sts=4:
+/*****************************************************************************
+ * See the file COPYING for full copying permissions. *
+ * *
+ * This program is free software: you can redistribute it and/or modify *
+ * it under the terms of the GNU General Public License as published by *
+ * the Free Software Foundation, either version 2 of the License, or *
+ * (at your option) any later version. *
+ * *
+ * This program is distributed in the hope that it will be useful, *
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of *
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
+ * GNU General Public License for more details. *
+ * *
+ * You should have received a copy of the GNU General Public License *
+ * along with this program. If not, see <http://www.gnu.org/licenses/>. *
+ *****************************************************************************/
+/*!
+ * \file
+ *
+ * \brief A simple Darcy test problem (cell-centered finite volume method) for the comparison of different diffusion laws.
+ */
+#ifndef DUMUX_DARCY_SUBPROBLEM_DIFFUSION_COMPARISON_HH
+#define DUMUX_DARCY_SUBPROBLEM_DIFFUSION_COMPARISON_HH
+
+#include <dune/grid/yaspgrid.hh>
+
+#include <dumux/discretization/cellcentered/tpfa/properties.hh>
+#include <dumux/flux/maxwellstefanslaw.hh>
+
+#include <dumux/porousmediumflow/1pnc/model.hh>
+#include <dumux/porousmediumflow/problem.hh>
+
+#include "./../spatialparams.hh"
+
+#include <dumux/material/fluidsystems/1padapter.hh>
+#include <dumux/material/fluidsystems/h2oair.hh>
+#include <dumux/material/fluidmatrixinteractions/diffusivityconstanttortuosity.hh>
+
+namespace Dumux
+{
+template <class TypeTag>
+class DarcySubProblem;
+
+namespace Properties
+{
+// Create new type tags
+namespace TTag {
+struct DarcyOnePTwoCTypeTag { using InheritsFrom = std::tuple<OnePNC, CCTpfaModel>; };
+} // end namespace TTag
+
+// Set the problem property
+template<class TypeTag>
+struct Problem<TypeTag, TTag::DarcyOnePTwoCTypeTag> { using type = Dumux::DarcySubProblem<TypeTag>; };
+
+// The fluid system
+template<class TypeTag>
+struct FluidSystem<TypeTag, TTag::DarcyOnePTwoCTypeTag>
+{
+ using H2OAir = FluidSystems::H2OAir<GetPropType<TypeTag, Properties::Scalar>>;
+ static constexpr auto phaseIdx = H2OAir::liquidPhaseIdx; // simulate the water phase
+ using type = FluidSystems::OnePAdapter<H2OAir, phaseIdx>;
+};
+
+// Use moles
+template<class TypeTag>
+struct UseMoles<TypeTag, TTag::DarcyOnePTwoCTypeTag> { static constexpr bool value = true; };
+
+// Do not replace one equation with a total mass balance
+template<class TypeTag>
+struct ReplaceCompEqIdx<TypeTag, TTag::DarcyOnePTwoCTypeTag> { static constexpr int value = 3; };
+
+//! Use a model with constant tortuosity for the effective diffusivity
+SET_TYPE_PROP(DarcyOnePTwoCTypeTag, EffectiveDiffusivityModel,
+ DiffusivityConstantTortuosity<GetPropType<TypeTag, Properties::Scalar>>);
+
+// Set the grid type
+template<class TypeTag>
+struct Grid<TypeTag, TTag::DarcyOnePTwoCTypeTag> { using type = Dune::YaspGrid<2>; };
+
+// Set the diffusion type
+template<class TypeTag>
+struct MolecularDiffusionType<TypeTag, TTag::DarcyOnePTwoCTypeTag> { using type = DIFFUSIONTYPE; };
+
+// Set the spatial paramaters type
+template<class TypeTag>
+struct SpatialParams<TypeTag, TTag::DarcyOnePTwoCTypeTag>
+{
+ using FVGridGeometry = GetPropType<TypeTag, Properties::FVGridGeometry>;
+ using Scalar = GetPropType<TypeTag, Properties::Scalar>;
+ using type = OnePSpatialParams<FVGridGeometry, Scalar>;
+};
+}
+
+template <class TypeTag>
+class DarcySubProblem : public PorousMediumFlowProblem<TypeTag>
+{
+ using ParentType = PorousMediumFlowProblem<TypeTag>;
+ using GridView = GetPropType<TypeTag, Properties::GridView>;
+ using Scalar = GetPropType<TypeTag, Properties::Scalar>;
+ using PrimaryVariables = GetPropType<TypeTag, Properties::PrimaryVariables>;
+ using FluidSystem = GetPropType<TypeTag, Properties::FluidSystem>;
+ using NumEqVector = GetPropType<TypeTag, Properties::NumEqVector>;
+ using BoundaryTypes = GetPropType<TypeTag, Properties::BoundaryTypes>;
+ using FVElementGeometry = typename GetPropType<TypeTag, Properties::FVGridGeometry>::LocalView;
+ using SubControlVolume = typename FVElementGeometry::SubControlVolume;
+ using SubControlVolumeFace = typename FVElementGeometry::SubControlVolumeFace;
+ using FVGridGeometry = GetPropType<TypeTag, Properties::FVGridGeometry>;
+ using DiffusionCoefficientAveragingType = typename StokesDarcyCouplingOptions::DiffusionCoefficientAveragingType;
+
+ // copy some indices for convenience
+ using Indices = typename GetPropType<TypeTag, Properties::ModelTraits>::Indices;
+ enum {
+ // grid and world dimension
+ dim = GridView::dimension,
+ dimworld = GridView::dimensionworld,
+
+ // primary variable indices
+ conti0EqIdx = Indices::conti0EqIdx,
+ pressureIdx = Indices::pressureIdx,
+ };
+
+ using Element = typename GridView::template Codim<0>::Entity;
+ using GlobalPosition = Dune::FieldVector<Scalar, dimworld>;
+
+ using CouplingManager = GetPropType<TypeTag, Properties::CouplingManager>;
+
+public:
+ DarcySubProblem(std::shared_ptr<const FVGridGeometry> fvGridGeometry,
+ std::shared_ptr<CouplingManager> couplingManager)
+ : ParentType(fvGridGeometry, "Darcy"), eps_(1e-7), couplingManager_(couplingManager)
+ {
+ pressure_ = getParamFromGroup<Scalar>(this->paramGroup(), "Problem.Pressure");
+ initialMoleFraction_ = getParamFromGroup<Scalar>(this->paramGroup(), "Problem.InitialMoleFraction");
+ }
+
+ /*!
+ * \name Simulation steering
+ */
+ // \{
+
+ /*!
+ * \brief Returns true if a restart file should be written to
+ * disk.
+ */
+ bool shouldWriteRestartFile() const
+ { return false; }
+
+ /*!
+ * \name Problem parameters
+ */
+ // \{
+
+ bool shouldWriteOutput() const // define output
+ { return true; }
+
+ /*!
+ * \brief Return the temperature within the domain in [K].
+ *
+ */
+ Scalar temperature() const
+ { return 273.15 + 10; } // 10°C
+ // \}
+
+ /*!
+ * \name Boundary conditions
+ */
+ // \{
+
+ /*!
+ * \brief Specifies which kind of boundary condition should be
+ * used for which equation on a given boundary control volume.
+ *
+ * \param element The element
+ * \param scvf The boundary sub control volume face
+ */
+ BoundaryTypes boundaryTypes(const Element& element, const SubControlVolumeFace& scvf) const
+ {
+ BoundaryTypes values;
+ values.setAllNeumann();
+
+ if (couplingManager().isCoupledEntity(CouplingManager::darcyIdx, scvf))
+ values.setAllCouplingNeumann();
+
+ return values;
+ }
+
+ /*!
+ * \brief Evaluate the boundary conditions for a Neumann control volume.
+ *
+ * \param element The element for which the Neumann boundary condition is set
+ * \param fvGeomentry The fvGeometry
+ * \param elemVolVars The element volume variables
+ * \param scvf The boundary sub control volume face
+ *
+ * For this method, the \a values variable stores primary variables.
+ */
+ template<class ElementVolumeVariables>
+ NumEqVector neumann(const Element& element,
+ const FVElementGeometry& fvGeometry,
+ const ElementVolumeVariables& elemVolVars,
+ const SubControlVolumeFace& scvf) const
+ {
+ NumEqVector values(0.0);
+
+ if (couplingManager().isCoupledEntity(CouplingManager::darcyIdx, scvf))
+ values = couplingManager().couplingData().massCouplingCondition(fvGeometry, elemVolVars, scvf, DiffusionCoefficientAveragingType::harmonic);
+
+ return values;
+ }
+
+ // \}
+
+ /*!
+ * \name Volume terms
+ */
+ // \{
+ /*!
+ * \brief Evaluate the source term for all phases within a given
+ * sub-control-volume.
+ *
+ * \param element The element for which the source term is set
+ * \param fvGeomentry The fvGeometry
+ * \param elemVolVars The element volume variables
+ * \param scv The subcontrolvolume
+ */
+ template<class ElementVolumeVariables>
+ NumEqVector source(const Element &element,
+ const FVElementGeometry& fvGeometry,
+ const ElementVolumeVariables& elemVolVars,
+ const SubControlVolume &scv) const
+ { return NumEqVector(0.0); }
+
+ // \}
+
+ /*!
+ * \brief Evaluate the initial value for a control volume.
+ *
+ * \param element The element
+ *
+ * For this method, the \a priVars parameter stores primary
+ * variables.
+ */
+ PrimaryVariables initial(const Element &element) const
+ {
+ PrimaryVariables values(0.0);
+ values[pressureIdx] = pressure_;
+ values[conti0EqIdx + 1] = initialMoleFraction_;
+
+ return values;
+ }
+
+ // \}
+
+ //! Get the coupling manager
+ const CouplingManager& couplingManager() const
+ { return *couplingManager_; }
+
+private:
+ bool onLeftBoundary_(const GlobalPosition &globalPos) const
+ { return globalPos[0] < this->fvGridGeometry().bBoxMin()[0] + eps_; }
+
+ bool onRightBoundary_(const GlobalPosition &globalPos) const
+ { return globalPos[0] > this->fvGridGeometry().bBoxMax()[0] - eps_; }
+
+ bool onLowerBoundary_(const GlobalPosition &globalPos) const
+ { return globalPos[1] < this->fvGridGeometry().bBoxMin()[1] + eps_; }
+
+ bool onUpperBoundary_(const GlobalPosition &globalPos) const
+ { return globalPos[1] > this->fvGridGeometry().bBoxMax()[1] - eps_; }
+
+ Scalar eps_;
+ Scalar pressure_;
+ Scalar initialMoleFraction_;
+
+ std::shared_ptr<CouplingManager> couplingManager_;
+};
+} //end namespace
+
+#endif //DUMUX_DARCY_SUBPROBLEM_HH
diff --git a/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/problem_stokes.hh b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/problem_stokes.hh
new file mode 100644
index 0000000000..8ff09445a9
--- /dev/null
+++ b/test/multidomain/boundary/stokesdarcy/1p2c_1p2c/diffusionlawcomparison/problem_stokes.hh
@@ -0,0 +1,322 @@
+// -*- mode: C++; tab-width: 4; indent-tabs-mode: nil; c-basic-offset: 4 -*-
+// vi: set et ts=4 sw=4 sts=4:
+/*****************************************************************************
+ * See the file COPYING for full copying permissions. *
+ * *
+ * This program is free software: you can redistribute it and/or modify *
+ * it under the terms of the GNU General Public License as published by *
+ * the Free Software Foundation, either version 2 of the License, or *
+ * (at your option) any later version. *
+ * *
+ * This program is distributed in the hope that it will be useful, *
+ * but WITHOUT ANY WARRANTY; without even the implied warranty of *
+ * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the *
+ * GNU General Public License for more details. *
+ * *
+ * You should have received a copy of the GNU General Public License *
+ * along with this program. If not, see <http://www.gnu.org/licenses/>. *
+ *****************************************************************************/
+/*!
+ * \file
+ * \ingroup NavierStokesTests
+ * \brief A simple Navier-Stokes test problem for the staggered grid (Navier-)Stokes model. This test uses the same set-up as the horizontal test but uses a harmonic averaging of the diffusion coefficients at the interface for ficks law and comares the solution with the maxwell stefan coupling
+ */
+#ifndef DUMUX_STOKES_SUBPROBLEM_DIFFUSION_COMPARISON_HH
+#define DUMUX_STOKES_SUBPROBLEM_DIFFUSION_COMPARISON_HH
+
+#include <dune/grid/yaspgrid.hh>
+
+#include <dumux/material/fluidsystems/1padapter.hh>
+#include <dumux/material/fluidsystems/h2oair.hh>
+#include <dumux/flux/maxwellstefanslaw.hh>
+
+#include <dumux/freeflow/navierstokes/problem.hh>
+#include <dumux/discretization/staggered/freeflow/properties.hh>
+#include <dumux/freeflow/compositional/navierstokesncmodel.hh>
+
+namespace Dumux
+{
+template <class TypeTag>
+class StokesSubProblem;
+
+namespace Properties
+{
+// Create new type tags
+namespace TTag {
+struct StokesOnePTwoCTypeTag { using InheritsFrom = std::tuple<NavierStokesNC, StaggeredFreeFlowModel>; };
+} // end namespace TTag
+
+// The fluid system
+template<class TypeTag>
+struct FluidSystem<TypeTag, TTag::StokesOnePTwoCTypeTag>
+{
+ using H2OAir = FluidSystems::H2OAir<GetPropType<TypeTag, Properties::Scalar>>;
+ static constexpr auto phaseIdx = H2OAir::liquidPhaseIdx; // simulate the water phase
+ using type = FluidSystems::OnePAdapter<H2OAir, phaseIdx>;
+};
+
+// Set the grid type
+template<class TypeTag>
+struct Grid<TypeTag, TTag::StokesOnePTwoCTypeTag> { using type = Dune::YaspGrid<2, Dune::EquidistantOffsetCoordinates<GetPropType<TypeTag, Properties::Scalar>, 2> >; };
+
+// Set the problem property
+template<class TypeTag>
+struct Problem<TypeTag, TTag::StokesOnePTwoCTypeTag> { using type = Dumux::StokesSubProblem<TypeTag> ; };
+
+template<class TypeTag>
+struct EnableFVGridGeometryCache<TypeTag, TTag::StokesOnePTwoCTypeTag> { static constexpr bool value = true; };
+template<class TypeTag>
+struct EnableGridFluxVariablesCache<TypeTag, TTag::StokesOnePTwoCTypeTag> { static constexpr bool value = true; };
+template<class TypeTag>
+struct EnableGridVolumeVariablesCache<TypeTag, TTag::StokesOnePTwoCTypeTag> { static constexpr bool value = true; };
+
+// Use moles
+template<class TypeTag>
+struct UseMoles<TypeTag, TTag::StokesOnePTwoCTypeTag> { static constexpr bool value = true; };
+
+// Set the grid type
+template<class TypeTag>
+struct MolecularDiffusionType<TypeTag, TTag::StokesOnePTwoCTypeTag> { using type = DIFFUSIONTYPE; };
+
+// Do not replace one equation with a total mass balance
+template<class TypeTag>
+struct ReplaceCompEqIdx<TypeTag, TTag::StokesOnePTwoCTypeTag> { static constexpr int value = 3; };
+}
+
+/*!
+ * \ingroup NavierStokesTests
+ * \brief Test problem for the 1pnc (Navier-) Stokes problem.
+ *
+ * Horizontal flow from left to right with a parabolic velocity profile.
+ */
+template <class TypeTag>
+class StokesSubProblem : public NavierStokesProblem<TypeTag>
+{
+ using ParentType = NavierStokesProblem<TypeTag>;
+ using GridView = GetPropType<TypeTag, Properties::GridView>;
+ using Scalar = GetPropType<TypeTag, Properties::Scalar>;
+ using Indices = typename GetPropType<TypeTag, Properties::ModelTraits>::Indices;
+ using FluidSystem = GetPropType<TypeTag, Properties::FluidSystem>;
+ using BoundaryTypes = GetPropType<TypeTag, Properties::BoundaryTypes>;
+ using FVGridGeometry = GetPropType<TypeTag, Properties::FVGridGeometry>;
+ using FVElementGeometry = typename FVGridGeometry::LocalView;
+ using SubControlVolumeFace = typename FVElementGeometry::SubControlVolumeFace;
+ using PrimaryVariables = GetPropType<TypeTag, Properties::PrimaryVariables>;
+ using NumEqVector = GetPropType<TypeTag, Properties::NumEqVector>;
+
+ using Element = typename GridView::template Codim<0>::Entity;
+ using GlobalPosition = typename Element::Geometry::GlobalCoordinate;
+
+ using CouplingManager = GetPropType<TypeTag, Properties::CouplingManager>;
+ using DiffusionCoefficientAveragingType = typename StokesDarcyCouplingOptions::DiffusionCoefficientAveragingType;
+
+public:
+ StokesSubProblem(std::shared_ptr<const FVGridGeometry> fvGridGeometry, std::shared_ptr<CouplingManager> couplingManager)
+ : ParentType(fvGridGeometry, "Stokes"), eps_(1e-6), injectionState_(false), couplingManager_(couplingManager)
+ {
+ inletVelocity_ = getParamFromGroup<Scalar>(this->paramGroup(), "Problem.Velocity");
+ pressure_ = getParamFromGroup<Scalar>(this->paramGroup(), "Problem.Pressure");
+ inletMoleFraction_ = getParamFromGroup<Scalar>(this->paramGroup(), "Problem.InletMoleFraction");
+ }
+
+ /*!
+ * \name Problem parameters
+ */
+ // \{
+
+ bool shouldWriteRestartFile() const
+ { return false; }
+
+ /*!
+ * \brief Return the temperature within the domain in [K].
+ *
+ * This problem assumes a temperature of 10 degrees Celsius.
+ */
+ Scalar temperature() const
+ { return 273.15 + 10; } // 10°C
+
+ /*!
+ * \brief Return the sources within the domain.
+ *
+ * \param globalPos The global position
+ */
+ NumEqVector sourceAtPos(const GlobalPosition &globalPos) const
+ { return NumEqVector(0.0); }
+ // \}
+
+ /*!
+ * \name Boundary conditions
+ */
+ // \{
+
+ /*!
+ * \brief Specifies which kind of boundary condition should be
+ * used for which equation on a given boundary segment.
+ *
+ * \param element The finite element
+ * \param scvf The sub control volume face
+ */
+ BoundaryTypes boundaryTypes(const Element& element,
+ const SubControlVolumeFace& scvf) const
+ {
+ BoundaryTypes values;
+
+ const auto& globalPos = scvf.dofPosition();
+
+ if(onLeftBoundary_(globalPos))
+ {
+ values.setDirichlet(Indices::conti0EqIdx + 1);
+ values.setDirichlet(Indices::velocityXIdx);
+ values.setDirichlet(Indices::velocityYIdx);
+ }
+ else if(onRightBoundary_(globalPos))
+ {
+ values.setDirichlet(Indices::pressureIdx);
+ values.setOutflow(Indices::conti0EqIdx + 1);
+ }
+ else
+ {
+ values.setDirichlet(Indices::velocityXIdx);
+ values.setDirichlet(Indices::velocityYIdx);
+ values.setNeumann(Indices::conti0EqIdx);
+ values.setNeumann(Indices::conti0EqIdx + 1);
+ }
+
+ if(couplingManager().isCoupledEntity(CouplingManager::stokesIdx, scvf))
+ {
+ values.setCouplingNeumann(Indices::conti0EqIdx);
+ values.setCouplingNeumann(Indices::conti0EqIdx + 1);
+ values.setCouplingNeumann(Indices::momentumYBalanceIdx);
+ values.setBJS(Indices::momentumXBalanceIdx);
+ }
+
+ return values;
+ }
+
+ /*!
+ * \brief Evaluate the boundary conditions for a Dirichlet control volume.
+ *
+ * \param element The element
+ * \param scvf The sub control volume face
+ */
+ PrimaryVariables dirichletAtPos(const GlobalPosition& globalPos) const
+ {
+ PrimaryVariables values(0.0);
+ values = initialAtPos(globalPos);
+
+ // start injecting after the velocity field had enough time to initialize
+ if(globalPos[0] < this->fvGridGeometry().bBoxMin()[0] + eps_ && isInjectionPeriod())
+ values[Indices::conti0EqIdx + 1] = inletMoleFraction_;
+
+ return values;
+ }
+
+ /*!
+ * \brief Evaluate the boundary conditions for a Neumann control volume.
+ *
+ * \param element The element for which the Neumann boundary condition is set
+ * \param fvGeomentry The fvGeometry
+ * \param elemVolVars The element volume variables
+ * \param elemFaceVars The element face variables
+ * \param scvf The boundary sub control volume face
+ */
+ template<class ElementVolumeVariables, class ElementFaceVariables>
+ NumEqVector neumann(const Element& element,
+ const FVElementGeometry& fvGeometry,
+ const ElementVolumeVariables& elemVolVars,
+ const ElementFaceVariables& elemFaceVars,
+ const SubControlVolumeFace& scvf) const
+ {
+ NumEqVector values(0.0);
+
+ if(couplingManager().isCoupledEntity(CouplingManager::stokesIdx, scvf))
+ {
+ values[Indices::momentumYBalanceIdx] = couplingManager().couplingData().momentumCouplingCondition(fvGeometry, elemVolVars, elemFaceVars, scvf);
+
+ const auto tmp = couplingManager().couplingData().massCouplingCondition(fvGeometry, elemVolVars, elemFaceVars, scvf, DiffusionCoefficientAveragingType::harmonic);
+ values[Indices::conti0EqIdx] = tmp[0];
+ values[Indices::conti0EqIdx + 1] = tmp[1];
+ }
+ return values;
+ }
+
+ // \}
+
+ //! Get the coupling manager
+ const CouplingManager& couplingManager() const
+ { return *couplingManager_; }
+
+ /*!
+ * \name Volume terms
+ */
+ // \{
+
+ /*!
+ * \brief Evaluate the initial value for a control volume.
+ *
+ * \param globalPos The global position
+ */
+ PrimaryVariables initialAtPos(const GlobalPosition &globalPos) const
+ {
+ PrimaryVariables values(0.0);
+ values[Indices::pressureIdx] = pressure_;
+ values[Indices::velocityXIdx] = inletVelocity_ * (globalPos[1] - this->fvGridGeometry().bBoxMin()[1])
+ * (this->fvGridGeometry().bBoxMax()[1] - globalPos[1])
+ / (0.25 * (this->fvGridGeometry().bBoxMax()[1] - this->fvGridGeometry().bBoxMin()[1])
+ * (this->fvGridGeometry().bBoxMax()[1] - this->fvGridGeometry().bBoxMin()[1]));
+
+ return values;
+ }
+
+ /*!
+ * \brief Returns the intrinsic permeability of required as input parameter for the Beavers-Joseph-Saffman boundary condition
+ */
+ Scalar permeability(const SubControlVolumeFace& scvf) const
+ {
+ return couplingManager().couplingData().darcyPermeability(scvf);
+ }
+
+ /*!
+ * \brief Returns the alpha value required as input parameter for the Beavers-Joseph-Saffman boundary condition
+ */
+ Scalar alphaBJ(const SubControlVolumeFace& scvf) const
+ {
+ return 1.0;
+ }
+
+ void setInjectionState(const bool yesNo)
+ {
+ injectionState_ = yesNo;
+ }
+
+ bool isInjectionPeriod() const
+ {
+ return injectionState_;
+ }
+
+ // \}
+
+private:
+ bool onLeftBoundary_(const GlobalPosition &globalPos) const
+ { return globalPos[0] < this->fvGridGeometry().bBoxMin()[0] + eps_; }
+
+ bool onRightBoundary_(const GlobalPosition &globalPos) const
+ { return globalPos[0] > this->fvGridGeometry().bBoxMax()[0] - eps_; }
+
+ bool onLowerBoundary_(const GlobalPosition &globalPos) const
+ { return globalPos[1] < this->fvGridGeometry().bBoxMin()[1] + eps_; }
+
+ bool onUpperBoundary_(const GlobalPosition &globalPos) const
+ { return globalPos[1] > this->fvGridGeometry().bBoxMax()[1] - eps_; }
+
+ Scalar eps_;
+ Scalar inletVelocity_;
+ Scalar pressure_;
+ Scalar inletMoleFraction_;
+ bool injectionState_;
+
+ std::shared_ptr<CouplingManager> couplingManager_;
+};
+} //end namespace
+
+#endif // DUMUX_STOKES_SUBPROBLEM_HH
diff --git a/test/references/test_md_darcy1p2c_stokes1p2c_comparison_darcy-reference.vtu b/test/references/test_md_darcy1p2c_stokes1p2c_comparison_darcy-reference.vtu
new file mode 100644
index 0000000000..c4f30a42ef
--- /dev/null
+++ b/test/references/test_md_darcy1p2c_stokes1p2c_comparison_darcy-reference.vtu
@@ -0,0 +1,658 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="400" NumberOfPoints="441">
+ <CellData Scalars="p">
+ <DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000
+ </DataArray>
+ <DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.69 999.69 999.69 999.689 999.689 999.689 999.689 999.689
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689
+ 999.692 999.692 999.691 999.691 999.69 999.69 999.69 999.69 999.69 999.69 999.69 999.69
+ 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.689 999.7 999.699 999.697 999.696
+ 999.695 999.694 999.693 999.692 999.692 999.692 999.691 999.691 999.691 999.691 999.69 999.69
+ 999.69 999.69 999.69 999.69 999.728 999.724 999.719 999.714 999.71 999.707 999.704 999.702
+ 999.701 999.699 999.698 999.697 999.696 999.695 999.694 999.694 999.693 999.693 999.693 999.692
+ 999.814 999.803 999.788 999.774 999.762 999.752 999.744 999.737 999.732 999.727 999.722 999.719
+ 999.715 999.713 999.71 999.708 999.706 999.705 999.703 999.702 1000.04 1000.01 999.979 999.948
+ 999.92 999.895 999.872 999.853 999.836 999.821 999.808 999.796 999.785 999.776 999.768 999.761
+ 999.754 999.748 999.743 999.74
+ </DataArray>
+ <DataArray type="Float32" Name="mu" NumberOfComponents="1" format="ascii">
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753 0.00130753
+ 0.00130753 0.00130753 0.00130753 0.00130753
+ </DataArray>
+ <DataArray type="Float32" Name="delp" NumberOfComponents="1" format="ascii">
+ 5.57338e-09 5.55883e-09 5.54428e-09 5.51518e-09 5.48607e-09 5.44242e-09 5.39876e-09 5.34055e-09 5.2969e-09 5.23869e-09 5.18048e-09 5.12227e-09
+ 5.06407e-09 5.02041e-09 4.9622e-09 4.9331e-09 4.88944e-09 4.86034e-09 4.84579e-09 4.84579e-09 5.57338e-09 5.57338e-09 5.54428e-09 5.52973e-09
+ 5.48607e-09 5.44242e-09 5.39876e-09 5.34055e-09 5.2969e-09 5.23869e-09 5.18048e-09 5.12227e-09 5.06407e-09 5.00586e-09 4.9622e-09 4.91855e-09
+ 4.88944e-09 4.86034e-09 4.84579e-09 4.83124e-09 5.60249e-09 5.58794e-09 5.57338e-09 5.54428e-09 5.50062e-09 5.45697e-09 5.41331e-09 5.3551e-09
+ 5.2969e-09 5.23869e-09 5.18048e-09 5.10772e-09 5.04951e-09 5.00586e-09 4.94765e-09 4.904e-09 4.87489e-09 4.84579e-09 4.81668e-09 4.81668e-09
+ 5.61704e-09 5.61704e-09 5.58794e-09 5.55883e-09 5.52973e-09 5.47152e-09 5.42786e-09 5.36966e-09 5.2969e-09 5.23869e-09 5.16593e-09 5.10772e-09
+ 5.04951e-09 4.99131e-09 4.9331e-09 4.88944e-09 4.84579e-09 4.81668e-09 4.80213e-09 4.78758e-09 5.6607e-09 5.64614e-09 5.63159e-09 5.60249e-09
+ 5.55883e-09 5.50062e-09 5.44242e-09 5.38421e-09 5.31145e-09 5.23869e-09 5.16593e-09 5.09317e-09 5.03496e-09 4.9622e-09 4.904e-09 4.86034e-09
+ 4.81668e-09 4.78758e-09 4.75848e-09 4.74392e-09 5.70435e-09 5.70435e-09 5.67525e-09 5.63159e-09 5.58794e-09 5.52973e-09 5.47152e-09 5.39876e-09
+ 5.326e-09 5.23869e-09 5.16593e-09 5.09317e-09 5.00586e-09 4.9331e-09 4.87489e-09 4.81668e-09 4.77303e-09 4.72937e-09 4.71482e-09 4.70027e-09
+ 5.77711e-09 5.76256e-09 5.73345e-09 5.6898e-09 5.63159e-09 5.57338e-09 5.50062e-09 5.42786e-09 5.34055e-09 5.25324e-09 5.16593e-09 5.07862e-09
+ 4.99131e-09 4.904e-09 4.83124e-09 4.77303e-09 4.71482e-09 4.67116e-09 4.64206e-09 4.64206e-09 5.84987e-09 5.83532e-09 5.80621e-09 5.74801e-09
+ 5.6898e-09 5.61704e-09 5.54428e-09 5.45697e-09 5.3551e-09 5.25324e-09 5.15138e-09 5.04951e-09 4.9622e-09 4.87489e-09 4.78758e-09 4.71482e-09
+ 4.65661e-09 4.61296e-09 4.5693e-09 4.55475e-09 5.93718e-09 5.92263e-09 5.89353e-09 5.83532e-09 5.76256e-09 5.67525e-09 5.58794e-09 5.48607e-09
+ 5.36966e-09 5.26779e-09 5.15138e-09 5.03496e-09 4.91855e-09 4.81668e-09 4.72937e-09 4.64206e-09 4.5693e-09 4.52565e-09 4.48199e-09 4.46744e-09
+ 6.0536e-09 6.02449e-09 5.99539e-09 5.92263e-09 5.84987e-09 5.74801e-09 5.64614e-09 5.52973e-09 5.39876e-09 5.26779e-09 5.13683e-09 5.00586e-09
+ 4.88944e-09 4.75848e-09 4.65661e-09 4.55475e-09 4.48199e-09 4.42378e-09 4.38013e-09 4.35102e-09 6.18456e-09 6.15546e-09 6.1118e-09 6.03904e-09
+ 5.95173e-09 5.83532e-09 5.70435e-09 5.57338e-09 5.42786e-09 5.28235e-09 5.12227e-09 4.97676e-09 4.83124e-09 4.70027e-09 4.5693e-09 4.46744e-09
+ 4.36557e-09 4.29281e-09 4.24916e-09 4.22006e-09 6.34464e-09 6.31553e-09 6.25732e-09 6.17001e-09 6.06815e-09 5.93718e-09 5.79166e-09 5.63159e-09
+ 5.47152e-09 5.2969e-09 5.12227e-09 4.94765e-09 4.77303e-09 4.62751e-09 4.48199e-09 4.35102e-09 4.23461e-09 4.1473e-09 4.08909e-09 4.05998e-09
+ 6.53381e-09 6.49015e-09 6.43195e-09 6.33008e-09 6.19912e-09 6.0536e-09 5.87897e-09 5.70435e-09 5.50062e-09 5.31145e-09 5.10772e-09 4.904e-09
+ 4.71482e-09 4.52565e-09 4.36557e-09 4.2055e-09 4.08909e-09 3.98722e-09 3.91447e-09 3.87081e-09 6.75209e-09 6.70843e-09 6.62112e-09 6.50471e-09
+ 6.35919e-09 6.18456e-09 5.98084e-09 5.77711e-09 5.54428e-09 5.326e-09 5.09317e-09 4.86034e-09 4.64206e-09 4.42378e-09 4.23461e-09 4.04543e-09
+ 3.89991e-09 3.7835e-09 3.69619e-09 3.65253e-09 7.01402e-09 6.97037e-09 6.8685e-09 6.72298e-09 6.54836e-09 6.33008e-09 6.1118e-09 5.86442e-09
+ 5.60249e-09 5.34055e-09 5.07862e-09 4.81668e-09 4.55475e-09 4.30737e-09 4.07454e-09 3.87081e-09 3.68163e-09 3.53612e-09 3.43425e-09 3.3906e-09
+ 7.33417e-09 7.27596e-09 7.14499e-09 6.97037e-09 6.75209e-09 6.50471e-09 6.24277e-09 5.95173e-09 5.6607e-09 5.3551e-09 5.04951e-09 4.75848e-09
+ 4.45289e-09 4.1764e-09 3.89991e-09 3.65253e-09 3.43425e-09 3.25963e-09 3.12866e-09 3.07045e-09 7.72707e-09 7.63976e-09 7.47968e-09 7.26141e-09
+ 6.99947e-09 6.70843e-09 6.38829e-09 6.06815e-09 5.7189e-09 5.38421e-09 5.03496e-09 4.68572e-09 4.35102e-09 4.01633e-09 3.69619e-09 3.40515e-09
+ 3.14321e-09 2.92493e-09 2.76486e-09 2.67755e-09 8.19273e-09 8.07631e-09 7.85803e-09 7.5961e-09 7.27596e-09 6.92671e-09 6.56291e-09 6.18456e-09
+ 5.79166e-09 5.39876e-09 5.00586e-09 4.61296e-09 4.23461e-09 3.85626e-09 3.47791e-09 3.14321e-09 2.82307e-09 2.53203e-09 2.32831e-09 2.19734e-09
+ 8.7748e-09 8.60018e-09 8.3237e-09 7.9599e-09 7.58155e-09 7.17409e-09 6.73754e-09 6.31553e-09 5.87897e-09 5.42786e-09 4.99131e-09 4.5402e-09
+ 4.10364e-09 3.66708e-09 3.24508e-09 2.83762e-09 2.44472e-09 2.08092e-09 1.78989e-09 1.60071e-09 9.54606e-09 9.24047e-09 8.83301e-09 8.3819e-09
+ 7.90169e-09 7.42148e-09 6.94126e-09 6.4465e-09 5.95173e-09 5.45697e-09 4.9622e-09 4.46744e-09 3.97267e-09 3.47791e-09 2.99769e-09 2.50293e-09
+ 2.03727e-09 1.57161e-09 1.1496e-09 8.14907e-10
+ </DataArray>
+ <DataArray type="Float32" Name="x^H2O_liquid" 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 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 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 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 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 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
+ 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 0.999999 0.999999 0.999999 0.999999 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 0.999995 0.999996 0.999996 0.999997 0.999998 0.999998 0.999998 0.999999 0.999999 0.999999 0.999999 0.999999
+ 0.999999 0.999999 0.999999 1 1 1 1 1 0.999982 0.999984 0.999986 0.999989
+ 0.999991 0.999992 0.999993 0.999994 0.999995 0.999996 0.999996 0.999997 0.999997 0.999997 0.999998 0.999998
+ 0.999998 0.999998 0.999999 0.999999 0.999936 0.999942 0.999951 0.999959 0.999965 0.99997 0.999975 0.999978
+ 0.999981 0.999983 0.999986 0.999987 0.999989 0.99999 0.999991 0.999992 0.999993 0.999994 0.999994 0.999994
+ 0.999794 0.999813 0.999837 0.99986 0.999879 0.999896 0.999909 0.99992 0.99993 0.999938 0.999945 0.999951
+ 0.999957 0.999961 0.999965 0.999969 0.999972 0.999974 0.999977 0.999978 0.999428 0.99947 0.999523 0.999574
+ 0.99962 0.999662 0.999698 0.99973 0.999758 0.999783 0.999805 0.999824 0.999841 0.999856 0.99987 0.999882
+ 0.999893 0.999903 0.999911 0.999916
+ </DataArray>
+ <DataArray type="Float32" Name="x^Air_liquid" NumberOfComponents="1" format="ascii">
+ 1.50841e-15 1.43417e-15 1.30183e-15 1.13667e-15 9.63829e-16 8.01879e-16 6.60899e-16 5.43893e-16 4.49462e-16 3.74218e-16 3.14368e-16 2.66518e-16
+ 2.27944e-16 1.96603e-16 1.71046e-16 1.50312e-16 1.33849e-16 1.21425e-16 1.1304e-16 1.088e-16 5.58393e-15 5.29696e-15 4.78887e-15 4.16097e-15
+ 3.51143e-15 2.90997e-15 2.39201e-15 1.96592e-15 1.62423e-15 1.35301e-15 1.13764e-15 9.65431e-16 8.26446e-16 7.13307e-16 6.20813e-16 5.45551e-16
+ 4.85579e-16 4.40151e-16 4.0938e-16 3.9378e-16 2.474e-14 2.34161e-14 2.10879e-14 1.82384e-14 1.53235e-14 1.26544e-14 1.03787e-14 8.52154e-15
+ 7.04045e-15 5.86845e-15 4.93876e-15 4.19517e-15 3.59425e-15 3.10418e-15 2.70261e-15 2.37493e-15 2.11295e-15 1.91378e-15 1.7784e-15 1.70959e-15
+ 1.1066e-13 1.04491e-13 9.37195e-14 8.06725e-14 6.74814e-14 5.55409e-14 4.5462e-14 3.72999e-14 3.08234e-14 2.5712e-14 2.166e-14 1.84169e-14
+ 1.57923e-14 1.36478e-14 1.18865e-14 1.04452e-14 9.28902e-15 8.40664e-15 7.80461e-15 7.49777e-15 4.94841e-13 4.66056e-13 4.16199e-13 3.56488e-13
+ 2.96871e-13 2.43552e-13 1.99003e-13 1.63197e-13 1.34919e-13 1.12649e-13 9.49997e-14 8.08599e-14 6.93982e-14 6.00141e-14 5.22888e-14 4.5949e-14
+ 4.0845e-14 3.69337e-14 3.42542e-14 3.28842e-14 2.2102e-12 2.07583e-12 1.84517e-12 1.57231e-12 1.30352e-12 1.06613e-12 8.69822e-13 7.13189e-13
+ 5.90011e-13 4.93162e-13 4.1639e-13 3.54806e-13 3.04797e-13 2.63767e-13 2.29908e-13 2.02039e-13 1.79519e-13 1.62184e-13 1.50255e-13 1.44134e-13
+ 9.85441e-12 9.22725e-12 8.16141e-12 6.91754e-12 5.70967e-12 4.6567e-12 3.79488e-12 3.112e-12 2.57694e-12 2.15668e-12 1.82332e-12 1.55549e-12
+ 1.33756e-12 1.15838e-12 1.01014e-12 8.87758e-13 7.88473e-13 7.11686e-13 6.58575e-13 6.31215e-13 4.38249e-11 4.09008e-11 3.59867e-11 3.03361e-11
+ 2.49321e-11 2.02839e-11 1.65179e-11 1.35529e-11 1.12365e-11 9.4179e-12 7.97361e-12 6.81098e-12 5.86292e-12 5.08162e-12 4.43365e-12 3.89701e-12
+ 3.45984e-12 3.11997e-12 2.88357e-12 2.76122e-12 1.94198e-10 1.80594e-10 1.5802e-10 1.32479e-10 1.08443e-10 8.80481e-11 7.16851e-11 5.88751e-11
+ 4.88889e-11 4.10458e-11 3.48063e-11 2.97722e-11 2.56575e-11 2.22584e-11 1.9432e-11 1.70836e-11 1.51622e-11 1.36599e-11 1.26083e-11 1.20612e-11
+ 8.56207e-10 7.93169e-10 6.90043e-10 5.75392e-10 4.69306e-10 3.80514e-10 3.09924e-10 2.54921e-10 2.12091e-10 1.78413e-10 1.51559e-10 1.29836e-10
+ 1.12034e-10 9.72904e-11 8.49983e-11 7.47507e-11 6.63273e-11 5.97004e-11 5.50281e-11 5.25819e-11 3.74879e-09 3.45849e-09 2.99117e-09 2.48142e-09
+ 2.01795e-09 1.63517e-09 1.33335e-09 1.099e-09 9.16503e-10 7.72706e-10 6.5771e-10 5.644e-10 4.87712e-10 4.2403e-10 3.70785e-10 3.26239e-10
+ 2.89443e-10 2.60297e-10 2.39579e-10 2.2865e-10 1.62575e-08 1.49333e-08 1.284e-08 1.06029e-08 8.60492e-09 6.97557e-09 5.69964e-09 4.7108e-09
+ 3.93972e-09 3.33026e-09 2.8411e-09 2.44279e-09 2.11438e-09 1.84086e-09 1.61147e-09 1.41885e-09 1.2589e-09 1.13125e-09 1.03966e-09 9.90909e-10
+ 6.95899e-08 6.36337e-08 5.44087e-08 4.47627e-08 3.62995e-08 2.9476e-08 2.41592e-08 2.00388e-08 1.68163e-08 1.42585e-08 1.21964e-08 1.05104e-08
+ 9.11527e-09 7.94943e-09 6.96843e-09 6.14139e-09 5.45076e-09 4.89502e-09 4.49189e-09 4.275e-09 2.92613e-07 2.66361e-07 2.26627e-07 1.86019e-07
+ 1.5099e-07 1.23012e-07 1.01265e-07 8.43739e-08 7.10998e-08 6.05049e-08 5.19177e-08 4.48628e-08 3.89999e-08 3.4082e-08 2.99283e-08 2.64113e-08
+ 2.34565e-08 2.10568e-08 1.9294e-08 1.83328e-08 1.20072e-06 1.08836e-06 9.22632e-07 7.57148e-07 6.16563e-07 5.04982e-07 4.18181e-07 3.50409e-07
+ 2.96779e-07 2.53668e-07 2.18497e-07 1.89434e-07 1.6516e-07 1.44705e-07 1.27355e-07 1.12592e-07 1.00106e-07 8.98623e-08 8.22245e-08 7.79894e-08
+ 4.76424e-06 4.3034e-06 3.64285e-06 2.99807e-06 2.4568e-06 2.02795e-06 1.6926e-06 1.42846e-06 1.21743e-06 1.04623e-06 9.0543e-07 7.88253e-07
+ 6.89776e-07 6.06336e-07 5.35185e-07 4.74305e-07 4.22435e-07 3.7939e-07 3.46729e-07 3.28224e-07 1.80392e-05 1.62643e-05 1.37987e-05 1.144e-05
+ 9.4721e-06 7.90528e-06 6.66625e-06 5.67728e-06 4.87675e-06 4.21965e-06 3.67362e-06 3.21512e-06 2.82675e-06 2.49537e-06 2.21094e-06 1.96594e-06
+ 1.75543e-06 1.57848e-06 1.44138e-06 1.36148e-06 6.39192e-05 5.77085e-05 4.93653e-05 4.14951e-05 3.4904e-05 2.9577e-05 2.52814e-05 2.17848e-05
+ 1.89033e-05 1.65009e-05 1.44773e-05 1.27579e-05 1.12864e-05 1.00191e-05 8.92208e-06 7.96916e-06 7.14234e-06 6.43694e-06 5.87657e-06 5.53727e-06
+ 0.000205656 0.000187039 0.0001629 0.000140141 0.00012066 0.000104414 9.08834e-05 7.95406e-05 6.99522e-05 6.17812e-05 5.47682e-05 4.87114e-05
+ 4.34524e-05 3.88648e-05 3.48472e-05 3.13184e-05 2.82192e-05 2.55291e-05 2.33256e-05 2.19184e-05 0.000571907 0.000529717 0.000477255 0.000426131
+ 0.000379662 0.000338394 0.000302037 0.000270076 0.000241964 0.000217194 0.000195322 0.000175962 0.000158787 0.000143517 0.000129913 0.000117774
+ 0.00010694 9.73373e-05 8.9165e-05 8.35209e-05
+ </DataArray>
+ <DataArray type="Float32" Name="X^H2O_liquid" 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 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 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 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 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 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
+ 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.999998 0.999998 0.999999 0.999999 0.999999 0.999999 0.999999 0.999999
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 0.999992 0.999993 0.999994 0.999995 0.999996 0.999997 0.999997 0.999998 0.999998 0.999998 0.999999 0.999999
+ 0.999999 0.999999 0.999999 0.999999 0.999999 0.999999 0.999999 0.999999 0.999971 0.999974 0.999978 0.999982
+ 0.999985 0.999987 0.999989 0.999991 0.999992 0.999993 0.999994 0.999995 0.999995 0.999996 0.999996 0.999997
+ 0.999997 0.999997 0.999998 0.999998 0.999897 0.999907 0.999921 0.999933 0.999944 0.999952 0.999959 0.999965
+ 0.99997 0.999973 0.999977 0.999979 0.999982 0.999984 0.999986 0.999987 0.999988 0.99999 0.999991 0.999991
+ 0.999669 0.999699 0.999738 0.999775 0.999806 0.999832 0.999854 0.999872 0.999888 0.999901 0.999912 0.999922
+ 0.99993 0.999938 0.999944 0.99995 0.999955 0.999959 0.999963 0.999965 0.999081 0.999149 0.999233 0.999315
+ 0.99939 0.999456 0.999515 0.999566 0.999611 0.999651 0.999686 0.999717 0.999745 0.999769 0.999791 0.999811
+ 0.999828 0.999844 0.999857 0.999866
+ </DataArray>
+ <DataArray type="Float32" Name="X^Air_liquid" NumberOfComponents="1" format="ascii">
+ 2.42482e-15 2.30548e-15 2.09274e-15 1.82723e-15 1.54939e-15 1.28905e-15 1.06242e-15 8.74326e-16 7.22526e-16 6.01568e-16 5.05357e-16 4.28436e-16
+ 3.66427e-16 3.16047e-16 2.74962e-16 2.41632e-16 2.15166e-16 1.95195e-16 1.81716e-16 1.749e-16 8.97636e-15 8.51504e-15 7.69826e-15 6.6889e-15
+ 5.64474e-15 4.67787e-15 3.84523e-15 3.16028e-15 2.61101e-15 2.17501e-15 1.82879e-15 1.55196e-15 1.32854e-15 1.14666e-15 9.97978e-16 8.76992e-16
+ 7.80584e-16 7.07558e-16 6.58091e-16 6.33014e-16 3.97703e-14 3.76422e-14 3.38995e-14 2.93188e-14 2.4633e-14 2.03423e-14 1.66841e-14 1.36987e-14
+ 1.13178e-14 9.43373e-15 7.93923e-15 6.74387e-15 5.77788e-15 4.99007e-15 4.34453e-15 3.81778e-15 3.39664e-15 3.07647e-15 2.85884e-15 2.74822e-15
+ 1.7789e-13 1.67973e-13 1.50657e-13 1.29684e-13 1.08479e-13 8.92838e-14 7.30817e-14 5.99608e-14 4.95497e-14 4.13329e-14 3.48191e-14 2.96058e-14
+ 2.53867e-14 2.19393e-14 1.91079e-14 1.67911e-14 1.49324e-14 1.3514e-14 1.25462e-14 1.20529e-14 7.95473e-13 7.492e-13 6.69053e-13 5.73066e-13
+ 4.7723e-13 3.91518e-13 3.19903e-13 2.62344e-13 2.16887e-13 1.81088e-13 1.52715e-13 1.29985e-13 1.1156e-13 9.64747e-14 8.40559e-14 7.38645e-14
+ 6.56596e-14 5.93721e-14 5.50647e-14 5.28624e-14 3.55297e-12 3.33697e-12 2.96617e-12 2.52755e-12 2.09546e-12 1.71385e-12 1.39827e-12 1.14647e-12
+ 9.48462e-13 7.92774e-13 6.6936e-13 5.70363e-13 4.89971e-13 4.24014e-13 3.69584e-13 3.24784e-13 2.88582e-13 2.60717e-13 2.4154e-13 2.31701e-13
+ 1.58413e-11 1.48331e-11 1.31197e-11 1.11202e-11 9.17848e-12 7.48581e-12 6.10039e-12 5.00264e-12 4.14251e-12 3.46694e-12 2.93105e-12 2.5005e-12
+ 2.15018e-12 1.86213e-12 1.62383e-12 1.4271e-12 1.2675e-12 1.14406e-12 1.05868e-12 1.0147e-12 7.045e-11 6.57493e-11 5.78499e-11 4.87664e-11
+ 4.00792e-11 3.26071e-11 2.65531e-11 2.17867e-11 1.80631e-11 1.51396e-11 1.28178e-11 1.09489e-11 9.42484e-12 8.16888e-12 7.12724e-12 6.26457e-12
+ 5.56181e-12 5.01546e-12 4.63543e-12 4.43876e-12 3.12179e-10 2.90311e-10 2.54023e-10 2.12964e-10 1.74325e-10 1.4154e-10 1.15236e-10 9.46436e-11
+ 7.85906e-11 6.59825e-11 5.59523e-11 4.78599e-11 4.12452e-11 3.5781e-11 3.12375e-11 2.74625e-11 2.43737e-11 2.19587e-11 2.02683e-11 1.93887e-11
+ 1.37638e-09 1.27504e-09 1.10927e-09 9.24961e-10 7.54426e-10 6.11688e-10 4.98213e-10 4.09793e-10 3.40943e-10 2.86805e-10 2.43636e-10 2.08716e-10
+ 1.80098e-10 1.56398e-10 1.36638e-10 1.20164e-10 1.06623e-10 9.59705e-11 8.84596e-11 8.45271e-11 6.0263e-09 5.55964e-09 4.80841e-09 3.98897e-09
+ 3.24391e-09 2.62859e-09 2.14341e-09 1.76668e-09 1.47331e-09 1.24215e-09 1.05729e-09 9.07291e-10 7.84013e-10 6.81643e-10 5.96049e-10 5.24441e-10
+ 4.6529e-10 4.18437e-10 3.85131e-10 3.67563e-10 2.61345e-08 2.40058e-08 2.06407e-08 1.70445e-08 1.38327e-08 1.12135e-08 9.16236e-09 7.57277e-09
+ 6.33323e-09 5.35351e-09 4.56717e-09 3.92687e-09 3.39894e-09 2.95925e-09 2.59049e-09 2.28085e-09 2.02372e-09 1.81853e-09 1.67128e-09 1.59292e-09
+ 1.11868e-07 1.02293e-07 8.74637e-08 7.19576e-08 5.83526e-08 4.73836e-08 3.88367e-08 3.22131e-08 2.70328e-08 2.2921e-08 1.96062e-08 1.68959e-08
+ 1.46531e-08 1.2779e-08 1.1202e-08 9.87248e-09 8.76228e-09 7.86891e-09 7.22085e-09 6.87221e-09 4.70385e-07 4.28184e-07 3.64311e-07 2.99031e-07
+ 2.42721e-07 1.97745e-07 1.62787e-07 1.35634e-07 1.14295e-07 9.72637e-08 8.34595e-08 7.21184e-08 6.26936e-08 5.47879e-08 4.81108e-08 4.2457e-08
+ 3.7707e-08 3.38496e-08 3.10157e-08 2.94706e-08 1.9302e-06 1.74958e-06 1.48316e-06 1.21714e-06 9.91146e-07 8.11774e-07 6.72239e-07 5.63294e-07
+ 4.77082e-07 4.0778e-07 3.51241e-07 3.04521e-07 2.65499e-07 2.32619e-07 2.04727e-07 1.80995e-07 1.60924e-07 1.44457e-07 1.32179e-07 1.25371e-07
+ 7.65866e-06 6.91783e-06 5.85599e-06 4.81948e-06 3.94938e-06 3.26e-06 2.7209e-06 2.2963e-06 1.95705e-06 1.68185e-06 1.45551e-06 1.26714e-06
+ 1.10884e-06 9.74704e-07 8.60327e-07 7.62461e-07 6.79079e-07 6.09882e-07 5.57379e-07 5.27631e-07 2.89983e-05 2.61451e-05 2.21817e-05 1.839e-05
+ 1.52266e-05 1.27079e-05 1.07162e-05 9.12638e-06 7.83951e-06 6.78322e-06 5.90545e-06 5.16841e-06 4.54409e-06 4.01138e-06 3.55416e-06 3.16031e-06
+ 2.82192e-06 2.53745e-06 2.31706e-06 2.18862e-06 0.000102748 9.2765e-05 7.9354e-05 6.67031e-05 5.61081e-05 4.75451e-05 4.06401e-05 3.50193e-05
+ 3.03873e-05 2.65254e-05 2.32725e-05 2.05086e-05 1.81431e-05 1.61059e-05 1.43425e-05 1.28106e-05 1.14815e-05 1.03476e-05 9.44674e-06 8.90131e-06
+ 0.000330557 0.000300637 0.000261841 0.000225262 0.00019395 0.000167838 0.00014609 0.000127858 0.000112446 9.93116e-05 8.80387e-05 7.83028e-05
+ 6.98493e-05 6.2475e-05 5.60169e-05 5.03444e-05 4.53625e-05 4.10383e-05 3.74961e-05 3.5234e-05 0.00091904 0.000851263 0.00076698 0.000684843
+ 0.000610178 0.000543867 0.000485446 0.000434085 0.000388907 0.000349101 0.000313949 0.000282835 0.00025523 0.000230688 0.000208823 0.000189312
+ 0.000171898 0.000156464 0.000143328 0.000134256
+ </DataArray>
+ <DataArray type="Float32" Name="process rank" 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 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 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 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 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 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 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 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>
+ </CellData>
+ <Points>
+ <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+ 0 0 0 0.05 0 0 0 0.05 0 0.05 0.05 0
+ 0.1 0 0 0.1 0.05 0 0.15 0 0 0.15 0.05 0
+ 0.2 0 0 0.2 0.05 0 0.25 0 0 0.25 0.05 0
+ 0.3 0 0 0.3 0.05 0 0.35 0 0 0.35 0.05 0
+ 0.4 0 0 0.4 0.05 0 0.45 0 0 0.45 0.05 0
+ 0.5 0 0 0.5 0.05 0 0.55 0 0 0.55 0.05 0
+ 0.6 0 0 0.6 0.05 0 0.65 0 0 0.65 0.05 0
+ 0.7 0 0 0.7 0.05 0 0.75 0 0 0.75 0.05 0
+ 0.8 0 0 0.8 0.05 0 0.85 0 0 0.85 0.05 0
+ 0.9 0 0 0.9 0.05 0 0.95 0 0 0.95 0.05 0
+ 1 0 0 1 0.05 0 0 0.1 0 0.05 0.1 0
+ 0.1 0.1 0 0.15 0.1 0 0.2 0.1 0 0.25 0.1 0
+ 0.3 0.1 0 0.35 0.1 0 0.4 0.1 0 0.45 0.1 0
+ 0.5 0.1 0 0.55 0.1 0 0.6 0.1 0 0.65 0.1 0
+ 0.7 0.1 0 0.75 0.1 0 0.8 0.1 0 0.85 0.1 0
+ 0.9 0.1 0 0.95 0.1 0 1 0.1 0 0 0.15 0
+ 0.05 0.15 0 0.1 0.15 0 0.15 0.15 0 0.2 0.15 0
+ 0.25 0.15 0 0.3 0.15 0 0.35 0.15 0 0.4 0.15 0
+ 0.45 0.15 0 0.5 0.15 0 0.55 0.15 0 0.6 0.15 0
+ 0.65 0.15 0 0.7 0.15 0 0.75 0.15 0 0.8 0.15 0
+ 0.85 0.15 0 0.9 0.15 0 0.95 0.15 0 1 0.15 0
+ 0 0.2 0 0.05 0.2 0 0.1 0.2 0 0.15 0.2 0
+ 0.2 0.2 0 0.25 0.2 0 0.3 0.2 0 0.35 0.2 0
+ 0.4 0.2 0 0.45 0.2 0 0.5 0.2 0 0.55 0.2 0
+ 0.6 0.2 0 0.65 0.2 0 0.7 0.2 0 0.75 0.2 0
+ 0.8 0.2 0 0.85 0.2 0 0.9 0.2 0 0.95 0.2 0
+ 1 0.2 0 0 0.25 0 0.05 0.25 0 0.1 0.25 0
+ 0.15 0.25 0 0.2 0.25 0 0.25 0.25 0 0.3 0.25 0
+ 0.35 0.25 0 0.4 0.25 0 0.45 0.25 0 0.5 0.25 0
+ 0.55 0.25 0 0.6 0.25 0 0.65 0.25 0 0.7 0.25 0
+ 0.75 0.25 0 0.8 0.25 0 0.85 0.25 0 0.9 0.25 0
+ 0.95 0.25 0 1 0.25 0 0 0.3 0 0.05 0.3 0
+ 0.1 0.3 0 0.15 0.3 0 0.2 0.3 0 0.25 0.3 0
+ 0.3 0.3 0 0.35 0.3 0 0.4 0.3 0 0.45 0.3 0
+ 0.5 0.3 0 0.55 0.3 0 0.6 0.3 0 0.65 0.3 0
+ 0.7 0.3 0 0.75 0.3 0 0.8 0.3 0 0.85 0.3 0
+ 0.9 0.3 0 0.95 0.3 0 1 0.3 0 0 0.35 0
+ 0.05 0.35 0 0.1 0.35 0 0.15 0.35 0 0.2 0.35 0
+ 0.25 0.35 0 0.3 0.35 0 0.35 0.35 0 0.4 0.35 0
+ 0.45 0.35 0 0.5 0.35 0 0.55 0.35 0 0.6 0.35 0
+ 0.65 0.35 0 0.7 0.35 0 0.75 0.35 0 0.8 0.35 0
+ 0.85 0.35 0 0.9 0.35 0 0.95 0.35 0 1 0.35 0
+ 0 0.4 0 0.05 0.4 0 0.1 0.4 0 0.15 0.4 0
+ 0.2 0.4 0 0.25 0.4 0 0.3 0.4 0 0.35 0.4 0
+ 0.4 0.4 0 0.45 0.4 0 0.5 0.4 0 0.55 0.4 0
+ 0.6 0.4 0 0.65 0.4 0 0.7 0.4 0 0.75 0.4 0
+ 0.8 0.4 0 0.85 0.4 0 0.9 0.4 0 0.95 0.4 0
+ 1 0.4 0 0 0.45 0 0.05 0.45 0 0.1 0.45 0
+ 0.15 0.45 0 0.2 0.45 0 0.25 0.45 0 0.3 0.45 0
+ 0.35 0.45 0 0.4 0.45 0 0.45 0.45 0 0.5 0.45 0
+ 0.55 0.45 0 0.6 0.45 0 0.65 0.45 0 0.7 0.45 0
+ 0.75 0.45 0 0.8 0.45 0 0.85 0.45 0 0.9 0.45 0
+ 0.95 0.45 0 1 0.45 0 0 0.5 0 0.05 0.5 0
+ 0.1 0.5 0 0.15 0.5 0 0.2 0.5 0 0.25 0.5 0
+ 0.3 0.5 0 0.35 0.5 0 0.4 0.5 0 0.45 0.5 0
+ 0.5 0.5 0 0.55 0.5 0 0.6 0.5 0 0.65 0.5 0
+ 0.7 0.5 0 0.75 0.5 0 0.8 0.5 0 0.85 0.5 0
+ 0.9 0.5 0 0.95 0.5 0 1 0.5 0 0 0.55 0
+ 0.05 0.55 0 0.1 0.55 0 0.15 0.55 0 0.2 0.55 0
+ 0.25 0.55 0 0.3 0.55 0 0.35 0.55 0 0.4 0.55 0
+ 0.45 0.55 0 0.5 0.55 0 0.55 0.55 0 0.6 0.55 0
+ 0.65 0.55 0 0.7 0.55 0 0.75 0.55 0 0.8 0.55 0
+ 0.85 0.55 0 0.9 0.55 0 0.95 0.55 0 1 0.55 0
+ 0 0.6 0 0.05 0.6 0 0.1 0.6 0 0.15 0.6 0
+ 0.2 0.6 0 0.25 0.6 0 0.3 0.6 0 0.35 0.6 0
+ 0.4 0.6 0 0.45 0.6 0 0.5 0.6 0 0.55 0.6 0
+ 0.6 0.6 0 0.65 0.6 0 0.7 0.6 0 0.75 0.6 0
+ 0.8 0.6 0 0.85 0.6 0 0.9 0.6 0 0.95 0.6 0
+ 1 0.6 0 0 0.65 0 0.05 0.65 0 0.1 0.65 0
+ 0.15 0.65 0 0.2 0.65 0 0.25 0.65 0 0.3 0.65 0
+ 0.35 0.65 0 0.4 0.65 0 0.45 0.65 0 0.5 0.65 0
+ 0.55 0.65 0 0.6 0.65 0 0.65 0.65 0 0.7 0.65 0
+ 0.75 0.65 0 0.8 0.65 0 0.85 0.65 0 0.9 0.65 0
+ 0.95 0.65 0 1 0.65 0 0 0.7 0 0.05 0.7 0
+ 0.1 0.7 0 0.15 0.7 0 0.2 0.7 0 0.25 0.7 0
+ 0.3 0.7 0 0.35 0.7 0 0.4 0.7 0 0.45 0.7 0
+ 0.5 0.7 0 0.55 0.7 0 0.6 0.7 0 0.65 0.7 0
+ 0.7 0.7 0 0.75 0.7 0 0.8 0.7 0 0.85 0.7 0
+ 0.9 0.7 0 0.95 0.7 0 1 0.7 0 0 0.75 0
+ 0.05 0.75 0 0.1 0.75 0 0.15 0.75 0 0.2 0.75 0
+ 0.25 0.75 0 0.3 0.75 0 0.35 0.75 0 0.4 0.75 0
+ 0.45 0.75 0 0.5 0.75 0 0.55 0.75 0 0.6 0.75 0
+ 0.65 0.75 0 0.7 0.75 0 0.75 0.75 0 0.8 0.75 0
+ 0.85 0.75 0 0.9 0.75 0 0.95 0.75 0 1 0.75 0
+ 0 0.8 0 0.05 0.8 0 0.1 0.8 0 0.15 0.8 0
+ 0.2 0.8 0 0.25 0.8 0 0.3 0.8 0 0.35 0.8 0
+ 0.4 0.8 0 0.45 0.8 0 0.5 0.8 0 0.55 0.8 0
+ 0.6 0.8 0 0.65 0.8 0 0.7 0.8 0 0.75 0.8 0
+ 0.8 0.8 0 0.85 0.8 0 0.9 0.8 0 0.95 0.8 0
+ 1 0.8 0 0 0.85 0 0.05 0.85 0 0.1 0.85 0
+ 0.15 0.85 0 0.2 0.85 0 0.25 0.85 0 0.3 0.85 0
+ 0.35 0.85 0 0.4 0.85 0 0.45 0.85 0 0.5 0.85 0
+ 0.55 0.85 0 0.6 0.85 0 0.65 0.85 0 0.7 0.85 0
+ 0.75 0.85 0 0.8 0.85 0 0.85 0.85 0 0.9 0.85 0
+ 0.95 0.85 0 1 0.85 0 0 0.9 0 0.05 0.9 0
+ 0.1 0.9 0 0.15 0.9 0 0.2 0.9 0 0.25 0.9 0
+ 0.3 0.9 0 0.35 0.9 0 0.4 0.9 0 0.45 0.9 0
+ 0.5 0.9 0 0.55 0.9 0 0.6 0.9 0 0.65 0.9 0
+ 0.7 0.9 0 0.75 0.9 0 0.8 0.9 0 0.85 0.9 0
+ 0.9 0.9 0 0.95 0.9 0 1 0.9 0 0 0.95 0
+ 0.05 0.95 0 0.1 0.95 0 0.15 0.95 0 0.2 0.95 0
+ 0.25 0.95 0 0.3 0.95 0 0.35 0.95 0 0.4 0.95 0
+ 0.45 0.95 0 0.5 0.95 0 0.55 0.95 0 0.6 0.95 0
+ 0.65 0.95 0 0.7 0.95 0 0.75 0.95 0 0.8 0.95 0
+ 0.85 0.95 0 0.9 0.95 0 0.95 0.95 0 1 0.95 0
+ 0 1 0 0.05 1 0 0.1 1 0 0.15 1 0
+ 0.2 1 0 0.25 1 0 0.3 1 0 0.35 1 0
+ 0.4 1 0 0.45 1 0 0.5 1 0 0.55 1 0
+ 0.6 1 0 0.65 1 0 0.7 1 0 0.75 1 0
+ 0.8 1 0 0.85 1 0 0.9 1 0 0.95 1 0
+ 1 1 0
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 1 4 5 3 4 6 7 5
+ 6 8 9 7 8 10 11 9 10 12 13 11
+ 12 14 15 13 14 16 17 15 16 18 19 17
+ 18 20 21 19 20 22 23 21 22 24 25 23
+ 24 26 27 25 26 28 29 27 28 30 31 29
+ 30 32 33 31 32 34 35 33 34 36 37 35
+ 36 38 39 37 38 40 41 39 2 3 43 42
+ 3 5 44 43 5 7 45 44 7 9 46 45
+ 9 11 47 46 11 13 48 47 13 15 49 48
+ 15 17 50 49 17 19 51 50 19 21 52 51
+ 21 23 53 52 23 25 54 53 25 27 55 54
+ 27 29 56 55 29 31 57 56 31 33 58 57
+ 33 35 59 58 35 37 60 59 37 39 61 60
+ 39 41 62 61 42 43 64 63 43 44 65 64
+ 44 45 66 65 45 46 67 66 46 47 68 67
+ 47 48 69 68 48 49 70 69 49 50 71 70
+ 50 51 72 71 51 52 73 72 52 53 74 73
+ 53 54 75 74 54 55 76 75 55 56 77 76
+ 56 57 78 77 57 58 79 78 58 59 80 79
+ 59 60 81 80 60 61 82 81 61 62 83 82
+ 63 64 85 84 64 65 86 85 65 66 87 86
+ 66 67 88 87 67 68 89 88 68 69 90 89
+ 69 70 91 90 70 71 92 91 71 72 93 92
+ 72 73 94 93 73 74 95 94 74 75 96 95
+ 75 76 97 96 76 77 98 97 77 78 99 98
+ 78 79 100 99 79 80 101 100 80 81 102 101
+ 81 82 103 102 82 83 104 103 84 85 106 105
+ 85 86 107 106 86 87 108 107 87 88 109 108
+ 88 89 110 109 89 90 111 110 90 91 112 111
+ 91 92 113 112 92 93 114 113 93 94 115 114
+ 94 95 116 115 95 96 117 116 96 97 118 117
+ 97 98 119 118 98 99 120 119 99 100 121 120
+ 100 101 122 121 101 102 123 122 102 103 124 123
+ 103 104 125 124 105 106 127 126 106 107 128 127
+ 107 108 129 128 108 109 130 129 109 110 131 130
+ 110 111 132 131 111 112 133 132 112 113 134 133
+ 113 114 135 134 114 115 136 135 115 116 137 136
+ 116 117 138 137 117 118 139 138 118 119 140 139
+ 119 120 141 140 120 121 142 141 121 122 143 142
+ 122 123 144 143 123 124 145 144 124 125 146 145
+ 126 127 148 147 127 128 149 148 128 129 150 149
+ 129 130 151 150 130 131 152 151 131 132 153 152
+ 132 133 154 153 133 134 155 154 134 135 156 155
+ 135 136 157 156 136 137 158 157 137 138 159 158
+ 138 139 160 159 139 140 161 160 140 141 162 161
+ 141 142 163 162 142 143 164 163 143 144 165 164
+ 144 145 166 165 145 146 167 166 147 148 169 168
+ 148 149 170 169 149 150 171 170 150 151 172 171
+ 151 152 173 172 152 153 174 173 153 154 175 174
+ 154 155 176 175 155 156 177 176 156 157 178 177
+ 157 158 179 178 158 159 180 179 159 160 181 180
+ 160 161 182 181 161 162 183 182 162 163 184 183
+ 163 164 185 184 164 165 186 185 165 166 187 186
+ 166 167 188 187 168 169 190 189 169 170 191 190
+ 170 171 192 191 171 172 193 192 172 173 194 193
+ 173 174 195 194 174 175 196 195 175 176 197 196
+ 176 177 198 197 177 178 199 198 178 179 200 199
+ 179 180 201 200 180 181 202 201 181 182 203 202
+ 182 183 204 203 183 184 205 204 184 185 206 205
+ 185 186 207 206 186 187 208 207 187 188 209 208
+ 189 190 211 210 190 191 212 211 191 192 213 212
+ 192 193 214 213 193 194 215 214 194 195 216 215
+ 195 196 217 216 196 197 218 217 197 198 219 218
+ 198 199 220 219 199 200 221 220 200 201 222 221
+ 201 202 223 222 202 203 224 223 203 204 225 224
+ 204 205 226 225 205 206 227 226 206 207 228 227
+ 207 208 229 228 208 209 230 229 210 211 232 231
+ 211 212 233 232 212 213 234 233 213 214 235 234
+ 214 215 236 235 215 216 237 236 216 217 238 237
+ 217 218 239 238 218 219 240 239 219 220 241 240
+ 220 221 242 241 221 222 243 242 222 223 244 243
+ 223 224 245 244 224 225 246 245 225 226 247 246
+ 226 227 248 247 227 228 249 248 228 229 250 249
+ 229 230 251 250 231 232 253 252 232 233 254 253
+ 233 234 255 254 234 235 256 255 235 236 257 256
+ 236 237 258 257 237 238 259 258 238 239 260 259
+ 239 240 261 260 240 241 262 261 241 242 263 262
+ 242 243 264 263 243 244 265 264 244 245 266 265
+ 245 246 267 266 246 247 268 267 247 248 269 268
+ 248 249 270 269 249 250 271 270 250 251 272 271
+ 252 253 274 273 253 254 275 274 254 255 276 275
+ 255 256 277 276 256 257 278 277 257 258 279 278
+ 258 259 280 279 259 260 281 280 260 261 282 281
+ 261 262 283 282 262 263 284 283 263 264 285 284
+ 264 265 286 285 265 266 287 286 266 267 288 287
+ 267 268 289 288 268 269 290 289 269 270 291 290
+ 270 271 292 291 271 272 293 292 273 274 295 294
+ 274 275 296 295 275 276 297 296 276 277 298 297
+ 277 278 299 298 278 279 300 299 279 280 301 300
+ 280 281 302 301 281 282 303 302 282 283 304 303
+ 283 284 305 304 284 285 306 305 285 286 307 306
+ 286 287 308 307 287 288 309 308 288 289 310 309
+ 289 290 311 310 290 291 312 311 291 292 313 312
+ 292 293 314 313 294 295 316 315 295 296 317 316
+ 296 297 318 317 297 298 319 318 298 299 320 319
+ 299 300 321 320 300 301 322 321 301 302 323 322
+ 302 303 324 323 303 304 325 324 304 305 326 325
+ 305 306 327 326 306 307 328 327 307 308 329 328
+ 308 309 330 329 309 310 331 330 310 311 332 331
+ 311 312 333 332 312 313 334 333 313 314 335 334
+ 315 316 337 336 316 317 338 337 317 318 339 338
+ 318 319 340 339 319 320 341 340 320 321 342 341
+ 321 322 343 342 322 323 344 343 323 324 345 344
+ 324 325 346 345 325 326 347 346 326 327 348 347
+ 327 328 349 348 328 329 350 349 329 330 351 350
+ 330 331 352 351 331 332 353 352 332 333 354 353
+ 333 334 355 354 334 335 356 355 336 337 358 357
+ 337 338 359 358 338 339 360 359 339 340 361 360
+ 340 341 362 361 341 342 363 362 342 343 364 363
+ 343 344 365 364 344 345 366 365 345 346 367 366
+ 346 347 368 367 347 348 369 368 348 349 370 369
+ 349 350 371 370 350 351 372 371 351 352 373 372
+ 352 353 374 373 353 354 375 374 354 355 376 375
+ 355 356 377 376 357 358 379 378 358 359 380 379
+ 359 360 381 380 360 361 382 381 361 362 383 382
+ 362 363 384 383 363 364 385 384 364 365 386 385
+ 365 366 387 386 366 367 388 387 367 368 389 388
+ 368 369 390 389 369 370 391 390 370 371 392 391
+ 371 372 393 392 372 373 394 393 373 374 395 394
+ 374 375 396 395 375 376 397 396 376 377 398 397
+ 378 379 400 399 379 380 401 400 380 381 402 401
+ 381 382 403 402 382 383 404 403 383 384 405 404
+ 384 385 406 405 385 386 407 406 386 387 408 407
+ 387 388 409 408 388 389 410 409 389 390 411 410
+ 390 391 412 411 391 392 413 412 392 393 414 413
+ 393 394 415 414 394 395 416 415 395 396 417 416
+ 396 397 418 417 397 398 419 418 399 400 421 420
+ 400 401 422 421 401 402 423 422 402 403 424 423
+ 403 404 425 424 404 405 426 425 405 406 427 426
+ 406 407 428 427 407 408 429 428 408 409 430 429
+ 409 410 431 430 410 411 432 431 411 412 433 432
+ 412 413 434 433 413 414 435 434 414 415 436 435
+ 415 416 437 436 416 417 438 437 417 418 439 438
+ 418 419 440 439
+ </DataArray>
+ <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+ 4 8 12 16 20 24 28 32 36 40 44 48
+ 52 56 60 64 68 72 76 80 84 88 92 96
+ 100 104 108 112 116 120 124 128 132 136 140 144
+ 148 152 156 160 164 168 172 176 180 184 188 192
+ 196 200 204 208 212 216 220 224 228 232 236 240
+ 244 248 252 256 260 264 268 272 276 280 284 288
+ 292 296 300 304 308 312 316 320 324 328 332 336
+ 340 344 348 352 356 360 364 368 372 376 380 384
+ 388 392 396 400 404 408 412 416 420 424 428 432
+ 436 440 444 448 452 456 460 464 468 472 476 480
+ 484 488 492 496 500 504 508 512 516 520 524 528
+ 532 536 540 544 548 552 556 560 564 568 572 576
+ 580 584 588 592 596 600 604 608 612 616 620 624
+ 628 632 636 640 644 648 652 656 660 664 668 672
+ 676 680 684 688 692 696 700 704 708 712 716 720
+ 724 728 732 736 740 744 748 752 756 760 764 768
+ 772 776 780 784 788 792 796 800 804 808 812 816
+ 820 824 828 832 836 840 844 848 852 856 860 864
+ 868 872 876 880 884 888 892 896 900 904 908 912
+ 916 920 924 928 932 936 940 944 948 952 956 960
+ 964 968 972 976 980 984 988 992 996 1000 1004 1008
+ 1012 1016 1020 1024 1028 1032 1036 1040 1044 1048 1052 1056
+ 1060 1064 1068 1072 1076 1080 1084 1088 1092 1096 1100 1104
+ 1108 1112 1116 1120 1124 1128 1132 1136 1140 1144 1148 1152
+ 1156 1160 1164 1168 1172 1176 1180 1184 1188 1192 1196 1200
+ 1204 1208 1212 1216 1220 1224 1228 1232 1236 1240 1244 1248
+ 1252 1256 1260 1264 1268 1272 1276 1280 1284 1288 1292 1296
+ 1300 1304 1308 1312 1316 1320 1324 1328 1332 1336 1340 1344
+ 1348 1352 1356 1360 1364 1368 1372 1376 1380 1384 1388 1392
+ 1396 1400 1404 1408 1412 1416 1420 1424 1428 1432 1436 1440
+ 1444 1448 1452 1456 1460 1464 1468 1472 1476 1480 1484 1488
+ 1492 1496 1500 1504 1508 1512 1516 1520 1524 1528 1532 1536
+ 1540 1544 1548 1552 1556 1560 1564 1568 1572 1576 1580 1584
+ 1588 1592 1596 1600
+ </DataArray>
+ <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
diff --git a/test/references/test_md_darcy1p2c_stokes1p2c_comparison_stokes-reference.vtu b/test/references/test_md_darcy1p2c_stokes1p2c_comparison_stokes-reference.vtu
new file mode 100644
index 0000000000..711dd62a31
--- /dev/null
+++ b/test/references/test_md_darcy1p2c_stokes1p2c_comparison_stokes-reference.vtu
@@ -0,0 +1,760 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="400" NumberOfPoints="441">
+ <CellData Scalars="p" Vectors="velocity_liquid (m/s)">
+ <DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000
+ </DataArray>
+ <DataArray type="Float32" Name="rhoMolar" NumberOfComponents="1" format="ascii">
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5 55491.5
+ 55491.5 55491.5 55491.5 55491.5
+ </DataArray>
+ <DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
+ 1000.28 1000.25 1000.22 1000.18 1000.15 1000.11 1000.08 1000.05 1000.02 999.987 999.962 999.939
+ 999.917 999.898 999.88 999.864 999.85 999.836 999.825 999.816 1000.29 1000.29 1000.28 1000.26
+ 1000.25 1000.22 1000.2 1000.18 1000.15 1000.12 1000.1 1000.07 1000.05 1000.03 1000 999.982
+ 999.962 999.944 999.926 999.912 1000.3 1000.3 1000.29 1000.29 1000.28 1000.28 1000.27 1000.26
+ 1000.25 1000.23 1000.21 1000.2 1000.18 1000.16 1000.14 1000.12 1000.1 1000.08 1000.06 1000.04
+ 1000.3 1000.3 1000.3 1000.3 1000.29 1000.29 1000.29 1000.29 1000.28 1000.28 1000.27 1000.26
+ 1000.25 1000.24 1000.23 1000.21 1000.2 1000.18 1000.17 1000.15 1000.3 1000.3 1000.3 1000.3
+ 1000.3 1000.3 1000.29 1000.29 1000.29 1000.29 1000.29 1000.28 1000.28 1000.27 1000.27 1000.26
+ 1000.25 1000.24 1000.23 1000.22 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3
+ 1000.29 1000.29 1000.29 1000.29 1000.29 1000.28 1000.28 1000.28 1000.27 1000.27 1000.26 1000.25
+ 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.29 1000.29 1000.29
+ 1000.29 1000.29 1000.29 1000.28 1000.28 1000.28 1000.27 1000.27 1000.3 1000.3 1000.3 1000.3
+ 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.29 1000.29 1000.29 1000.29 1000.29 1000.29
+ 1000.29 1000.28 1000.28 1000.28 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3
+ 1000.3 1000.3 1000.3 1000.29 1000.29 1000.29 1000.29 1000.29 1000.29 1000.29 1000.28 1000.28
+ 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.29
+ 1000.29 1000.29 1000.29 1000.29 1000.29 1000.29 1000.28 1000.28 1000.3 1000.3 1000.3 1000.3
+ 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.29 1000.29 1000.29 1000.29 1000.29
+ 1000.29 1000.29 1000.28 1000.28 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3
+ 1000.3 1000.3 1000.3 1000.29 1000.29 1000.29 1000.29 1000.29 1000.29 1000.29 1000.28 1000.28
+ 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.29 1000.29
+ 1000.29 1000.29 1000.29 1000.29 1000.29 1000.28 1000.28 1000.28 1000.3 1000.3 1000.3 1000.3
+ 1000.3 1000.3 1000.3 1000.3 1000.3 1000.29 1000.29 1000.29 1000.29 1000.29 1000.29 1000.28
+ 1000.28 1000.28 1000.27 1000.27 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3
+ 1000.29 1000.29 1000.29 1000.29 1000.29 1000.28 1000.28 1000.28 1000.27 1000.27 1000.26 1000.25
+ 1000.3 1000.3 1000.3 1000.3 1000.3 1000.3 1000.29 1000.29 1000.29 1000.29 1000.29 1000.28
+ 1000.28 1000.27 1000.27 1000.26 1000.25 1000.24 1000.23 1000.22 1000.3 1000.3 1000.3 1000.3
+ 1000.29 1000.29 1000.29 1000.29 1000.28 1000.28 1000.27 1000.26 1000.25 1000.24 1000.23 1000.22
+ 1000.2 1000.19 1000.17 1000.16 1000.3 1000.3 1000.29 1000.29 1000.29 1000.28 1000.27 1000.26
+ 1000.25 1000.24 1000.22 1000.2 1000.18 1000.17 1000.15 1000.13 1000.11 1000.09 1000.07 1000.05
+ 1000.3 1000.29 1000.28 1000.27 1000.26 1000.24 1000.22 1000.19 1000.17 1000.14 1000.12 1000.09
+ 1000.07 1000.04 1000.02 999.996 999.975 999.955 999.937 999.921 1000.29 1000.28 1000.26 1000.23
+ 1000.2 1000.16 1000.13 1000.09 1000.06 1000.03 1000 999.974 999.95 999.928 999.908 999.889
+ 999.872 999.857 999.844 999.833
+ </DataArray>
+ <DataArray type="Float32" Name="X^H2O_liquid" NumberOfComponents="1" format="ascii">
+ 0.998445 0.998519 0.998603 0.998694 0.998789 0.998883 0.998973 0.999058 0.999137 0.999211 0.999278 0.999339
+ 0.999395 0.999447 0.999493 0.999536 0.999575 0.99961 0.999641 0.999665 0.9984 0.998417 0.998444 0.998482
+ 0.998529 0.998584 0.998645 0.998711 0.998779 0.998847 0.998915 0.998982 0.999047 0.999109 0.999168 0.999224
+ 0.999277 0.999326 0.999372 0.999411 0.998394 0.998397 0.998402 0.99841 0.998424 0.998442 0.998466 0.998495
+ 0.998529 0.998568 0.99861 0.998656 0.998705 0.998755 0.998807 0.998859 0.998912 0.998964 0.999015 0.999061
+ 0.998394 0.998394 0.998395 0.998396 0.998399 0.998404 0.998411 0.99842 0.998433 0.998449 0.998468 0.998491
+ 0.998517 0.998546 0.998578 0.998613 0.99865 0.998688 0.998729 0.998768 0.998393 0.998393 0.998394 0.998394
+ 0.998395 0.998396 0.998398 0.998401 0.998406 0.998412 0.99842 0.99843 0.998443 0.998457 0.998474 0.998494
+ 0.998516 0.99854 0.998567 0.998593 0.998393 0.998393 0.998393 0.998394 0.998394 0.998394 0.998395 0.998397
+ 0.998398 0.998401 0.998405 0.99841 0.998416 0.998424 0.998434 0.998445 0.998458 0.998473 0.99849 0.998508
+ 0.998393 0.998393 0.998393 0.998394 0.998394 0.998394 0.998394 0.998395 0.998396 0.998398 0.9984 0.998403
+ 0.998407 0.998412 0.998418 0.998425 0.998434 0.998444 0.998455 0.998467 0.998393 0.998393 0.998393 0.998393
+ 0.998394 0.998394 0.998394 0.998394 0.998395 0.998396 0.998398 0.9984 0.998403 0.998406 0.99841 0.998416
+ 0.998422 0.99843 0.998438 0.998448 0.998393 0.998393 0.998393 0.998393 0.998394 0.998394 0.998394 0.998394
+ 0.998395 0.998396 0.998397 0.998399 0.998401 0.998404 0.998407 0.998411 0.998417 0.998423 0.99843 0.998438
+ 0.998393 0.998393 0.998393 0.998393 0.998394 0.998394 0.998394 0.998394 0.998395 0.998395 0.998397 0.998398
+ 0.9984 0.998402 0.998406 0.99841 0.998414 0.99842 0.998427 0.998434 0.998393 0.998393 0.998393 0.998393
+ 0.998394 0.998394 0.998394 0.998394 0.998395 0.998395 0.998397 0.998398 0.9984 0.998402 0.998406 0.99841
+ 0.998414 0.99842 0.998427 0.998434 0.998393 0.998393 0.998393 0.998393 0.998394 0.998394 0.998394 0.998394
+ 0.998395 0.998396 0.998397 0.998399 0.998401 0.998404 0.998407 0.998411 0.998417 0.998423 0.99843 0.998438
+ 0.998393 0.998393 0.998393 0.998393 0.998394 0.998394 0.998394 0.998394 0.998395 0.998396 0.998398 0.9984
+ 0.998403 0.998406 0.99841 0.998416 0.998422 0.99843 0.998438 0.998448 0.998393 0.998393 0.998393 0.998394
+ 0.998394 0.998394 0.998394 0.998395 0.998396 0.998398 0.9984 0.998403 0.998407 0.998412 0.998418 0.998425
+ 0.998433 0.998443 0.998455 0.998467 0.998393 0.998393 0.998393 0.998394 0.998394 0.998394 0.998395 0.998396
+ 0.998398 0.998401 0.998405 0.99841 0.998416 0.998424 0.998434 0.998445 0.998458 0.998473 0.99849 0.998507
+ 0.998393 0.998393 0.998394 0.998394 0.998395 0.998396 0.998398 0.998401 0.998405 0.998412 0.998419 0.998429
+ 0.998442 0.998456 0.998473 0.998493 0.998515 0.998539 0.998565 0.998592 0.998393 0.998394 0.998394 0.998396
+ 0.998398 0.998403 0.998409 0.998418 0.99843 0.998446 0.998465 0.998487 0.998513 0.998541 0.998573 0.998608
+ 0.998645 0.998683 0.998724 0.998763 0.998394 0.998395 0.998399 0.998406 0.998417 0.998434 0.998455 0.998483
+ 0.998515 0.998553 0.998595 0.998641 0.998689 0.998739 0.998791 0.998844 0.998897 0.998949 0.999001 0.999048
+ 0.998396 0.998406 0.998424 0.998453 0.998493 0.998543 0.9986 0.998664 0.998731 0.998799 0.998869 0.998937
+ 0.999004 0.999069 0.99913 0.999189 0.999244 0.999296 0.999345 0.999386 0.998406 0.998443 0.998501 0.998575
+ 0.998661 0.998752 0.998844 0.998934 0.999019 0.9991 0.999176 0.999245 0.999309 0.999368 0.999421 0.99947
+ 0.999515 0.999555 0.999591 0.999618
+ </DataArray>
+ <DataArray type="Float32" Name="x^H2O_liquid" NumberOfComponents="1" format="ascii">
+ 0.999032 0.999078 0.99913 0.999187 0.999246 0.999305 0.999361 0.999414 0.999463 0.999509 0.999551 0.999589
+ 0.999624 0.999656 0.999685 0.999711 0.999735 0.999757 0.999777 0.999791 0.999004 0.999015 0.999032 0.999055
+ 0.999084 0.999119 0.999157 0.999198 0.99924 0.999283 0.999325 0.999367 0.999407 0.999445 0.999482 0.999517
+ 0.99955 0.999581 0.999609 0.999634 0.999 0.999002 0.999005 0.999011 0.999019 0.99903 0.999045 0.999063
+ 0.999084 0.999108 0.999135 0.999164 0.999194 0.999225 0.999258 0.99929 0.999323 0.999355 0.999387 0.999416
+ 0.999 0.999 0.999001 0.999002 0.999004 0.999006 0.999011 0.999017 0.999025 0.999035 0.999047 0.999061
+ 0.999077 0.999095 0.999115 0.999137 0.99916 0.999184 0.999209 0.999233 0.999 0.999 0.999 0.999
+ 0.999001 0.999002 0.999003 0.999005 0.999008 0.999012 0.999017 0.999023 0.999031 0.99904 0.99905 0.999063
+ 0.999076 0.999091 0.999108 0.999125 0.999 0.999 0.999 0.999 0.999 0.999001 0.999001 0.999002
+ 0.999003 0.999005 0.999007 0.99901 0.999014 0.999019 0.999025 0.999032 0.99904 0.99905 0.99906 0.999071
+ 0.999 0.999 0.999 0.999 0.999 0.999 0.999001 0.999001 0.999002 0.999003 0.999004 0.999006
+ 0.999008 0.999011 0.999015 0.99902 0.999025 0.999031 0.999038 0.999046 0.999 0.999 0.999 0.999
+ 0.999 0.999 0.999 0.999001 0.999001 0.999002 0.999003 0.999004 0.999006 0.999008 0.999011 0.999014
+ 0.999018 0.999023 0.999028 0.999034 0.999 0.999 0.999 0.999 0.999 0.999 0.999 0.999
+ 0.999001 0.999001 0.999002 0.999003 0.999005 0.999006 0.999008 0.999011 0.999014 0.999018 0.999023 0.999028
+ 0.999 0.999 0.999 0.999 0.999 0.999 0.999 0.999 0.999001 0.999001 0.999002 0.999003
+ 0.999004 0.999006 0.999008 0.99901 0.999013 0.999017 0.999021 0.999025 0.999 0.999 0.999 0.999
+ 0.999 0.999 0.999 0.999 0.999001 0.999001 0.999002 0.999003 0.999004 0.999006 0.999008 0.99901
+ 0.999013 0.999017 0.999021 0.999025 0.999 0.999 0.999 0.999 0.999 0.999 0.999 0.999
+ 0.999001 0.999001 0.999002 0.999003 0.999005 0.999006 0.999008 0.999011 0.999014 0.999018 0.999023 0.999028
+ 0.999 0.999 0.999 0.999 0.999 0.999 0.999 0.999001 0.999001 0.999002 0.999003 0.999004
+ 0.999006 0.999008 0.999011 0.999014 0.999018 0.999023 0.999028 0.999034 0.999 0.999 0.999 0.999
+ 0.999 0.999 0.999001 0.999001 0.999002 0.999003 0.999004 0.999006 0.999008 0.999011 0.999015 0.99902
+ 0.999025 0.999031 0.999038 0.999046 0.999 0.999 0.999 0.999 0.999 0.999001 0.999001 0.999002
+ 0.999003 0.999005 0.999007 0.99901 0.999014 0.999019 0.999025 0.999032 0.99904 0.99905 0.99906 0.999071
+ 0.999 0.999 0.999 0.999 0.999001 0.999002 0.999003 0.999005 0.999007 0.999011 0.999016 0.999022
+ 0.99903 0.999039 0.99905 0.999062 0.999076 0.999091 0.999107 0.999124 0.999 0.999 0.999001 0.999001
+ 0.999003 0.999006 0.99901 0.999015 0.999023 0.999033 0.999044 0.999058 0.999074 0.999092 0.999112 0.999133
+ 0.999156 0.999181 0.999206 0.99923 0.999 0.999001 0.999003 0.999008 0.999015 0.999025 0.999039 0.999056
+ 0.999076 0.999099 0.999126 0.999154 0.999184 0.999215 0.999248 0.999281 0.999313 0.999346 0.999378 0.999408
+ 0.999002 0.999008 0.999019 0.999037 0.999062 0.999093 0.999129 0.999168 0.99921 0.999253 0.999296 0.999339
+ 0.99938 0.99942 0.999459 0.999495 0.99953 0.999562 0.999592 0.999618 0.999008 0.999031 0.999067 0.999113
+ 0.999166 0.999223 0.99928 0.999336 0.99939 0.99944 0.999487 0.99953 0.99957 0.999607 0.99964 0.99967
+ 0.999698 0.999723 0.999745 0.999762
+ </DataArray>
+ <DataArray type="Float32" Name="X^Air_liquid" NumberOfComponents="1" format="ascii">
+ 0.00155497 0.00148115 0.0013971 0.00130564 0.00121111 0.00111736 0.00102712 0.000941969 0.000862679 0.000789466 0.000722228 0.000660689
+ 0.000604492 0.000553245 0.000506552 0.000464033 0.000425346 0.000390278 0.00035921 0.000335428 0.00159967 0.00158271 0.00155555 0.00151804
+ 0.00147096 0.00141585 0.0013546 0.00128918 0.0012214 0.00115281 0.00108467 0.00101793 0.000953322 0.000891316 0.000832231 0.000776242
+ 0.000723431 0.000673845 0.000627848 0.000588945 0.00160586 0.00160346 0.00159841 0.00158971 0.00157649 0.00155817 0.00153446 0.00150538
+ 0.00147121 0.00143243 0.00138968 0.00134364 0.00129503 0.00124456 0.00119289 0.00114061 0.00108826 0.00103632 0.000985416 0.000939129
+ 0.00160649 0.00160619 0.00160542 0.00160379 0.00160088 0.00159619 0.00158926 0.00157966 0.00156703 0.00155115 0.00153188 0.0015092
+ 0.0014832 0.00145405 0.001422 0.00138734 0.0013504 0.00131153 0.00127126 0.0012323 0.00160655 0.0016065 0.00160635 0.001606
+ 0.00160529 0.00160401 0.00160191 0.00159874 0.0015942 0.00158803 0.00157998 0.00156984 0.00155746 0.00154271 0.00152554 0.00150595
+ 0.00148397 0.00145969 0.00143333 0.00140653 0.00160655 0.00160654 0.0016065 0.00160638 0.00160613 0.00160565 0.00160482 0.0016035
+ 0.00160153 0.00159872 0.0015949 0.00158988 0.0015835 0.00157561 0.00156608 0.00155481 0.00154172 0.00152678 0.00151002 0.00149237
+ 0.00160656 0.00160655 0.00160653 0.00160647 0.00160635 0.0016061 0.00160566 0.00160495 0.00160385 0.00160226 0.00160004 0.00159707
+ 0.00159321 0.00158834 0.00158232 0.00157507 0.00156648 0.00155648 0.00154503 0.00153272 0.00160656 0.00160655 0.00160654 0.0016065
+ 0.00160642 0.00160626 0.00160598 0.00160551 0.00160477 0.00160369 0.00160216 0.00160008 0.00159736 0.00159387 0.00158953 0.00158422
+ 0.00157787 0.00157039 0.00156174 0.0015523 0.00160656 0.00160655 0.00160654 0.00160652 0.00160646 0.00160634 0.00160612 0.00160575
+ 0.00160518 0.00160432 0.00160311 0.00160146 0.00159927 0.00159645 0.00159291 0.00158857 0.00158333 0.00157712 0.00156989 0.00156195
+ 0.00160656 0.00160655 0.00160655 0.00160652 0.00160647 0.00160636 0.00160617 0.00160585 0.00160534 0.00160458 0.0016035 0.00160202
+ 0.00160005 0.00159751 0.00159431 0.00159036 0.00158559 0.00157992 0.0015733 0.001566 0.00160656 0.00160655 0.00160655 0.00160652
+ 0.00160647 0.00160636 0.00160617 0.00160585 0.00160534 0.00160458 0.0016035 0.00160202 0.00160005 0.00159751 0.0015943 0.00159036
+ 0.00158559 0.00157992 0.0015733 0.001566 0.00160656 0.00160655 0.00160654 0.00160652 0.00160646 0.00160634 0.00160612 0.00160575
+ 0.00160518 0.00160432 0.00160311 0.00160146 0.00159927 0.00159645 0.00159291 0.00158857 0.00158333 0.00157712 0.00156989 0.00156195
+ 0.00160656 0.00160655 0.00160654 0.0016065 0.00160642 0.00160626 0.00160598 0.00160551 0.00160477 0.00160369 0.00160216 0.00160008
+ 0.00159736 0.00159387 0.00158953 0.00158423 0.00157787 0.00157039 0.00156174 0.00155231 0.00160656 0.00160655 0.00160653 0.00160647
+ 0.00160635 0.0016061 0.00160566 0.00160495 0.00160386 0.00160227 0.00160005 0.00159708 0.00159323 0.00158835 0.00158235 0.0015751
+ 0.00156651 0.00155652 0.00154508 0.00153277 0.00160655 0.00160654 0.0016065 0.00160638 0.00160613 0.00160566 0.00160484 0.00160353
+ 0.00160157 0.00159878 0.00159497 0.00158997 0.00158362 0.00157575 0.00156625 0.00155501 0.00154195 0.00152704 0.00151031 0.00149269
+ 0.00160655 0.00160651 0.00160637 0.00160604 0.00160536 0.00160413 0.0016021 0.001599 0.00159456 0.00158849 0.00158055 0.00157053
+ 0.00155827 0.00154364 0.00152658 0.00150709 0.0014852 0.00146102 0.00143473 0.00140798 0.00160651 0.00160629 0.00160565 0.00160425
+ 0.00160164 0.00159734 0.00159084 0.00158171 0.00156957 0.00155415 0.00153533 0.00151305 0.00148741 0.00145856 0.00142674 0.00139225
+ 0.00135542 0.00131662 0.00127636 0.00123737 0.00160619 0.00160465 0.00160096 0.00159399 0.00158272 0.0015664 0.0015446 0.00151723
+ 0.0014845 0.00144685 0.00140491 0.00135938 0.001311 0.00126052 0.00120863 0.00115598 0.00110312 0.00105058 0.000999016 0.000952087
+ 0.00160383 0.00159439 0.0015758 0.00154657 0.00150662 0.00145705 0.00139963 0.00133643 0.00126946 0.00120054 0.00113121 0.00106267
+ 0.000995819 0.000931326 0.000869628 0.000810998 0.000755587 0.000703496 0.000655149 0.000614271 0.00159351 0.00155675 0.00149926 0.0014249
+ 0.00133943 0.00124845 0.00115643 0.00106648 0.00098053 0.000899678 0.000824417 0.000754866 0.000690913 0.000632312 0.000578744 0.000529861
+ 0.000485326 0.000444932 0.000409141 0.00038176
+ </DataArray>
+ <DataArray type="Float32" Name="x^Air_liquid" NumberOfComponents="1" format="ascii">
+ 0.000967873 0.000921896 0.000869556 0.000812603 0.000753741 0.000695373 0.000639188 0.00058618 0.000536823 0.000491251 0.0004494 0.000411098
+ 0.000376123 0.00034423 0.000315172 0.000288712 0.000264638 0.000242816 0.000223484 0.000208687 0.000995711 0.000985147 0.000968234 0.000944869
+ 0.000915551 0.00088123 0.00084309 0.000802353 0.000760148 0.000717441 0.000675016 0.000633471 0.000593248 0.000554649 0.00051787 0.00048302
+ 0.000450148 0.000419286 0.000390659 0.000366447 0.000999565 0.000998072 0.000994928 0.000989506 0.000981275 0.000969864 0.000955096 0.000936983
+ 0.000915703 0.000891558 0.000864933 0.000836263 0.000805996 0.00077457 0.000742396 0.000709848 0.000677256 0.000644915 0.000613227 0.000584413
+ 0.000999958 0.000999774 0.000999289 0.000998277 0.000996463 0.000993546 0.000989229 0.000983247 0.000975384 0.000965491 0.000953488 0.000939364
+ 0.000923173 0.000905021 0.000885061 0.000863475 0.000840471 0.000816272 0.000791195 0.000766937 0.000999994 0.000999965 0.000999874 0.000999655
+ 0.00099921 0.000998412 0.000997108 0.00099513 0.000992303 0.00098846 0.000983447 0.000977135 0.000969421 0.000960236 0.000949544 0.00093734
+ 0.000923652 0.000908534 0.000892117 0.000875426 0.000999999 0.00099999 0.000999962 0.00099989 0.000999733 0.000999434 0.000998919 0.000998098
+ 0.000996868 0.000995119 0.000992738 0.000989613 0.00098564 0.000980725 0.000974789 0.000967769 0.000959619 0.000950315 0.000939876 0.000928887
+ 0.000999999 0.000999996 0.000999982 0.000999947 0.000999868 0.000999714 0.000999442 0.000998997 0.000998316 0.000997324 0.000995944 0.000994092
+ 0.000991688 0.000988651 0.000984907 0.000980388 0.000975038 0.000968808 0.000961681 0.000954013 0.001 0.000999997 0.000999989 0.000999967
+ 0.000999917 0.000999818 0.00099964 0.000999345 0.000998887 0.000998211 0.00099726 0.000995968 0.00099427 0.0009921 0.000989393 0.00098609
+ 0.000982133 0.000977475 0.000972085 0.00096621 0.001 0.000999998 0.000999992 0.000999975 0.000999938 0.000999862 0.000999725 0.000999497
+ 0.00099914 0.000998609 0.000997855 0.000996825 0.000995461 0.000993706 0.000991503 0.000988796 0.000985532 0.000981666 0.000977163 0.000972219
+ 0.001 0.000999998 0.000999993 0.000999979 0.000999946 0.00099988 0.000999759 0.000999558 0.000999241 0.000998768 0.000998095 0.000997172
+ 0.000995946 0.000994363 0.000992369 0.000989912 0.000986941 0.00098341 0.000979287 0.000974743 0.001 0.000999998 0.000999993 0.000999979
+ 0.000999946 0.00099988 0.000999759 0.000999558 0.000999241 0.000998768 0.000998095 0.000997172 0.000995946 0.000994363 0.000992369 0.000989912
+ 0.000986941 0.00098341 0.000979287 0.000974743 0.001 0.000999998 0.000999992 0.000999975 0.000999938 0.000999862 0.000999725 0.000999497
+ 0.00099914 0.000998609 0.000997855 0.000996825 0.000995461 0.000993706 0.000991503 0.000988796 0.000985532 0.000981666 0.000977163 0.000972219
+ 0.001 0.000999997 0.000999989 0.000999967 0.000999917 0.000999818 0.00099964 0.000999345 0.000998887 0.000998212 0.00099726 0.000995968
+ 0.000994271 0.000992101 0.000989395 0.000986091 0.000982135 0.000977477 0.000972088 0.000966214 0.000999999 0.000999996 0.000999982 0.000999947
+ 0.000999868 0.000999715 0.000999443 0.000998999 0.000998318 0.000997328 0.000995949 0.000994099 0.000991697 0.000988662 0.000984921 0.000980406
+ 0.000975059 0.000968833 0.00096171 0.000954046 0.000999999 0.00099999 0.000999963 0.000999891 0.000999736 0.000999441 0.00099893 0.000998115
+ 0.000996893 0.000995154 0.000992784 0.000989672 0.000985714 0.000980815 0.000974896 0.000967894 0.000959763 0.000950477 0.000940056 0.000929083
+ 0.000999995 0.000999969 0.000999885 0.000999679 0.000999256 0.000998488 0.000997224 0.000995294 0.000992525 0.000988747 0.000983804 0.000977564
+ 0.000969924 0.000960812 0.000950189 0.000938051 0.000924422 0.000909358 0.000892986 0.000876333 0.000999971 0.000999831 0.000999435 0.000998563
+ 0.00099694 0.000994259 0.000990214 0.000984523 0.000976961 0.000967362 0.000955636 0.000941765 0.000925793 0.000907826 0.000888011 0.000866533
+ 0.000843601 0.000819439 0.000794369 0.000770094 0.00099977 0.000998813 0.000996512 0.00099217 0.000985154 0.000974992 0.000961414 0.000944367
+ 0.000923982 0.000900539 0.000874419 0.000846064 0.00081594 0.000784507 0.0007522 0.000719416 0.000686508 0.000653794 0.000621694 0.000592479
+ 0.000998302 0.000992424 0.000980846 0.000962639 0.000937761 0.00090689 0.000871133 0.000831772 0.000790071 0.000747163 0.000703995 0.000661319
+ 0.000619704 0.000579555 0.000541148 0.000504653 0.000470163 0.000437741 0.00040765 0.000382209 0.000991872 0.000968977 0.000933178 0.000886869
+ 0.000833643 0.000776994 0.000719698 0.000663691 0.000610185 0.000559854 0.000513006 0.000469714 0.000429909 0.000393437 0.000360099 0.000329677
+ 0.000301963 0.000276826 0.000254554 0.000237516
+ </DataArray>
+ <DataArray type="Float32" Name="D^Air_liquid" NumberOfComponents="1" format="ascii">
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ 1.90888e-09 1.90888e-09 1.90888e-09 1.90888e-09
+ </DataArray>
+ <DataArray type="Float32" Name="velocity_liquid (m/s)" NumberOfComponents="3" format="ascii">
+ 9.78509e-08 -3.50913e-10 0 9.84899e-08 -2.88089e-10 0 9.89534e-08 -1.75455e-10 0 9.92314e-08 -1.02466e-10 0
+ 9.93949e-08 -6.10352e-11 0 9.94935e-08 -3.75577e-11 0 9.95548e-08 -2.37563e-11 0 9.95938e-08 -1.52926e-11 0
+ 9.9619e-08 -9.92122e-12 0 9.96354e-08 -6.43426e-12 0 9.9646e-08 -4.14369e-12 0 9.96528e-08 -2.63432e-12 0
+ 9.9657e-08 -1.64362e-12 0 9.96597e-08 -1.00009e-12 0 9.96613e-08 -5.90025e-13 0 9.96622e-08 -3.3812e-13 0
+ 9.96627e-08 -1.962e-13 0 9.96631e-08 -1.35822e-13 0 9.96633e-08 -1.41669e-13 0 9.96637e-08 -1.93741e-13 0
+ 2.7748e-07 -6.81674e-10 0 2.7759e-07 -7.06311e-10 0 2.77881e-07 -5.12214e-10 0 2.7818e-07 -3.42244e-10 0
+ 2.7842e-07 -2.25104e-10 0 2.78597e-07 -1.4859e-10 0 2.78721e-07 -9.87399e-11 0 2.78808e-07 -6.5875e-11 0
+ 2.78867e-07 -4.39348e-11 0 2.78908e-07 -2.91654e-11 0 2.78935e-07 -1.9193e-11 0 2.78953e-07 -1.24724e-11 0
+ 2.78965e-07 -7.97229e-12 0 2.78973e-07 -4.99212e-12 0 2.78977e-07 -3.05264e-12 0 2.78981e-07 -1.82753e-12 0
+ 2.78982e-07 -1.09917e-12 0 2.78984e-07 -7.27468e-13 0 2.78984e-07 -6.20885e-13 0 2.78985e-07 -6.85268e-13 0
+ 4.3743e-07 -5.91036e-10 0 4.37365e-07 -8.42511e-10 0 4.37442e-07 -7.44241e-10 0 4.37604e-07 -5.7101e-10 0
+ 4.37782e-07 -4.14154e-10 0 4.37938e-07 -2.92937e-10 0 4.38064e-07 -2.04502e-10 0 4.38159e-07 -1.41469e-10 0
+ 4.38228e-07 -9.7022e-11 0 4.38277e-07 -6.58985e-11 0 4.38312e-07 -4.42512e-11 0 4.38336e-07 -2.93156e-11 0
+ 4.38352e-07 -1.91133e-11 0 4.38363e-07 -1.22323e-11 0 4.3837e-07 -7.66974e-12 0 4.38374e-07 -4.72081e-12 0
+ 4.38377e-07 -2.89853e-12 0 4.38379e-07 -1.87251e-12 0 4.3838e-07 -1.41337e-12 0 4.38381e-07 -1.32249e-12 0
+ 5.77438e-07 -4.58939e-10 0 5.7733e-07 -8.0211e-10 0 5.77284e-07 -8.15326e-10 0 5.77318e-07 -6.95518e-10 0
+ 5.77397e-07 -5.46385e-10 0 5.7749e-07 -4.10259e-10 0 5.77576e-07 -2.99536e-10 0 5.7765e-07 -2.14403e-10 0
+ 5.77707e-07 -1.51022e-10 0 5.77751e-07 -1.04834e-10 0 5.77783e-07 -7.17253e-11 0 5.77805e-07 -4.83343e-11 0
+ 5.77821e-07 -3.20412e-11 0 5.77832e-07 -2.08617e-11 0 5.77839e-07 -1.33246e-11 0 5.77844e-07 -8.35998e-12 0
+ 5.77847e-07 -5.20533e-12 0 5.77849e-07 -3.32738e-12 0 5.7785e-07 -2.35268e-12 0 5.77851e-07 -1.9983e-12 0
+ 6.97452e-07 -3.4938e-10 0 6.97341e-07 -6.92649e-10 0 6.97238e-07 -7.75302e-10 0 6.97184e-07 -7.15354e-10 0
+ 6.97177e-07 -5.98348e-10 0 6.97198e-07 -4.72047e-10 0 6.97231e-07 -3.58267e-10 0 6.97267e-07 -2.64393e-10 0
+ 6.97299e-07 -1.9084e-10 0 6.97326e-07 -1.3516e-10 0 6.97347e-07 -9.40702e-11 0 6.97362e-07 -6.43703e-11 0
+ 6.97374e-07 -4.32933e-11 0 6.97382e-07 -2.85971e-11 0 6.97387e-07 -1.85388e-11 0 6.97391e-07 -1.18057e-11 0
+ 6.97394e-07 -7.43557e-12 0 6.97396e-07 -4.74036e-12 0 6.97397e-07 -3.239e-12 0 6.97397e-07 -2.60004e-12 0
+ 7.97462e-07 -2.63303e-10 0 7.97358e-07 -5.6389e-10 0 7.97231e-07 -6.73666e-10 0 7.97125e-07 -6.57868e-10 0
+ 7.97054e-07 -5.77226e-10 0 7.97014e-07 -4.73685e-10 0 7.96995e-07 -3.71227e-10 0 7.96989e-07 -2.81186e-10 0
+ 7.9699e-07 -2.07335e-10 0 7.96994e-07 -1.49475e-10 0 7.96999e-07 -1.05631e-10 0 7.97004e-07 -7.327e-11 0
+ 7.97008e-07 -4.99082e-11 0 7.97011e-07 -3.33787e-11 0 7.97014e-07 -2.19126e-11 0 7.97016e-07 -1.41303e-11 0
+ 7.97017e-07 -8.99384e-12 0 7.97018e-07 -5.7481e-12 0 7.97019e-07 -3.86594e-12 0 7.97019e-07 -3.00792e-12 0
+ 8.77468e-07 -1.93209e-10 0 8.77372e-07 -4.33796e-10 0 8.77236e-07 -5.41222e-10 0 8.77102e-07 -5.50129e-10 0
+ 8.76988e-07 -5.00054e-10 0 8.769e-07 -4.22952e-10 0 8.76837e-07 -3.40001e-10 0 8.76794e-07 -2.63062e-10 0
+ 8.76765e-07 -1.97457e-10 0 8.76746e-07 -1.44526e-10 0 8.76734e-07 -1.03492e-10 0 8.76726e-07 -7.26471e-11 0
+ 8.76722e-07 -5.00425e-11 0 8.76719e-07 -3.38406e-11 0 8.76717e-07 -2.24679e-11 0 8.76717e-07 -1.46566e-11 0
+ 8.76716e-07 -9.4306e-12 0 8.76716e-07 -6.06906e-12 0 8.76716e-07 -4.06937e-12 0 8.76716e-07 -3.12354e-12 0
+ 9.37472e-07 -1.32976e-10 0 9.37381e-07 -3.0728e-10 0 9.37243e-07 -3.94342e-10 0 9.37093e-07 -4.11931e-10 0
+ 9.36952e-07 -3.83937e-10 0 9.36833e-07 -3.32026e-10 0 9.36737e-07 -2.72102e-10 0 9.36664e-07 -2.14055e-10 0
+ 9.36609e-07 -1.63001e-10 0 9.3657e-07 -1.20827e-10 0 9.36542e-07 -8.75172e-11 0 9.36523e-07 -6.20957e-11 0
+ 9.36511e-07 -4.32253e-11 0 9.36502e-07 -2.9546e-11 0 9.36497e-07 -1.9842e-11 0 9.36493e-07 -1.31056e-11 0
+ 9.36491e-07 -8.54568e-12 0 9.3649e-07 -5.57112e-12 0 9.36489e-07 -3.76987e-12 0 9.36489e-07 -2.89869e-12 0
+ 9.77474e-07 -7.85253e-11 0 9.77387e-07 -1.84591e-10 0 9.77249e-07 -2.41092e-10 0 9.77091e-07 -2.56358e-10 0
+ 9.76935e-07 -2.43029e-10 0 9.76795e-07 -2.13501e-10 0 9.76678e-07 -1.77498e-10 0 9.76584e-07 -1.41476e-10 0
+ 9.76512e-07 -1.09051e-10 0 9.76458e-07 -8.17766e-11 0 9.76419e-07 -5.99103e-11 0 9.76391e-07 -4.30054e-11 0
+ 9.76372e-07 -3.03092e-11 0 9.76359e-07 -2.10026e-11 0 9.7635e-07 -1.43267e-11 0 9.76345e-07 -9.63838e-12 0
+ 9.76341e-07 -6.4248e-12 0 9.76339e-07 -4.29935e-12 0 9.76338e-07 -2.99277e-12 0 9.76337e-07 -2.35105e-12 0
+ 9.97475e-07 -2.71203e-11 0 9.9739e-07 -6.45007e-11 0 9.97252e-07 -8.5312e-11 0 9.9709e-07 -9.19871e-11 0
+ 9.96927e-07 -8.85207e-11 0 9.96778e-07 -7.90138e-11 0 9.9665e-07 -6.68153e-11 0 9.96546e-07 -5.42414e-11 0
+ 9.96465e-07 -4.26573e-11 0 9.96403e-07 -3.27089e-11 0 9.96358e-07 -2.45702e-11 0 9.96325e-07 -1.81456e-11 0
+ 9.96303e-07 -1.3212e-11 0 9.96287e-07 -9.50698e-12 0 9.96277e-07 -6.77792e-12 0 9.9627e-07 -4.80556e-12 0
+ 9.96265e-07 -3.41212e-12 0 9.96263e-07 -2.46235e-12 0 9.96261e-07 -1.86256e-12 0 9.9626e-07 -1.56223e-12 0
+ 9.97475e-07 2.32996e-11 0 9.9739e-07 5.46312e-11 0 9.97252e-07 7.10773e-11 0 9.97091e-07 7.50504e-11 0
+ 9.96928e-07 7.03423e-11 0 9.96778e-07 6.07571e-11 0 9.9665e-07 4.93264e-11 0 9.96546e-07 3.80768e-11 0
+ 9.96464e-07 2.81337e-11 0 9.96402e-07 1.99592e-11 0 9.96357e-07 1.3595e-11 0 9.96324e-07 8.85632e-12 0
+ 9.96301e-07 5.46391e-12 0 9.96286e-07 3.12378e-12 0 9.96275e-07 1.56871e-12 0 9.96268e-07 5.7554e-13 0
+ 9.96264e-07 -3.14821e-14 0 9.96261e-07 -3.84003e-13 0 9.96259e-07 -5.74673e-13 0 9.96258e-07 -6.60433e-13 0
+ 9.77474e-07 7.45929e-11 0 9.77387e-07 1.74525e-10 0 9.7725e-07 2.26708e-10 0 9.77091e-07 2.39401e-10 0
+ 9.76935e-07 2.24986e-10 0 9.76795e-07 1.95517e-10 0 9.76677e-07 1.60381e-10 0 9.76583e-07 1.2574e-10 0
+ 9.7651e-07 9.49736e-11 0 9.76456e-07 6.94615e-11 0 9.76416e-07 4.93388e-11 0 9.76388e-07 3.40776e-11 0
+ 9.76368e-07 2.28754e-11 0 9.76355e-07 1.48862e-11 0 9.76346e-07 9.33915e-12 0 9.7634e-07 5.58893e-12 0
+ 9.76336e-07 3.12583e-12 0 9.76334e-07 1.56777e-12 0 9.76332e-07 6.47776e-13 0 9.76331e-07 2.07589e-13 0
+ 9.37472e-07 1.28806e-10 0 9.37382e-07 2.96814e-10 0 9.37244e-07 3.79678e-10 0 9.37094e-07 3.94975e-10 0
+ 9.36953e-07 3.66211e-10 0 9.36833e-07 3.1463e-10 0 9.36736e-07 2.55759e-10 0 9.36662e-07 1.9919e-10 0
+ 9.36606e-07 1.49817e-10 0 9.36566e-07 1.09373e-10 0 9.36538e-07 7.77383e-11 0 9.36518e-07 5.3873e-11 0
+ 9.36504e-07 3.64019e-11 0 9.36495e-07 2.39459e-11 0 9.36489e-07 1.52823e-11 0 9.36485e-07 9.40429e-12 0
+ 9.36483e-07 5.52553e-12 0 9.36481e-07 3.06095e-12 0 9.3648e-07 1.60303e-12 0 9.36479e-07 9.08367e-13 0
+ 8.77468e-07 1.88642e-10 0 8.77373e-07 4.22716e-10 0 8.77238e-07 5.26194e-10 0 8.77103e-07 5.33284e-10 0
+ 8.76989e-07 4.82929e-10 0 8.769e-07 4.06542e-10 0 8.76836e-07 3.24887e-10 0 8.76791e-07 2.49532e-10 0
+ 8.7676e-07 1.85609e-10 0 8.7674e-07 1.34336e-10 0 8.76726e-07 9.48617e-11 0 8.76718e-07 6.54358e-11 0
+ 8.76712e-07 4.40878e-11 0 8.76708e-07 2.89711e-11 0 8.76706e-07 1.85115e-11 0 8.76705e-07 1.14454e-11 0
+ 8.76704e-07 6.8031e-12 0 8.76703e-07 3.87126e-12 0 8.76703e-07 2.15451e-12 0 8.76703e-07 1.35012e-12 0
+ 7.97462e-07 2.58118e-10 0 7.97359e-07 5.51981e-10 0 7.97233e-07 6.58298e-10 0 7.97127e-07 6.4141e-10 0
+ 7.97055e-07 5.61144e-10 0 7.97013e-07 4.58775e-10 0 7.96992e-07 3.57855e-10 0 7.96984e-07 2.69467e-10 0
+ 7.96983e-07 1.97241e-10 0 7.96986e-07 1.40905e-10 0 7.96989e-07 9.84465e-11 0 7.96992e-07 6.73147e-11 0
+ 7.96995e-07 4.50212e-11 0 7.96997e-07 2.94006e-11 0 7.96999e-07 1.8689e-11 0 7.97e-07 1.15134e-11 0
+ 7.97001e-07 6.84307e-12 0 7.97001e-07 3.93103e-12 0 7.97001e-07 2.25964e-12 0 7.97002e-07 1.50091e-12 0
+ 6.97453e-07 3.43259e-10 0 6.97343e-07 6.79763e-10 0 6.9724e-07 7.59851e-10 0 6.97185e-07 6.99832e-10 0
+ 6.97176e-07 5.83966e-10 0 6.97195e-07 4.5927e-10 0 6.97226e-07 3.47186e-10 0 6.97259e-07 2.54929e-10 0
+ 6.97289e-07 1.82849e-10 0 6.97314e-07 1.28478e-10 0 6.97333e-07 8.85348e-11 0 6.97347e-07 5.98246e-11 0
+ 6.97357e-07 3.95904e-11 0 6.97364e-07 2.55992e-11 0 6.97369e-07 1.6117e-11 0 6.97372e-07 9.83825e-12 0
+ 6.97374e-07 5.80745e-12 0 6.97375e-07 3.34386e-12 0 6.97375e-07 1.97668e-12 0 6.97376e-07 1.3916e-12 0
+ 5.77439e-07 4.51437e-10 0 5.77332e-07 7.88359e-10 0 5.77286e-07 8.00502e-10 0 5.77318e-07 6.81862e-10 0
+ 5.77394e-07 5.34571e-10 0 5.77484e-07 4.00304e-10 0 5.77568e-07 2.91239e-10 0 5.77638e-07 2.07525e-10 0
+ 5.77693e-07 1.45342e-10 0 5.77734e-07 1.00164e-10 0 5.77764e-07 6.79062e-11 0 5.77785e-07 4.52299e-11 0
+ 5.778e-07 2.95329e-11 0 5.7781e-07 1.88436e-11 0 5.77816e-07 1.16999e-11 0 5.7782e-07 7.03817e-12 0
+ 5.77823e-07 4.1005e-12 0 5.77824e-07 2.358e-12 0 5.77825e-07 1.44581e-12 0 5.77826e-07 1.10166e-12 0
+ 4.37431e-07 5.81715e-10 0 4.37367e-07 8.2876e-10 0 4.37442e-07 7.31495e-10 0 4.376e-07 5.60539e-10 0
+ 4.37774e-07 4.05838e-10 0 4.37927e-07 2.86354e-10 0 4.38049e-07 1.99256e-10 0 4.38141e-07 1.37257e-10 0
+ 4.38208e-07 9.36243e-11 0 4.38255e-07 6.31528e-11 0 4.38288e-07 4.20354e-11 0 4.38311e-07 2.75336e-11 0
+ 4.38326e-07 1.76861e-11 0 4.38336e-07 1.1092e-11 0 4.38342e-07 6.75543e-12 0 4.38346e-07 3.97557e-12 0
+ 4.38348e-07 2.26684e-12 0 4.38349e-07 1.29931e-12 0 4.3835e-07 8.47472e-13 0 4.38351e-07 7.31554e-13 0
+ 2.7748e-07 6.71079e-10 0 2.77588e-07 6.95079e-10 0 2.77875e-07 5.0375e-10 0 2.78168e-07 3.36235e-10 0
+ 2.78404e-07 2.20795e-10 0 2.78577e-07 1.4541e-10 0 2.78699e-07 9.63235e-11 0 2.78783e-07 6.39968e-11 0
+ 2.7884e-07 4.24535e-11 0 2.78879e-07 2.79881e-11 0 2.78905e-07 1.82551e-11 0 2.78922e-07 1.17263e-11 0
+ 2.78933e-07 7.38038e-12 0 2.78941e-07 4.5231e-12 0 2.78945e-07 2.67871e-12 0 2.78948e-07 1.52234e-12 0
+ 2.78949e-07 8.35727e-13 0 2.7895e-07 4.7642e-13 0 2.7895e-07 3.50861e-13 0 2.78951e-07 3.74146e-13 0
+ 9.78455e-08 3.45479e-10 0 9.84745e-08 2.83556e-10 0 9.89307e-08 1.72607e-10 0 9.9204e-08 1.0071e-10 0
+ 9.93646e-08 5.98983e-11 0 9.94613e-08 3.67748e-11 0 9.95212e-08 2.31878e-11 0 9.95593e-08 1.48637e-11 0
+ 9.95837e-08 9.58982e-12 0 9.95995e-08 6.17484e-12 0 9.96096e-08 3.93957e-12 0 9.9616e-08 2.47374e-12 0
+ 9.962e-08 1.51758e-12 0 9.96224e-08 9.01291e-13 0 9.96239e-08 5.11962e-13 0 9.96246e-08 2.74539e-13 0
+ 9.96251e-08 1.40239e-13 0 9.96253e-08 7.90896e-14 0 9.96254e-08 7.34092e-14 0 9.96256e-08 1.04479e-13 0
+ </DataArray>
+ <DataArray type="Float32" Name="process rank" 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 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 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 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 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 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 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 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>
+ </CellData>
+ <Points>
+ <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii">
+ 0 1 0 0.05 1 0 0 1.05 0 0.05 1.05 0
+ 0.1 1 0 0.1 1.05 0 0.15 1 0 0.15 1.05 0
+ 0.2 1 0 0.2 1.05 0 0.25 1 0 0.25 1.05 0
+ 0.3 1 0 0.3 1.05 0 0.35 1 0 0.35 1.05 0
+ 0.4 1 0 0.4 1.05 0 0.45 1 0 0.45 1.05 0
+ 0.5 1 0 0.5 1.05 0 0.55 1 0 0.55 1.05 0
+ 0.6 1 0 0.6 1.05 0 0.65 1 0 0.65 1.05 0
+ 0.7 1 0 0.7 1.05 0 0.75 1 0 0.75 1.05 0
+ 0.8 1 0 0.8 1.05 0 0.85 1 0 0.85 1.05 0
+ 0.9 1 0 0.9 1.05 0 0.95 1 0 0.95 1.05 0
+ 1 1 0 1 1.05 0 0 1.1 0 0.05 1.1 0
+ 0.1 1.1 0 0.15 1.1 0 0.2 1.1 0 0.25 1.1 0
+ 0.3 1.1 0 0.35 1.1 0 0.4 1.1 0 0.45 1.1 0
+ 0.5 1.1 0 0.55 1.1 0 0.6 1.1 0 0.65 1.1 0
+ 0.7 1.1 0 0.75 1.1 0 0.8 1.1 0 0.85 1.1 0
+ 0.9 1.1 0 0.95 1.1 0 1 1.1 0 0 1.15 0
+ 0.05 1.15 0 0.1 1.15 0 0.15 1.15 0 0.2 1.15 0
+ 0.25 1.15 0 0.3 1.15 0 0.35 1.15 0 0.4 1.15 0
+ 0.45 1.15 0 0.5 1.15 0 0.55 1.15 0 0.6 1.15 0
+ 0.65 1.15 0 0.7 1.15 0 0.75 1.15 0 0.8 1.15 0
+ 0.85 1.15 0 0.9 1.15 0 0.95 1.15 0 1 1.15 0
+ 0 1.2 0 0.05 1.2 0 0.1 1.2 0 0.15 1.2 0
+ 0.2 1.2 0 0.25 1.2 0 0.3 1.2 0 0.35 1.2 0
+ 0.4 1.2 0 0.45 1.2 0 0.5 1.2 0 0.55 1.2 0
+ 0.6 1.2 0 0.65 1.2 0 0.7 1.2 0 0.75 1.2 0
+ 0.8 1.2 0 0.85 1.2 0 0.9 1.2 0 0.95 1.2 0
+ 1 1.2 0 0 1.25 0 0.05 1.25 0 0.1 1.25 0
+ 0.15 1.25 0 0.2 1.25 0 0.25 1.25 0 0.3 1.25 0
+ 0.35 1.25 0 0.4 1.25 0 0.45 1.25 0 0.5 1.25 0
+ 0.55 1.25 0 0.6 1.25 0 0.65 1.25 0 0.7 1.25 0
+ 0.75 1.25 0 0.8 1.25 0 0.85 1.25 0 0.9 1.25 0
+ 0.95 1.25 0 1 1.25 0 0 1.3 0 0.05 1.3 0
+ 0.1 1.3 0 0.15 1.3 0 0.2 1.3 0 0.25 1.3 0
+ 0.3 1.3 0 0.35 1.3 0 0.4 1.3 0 0.45 1.3 0
+ 0.5 1.3 0 0.55 1.3 0 0.6 1.3 0 0.65 1.3 0
+ 0.7 1.3 0 0.75 1.3 0 0.8 1.3 0 0.85 1.3 0
+ 0.9 1.3 0 0.95 1.3 0 1 1.3 0 0 1.35 0
+ 0.05 1.35 0 0.1 1.35 0 0.15 1.35 0 0.2 1.35 0
+ 0.25 1.35 0 0.3 1.35 0 0.35 1.35 0 0.4 1.35 0
+ 0.45 1.35 0 0.5 1.35 0 0.55 1.35 0 0.6 1.35 0
+ 0.65 1.35 0 0.7 1.35 0 0.75 1.35 0 0.8 1.35 0
+ 0.85 1.35 0 0.9 1.35 0 0.95 1.35 0 1 1.35 0
+ 0 1.4 0 0.05 1.4 0 0.1 1.4 0 0.15 1.4 0
+ 0.2 1.4 0 0.25 1.4 0 0.3 1.4 0 0.35 1.4 0
+ 0.4 1.4 0 0.45 1.4 0 0.5 1.4 0 0.55 1.4 0
+ 0.6 1.4 0 0.65 1.4 0 0.7 1.4 0 0.75 1.4 0
+ 0.8 1.4 0 0.85 1.4 0 0.9 1.4 0 0.95 1.4 0
+ 1 1.4 0 0 1.45 0 0.05 1.45 0 0.1 1.45 0
+ 0.15 1.45 0 0.2 1.45 0 0.25 1.45 0 0.3 1.45 0
+ 0.35 1.45 0 0.4 1.45 0 0.45 1.45 0 0.5 1.45 0
+ 0.55 1.45 0 0.6 1.45 0 0.65 1.45 0 0.7 1.45 0
+ 0.75 1.45 0 0.8 1.45 0 0.85 1.45 0 0.9 1.45 0
+ 0.95 1.45 0 1 1.45 0 0 1.5 0 0.05 1.5 0
+ 0.1 1.5 0 0.15 1.5 0 0.2 1.5 0 0.25 1.5 0
+ 0.3 1.5 0 0.35 1.5 0 0.4 1.5 0 0.45 1.5 0
+ 0.5 1.5 0 0.55 1.5 0 0.6 1.5 0 0.65 1.5 0
+ 0.7 1.5 0 0.75 1.5 0 0.8 1.5 0 0.85 1.5 0
+ 0.9 1.5 0 0.95 1.5 0 1 1.5 0 0 1.55 0
+ 0.05 1.55 0 0.1 1.55 0 0.15 1.55 0 0.2 1.55 0
+ 0.25 1.55 0 0.3 1.55 0 0.35 1.55 0 0.4 1.55 0
+ 0.45 1.55 0 0.5 1.55 0 0.55 1.55 0 0.6 1.55 0
+ 0.65 1.55 0 0.7 1.55 0 0.75 1.55 0 0.8 1.55 0
+ 0.85 1.55 0 0.9 1.55 0 0.95 1.55 0 1 1.55 0
+ 0 1.6 0 0.05 1.6 0 0.1 1.6 0 0.15 1.6 0
+ 0.2 1.6 0 0.25 1.6 0 0.3 1.6 0 0.35 1.6 0
+ 0.4 1.6 0 0.45 1.6 0 0.5 1.6 0 0.55 1.6 0
+ 0.6 1.6 0 0.65 1.6 0 0.7 1.6 0 0.75 1.6 0
+ 0.8 1.6 0 0.85 1.6 0 0.9 1.6 0 0.95 1.6 0
+ 1 1.6 0 0 1.65 0 0.05 1.65 0 0.1 1.65 0
+ 0.15 1.65 0 0.2 1.65 0 0.25 1.65 0 0.3 1.65 0
+ 0.35 1.65 0 0.4 1.65 0 0.45 1.65 0 0.5 1.65 0
+ 0.55 1.65 0 0.6 1.65 0 0.65 1.65 0 0.7 1.65 0
+ 0.75 1.65 0 0.8 1.65 0 0.85 1.65 0 0.9 1.65 0
+ 0.95 1.65 0 1 1.65 0 0 1.7 0 0.05 1.7 0
+ 0.1 1.7 0 0.15 1.7 0 0.2 1.7 0 0.25 1.7 0
+ 0.3 1.7 0 0.35 1.7 0 0.4 1.7 0 0.45 1.7 0
+ 0.5 1.7 0 0.55 1.7 0 0.6 1.7 0 0.65 1.7 0
+ 0.7 1.7 0 0.75 1.7 0 0.8 1.7 0 0.85 1.7 0
+ 0.9 1.7 0 0.95 1.7 0 1 1.7 0 0 1.75 0
+ 0.05 1.75 0 0.1 1.75 0 0.15 1.75 0 0.2 1.75 0
+ 0.25 1.75 0 0.3 1.75 0 0.35 1.75 0 0.4 1.75 0
+ 0.45 1.75 0 0.5 1.75 0 0.55 1.75 0 0.6 1.75 0
+ 0.65 1.75 0 0.7 1.75 0 0.75 1.75 0 0.8 1.75 0
+ 0.85 1.75 0 0.9 1.75 0 0.95 1.75 0 1 1.75 0
+ 0 1.8 0 0.05 1.8 0 0.1 1.8 0 0.15 1.8 0
+ 0.2 1.8 0 0.25 1.8 0 0.3 1.8 0 0.35 1.8 0
+ 0.4 1.8 0 0.45 1.8 0 0.5 1.8 0 0.55 1.8 0
+ 0.6 1.8 0 0.65 1.8 0 0.7 1.8 0 0.75 1.8 0
+ 0.8 1.8 0 0.85 1.8 0 0.9 1.8 0 0.95 1.8 0
+ 1 1.8 0 0 1.85 0 0.05 1.85 0 0.1 1.85 0
+ 0.15 1.85 0 0.2 1.85 0 0.25 1.85 0 0.3 1.85 0
+ 0.35 1.85 0 0.4 1.85 0 0.45 1.85 0 0.5 1.85 0
+ 0.55 1.85 0 0.6 1.85 0 0.65 1.85 0 0.7 1.85 0
+ 0.75 1.85 0 0.8 1.85 0 0.85 1.85 0 0.9 1.85 0
+ 0.95 1.85 0 1 1.85 0 0 1.9 0 0.05 1.9 0
+ 0.1 1.9 0 0.15 1.9 0 0.2 1.9 0 0.25 1.9 0
+ 0.3 1.9 0 0.35 1.9 0 0.4 1.9 0 0.45 1.9 0
+ 0.5 1.9 0 0.55 1.9 0 0.6 1.9 0 0.65 1.9 0
+ 0.7 1.9 0 0.75 1.9 0 0.8 1.9 0 0.85 1.9 0
+ 0.9 1.9 0 0.95 1.9 0 1 1.9 0 0 1.95 0
+ 0.05 1.95 0 0.1 1.95 0 0.15 1.95 0 0.2 1.95 0
+ 0.25 1.95 0 0.3 1.95 0 0.35 1.95 0 0.4 1.95 0
+ 0.45 1.95 0 0.5 1.95 0 0.55 1.95 0 0.6 1.95 0
+ 0.65 1.95 0 0.7 1.95 0 0.75 1.95 0 0.8 1.95 0
+ 0.85 1.95 0 0.9 1.95 0 0.95 1.95 0 1 1.95 0
+ 0 2 0 0.05 2 0 0.1 2 0 0.15 2 0
+ 0.2 2 0 0.25 2 0 0.3 2 0 0.35 2 0
+ 0.4 2 0 0.45 2 0 0.5 2 0 0.55 2 0
+ 0.6 2 0 0.65 2 0 0.7 2 0 0.75 2 0
+ 0.8 2 0 0.85 2 0 0.9 2 0 0.95 2 0
+ 1 2 0
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 1 4 5 3 4 6 7 5
+ 6 8 9 7 8 10 11 9 10 12 13 11
+ 12 14 15 13 14 16 17 15 16 18 19 17
+ 18 20 21 19 20 22 23 21 22 24 25 23
+ 24 26 27 25 26 28 29 27 28 30 31 29
+ 30 32 33 31 32 34 35 33 34 36 37 35
+ 36 38 39 37 38 40 41 39 2 3 43 42
+ 3 5 44 43 5 7 45 44 7 9 46 45
+ 9 11 47 46 11 13 48 47 13 15 49 48
+ 15 17 50 49 17 19 51 50 19 21 52 51
+ 21 23 53 52 23 25 54 53 25 27 55 54
+ 27 29 56 55 29 31 57 56 31 33 58 57
+ 33 35 59 58 35 37 60 59 37 39 61 60
+ 39 41 62 61 42 43 64 63 43 44 65 64
+ 44 45 66 65 45 46 67 66 46 47 68 67
+ 47 48 69 68 48 49 70 69 49 50 71 70
+ 50 51 72 71 51 52 73 72 52 53 74 73
+ 53 54 75 74 54 55 76 75 55 56 77 76
+ 56 57 78 77 57 58 79 78 58 59 80 79
+ 59 60 81 80 60 61 82 81 61 62 83 82
+ 63 64 85 84 64 65 86 85 65 66 87 86
+ 66 67 88 87 67 68 89 88 68 69 90 89
+ 69 70 91 90 70 71 92 91 71 72 93 92
+ 72 73 94 93 73 74 95 94 74 75 96 95
+ 75 76 97 96 76 77 98 97 77 78 99 98
+ 78 79 100 99 79 80 101 100 80 81 102 101
+ 81 82 103 102 82 83 104 103 84 85 106 105
+ 85 86 107 106 86 87 108 107 87 88 109 108
+ 88 89 110 109 89 90 111 110 90 91 112 111
+ 91 92 113 112 92 93 114 113 93 94 115 114
+ 94 95 116 115 95 96 117 116 96 97 118 117
+ 97 98 119 118 98 99 120 119 99 100 121 120
+ 100 101 122 121 101 102 123 122 102 103 124 123
+ 103 104 125 124 105 106 127 126 106 107 128 127
+ 107 108 129 128 108 109 130 129 109 110 131 130
+ 110 111 132 131 111 112 133 132 112 113 134 133
+ 113 114 135 134 114 115 136 135 115 116 137 136
+ 116 117 138 137 117 118 139 138 118 119 140 139
+ 119 120 141 140 120 121 142 141 121 122 143 142
+ 122 123 144 143 123 124 145 144 124 125 146 145
+ 126 127 148 147 127 128 149 148 128 129 150 149
+ 129 130 151 150 130 131 152 151 131 132 153 152
+ 132 133 154 153 133 134 155 154 134 135 156 155
+ 135 136 157 156 136 137 158 157 137 138 159 158
+ 138 139 160 159 139 140 161 160 140 141 162 161
+ 141 142 163 162 142 143 164 163 143 144 165 164
+ 144 145 166 165 145 146 167 166 147 148 169 168
+ 148 149 170 169 149 150 171 170 150 151 172 171
+ 151 152 173 172 152 153 174 173 153 154 175 174
+ 154 155 176 175 155 156 177 176 156 157 178 177
+ 157 158 179 178 158 159 180 179 159 160 181 180
+ 160 161 182 181 161 162 183 182 162 163 184 183
+ 163 164 185 184 164 165 186 185 165 166 187 186
+ 166 167 188 187 168 169 190 189 169 170 191 190
+ 170 171 192 191 171 172 193 192 172 173 194 193
+ 173 174 195 194 174 175 196 195 175 176 197 196
+ 176 177 198 197 177 178 199 198 178 179 200 199
+ 179 180 201 200 180 181 202 201 181 182 203 202
+ 182 183 204 203 183 184 205 204 184 185 206 205
+ 185 186 207 206 186 187 208 207 187 188 209 208
+ 189 190 211 210 190 191 212 211 191 192 213 212
+ 192 193 214 213 193 194 215 214 194 195 216 215
+ 195 196 217 216 196 197 218 217 197 198 219 218
+ 198 199 220 219 199 200 221 220 200 201 222 221
+ 201 202 223 222 202 203 224 223 203 204 225 224
+ 204 205 226 225 205 206 227 226 206 207 228 227
+ 207 208 229 228 208 209 230 229 210 211 232 231
+ 211 212 233 232 212 213 234 233 213 214 235 234
+ 214 215 236 235 215 216 237 236 216 217 238 237
+ 217 218 239 238 218 219 240 239 219 220 241 240
+ 220 221 242 241 221 222 243 242 222 223 244 243
+ 223 224 245 244 224 225 246 245 225 226 247 246
+ 226 227 248 247 227 228 249 248 228 229 250 249
+ 229 230 251 250 231 232 253 252 232 233 254 253
+ 233 234 255 254 234 235 256 255 235 236 257 256
+ 236 237 258 257 237 238 259 258 238 239 260 259
+ 239 240 261 260 240 241 262 261 241 242 263 262
+ 242 243 264 263 243 244 265 264 244 245 266 265
+ 245 246 267 266 246 247 268 267 247 248 269 268
+ 248 249 270 269 249 250 271 270 250 251 272 271
+ 252 253 274 273 253 254 275 274 254 255 276 275
+ 255 256 277 276 256 257 278 277 257 258 279 278
+ 258 259 280 279 259 260 281 280 260 261 282 281
+ 261 262 283 282 262 263 284 283 263 264 285 284
+ 264 265 286 285 265 266 287 286 266 267 288 287
+ 267 268 289 288 268 269 290 289 269 270 291 290
+ 270 271 292 291 271 272 293 292 273 274 295 294
+ 274 275 296 295 275 276 297 296 276 277 298 297
+ 277 278 299 298 278 279 300 299 279 280 301 300
+ 280 281 302 301 281 282 303 302 282 283 304 303
+ 283 284 305 304 284 285 306 305 285 286 307 306
+ 286 287 308 307 287 288 309 308 288 289 310 309
+ 289 290 311 310 290 291 312 311 291 292 313 312
+ 292 293 314 313 294 295 316 315 295 296 317 316
+ 296 297 318 317 297 298 319 318 298 299 320 319
+ 299 300 321 320 300 301 322 321 301 302 323 322
+ 302 303 324 323 303 304 325 324 304 305 326 325
+ 305 306 327 326 306 307 328 327 307 308 329 328
+ 308 309 330 329 309 310 331 330 310 311 332 331
+ 311 312 333 332 312 313 334 333 313 314 335 334
+ 315 316 337 336 316 317 338 337 317 318 339 338
+ 318 319 340 339 319 320 341 340 320 321 342 341
+ 321 322 343 342 322 323 344 343 323 324 345 344
+ 324 325 346 345 325 326 347 346 326 327 348 347
+ 327 328 349 348 328 329 350 349 329 330 351 350
+ 330 331 352 351 331 332 353 352 332 333 354 353
+ 333 334 355 354 334 335 356 355 336 337 358 357
+ 337 338 359 358 338 339 360 359 339 340 361 360
+ 340 341 362 361 341 342 363 362 342 343 364 363
+ 343 344 365 364 344 345 366 365 345 346 367 366
+ 346 347 368 367 347 348 369 368 348 349 370 369
+ 349 350 371 370 350 351 372 371 351 352 373 372
+ 352 353 374 373 353 354 375 374 354 355 376 375
+ 355 356 377 376 357 358 379 378 358 359 380 379
+ 359 360 381 380 360 361 382 381 361 362 383 382
+ 362 363 384 383 363 364 385 384 364 365 386 385
+ 365 366 387 386 366 367 388 387 367 368 389 388
+ 368 369 390 389 369 370 391 390 370 371 392 391
+ 371 372 393 392 372 373 394 393 373 374 395 394
+ 374 375 396 395 375 376 397 396 376 377 398 397
+ 378 379 400 399 379 380 401 400 380 381 402 401
+ 381 382 403 402 382 383 404 403 383 384 405 404
+ 384 385 406 405 385 386 407 406 386 387 408 407
+ 387 388 409 408 388 389 410 409 389 390 411 410
+ 390 391 412 411 391 392 413 412 392 393 414 413
+ 393 394 415 414 394 395 416 415 395 396 417 416
+ 396 397 418 417 397 398 419 418 399 400 421 420
+ 400 401 422 421 401 402 423 422 402 403 424 423
+ 403 404 425 424 404 405 426 425 405 406 427 426
+ 406 407 428 427 407 408 429 428 408 409 430 429
+ 409 410 431 430 410 411 432 431 411 412 433 432
+ 412 413 434 433 413 414 435 434 414 415 436 435
+ 415 416 437 436 416 417 438 437 417 418 439 438
+ 418 419 440 439
+ </DataArray>
+ <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+ 4 8 12 16 20 24 28 32 36 40 44 48
+ 52 56 60 64 68 72 76 80 84 88 92 96
+ 100 104 108 112 116 120 124 128 132 136 140 144
+ 148 152 156 160 164 168 172 176 180 184 188 192
+ 196 200 204 208 212 216 220 224 228 232 236 240
+ 244 248 252 256 260 264 268 272 276 280 284 288
+ 292 296 300 304 308 312 316 320 324 328 332 336
+ 340 344 348 352 356 360 364 368 372 376 380 384
+ 388 392 396 400 404 408 412 416 420 424 428 432
+ 436 440 444 448 452 456 460 464 468 472 476 480
+ 484 488 492 496 500 504 508 512 516 520 524 528
+ 532 536 540 544 548 552 556 560 564 568 572 576
+ 580 584 588 592 596 600 604 608 612 616 620 624
+ 628 632 636 640 644 648 652 656 660 664 668 672
+ 676 680 684 688 692 696 700 704 708 712 716 720
+ 724 728 732 736 740 744 748 752 756 760 764 768
+ 772 776 780 784 788 792 796 800 804 808 812 816
+ 820 824 828 832 836 840 844 848 852 856 860 864
+ 868 872 876 880 884 888 892 896 900 904 908 912
+ 916 920 924 928 932 936 940 944 948 952 956 960
+ 964 968 972 976 980 984 988 992 996 1000 1004 1008
+ 1012 1016 1020 1024 1028 1032 1036 1040 1044 1048 1052 1056
+ 1060 1064 1068 1072 1076 1080 1084 1088 1092 1096 1100 1104
+ 1108 1112 1116 1120 1124 1128 1132 1136 1140 1144 1148 1152
+ 1156 1160 1164 1168 1172 1176 1180 1184 1188 1192 1196 1200
+ 1204 1208 1212 1216 1220 1224 1228 1232 1236 1240 1244 1248
+ 1252 1256 1260 1264 1268 1272 1276 1280 1284 1288 1292 1296
+ 1300 1304 1308 1312 1316 1320 1324 1328 1332 1336 1340 1344
+ 1348 1352 1356 1360 1364 1368 1372 1376 1380 1384 1388 1392
+ 1396 1400 1404 1408 1412 1416 1420 1424 1428 1432 1436 1440
+ 1444 1448 1452 1456 1460 1464 1468 1472 1476 1480 1484 1488
+ 1492 1496 1500 1504 1508 1512 1516 1520 1524 1528 1532 1536
+ 1540 1544 1548 1552 1556 1560 1564 1568 1572 1576 1580 1584
+ 1588 1592 1596 1600
+ </DataArray>
+ <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>
--
GitLab