diff --git a/test/implicit/richards/CMakeLists.txt b/test/implicit/richards/CMakeLists.txt
index 0d4bbe12a5a3a025a32c77cd049a7e775664c691..50aacbfa60ac27f1f493af3af478f4a520169eb0 100644
--- a/test/implicit/richards/CMakeLists.txt
+++ b/test/implicit/richards/CMakeLists.txt
@@ -24,6 +24,14 @@ add_dumux_test(test_ccrichards test_ccrichards test_ccrichards.cc
${CMAKE_CURRENT_BINARY_DIR}/richardslenscc-00008.vtu
--command "${CMAKE_CURRENT_BINARY_DIR}/test_ccrichards")
+# comparison to analytical solution - only with cc
+add_dumux_test(test_ccrichardsanalytical test_ccrichardsanalytical test_ccrichardsanalytical.cc
+ python ${CMAKE_SOURCE_DIR}/bin/runtest.py
+ --script fuzzy
+ --files ${CMAKE_SOURCE_DIR}/test/references/richardsanalyticalcc-reference.vtu
+ ${CMAKE_CURRENT_BINARY_DIR}/richardsanalyticalcc-01002.vtu
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_ccrichardsanalytical")
+
# non-isothermal tests
add_dumux_test(test_boxrichardsniconvection test_boxrichardsniconvection test_boxrichardsniconvection.cc
python ${CMAKE_SOURCE_DIR}/bin/runtest.py
@@ -51,4 +59,4 @@ add_dumux_test(test_ccrichardsniconduction test_ccrichardsniconduction test_ccri
--script fuzzy
--files ${CMAKE_SOURCE_DIR}/test/references/richardsniconductioncc-reference.vtu
${CMAKE_CURRENT_BINARY_DIR}/test_ccrichardsniconduction-00007.vtu
- --command "${CMAKE_CURRENT_BINARY_DIR}/test_ccrichardsniconduction")
\ No newline at end of file
+ --command "${CMAKE_CURRENT_BINARY_DIR}/test_ccrichardsniconduction")
diff --git a/test/implicit/richards/RichardsAnalytical.m b/test/implicit/richards/RichardsAnalytical.m
new file mode 100644
index 0000000000000000000000000000000000000000..e32f82e26226171cecf403118ed2aa405778fd53
--- /dev/null
+++ b/test/implicit/richards/RichardsAnalytical.m
@@ -0,0 +1,74 @@
+%% Konstanten
+phi=0.4; %Porosität
+rho_w=1000; %Wasserdichte
+mu_w=10^(-3); %dynamische Viskostät des Wassers bei T=285.15k
+K_int=5*10^(-12); %intrinsische Permeabilität
+g=[0;0;-9.81]; %Gravitationsbeschleunigung
+p_n=10^5; %Luftdruck
+
+pcsweLow=10^4;
+pcsweHigh=1000;
+m=-5*10^5;
+sThres=0.01;
+
+%% Linearmodell-Parameter
+EntryPc=0; %[Pa]
+MaxPc=10^10; %[Pa]
+
+
+%% Brooks-corey-Modell-Parameter
+p_d=1000; %Eindringdruck [Pa]
+lambda=2; %Verteilung der Porengröße
+
+%% Variablen
+x=sym('x'); %Raum
+globalPos1=sym('globalPos1'); %Raum y
+z=sym('z'); %Raum
+time=sym('time'); %Zeit
+pwBottom=sym('pwBottom');
+pwTop=sym('pwTop');
+
+
+%syms x y z t
+
+%% Lösung Wasserdruck
+%p_w=pwBottom+(pwTop-pwBottom)*(globalPos1*globalPos1*globalPos1*time/100)/64; %einfach
+p_w=pwBottom+0.5*(pwTop-pwBottom)*(1+tanh(5*globalPos1-15+time/10)); %komplex
+
+%% Sättigung
+%S_w=0.7;
+%S_w=(p_n-p_w-MaxPc)/(EntryPc-MaxPc); %Lineares Modell
+%S_w=((p_n-p_w)/p_d)^(-lambda)); %Brooks-Corey-Modell nicht regularisiert
+
+% S_w=@(p_w) ( ...
+% (p_w<=10000) .* (p_w) + ... % regularization towards zero
+% and((p_w>10000),(p_w < 30000)) .* (((p_n-p_w)/p_d)^(-lambda)) + ... % Brooks-Corey part
+% (p_w >= 30000) .* (p_w) ... % regularization towards infinity
+% ); % regularisiertes Brooks-Corey-Modell
+
+
+% p_c=p_n-p_w;
+% S_w=@(p_c) ( ...
+% (sThres+(p_c-pcsweLow)/m) .* (p_c>=pcsweLow) + ... % regularization towards zero
+% (((p_c)/p_d)^(-lambda)) .* ((p_c < pcsweLow) && (p_c > pcsweHigh)) +... % Brooks-Corey part
+% (1+(p_c-pcsweHigh)/m) .* (p_c <=pcsweHigh) ... % regularization towards infinity
+% );
+
+p_c=p_n-p_w;
+S_w=@(p_c) ( ...
+ (sThres+(p_c-pcsweLow)/m)*heaviside(p_c-pcsweLow) + ... % regularization towards zero
+ (((p_c)/p_d)^(-lambda))*(heaviside(p_c-pcsweHigh)-heaviside(p_c-pcsweLow))+... % Brooks-Corey part
+ (1+(p_c-pcsweHigh)/m) .* heaviside((-p_c)+pcsweHigh) ... % regularization towards infinity
+ );
+
+
+
+
+ezplot(S_w,[0,5000]);
+%% Permeabilität
+%Krw=max(min(S_w,1.0),0.0);
+%K=K_int*S_w; %Lineares Modell
+K=K_int*S_w(p_c)^(2/lambda+3); %Brooks-Corey-Modell
+
+%% Richards-Quellterm
+q_w=phi*diff((rho_w*S_w(p_c)),time)-divergence((rho_w/mu_w*K*([diff(p_w,x);diff(p_w,globalPos1);diff(p_w,z)]-rho_w*g)),[x,globalPos1,z]);
diff --git a/test/implicit/richards/grids/richardsanalytical-1x64.dgf b/test/implicit/richards/grids/richardsanalytical-1x64.dgf
new file mode 100644
index 0000000000000000000000000000000000000000..1d03a21af4248a436d3f38677be49cbb9280b697
--- /dev/null
+++ b/test/implicit/richards/grids/richardsanalytical-1x64.dgf
@@ -0,0 +1,6 @@
+DGF
+Interval
+0 0 % first corner
+6 4 % second corner
+1 512 % cells in x and in y direction
+#
diff --git a/test/implicit/richards/richardsanalyticalproblem.hh b/test/implicit/richards/richardsanalyticalproblem.hh
new file mode 100644
index 0000000000000000000000000000000000000000..018389930f0485618e98b10f7f880b2d98baeac2
--- /dev/null
+++ b/test/implicit/richards/richardsanalyticalproblem.hh
@@ -0,0 +1,418 @@
+// -*- 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 one-dimensional infiltration problem with a smoth, given solution.
+ *
+ * The source term is calculated analytically. Thus, this example can be used
+ * to calclate the L2 error and to show convergence for grid and time-step
+ * refinement.
+ */
+#ifndef DUMUX_RICHARDS_ANALYTICALPROBLEM_HH
+#define DUMUX_RICHARDS_ANALYTICALPROBLEM_HH
+
+#include <cmath>
+#include <dune/geometry/type.hh>
+#include <dune/geometry/quadraturerules.hh>
+#include <dune/grid/io/file/dgfparser.hh>
+
+#include <dumux/implicit/richards/richardsmodel.hh>
+#include <dumux/material/components/simpleh2o.hh>
+#include <dumux/material/fluidsystems/liquidphase.hh>
+
+#include "richardsanalyticalspatialparams.hh"
+
+namespace Dumux
+{
+
+template <class TypeTag>
+class RichardsAnalyticalProblem;
+
+//////////
+// Specify the properties for the analytical problem
+//////////
+namespace Properties
+{
+NEW_TYPE_TAG(RichardsAnalyticalProblem, INHERITS_FROM(Richards, RichardsAnalyticalSpatialParams));
+NEW_TYPE_TAG(RichardsAnalyticalBoxProblem, INHERITS_FROM(BoxModel, RichardsAnalyticalProblem));
+NEW_TYPE_TAG(RichardsAnalyticalCCProblem, INHERITS_FROM(CCModel, RichardsAnalyticalProblem));
+
+// Use 2d YaspGrid
+SET_TYPE_PROP(RichardsAnalyticalProblem, Grid, Dune::YaspGrid<2>);
+
+// Set the physical problem to be solved
+SET_PROP(RichardsAnalyticalProblem, Problem)
+{ typedef Dumux::RichardsAnalyticalProblem<TypeTag> type; };
+
+// Set the wetting phase
+SET_PROP(RichardsAnalyticalProblem, WettingPhase)
+{
+private:
+ typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar;
+public:
+ typedef Dumux::LiquidPhase<Scalar, Dumux::SimpleH2O<Scalar> > type;
+};
+
+// Enable gravity
+SET_BOOL_PROP(RichardsAnalyticalProblem, ProblemEnableGravity, true);
+
+// Enable partial reassembly of the Jacobian matrix
+SET_BOOL_PROP(RichardsAnalyticalProblem, ImplicitEnablePartialReassemble, true);
+
+// Enable re-use of the Jacobian matrix for the first iteration of a time step
+SET_BOOL_PROP(RichardsAnalyticalProblem, ImplicitEnableJacobianRecycling, true);
+
+// Use forward differences to approximate the Jacobian matrix
+SET_INT_PROP(RichardsAnalyticalProblem, ImplicitNumericDifferenceMethod, +1);
+
+// Set the maximum number of newton iterations of a time step
+SET_INT_PROP(RichardsAnalyticalProblem, NewtonMaxSteps, 28);
+
+// Set the "desireable" number of newton iterations of a time step
+SET_INT_PROP(RichardsAnalyticalProblem, NewtonTargetSteps, 18);
+
+// Do not write the intermediate results of the newton method
+SET_BOOL_PROP(RichardsAnalyticalProblem, NewtonWriteConvergence, false);
+}
+
+/*!
+ * \ingroup RichardsModel
+ * \ingroup ImplicitTestProblems
+ *
+ *\brief A water infiltration problem using Richards model and comparing
+ * to an analytical solution. Implemented by using the source term
+ * defined as the analytical solution.
+ *
+ * The domain is box shaped. Top and bottom boundaries are Dirichlet
+ * boundaries with fixed water pressure (fixed Saturation \f$S_w = 0\f$),
+ * left and right boundary are closed (Neumann 0 boundary).
+ * This problem uses the \ref RichardsModel
+ *
+ * The L2 error is decreasing with decreasing time and space discretization.
+ */
+template <class TypeTag>
+class RichardsAnalyticalProblem : public RichardsProblem<TypeTag>
+{
+ typedef RichardsProblem<TypeTag> ParentType;
+
+ typedef typename GET_PROP_TYPE(TypeTag, GridView) GridView;
+ typedef typename GridView::template Codim<0>::Entity::Geometry Geometry;
+ typedef typename GET_PROP_TYPE(TypeTag, PrimaryVariables) PrimaryVariables;
+ typedef typename GET_PROP_TYPE(TypeTag, BoundaryTypes) BoundaryTypes;
+ typedef typename GET_PROP_TYPE(TypeTag, TimeManager) TimeManager;
+ typedef typename GET_PROP_TYPE(TypeTag, FVElementGeometry) FVElementGeometry;
+ typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar;
+
+ typedef typename GET_PROP_TYPE(TypeTag, Indices) Indices;
+ enum {
+ // copy some indices for convenience
+ pwIdx = Indices::pwIdx,
+ contiEqIdx = Indices::contiEqIdx
+ };
+ enum {
+ // Grid and world dimension
+ dim = GridView::dimension,
+ dimWorld = GridView::dimensionworld
+ };
+
+ typedef typename GridView::template Codim<0>::Entity Element;
+ typedef Dune::FieldVector<Scalar, dimWorld> GlobalPosition;
+
+public:
+ /*!
+ * \brief Constructor
+ *
+ * \param timeManager The Dumux TimeManager for simulation management.
+ * \param gridView The grid view on the spatial domain of the problem
+ */
+ RichardsAnalyticalProblem(TimeManager &timeManager,
+ const GridView &gridView)
+ : ParentType(timeManager, gridView)
+ {
+ eps_ = 3e-6;
+ pnRef_ = 1e5; // air pressure
+ name_ = GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, std::string, Problem, Name);
+ }
+
+ /*!
+ * \name Problem parameters
+ */
+ // \{
+
+ /*!
+ * \brief The problem name
+ *
+ * This is used as a prefix for files generated by the simulation.
+ */
+ const std::string name() const
+ { return name_; }
+
+ /*!
+ * \brief Returns the temperature [K] within a finite volume
+ *
+ * This problem assumes a temperature of 10 degrees Celsius.
+ */
+ Scalar temperature() const
+ { return 273.15 + 10; } // -> 10°C
+
+ /*!
+ * \brief Returns the reference pressure [Pa] of the non-wetting
+ * fluid phase within a finite volume
+ *
+ * This problem assumes a constant reference pressure of 1 bar.
+ *
+ * \param element The DUNE Codim<0> entity which intersects with
+ * the finite volume in question
+ * \param fvGeometry The finite volume geometry of the element
+ * \param scvIdx The sub control volume index inside the finite
+ * volume geometry
+ */
+ Scalar referencePressure(const Element &element,
+ const FVElementGeometry &fvGeometry,
+ const int scvIdx) const
+ { return pnRef_; }
+
+ /*!
+ * \brief Evaluate the source values for a control volume.
+ *
+ * For this method, the \a values parameter stores primary
+ * variables. For this test case, the analytical solution is
+ * used to calculate the source term. See the Matlab script
+ * Richards.m which uses Matlab's Symbolic Toolbox to calclate
+ * the source term.
+ *
+ * \param values Storage for all primary variables of the source term
+ * \param globalPos The position for which the source term is set
+ */
+ void sourceAtPos(PrimaryVariables &values,
+ const GlobalPosition &globalPos) const
+ {
+ const Scalar time = this->timeManager().time() + this->timeManager().timeStepSize();
+ const Scalar pwTop = 98942.8;
+ const Scalar pwBottom = 95641.1;
+
+ // linear model with complex solution
+ // calcluated with Matlab script "Richards.m"
+ values = (std::pow(std::tanh(globalPos[1]*5.0+time*(1.0/1.0E1)-1.5E1),2.0)*(1.0/1.0E1)
+ -1.0/1.0E1)*(pwBottom*(1.0/2.0)-pwTop*(1.0/2.0))*4.0E-8-((std::pow(std::tanh(globalPos[1]
+ *5.0+time*(1.0/1.0E1)-1.5E1),2.0)*5.0-5.0)*(pwBottom*(1.0/2.0)-pwTop*(1.0/2.0))-1.0E3)
+ *(std::pow(std::tanh(globalPos[1]*5.0+time*(1.0/1.0E1)-1.5E1),2.0)*5.0-5.0)*(pwBottom
+ *(1.0/2.0)-pwTop*(1.0/2.0))*5.0E-16+std::tanh(globalPos[1]*5.0+time*(1.0/1.0E1)-1.5E1)
+ *(std::pow(std::tanh(globalPos[1]*5.0+time*(1.0/1.0E1)-1.5E1),2.0)*5.0-5.0)*(pwBottom
+ *(1.0/2.0)-pwTop*(1.0/2.0))*(pwBottom*5.0E-16-(std::tanh(globalPos[1]*5.0+time*(1.0/1.0E1)
+ -1.5E1)+1.0)*(pwBottom*(1.0/2.0)-pwTop*(1.0/2.0))*5.0E-16+4.99995E-6)*1.0E1;
+ }
+
+ // \}
+
+ /*!
+ * \name Boundary conditions
+ */
+ // \{
+
+ /*!
+ * \brief Specifies which kind of boundary condition should be
+ * used for which equation on a given boundary segment.
+ *
+ * \param values The boundary types for the conservation equations
+ * \param globalPos The position for which the boundary type is set
+ */
+ void boundaryTypesAtPos(BoundaryTypes &values,
+ const GlobalPosition &globalPos) const
+ {
+ if (onLowerBoundary_(globalPos) ||
+ onUpperBoundary_(globalPos))
+ {
+ values.setAllDirichlet();
+ }
+ else
+ values.setAllNeumann();
+ }
+
+ /*!
+ * \brief Evaluate the boundary conditions for a dirichlet
+ * boundary segment.
+ *
+ * \param values The dirichlet values for the primary variables
+ * \param globalPos The position for which the Dirichlet value is set
+ *
+ * For this method, the \a values parameter stores primary variables.
+ */
+ void dirichletAtPos(PrimaryVariables &values,
+ const GlobalPosition &globalPos) const
+ {
+ const Scalar time = this->timeManager().time() + this->timeManager().timeStepSize();
+ const Scalar pwTop = 98942.8;
+ const Scalar pwBottom = 95641.1;
+ Scalar pw = pwBottom
+ + 0.5 * (std::tanh( (5.0 * globalPos[1]) - 15.0 + time/10.0) + 1.0) * (pwTop - pwBottom);
+
+ values[pwIdx] = pw;
+ }
+
+ /*!
+ * \brief Evaluate the boundary conditions for a neumann
+ * boundary segment.
+ *
+ * For this method, the \a values parameter stores the mass flux
+ * in normal direction of each phase. Negative values mean influx.
+ *
+ * \param values The neumann values for the conservation equations
+ * \param globalPos The position for which the Neumann value is set
+ */
+ void neumannAtPos(PrimaryVariables &values,
+ const GlobalPosition &globalPos) const
+ {
+ values = 0.0;
+ }
+
+ /*!
+ * \name Volume terms
+ */
+ // \{
+
+ /*!
+ * \brief Evaluate the initial values for a control volume.
+ *
+ * For this method, the \a values parameter stores primary
+ * variables.
+ *
+ * \param values Storage for all primary variables of the initial condition
+ * \param globalPos The position for which the boundary type is set
+ */
+ void initialAtPos(PrimaryVariables &values,
+ const GlobalPosition &globalPos) const
+ {
+ const Scalar time = this->timeManager().time();
+ analyticalSolution(values, time, globalPos);
+ }
+
+ // \}
+
+ /*!
+ * \brief Evaluate the analytical solution.
+ *
+ * \param values The dirichlet values for the primary variables
+ * \param time The time at wich the solution should be evaluated
+ * \param globalPos The position for which the Dirichlet value is set
+ *
+ * For this method, the \a values parameter stores primary variables.
+ */
+ void analyticalSolution(PrimaryVariables &values,
+ const Scalar time,
+ const GlobalPosition &globalPos) const
+ {
+ const Scalar pwTop = 98942.8;
+ const Scalar pwBottom = 95641.1;
+ Scalar pw = pwBottom
+ + 0.5 * (std::tanh( (5.0 * globalPos[1]) - 15.0 + time/10.0) + 1.0) * (pwTop - pwBottom);
+
+ values[pwIdx] = pw;
+ }
+
+ /*!
+ * \brief Calculate the L2 error between the solution given by
+ * dirichletAtPos and the numerical approximation.
+ * \note Works for cell-centered FV only because the numerical
+ * approximation is only evaluated in the cell center (once).
+ * To extend this function to the box method the evaluation
+ * has to be exted to box' subvolumes.
+ */
+ Scalar calculateL2Error()
+ {
+ const unsigned int qOrder = 4;
+ Scalar l2error = 0.0;
+ Scalar l2analytic = 0.0;
+ const Scalar time = this->timeManager().time() + this->timeManager().timeStepSize();
+
+ for (auto elementIt = this->gridView().template begin<0>();
+ elementIt != this->gridView().template end<0>(); ++elementIt)
+ {
+ // value from numerical approximation
+#if DUNE_VERSION_NEWER(DUNE_COMMON, 2, 4)
+ Scalar numericalSolution =
+ this->model().curSol()[this->model().dofMapper().subIndex(*elementIt, 0, 0)];
+#else
+ Scalar numericalSolution =
+ this->model().curSol()[this->model().dofMapper().map(*elementIt, 0, 0)];
+#endif
+ // integrate over element using a quadrature rule
+ Geometry geometry = elementIt->geometry();
+ Dune::GeometryType gt = geometry.type();
+ Dune::QuadratureRule<Scalar, dim> rule =
+ Dune::QuadratureRules<Scalar, dim>::rule(gt, qOrder);
+
+ for (auto qIt = rule.begin(); qIt != rule.end(); ++qIt)
+ {
+ // evaluate analytical solution
+ Dune::FieldVector<Scalar, dim> globalPos = geometry.global(qIt->position());
+ Dune::FieldVector<Scalar, 1> values(0.0);
+ analyticalSolution(values, time, globalPos);
+ // add contributino of current quadrature point
+ l2error += (numericalSolution - values[0]) * (numericalSolution - values[0]) *
+ qIt->weight() * geometry.integrationElement(qIt->position());
+ l2analytic += values[0] * values[0] *
+ qIt->weight() * geometry.integrationElement(qIt->position());
+ }
+ }
+ return std::sqrt(l2error/l2analytic);
+ }
+
+ /*!
+ * \brief Write the relevant secondary variables of the current
+ * solution into an VTK output file.
+ */
+ void writeOutput(const bool verbose = true)
+ {
+ ParentType::writeOutput(verbose);
+
+ Scalar l2error = calculateL2Error();
+
+ // compute L2 error if analytical solution is available
+ std::cout.precision(8);
+ std::cout << "L2 error for "
+ << std::setw(6) << this->gridView().size(0)
+ << " elements: "
+ << std::scientific
+ << l2error
+ << std::endl;
+ }
+
+private:
+
+ // evalutates if global position is at lower boundary
+ bool onLowerBoundary_(const GlobalPosition &globalPos) const
+ {
+ return globalPos[1] < this->bBoxMin()[1] + eps_;
+ }
+
+ // evalutates if global position is at upper boundary
+ bool onUpperBoundary_(const GlobalPosition &globalPos) const
+ {
+ return globalPos[1] > this->bBoxMax()[1] - eps_;
+ }
+
+ Scalar eps_;
+ Scalar pnRef_;
+ std::string name_;
+};
+} //end namespace
+
+#endif
diff --git a/test/implicit/richards/richardsanalyticalspatialparams.hh b/test/implicit/richards/richardsanalyticalspatialparams.hh
new file mode 100644
index 0000000000000000000000000000000000000000..45562277a3a5e77837cd46ed45cc5abed8bf3631
--- /dev/null
+++ b/test/implicit/richards/richardsanalyticalspatialparams.hh
@@ -0,0 +1,168 @@
+// -*- 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 spatial parameters for the RichardsAnalyticalProblem
+ */
+#ifndef DUMUX_RICHARDS_ANALYTICAL_SPATIAL_PARAMETERS_HH
+#define DUMUX_RICHARDS_ANALYTICAL_SPATIAL_PARAMETERS_HH
+
+#include <dumux/material/spatialparams/implicitspatialparams.hh>
+#include <dumux/material/fluidmatrixinteractions/2p/regularizedvangenuchten.hh>
+#include <dumux/material/fluidmatrixinteractions/2p/regularizedbrookscorey.hh>
+#include <dumux/material/fluidmatrixinteractions/2p/brookscorey.hh>
+#include <dumux/material/fluidmatrixinteractions/2p/linearmaterial.hh>
+#include <dumux/material/fluidmatrixinteractions/2p/efftoabslaw.hh>
+
+#include <dumux/implicit/richards/richardsmodel.hh>
+
+namespace Dumux
+{
+
+// forward declaration
+template<class TypeTag>
+class RichardsAnalyticalSpatialParams;
+
+namespace Properties
+{
+// The spatial parameters TypeTag
+NEW_TYPE_TAG(RichardsAnalyticalSpatialParams);
+
+// Set the spatial parameters
+SET_TYPE_PROP(RichardsAnalyticalSpatialParams, SpatialParams, Dumux::RichardsAnalyticalSpatialParams<TypeTag>);
+
+// Set the material law
+SET_PROP(RichardsAnalyticalSpatialParams, MaterialLaw)
+{
+private:
+ typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar;
+ typedef LinearMaterial<Scalar> EffectiveLaw;
+public:
+ typedef EffToAbsLaw<EffectiveLaw> type;
+};
+}
+
+/*!
+ * \ingroup RichardsModel
+ * \ingroup ImplicitTestProblems
+ * \brief The spatial parameters for the RichardsAnalyticalProblem
+ */
+template<class TypeTag>
+class RichardsAnalyticalSpatialParams : public ImplicitSpatialParams<TypeTag>
+{
+ typedef ImplicitSpatialParams<TypeTag> ParentType;
+ typedef typename GET_PROP_TYPE(TypeTag, Grid) Grid;
+ typedef typename GET_PROP_TYPE(TypeTag, GridView) GridView;
+ typedef typename GET_PROP_TYPE(TypeTag, Scalar) Scalar;
+ typedef typename Grid::ctype CoordScalar;
+
+ enum {
+ dim=GridView::dimension,
+ dimWorld=GridView::dimensionworld
+ };
+
+ typedef Dune::FieldVector<CoordScalar,dimWorld> GlobalPosition;
+
+ typedef typename GridView::template Codim<0>::Entity Element;
+ typedef typename GET_PROP_TYPE(TypeTag, FVElementGeometry) FVElementGeometry;
+
+public:
+ typedef typename GET_PROP_TYPE(TypeTag, MaterialLaw) MaterialLaw;
+ //! The parameters of the material law to be used
+ typedef typename MaterialLaw::Params MaterialLawParams;
+
+ /*!
+ * \brief Constructor
+ *
+ * \param gridView The DUNE GridView representing the spatial
+ * domain of the problem.
+ */
+ RichardsAnalyticalSpatialParams(const GridView& gridView)
+ : ParentType(gridView)
+ {
+ K_ = 5e-12;
+ materialParams_.setSwr(0.0);
+ materialParams_.setSnr(0.0);
+ materialParams_.setEntryPc(0);
+ materialParams_.setMaxPc(1e10);
+ }
+
+ /*!
+ * \brief Returns the intrinsic permeability tensor [m^2] at a given location
+ *
+ * \param element An arbitrary DUNE Codim<0> entity of the grid view
+ * \param fvGeometry The current finite volume geometry of the element
+ * \param scvIdx The index of the sub-control volume
+ */
+ Scalar intrinsicPermeability(const Element &element,
+ const FVElementGeometry &fvGeometry,
+ int scvIdx) const
+ {
+ return K_;
+ }
+
+ /*!
+ * \brief Returns the porosity [] at a given location
+ *
+ * \param element An arbitrary DUNE Codim<0> entity of the grid view
+ * \param fvGeometry The current finite volume geometry of the element
+ * \param scvIdx The index of the sub-control volume
+ */
+ Scalar porosity(const Element &element,
+ const FVElementGeometry &fvGeometry,
+ int scvIdx) const
+ { return 0.4; }
+
+ /*!
+ * \brief Returns the parameters for the material law at a given location
+ *
+ * \param element An arbitrary DUNE Codim<0> entity of the grid view
+ * \param fvGeometry The current finite volume geometry of the element
+ * \param scvIdx The index of the sub-control volume
+ */
+ const MaterialLawParams& materialLawParams(const Element &element,
+ const FVElementGeometry &fvGeometry,
+ int scvIdx) const
+ {
+ return materialLawParams(fvGeometry.subContVol[scvIdx].global);
+ }
+
+ /*!
+ * \brief Returns the parameters for the material law at a given location
+ *
+ * This method is not actually required by the Richards model, but provided
+ * for the convenience of the RichardsAnalyticalProblem
+ *
+ * \param globalPos A global coordinate vector
+ */
+ const MaterialLawParams& materialLawParams(const GlobalPosition &globalPos) const
+ {
+ return materialParams_;
+ }
+
+private:
+ // intrinsic permeability
+ Scalar K_;
+
+ MaterialLawParams materialParams_;
+};
+} // end namespace
+
+#endif
diff --git a/test/implicit/richards/test_ccrichardsanalytical.cc b/test/implicit/richards/test_ccrichardsanalytical.cc
new file mode 100644
index 0000000000000000000000000000000000000000..5a876d8ba85e54c428713b5b07231abe1bd87468
--- /dev/null
+++ b/test/implicit/richards/test_ccrichardsanalytical.cc
@@ -0,0 +1,58 @@
+// -*- 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 Test for the Richards CC analytical model.
+ */
+#include "config.h"
+#include "richardsanalyticalproblem.hh"
+#include <dumux/common/start.hh>
+
+/*!
+ * \brief Provides an interface for customizing error messages associated with
+ * reading in parameters.
+ *
+ * \param progName The name of the program, that was tried to be started.
+ * \param errorMsg The error message that was issued by the start function.
+ * Comprises the thing that went wrong and a general help message.
+ */
+void usage(const char *progName, const std::string &errorMsg)
+{
+ if (errorMsg.size() > 0) {
+ std::string errorMessageOut = "\nUsage: ";
+ errorMessageOut += progName;
+ errorMessageOut += " [options]\n";
+ errorMessageOut += errorMsg;
+ errorMessageOut += "\n\nThe list of mandatory options for this program is:\n"
+ "\t-TimeManager.TEnd End of the simulation [s] \n"
+ "\t-TimeManager.DtInitial Initial timestep size [s] \n"
+ "\t-Grid.File Name of the file containing the grid \n"
+ "\t definition in DGF format\n";
+
+ std::cout << errorMessageOut
+ << "\n";
+ }
+}
+
+int main(int argc, char** argv)
+{
+ typedef TTAG(RichardsAnalyticalCCProblem) ProblemTypeTag;
+ return Dumux::start<ProblemTypeTag>(argc, argv, usage);
+}
diff --git a/test/implicit/richards/test_ccrichardsanalytical.input b/test/implicit/richards/test_ccrichardsanalytical.input
new file mode 100644
index 0000000000000000000000000000000000000000..73435ae63ab9999e29044b950f7b8d181058c229
--- /dev/null
+++ b/test/implicit/richards/test_ccrichardsanalytical.input
@@ -0,0 +1,11 @@
+[TimeManager]
+DtInitial = 1e-1 # [s]
+TEnd = 100 # [s]
+MaxTimeStepSize = 1e-1 # [s]
+
+[Grid]
+File = ./grids/richardsanalytical-1x64.dgf
+
+[Problem]
+Name = richardsanalyticalcc
+EnableGravity = false
diff --git a/test/references/richardsanalyticalcc-reference.vtu b/test/references/richardsanalyticalcc-reference.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..203363fcac3cc3fad3620a215607766e1f98028b
--- /dev/null
+++ b/test/references/richardsanalyticalcc-reference.vtu
@@ -0,0 +1,1165 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="512" NumberOfPoints="1026">
+ <CellData Scalars="Sn">
+ <DataArray type="Float32" Name="Sn" NumberOfComponents="1" format="ascii">
+ 4.35874e-07 4.35873e-07 4.35872e-07 4.3587e-07 4.35869e-07 4.35867e-07 4.35865e-07 4.35863e-07 4.35861e-07 4.35859e-07 4.35856e-07 4.35853e-07
+ 4.3585e-07 4.35847e-07 4.35844e-07 4.3584e-07 4.35836e-07 4.35831e-07 4.35826e-07 4.35821e-07 4.35816e-07 4.3581e-07 4.35803e-07 4.35796e-07
+ 4.35788e-07 4.3578e-07 4.35771e-07 4.35762e-07 4.35751e-07 4.3574e-07 4.35728e-07 4.35715e-07 4.357e-07 4.35685e-07 4.35668e-07 4.3565e-07
+ 4.35631e-07 4.3561e-07 4.35587e-07 4.35562e-07 4.35536e-07 4.35507e-07 4.35476e-07 4.35442e-07 4.35406e-07 4.35367e-07 4.35324e-07 4.35278e-07
+ 4.35229e-07 4.35175e-07 4.35117e-07 4.35055e-07 4.34987e-07 4.34914e-07 4.34835e-07 4.34749e-07 4.34657e-07 4.34557e-07 4.34449e-07 4.34333e-07
+ 4.34207e-07 4.34071e-07 4.33924e-07 4.33765e-07 4.33594e-07 4.33409e-07 4.33209e-07 4.32993e-07 4.3276e-07 4.32508e-07 4.32236e-07 4.31943e-07
+ 4.31626e-07 4.31284e-07 4.30916e-07 4.30518e-07 4.30089e-07 4.29627e-07 4.29129e-07 4.28591e-07 4.28012e-07 4.27389e-07 4.26717e-07 4.25994e-07
+ 4.25216e-07 4.24379e-07 4.23478e-07 4.22511e-07 4.21471e-07 4.20354e-07 4.19156e-07 4.1787e-07 4.16492e-07 4.15015e-07 4.13434e-07 4.11743e-07
+ 4.09935e-07 4.08004e-07 4.05943e-07 4.03746e-07 4.01407e-07 3.98919e-07 3.96275e-07 3.93469e-07 3.90495e-07 3.87349e-07 3.84023e-07 3.80515e-07
+ 3.7682e-07 3.72935e-07 3.68857e-07 3.64586e-07 3.6012e-07 3.55462e-07 3.50614e-07 3.45579e-07 3.40363e-07 3.34972e-07 3.29415e-07 3.23701e-07
+ 3.17843e-07 3.11853e-07 3.05745e-07 2.99536e-07 2.93241e-07 2.8688e-07 2.80469e-07 2.7403e-07 2.6758e-07 2.61141e-07 2.5473e-07 2.48369e-07
+ 2.42074e-07 2.35865e-07 2.29757e-07 2.23767e-07 2.17909e-07 2.12195e-07 2.06638e-07 2.01247e-07 1.96031e-07 1.90996e-07 1.86148e-07 1.8149e-07
+ 1.77024e-07 1.72753e-07 1.68675e-07 1.6479e-07 1.61095e-07 1.57587e-07 1.54261e-07 1.51115e-07 1.48141e-07 1.45335e-07 1.42691e-07 1.40203e-07
+ 1.37864e-07 1.35667e-07 1.33606e-07 1.31675e-07 1.29867e-07 1.28176e-07 1.26595e-07 1.25118e-07 1.2374e-07 1.22454e-07 1.21256e-07 1.20139e-07
+ 1.19099e-07 1.18132e-07 1.17231e-07 1.16394e-07 1.15616e-07 1.14893e-07 1.14221e-07 1.13598e-07 1.13019e-07 1.12481e-07 1.11983e-07 1.11521e-07
+ 1.11092e-07 1.10694e-07 1.10326e-07 1.09984e-07 1.09667e-07 1.09374e-07 1.09102e-07 1.0885e-07 1.08617e-07 1.08401e-07 1.08201e-07 1.08016e-07
+ 1.07845e-07 1.07686e-07 1.07539e-07 1.07403e-07 1.07277e-07 1.07161e-07 1.07053e-07 1.06953e-07 1.06861e-07 1.06775e-07 1.06696e-07 1.06623e-07
+ 1.06555e-07 1.06493e-07 1.06435e-07 1.06381e-07 1.06332e-07 1.06286e-07 1.06243e-07 1.06204e-07 1.06168e-07 1.06134e-07 1.06103e-07 1.06074e-07
+ 1.06048e-07 1.06023e-07 1.06e-07 1.05979e-07 1.0596e-07 1.05942e-07 1.05925e-07 1.0591e-07 1.05895e-07 1.05882e-07 1.0587e-07 1.05859e-07
+ 1.05848e-07 1.05839e-07 1.0583e-07 1.05822e-07 1.05814e-07 1.05807e-07 1.058e-07 1.05794e-07 1.05789e-07 1.05784e-07 1.05779e-07 1.05774e-07
+ 1.0577e-07 1.05767e-07 1.05763e-07 1.0576e-07 1.05757e-07 1.05754e-07 1.05751e-07 1.05749e-07 1.05747e-07 1.05745e-07 1.05743e-07 1.05741e-07
+ 1.0574e-07 1.05738e-07 1.05737e-07 1.05736e-07 1.05734e-07 1.05733e-07 1.05732e-07 1.05731e-07 1.05731e-07 1.0573e-07 1.05729e-07 1.05728e-07
+ 1.05728e-07 1.05727e-07 1.05727e-07 1.05726e-07 1.05726e-07 1.05725e-07 1.05725e-07 1.05724e-07 1.05724e-07 1.05724e-07 1.05724e-07 1.05723e-07
+ 1.05723e-07 1.05723e-07 1.05723e-07 1.05722e-07 1.05722e-07 1.05722e-07 1.05722e-07 1.05722e-07 1.05722e-07 1.05722e-07 1.05721e-07 1.05721e-07
+ 1.05721e-07 1.05721e-07 1.05721e-07 1.05721e-07 1.05721e-07 1.05721e-07 1.05721e-07 1.05721e-07 1.05721e-07 1.05721e-07 1.05721e-07 1.05721e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07 1.0572e-07
+ </DataArray>
+ <DataArray type="Float32" Name="Sw" NumberOfComponents="1" format="ascii">
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1 1 1 1 1
+ 1 1 1 1 1 1 1 1
+ </DataArray>
+ <DataArray type="Float32" Name="pn" 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 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 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="pw" NumberOfComponents="1" format="ascii">
+ 95641.3 95641.3 95641.3 95641.3 95641.3 95641.3 95641.4 95641.4 95641.4 95641.4 95641.4 95641.5
+ 95641.5 95641.5 95641.6 95641.6 95641.6 95641.7 95641.7 95641.8 95641.8 95641.9 95642 95642
+ 95642.1 95642.2 95642.3 95642.4 95642.5 95642.6 95642.7 95642.9 95643 95643.1 95643.3 95643.5
+ 95643.7 95643.9 95644.1 95644.4 95644.6 95644.9 95645.2 95645.6 95645.9 95646.3 95646.8 95647.2
+ 95647.7 95648.2 95648.8 95649.5 95650.1 95650.9 95651.6 95652.5 95653.4 95654.4 95655.5 95656.7
+ 95657.9 95659.3 95660.8 95662.3 95664.1 95665.9 95667.9 95670.1 95672.4 95674.9 95677.6 95680.6
+ 95683.7 95687.2 95690.8 95694.8 95699.1 95703.7 95708.7 95714.1 95719.9 95726.1 95732.8 95740.1
+ 95747.8 95756.2 95765.2 95774.9 95785.3 95796.5 95808.4 95821.3 95835.1 95849.8 95865.7 95882.6
+ 95900.6 95920 95940.6 95962.5 95985.9 96010.8 96037.3 96065.3 96095 96126.5 96159.8 96194.9
+ 96231.8 96270.7 96311.4 96354.1 96398.8 96445.4 96493.9 96544.2 96596.4 96650.3 96705.9 96763
+ 96821.6 96881.5 96942.5 97004.6 97067.6 97131.2 97195.3 97259.7 97324.2 97388.6 97452.7 97516.3
+ 97579.3 97641.4 97702.4 97762.3 97820.9 97878 97933.6 97987.5 98039.7 98090 98138.5 98185.1
+ 98229.8 98272.5 98313.2 98352.1 98389.1 98424.1 98457.4 98488.9 98518.6 98546.6 98573.1 98598
+ 98621.4 98643.3 98663.9 98683.2 98701.3 98718.2 98734.1 98748.8 98762.6 98775.5 98787.4 98798.6
+ 98809 98818.7 98827.7 98836.1 98843.8 98851.1 98857.8 98864 98869.8 98875.2 98880.2 98884.8
+ 98889.1 98893.1 98896.7 98900.2 98903.3 98906.3 98909 98911.5 98913.8 98916 98918 98919.8
+ 98921.6 98923.1 98924.6 98926 98927.2 98928.4 98929.5 98930.5 98931.4 98932.2 98933 98933.8
+ 98934.4 98935.1 98935.7 98936.2 98936.7 98937.1 98937.6 98938 98938.3 98938.7 98939 98939.3
+ 98939.5 98939.8 98940 98940.2 98940.4 98940.6 98940.8 98940.9 98941 98941.2 98941.3 98941.4
+ 98941.5 98941.6 98941.7 98941.8 98941.9 98941.9 98942 98942.1 98942.1 98942.2 98942.2 98942.3
+ 98942.3 98942.3 98942.4 98942.4 98942.4 98942.5 98942.5 98942.5 98942.5 98942.6 98942.6 98942.6
+ 98942.6 98942.6 98942.6 98942.6 98942.7 98942.7 98942.7 98942.7 98942.7 98942.7 98942.7 98942.7
+ 98942.7 98942.7 98942.7 98942.7 98942.7 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8 98942.8
+ </DataArray>
+ <DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii">
+ 4358.74 4358.73 4358.72 4358.7 4358.69 4358.67 4358.65 4358.63 4358.61 4358.59 4358.56 4358.53
+ 4358.5 4358.47 4358.44 4358.4 4358.36 4358.31 4358.26 4358.21 4358.16 4358.1 4358.03 4357.96
+ 4357.88 4357.8 4357.71 4357.62 4357.51 4357.4 4357.28 4357.15 4357 4356.85 4356.68 4356.5
+ 4356.31 4356.1 4355.87 4355.62 4355.36 4355.07 4354.76 4354.42 4354.06 4353.67 4353.24 4352.78
+ 4352.29 4351.75 4351.17 4350.55 4349.87 4349.14 4348.35 4347.49 4346.57 4345.57 4344.49 4343.33
+ 4342.07 4340.71 4339.24 4337.65 4335.94 4334.09 4332.09 4329.93 4327.6 4325.08 4322.36 4319.43
+ 4316.26 4312.84 4309.16 4305.18 4300.89 4296.27 4291.29 4285.91 4280.12 4273.89 4267.17 4259.94
+ 4252.16 4243.79 4234.78 4225.11 4214.71 4203.54 4191.56 4178.7 4164.92 4150.15 4134.34 4117.43
+ 4099.35 4080.04 4059.43 4037.46 4014.07 3989.19 3962.75 3934.69 3904.95 3873.49 3840.23 3805.15
+ 3768.2 3729.35 3688.57 3645.86 3601.2 3554.62 3506.14 3455.79 3403.63 3349.72 3294.15 3237.01
+ 3178.43 3118.53 3057.45 2995.36 2932.41 2868.8 2804.69 2740.3 2675.8 2611.41 2547.3 2483.69
+ 2420.74 2358.65 2297.57 2237.67 2179.09 2121.95 2066.38 2012.47 1960.31 1909.96 1861.48 1814.9
+ 1770.24 1727.53 1686.75 1647.9 1610.95 1575.87 1542.61 1511.15 1481.41 1453.35 1426.91 1402.03
+ 1378.64 1356.67 1336.06 1316.75 1298.67 1281.76 1265.95 1251.18 1237.4 1224.54 1212.56 1201.39
+ 1190.99 1181.32 1172.31 1163.94 1156.16 1148.93 1142.21 1135.98 1130.19 1124.81 1119.83 1115.21
+ 1110.92 1106.94 1103.26 1099.84 1096.67 1093.74 1091.02 1088.5 1086.17 1084.01 1082.01 1080.16
+ 1078.45 1076.86 1075.39 1074.03 1072.77 1071.61 1070.53 1069.53 1068.61 1067.75 1066.96 1066.23
+ 1065.55 1064.93 1064.35 1063.81 1063.32 1062.86 1062.43 1062.04 1061.68 1061.34 1061.03 1060.74
+ 1060.48 1060.23 1060 1059.79 1059.6 1059.42 1059.25 1059.1 1058.95 1058.82 1058.7 1058.59
+ 1058.48 1058.39 1058.3 1058.22 1058.14 1058.07 1058 1057.94 1057.89 1057.84 1057.79 1057.74
+ 1057.7 1057.67 1057.63 1057.6 1057.57 1057.54 1057.51 1057.49 1057.47 1057.45 1057.43 1057.41
+ 1057.4 1057.38 1057.37 1057.36 1057.34 1057.33 1057.32 1057.31 1057.31 1057.3 1057.29 1057.28
+ 1057.28 1057.27 1057.27 1057.26 1057.26 1057.25 1057.25 1057.24 1057.24 1057.24 1057.24 1057.23
+ 1057.23 1057.23 1057.23 1057.22 1057.22 1057.22 1057.22 1057.22 1057.22 1057.22 1057.21 1057.21
+ 1057.21 1057.21 1057.21 1057.21 1057.21 1057.21 1057.21 1057.21 1057.21 1057.21 1057.21 1057.21
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2 1057.2
+ </DataArray>
+ <DataArray type="Float32" Name="rhoW" NumberOfComponents="1" format="ascii">
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000
+ </DataArray>
+ <DataArray type="Float32" Name="rhoN" NumberOfComponents="1" format="ascii">
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
+ </DataArray>
+ <DataArray type="Float32" Name="mobW" NumberOfComponents="1" format="ascii">
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
+ 1000 1000 1000 1000 1000 1000 1000 1000
+ </DataArray>
+ <DataArray type="Float32" Name="mobN" NumberOfComponents="1" format="ascii">
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10 1e+10
+ </DataArray>
+ <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ </DataArray>
+ <DataArray type="Float32" Name="temperature" NumberOfComponents="1" format="ascii">
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
+ </DataArray>
+ <DataArray type="Float32" Name="water content" NumberOfComponents="1" format="ascii">
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
+ </DataArray>
+ <DataArray type="Float32" Name="source" NumberOfComponents="1" format="ascii">
+ -7.79181e-05 -8.42482e-05 -9.10926e-05 -9.84928e-05 -0.000106494 -0.000115145 -0.000124499 -0.000134612 -0.000145546 -0.000157369 -0.000170151 -0.000183971
+ -0.000198912 -0.000215067 -0.000232533 -0.000251417 -0.000271833 -0.000293906 -0.00031777 -0.000343569 -0.000371462 -0.000401617 -0.000434217 -0.00046946
+ -0.00050756 -0.000548749 -0.000593275 -0.000641408 -0.00069344 -0.000749684 -0.000810482 -0.0008762 -0.000947235 -0.00102401 -0.001107 -0.00119669
+ -0.00129362 -0.00139838 -0.0015116 -0.00163394 -0.00176614 -0.00190898 -0.00206332 -0.00223008 -0.00241022 -0.00260483 -0.00281504 -0.00304209
+ -0.0032873 -0.0035521 -0.00383803 -0.00414674 -0.00448 -0.00483973 -0.00522796 -0.00564689 -0.00609888 -0.00658644 -0.00711228 -0.00767928
+ -0.00829053 -0.00894932 -0.00965915 -0.0104238 -0.0112471 -0.0121335 -0.0130873 -0.0141133 -0.0152164 -0.016402 -0.0176755 -0.0190428
+ -0.0205098 -0.0220828 -0.0237684 -0.0255732 -0.0275042 -0.0295682 -0.0317724 -0.0341238 -0.0366294 -0.039296 -0.0421303 -0.0451382
+ -0.0483254 -0.0516967 -0.0552559 -0.0590058 -0.0629477 -0.0670809 -0.0714029 -0.0759088 -0.0805907 -0.0854377 -0.090435 -0.0955639
+ -0.100801 -0.106118 -0.111481 -0.11685 -0.12218 -0.127417 -0.132503 -0.137373 -0.141953 -0.146166 -0.149928 -0.153151
+ -0.155741 -0.157607 -0.158652 -0.158786 -0.15792 -0.155973 -0.152874 -0.148564 -0.143 -0.136157 -0.12803 -0.118637
+ -0.108021 -0.0962489 -0.0834128 -0.0696287 -0.0550348 -0.0397891 -0.0240654 -0.00805012 0.00806326 0.0240785 0.0398021 0.0550478
+ 0.0696415 0.0834254 0.0962612 0.108033 0.118649 0.128041 0.136168 0.143011 0.148575 0.152884 0.155983 0.157929
+ 0.158795 0.158661 0.157615 0.155749 0.153158 0.149935 0.146173 0.141959 0.137379 0.132509 0.127422 0.122185
+ 0.116855 0.111485 0.106122 0.100805 0.0955675 0.0904384 0.0854408 0.0805937 0.0759115 0.0714055 0.0670833 0.0629499
+ 0.0590079 0.0552579 0.0516985 0.048327 0.0451397 0.0421317 0.0392974 0.0366306 0.0341249 0.0317735 0.0295692 0.0275051
+ 0.0255741 0.0237692 0.0220836 0.0205104 0.0190434 0.0176761 0.0164025 0.0152169 0.0141137 0.0130877 0.0121339 0.0112475
+ 0.0104241 0.00965947 0.00894961 0.0082908 0.00767953 0.00711251 0.00658665 0.00609908 0.00564707 0.00522813 0.00483989 0.00448015
+ 0.00414688 0.00383816 0.00355222 0.00328741 0.00304219 0.00281513 0.00260491 0.0024103 0.00223015 0.00206339 0.00190905 0.00176619
+ 0.00163399 0.00151165 0.00139843 0.00129367 0.00119673 0.00110703 0.00102405 0.000947265 0.000876229 0.000810508 0.000749709 0.000693462
+ 0.000641429 0.000593294 0.000548766 0.000507577 0.000469475 0.000434231 0.00040163 0.000371474 0.00034358 0.00031778 0.000293915 0.000271842
+ 0.000251425 0.000232541 0.000215074 0.000198919 0.000183977 0.000170156 0.000157374 0.000145551 0.000134616 0.000124503 0.000115149 0.000106497
+ 9.8496e-05 9.10955e-05 8.4251e-05 7.79206e-05 7.20658e-05 6.66508e-05 6.16427e-05 5.70108e-05 5.27269e-05 4.87649e-05 4.51006e-05 4.17116e-05
+ 3.85772e-05 3.56783e-05 3.29973e-05 3.05177e-05 2.82244e-05 2.61034e-05 2.41418e-05 2.23276e-05 2.06498e-05 1.9098e-05 1.76628e-05 1.63355e-05
+ 1.51079e-05 1.39726e-05 1.29225e-05 1.19514e-05 1.10533e-05 1.02226e-05 9.45441e-06 8.74392e-06 8.08681e-06 7.47909e-06 6.91704e-06 6.39722e-06
+ 5.91647e-06 5.47185e-06 5.06063e-06 4.68033e-06 4.3286e-06 4.0033e-06 3.70245e-06 3.42421e-06 3.16688e-06 2.92889e-06 2.70878e-06 2.50521e-06
+ 2.31694e-06 2.14282e-06 1.98179e-06 1.83285e-06 1.69511e-06 1.56772e-06 1.44991e-06 1.34095e-06 1.24017e-06 1.14697e-06 1.06078e-06 9.81059e-07
+ 9.07331e-07 8.39144e-07 7.76082e-07 7.17758e-07 6.63818e-07 6.13931e-07 5.67794e-07 5.25123e-07 4.8566e-07 4.49162e-07 4.15407e-07 3.84189e-07
+ 3.55316e-07 3.28614e-07 3.03918e-07 2.81078e-07 2.59955e-07 2.40419e-07 2.22351e-07 2.05641e-07 1.90187e-07 1.75894e-07 1.62676e-07 1.50451e-07
+ 1.39144e-07 1.28687e-07 1.19016e-07 1.10072e-07 1.018e-07 9.41495e-08 8.70741e-08 8.05304e-08 7.44784e-08 6.88813e-08 6.37048e-08 5.89173e-08
+ 5.44896e-08 5.03946e-08 4.66074e-08 4.31048e-08 3.98654e-08 3.68695e-08 3.40987e-08 3.15361e-08 2.91662e-08 2.69743e-08 2.49471e-08 2.30723e-08
+ 2.13384e-08 1.97348e-08 1.82517e-08 1.68801e-08 1.56115e-08 1.44383e-08 1.33532e-08 1.23497e-08 1.14216e-08 1.05633e-08 9.76944e-09 9.03526e-09
+ 8.35625e-09 7.72827e-09 7.14748e-09 6.61034e-09 6.11356e-09 5.65412e-09 5.22921e-09 4.83623e-09 4.47278e-09 4.13664e-09 3.82577e-09 3.53826e-09
+ 3.27235e-09 3.02643e-09 2.79899e-09 2.58864e-09 2.39411e-09 2.21419e-09 2.04779e-09 1.89389e-09 1.75156e-09 1.61993e-09 1.49819e-09 1.3856e-09
+ 1.28147e-09 1.18517e-09 1.0961e-09 1.01373e-09 9.37545e-10 8.67087e-10 8.01925e-10 7.41659e-10 6.85923e-10 6.34375e-10 5.86701e-10 5.42609e-10
+ 5.01832e-10 4.64118e-10 4.29239e-10 3.96982e-10 3.67148e-10 3.39556e-10 3.14038e-10 2.90438e-10 2.68611e-10 2.48425e-10 2.29755e-10 2.12489e-10
+ 1.9652e-10 1.81751e-10 1.68093e-10 1.5546e-10 1.43777e-10 1.32972e-10 1.22979e-10 1.13737e-10 1.0519e-10 9.72845e-11 8.99734e-11 8.32119e-11
+ 7.69584e-11 7.11749e-11 6.5826e-11 6.08791e-11 5.6304e-11 5.20726e-11 4.81594e-11 4.45402e-11 4.11928e-11 3.80971e-11 3.52341e-11 3.25863e-11
+ 3.01373e-11 2.78725e-11 2.57778e-11 2.38406e-11 2.2049e-11 2.03919e-11 1.88595e-11 1.74421e-11 1.61313e-11 1.4919e-11 1.37979e-11 1.2761e-11
+ 1.1802e-11 1.0915e-11 1.00947e-11 9.33612e-12 8.63451e-12 7.98561e-12 7.38547e-12 6.83048e-12 6.31714e-12 5.84236e-12 5.40329e-12 4.9972e-12
+ 4.62169e-12 4.27433e-12 3.95314e-12 3.65608e-12 3.3813e-12 3.12719e-12 2.89222e-12 2.67485e-12 2.47375e-12 2.28797e-12 2.11598e-12 1.95689e-12
+ 1.80997e-12 1.67382e-12 1.54809e-12 1.43174e-12 1.32412e-12 1.22463e-12 1.13262e-12 1.04751e-12 9.68766e-13 8.95966e-13 8.28591e-13 7.66274e-13
+ 7.08796e-13 6.55497e-13 6.0623e-13 5.60703e-13 5.18474e-13 4.79544e-13 4.43547e-13 4.1019e-13 3.79398e-13 3.50879e-13 3.24486e-13 3.00146e-13
+ 2.77565e-13 2.56671e-13 2.37463e-13 2.19574e-13 2.03079e-13 1.87756e-13 1.7368e-13 1.6063e-13
+ </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 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 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 6 0 0 0 0.0078125 0 6 0.0078125 0
+ 0 0.015625 0 6 0.015625 0 0 0.0234375 0 6 0.0234375 0
+ 0 0.03125 0 6 0.03125 0 0 0.0390625 0 6 0.0390625 0
+ 0 0.046875 0 6 0.046875 0 0 0.0546875 0 6 0.0546875 0
+ 0 0.0625 0 6 0.0625 0 0 0.0703125 0 6 0.0703125 0
+ 0 0.078125 0 6 0.078125 0 0 0.0859375 0 6 0.0859375 0
+ 0 0.09375 0 6 0.09375 0 0 0.101562 0 6 0.101562 0
+ 0 0.109375 0 6 0.109375 0 0 0.117188 0 6 0.117188 0
+ 0 0.125 0 6 0.125 0 0 0.132812 0 6 0.132812 0
+ 0 0.140625 0 6 0.140625 0 0 0.148438 0 6 0.148438 0
+ 0 0.15625 0 6 0.15625 0 0 0.164062 0 6 0.164062 0
+ 0 0.171875 0 6 0.171875 0 0 0.179688 0 6 0.179688 0
+ 0 0.1875 0 6 0.1875 0 0 0.195312 0 6 0.195312 0
+ 0 0.203125 0 6 0.203125 0 0 0.210938 0 6 0.210938 0
+ 0 0.21875 0 6 0.21875 0 0 0.226562 0 6 0.226562 0
+ 0 0.234375 0 6 0.234375 0 0 0.242188 0 6 0.242188 0
+ 0 0.25 0 6 0.25 0 0 0.257812 0 6 0.257812 0
+ 0 0.265625 0 6 0.265625 0 0 0.273438 0 6 0.273438 0
+ 0 0.28125 0 6 0.28125 0 0 0.289062 0 6 0.289062 0
+ 0 0.296875 0 6 0.296875 0 0 0.304688 0 6 0.304688 0
+ 0 0.3125 0 6 0.3125 0 0 0.320312 0 6 0.320312 0
+ 0 0.328125 0 6 0.328125 0 0 0.335938 0 6 0.335938 0
+ 0 0.34375 0 6 0.34375 0 0 0.351562 0 6 0.351562 0
+ 0 0.359375 0 6 0.359375 0 0 0.367188 0 6 0.367188 0
+ 0 0.375 0 6 0.375 0 0 0.382812 0 6 0.382812 0
+ 0 0.390625 0 6 0.390625 0 0 0.398438 0 6 0.398438 0
+ 0 0.40625 0 6 0.40625 0 0 0.414062 0 6 0.414062 0
+ 0 0.421875 0 6 0.421875 0 0 0.429688 0 6 0.429688 0
+ 0 0.4375 0 6 0.4375 0 0 0.445312 0 6 0.445312 0
+ 0 0.453125 0 6 0.453125 0 0 0.460938 0 6 0.460938 0
+ 0 0.46875 0 6 0.46875 0 0 0.476562 0 6 0.476562 0
+ 0 0.484375 0 6 0.484375 0 0 0.492188 0 6 0.492188 0
+ 0 0.5 0 6 0.5 0 0 0.507812 0 6 0.507812 0
+ 0 0.515625 0 6 0.515625 0 0 0.523438 0 6 0.523438 0
+ 0 0.53125 0 6 0.53125 0 0 0.539062 0 6 0.539062 0
+ 0 0.546875 0 6 0.546875 0 0 0.554688 0 6 0.554688 0
+ 0 0.5625 0 6 0.5625 0 0 0.570312 0 6 0.570312 0
+ 0 0.578125 0 6 0.578125 0 0 0.585938 0 6 0.585938 0
+ 0 0.59375 0 6 0.59375 0 0 0.601562 0 6 0.601562 0
+ 0 0.609375 0 6 0.609375 0 0 0.617188 0 6 0.617188 0
+ 0 0.625 0 6 0.625 0 0 0.632812 0 6 0.632812 0
+ 0 0.640625 0 6 0.640625 0 0 0.648438 0 6 0.648438 0
+ 0 0.65625 0 6 0.65625 0 0 0.664062 0 6 0.664062 0
+ 0 0.671875 0 6 0.671875 0 0 0.679688 0 6 0.679688 0
+ 0 0.6875 0 6 0.6875 0 0 0.695312 0 6 0.695312 0
+ 0 0.703125 0 6 0.703125 0 0 0.710938 0 6 0.710938 0
+ 0 0.71875 0 6 0.71875 0 0 0.726562 0 6 0.726562 0
+ 0 0.734375 0 6 0.734375 0 0 0.742188 0 6 0.742188 0
+ 0 0.75 0 6 0.75 0 0 0.757812 0 6 0.757812 0
+ 0 0.765625 0 6 0.765625 0 0 0.773438 0 6 0.773438 0
+ 0 0.78125 0 6 0.78125 0 0 0.789062 0 6 0.789062 0
+ 0 0.796875 0 6 0.796875 0 0 0.804688 0 6 0.804688 0
+ 0 0.8125 0 6 0.8125 0 0 0.820312 0 6 0.820312 0
+ 0 0.828125 0 6 0.828125 0 0 0.835938 0 6 0.835938 0
+ 0 0.84375 0 6 0.84375 0 0 0.851562 0 6 0.851562 0
+ 0 0.859375 0 6 0.859375 0 0 0.867188 0 6 0.867188 0
+ 0 0.875 0 6 0.875 0 0 0.882812 0 6 0.882812 0
+ 0 0.890625 0 6 0.890625 0 0 0.898438 0 6 0.898438 0
+ 0 0.90625 0 6 0.90625 0 0 0.914062 0 6 0.914062 0
+ 0 0.921875 0 6 0.921875 0 0 0.929688 0 6 0.929688 0
+ 0 0.9375 0 6 0.9375 0 0 0.945312 0 6 0.945312 0
+ 0 0.953125 0 6 0.953125 0 0 0.960938 0 6 0.960938 0
+ 0 0.96875 0 6 0.96875 0 0 0.976562 0 6 0.976562 0
+ 0 0.984375 0 6 0.984375 0 0 0.992188 0 6 0.992188 0
+ 0 1 0 6 1 0 0 1.00781 0 6 1.00781 0
+ 0 1.01562 0 6 1.01562 0 0 1.02344 0 6 1.02344 0
+ 0 1.03125 0 6 1.03125 0 0 1.03906 0 6 1.03906 0
+ 0 1.04688 0 6 1.04688 0 0 1.05469 0 6 1.05469 0
+ 0 1.0625 0 6 1.0625 0 0 1.07031 0 6 1.07031 0
+ 0 1.07812 0 6 1.07812 0 0 1.08594 0 6 1.08594 0
+ 0 1.09375 0 6 1.09375 0 0 1.10156 0 6 1.10156 0
+ 0 1.10938 0 6 1.10938 0 0 1.11719 0 6 1.11719 0
+ 0 1.125 0 6 1.125 0 0 1.13281 0 6 1.13281 0
+ 0 1.14062 0 6 1.14062 0 0 1.14844 0 6 1.14844 0
+ 0 1.15625 0 6 1.15625 0 0 1.16406 0 6 1.16406 0
+ 0 1.17188 0 6 1.17188 0 0 1.17969 0 6 1.17969 0
+ 0 1.1875 0 6 1.1875 0 0 1.19531 0 6 1.19531 0
+ 0 1.20312 0 6 1.20312 0 0 1.21094 0 6 1.21094 0
+ 0 1.21875 0 6 1.21875 0 0 1.22656 0 6 1.22656 0
+ 0 1.23438 0 6 1.23438 0 0 1.24219 0 6 1.24219 0
+ 0 1.25 0 6 1.25 0 0 1.25781 0 6 1.25781 0
+ 0 1.26562 0 6 1.26562 0 0 1.27344 0 6 1.27344 0
+ 0 1.28125 0 6 1.28125 0 0 1.28906 0 6 1.28906 0
+ 0 1.29688 0 6 1.29688 0 0 1.30469 0 6 1.30469 0
+ 0 1.3125 0 6 1.3125 0 0 1.32031 0 6 1.32031 0
+ 0 1.32812 0 6 1.32812 0 0 1.33594 0 6 1.33594 0
+ 0 1.34375 0 6 1.34375 0 0 1.35156 0 6 1.35156 0
+ 0 1.35938 0 6 1.35938 0 0 1.36719 0 6 1.36719 0
+ 0 1.375 0 6 1.375 0 0 1.38281 0 6 1.38281 0
+ 0 1.39062 0 6 1.39062 0 0 1.39844 0 6 1.39844 0
+ 0 1.40625 0 6 1.40625 0 0 1.41406 0 6 1.41406 0
+ 0 1.42188 0 6 1.42188 0 0 1.42969 0 6 1.42969 0
+ 0 1.4375 0 6 1.4375 0 0 1.44531 0 6 1.44531 0
+ 0 1.45312 0 6 1.45312 0 0 1.46094 0 6 1.46094 0
+ 0 1.46875 0 6 1.46875 0 0 1.47656 0 6 1.47656 0
+ 0 1.48438 0 6 1.48438 0 0 1.49219 0 6 1.49219 0
+ 0 1.5 0 6 1.5 0 0 1.50781 0 6 1.50781 0
+ 0 1.51562 0 6 1.51562 0 0 1.52344 0 6 1.52344 0
+ 0 1.53125 0 6 1.53125 0 0 1.53906 0 6 1.53906 0
+ 0 1.54688 0 6 1.54688 0 0 1.55469 0 6 1.55469 0
+ 0 1.5625 0 6 1.5625 0 0 1.57031 0 6 1.57031 0
+ 0 1.57812 0 6 1.57812 0 0 1.58594 0 6 1.58594 0
+ 0 1.59375 0 6 1.59375 0 0 1.60156 0 6 1.60156 0
+ 0 1.60938 0 6 1.60938 0 0 1.61719 0 6 1.61719 0
+ 0 1.625 0 6 1.625 0 0 1.63281 0 6 1.63281 0
+ 0 1.64062 0 6 1.64062 0 0 1.64844 0 6 1.64844 0
+ 0 1.65625 0 6 1.65625 0 0 1.66406 0 6 1.66406 0
+ 0 1.67188 0 6 1.67188 0 0 1.67969 0 6 1.67969 0
+ 0 1.6875 0 6 1.6875 0 0 1.69531 0 6 1.69531 0
+ 0 1.70312 0 6 1.70312 0 0 1.71094 0 6 1.71094 0
+ 0 1.71875 0 6 1.71875 0 0 1.72656 0 6 1.72656 0
+ 0 1.73438 0 6 1.73438 0 0 1.74219 0 6 1.74219 0
+ 0 1.75 0 6 1.75 0 0 1.75781 0 6 1.75781 0
+ 0 1.76562 0 6 1.76562 0 0 1.77344 0 6 1.77344 0
+ 0 1.78125 0 6 1.78125 0 0 1.78906 0 6 1.78906 0
+ 0 1.79688 0 6 1.79688 0 0 1.80469 0 6 1.80469 0
+ 0 1.8125 0 6 1.8125 0 0 1.82031 0 6 1.82031 0
+ 0 1.82812 0 6 1.82812 0 0 1.83594 0 6 1.83594 0
+ 0 1.84375 0 6 1.84375 0 0 1.85156 0 6 1.85156 0
+ 0 1.85938 0 6 1.85938 0 0 1.86719 0 6 1.86719 0
+ 0 1.875 0 6 1.875 0 0 1.88281 0 6 1.88281 0
+ 0 1.89062 0 6 1.89062 0 0 1.89844 0 6 1.89844 0
+ 0 1.90625 0 6 1.90625 0 0 1.91406 0 6 1.91406 0
+ 0 1.92188 0 6 1.92188 0 0 1.92969 0 6 1.92969 0
+ 0 1.9375 0 6 1.9375 0 0 1.94531 0 6 1.94531 0
+ 0 1.95312 0 6 1.95312 0 0 1.96094 0 6 1.96094 0
+ 0 1.96875 0 6 1.96875 0 0 1.97656 0 6 1.97656 0
+ 0 1.98438 0 6 1.98438 0 0 1.99219 0 6 1.99219 0
+ 0 2 0 6 2 0 0 2.00781 0 6 2.00781 0
+ 0 2.01562 0 6 2.01562 0 0 2.02344 0 6 2.02344 0
+ 0 2.03125 0 6 2.03125 0 0 2.03906 0 6 2.03906 0
+ 0 2.04688 0 6 2.04688 0 0 2.05469 0 6 2.05469 0
+ 0 2.0625 0 6 2.0625 0 0 2.07031 0 6 2.07031 0
+ 0 2.07812 0 6 2.07812 0 0 2.08594 0 6 2.08594 0
+ 0 2.09375 0 6 2.09375 0 0 2.10156 0 6 2.10156 0
+ 0 2.10938 0 6 2.10938 0 0 2.11719 0 6 2.11719 0
+ 0 2.125 0 6 2.125 0 0 2.13281 0 6 2.13281 0
+ 0 2.14062 0 6 2.14062 0 0 2.14844 0 6 2.14844 0
+ 0 2.15625 0 6 2.15625 0 0 2.16406 0 6 2.16406 0
+ 0 2.17188 0 6 2.17188 0 0 2.17969 0 6 2.17969 0
+ 0 2.1875 0 6 2.1875 0 0 2.19531 0 6 2.19531 0
+ 0 2.20312 0 6 2.20312 0 0 2.21094 0 6 2.21094 0
+ 0 2.21875 0 6 2.21875 0 0 2.22656 0 6 2.22656 0
+ 0 2.23438 0 6 2.23438 0 0 2.24219 0 6 2.24219 0
+ 0 2.25 0 6 2.25 0 0 2.25781 0 6 2.25781 0
+ 0 2.26562 0 6 2.26562 0 0 2.27344 0 6 2.27344 0
+ 0 2.28125 0 6 2.28125 0 0 2.28906 0 6 2.28906 0
+ 0 2.29688 0 6 2.29688 0 0 2.30469 0 6 2.30469 0
+ 0 2.3125 0 6 2.3125 0 0 2.32031 0 6 2.32031 0
+ 0 2.32812 0 6 2.32812 0 0 2.33594 0 6 2.33594 0
+ 0 2.34375 0 6 2.34375 0 0 2.35156 0 6 2.35156 0
+ 0 2.35938 0 6 2.35938 0 0 2.36719 0 6 2.36719 0
+ 0 2.375 0 6 2.375 0 0 2.38281 0 6 2.38281 0
+ 0 2.39062 0 6 2.39062 0 0 2.39844 0 6 2.39844 0
+ 0 2.40625 0 6 2.40625 0 0 2.41406 0 6 2.41406 0
+ 0 2.42188 0 6 2.42188 0 0 2.42969 0 6 2.42969 0
+ 0 2.4375 0 6 2.4375 0 0 2.44531 0 6 2.44531 0
+ 0 2.45312 0 6 2.45312 0 0 2.46094 0 6 2.46094 0
+ 0 2.46875 0 6 2.46875 0 0 2.47656 0 6 2.47656 0
+ 0 2.48438 0 6 2.48438 0 0 2.49219 0 6 2.49219 0
+ 0 2.5 0 6 2.5 0 0 2.50781 0 6 2.50781 0
+ 0 2.51562 0 6 2.51562 0 0 2.52344 0 6 2.52344 0
+ 0 2.53125 0 6 2.53125 0 0 2.53906 0 6 2.53906 0
+ 0 2.54688 0 6 2.54688 0 0 2.55469 0 6 2.55469 0
+ 0 2.5625 0 6 2.5625 0 0 2.57031 0 6 2.57031 0
+ 0 2.57812 0 6 2.57812 0 0 2.58594 0 6 2.58594 0
+ 0 2.59375 0 6 2.59375 0 0 2.60156 0 6 2.60156 0
+ 0 2.60938 0 6 2.60938 0 0 2.61719 0 6 2.61719 0
+ 0 2.625 0 6 2.625 0 0 2.63281 0 6 2.63281 0
+ 0 2.64062 0 6 2.64062 0 0 2.64844 0 6 2.64844 0
+ 0 2.65625 0 6 2.65625 0 0 2.66406 0 6 2.66406 0
+ 0 2.67188 0 6 2.67188 0 0 2.67969 0 6 2.67969 0
+ 0 2.6875 0 6 2.6875 0 0 2.69531 0 6 2.69531 0
+ 0 2.70312 0 6 2.70312 0 0 2.71094 0 6 2.71094 0
+ 0 2.71875 0 6 2.71875 0 0 2.72656 0 6 2.72656 0
+ 0 2.73438 0 6 2.73438 0 0 2.74219 0 6 2.74219 0
+ 0 2.75 0 6 2.75 0 0 2.75781 0 6 2.75781 0
+ 0 2.76562 0 6 2.76562 0 0 2.77344 0 6 2.77344 0
+ 0 2.78125 0 6 2.78125 0 0 2.78906 0 6 2.78906 0
+ 0 2.79688 0 6 2.79688 0 0 2.80469 0 6 2.80469 0
+ 0 2.8125 0 6 2.8125 0 0 2.82031 0 6 2.82031 0
+ 0 2.82812 0 6 2.82812 0 0 2.83594 0 6 2.83594 0
+ 0 2.84375 0 6 2.84375 0 0 2.85156 0 6 2.85156 0
+ 0 2.85938 0 6 2.85938 0 0 2.86719 0 6 2.86719 0
+ 0 2.875 0 6 2.875 0 0 2.88281 0 6 2.88281 0
+ 0 2.89062 0 6 2.89062 0 0 2.89844 0 6 2.89844 0
+ 0 2.90625 0 6 2.90625 0 0 2.91406 0 6 2.91406 0
+ 0 2.92188 0 6 2.92188 0 0 2.92969 0 6 2.92969 0
+ 0 2.9375 0 6 2.9375 0 0 2.94531 0 6 2.94531 0
+ 0 2.95312 0 6 2.95312 0 0 2.96094 0 6 2.96094 0
+ 0 2.96875 0 6 2.96875 0 0 2.97656 0 6 2.97656 0
+ 0 2.98438 0 6 2.98438 0 0 2.99219 0 6 2.99219 0
+ 0 3 0 6 3 0 0 3.00781 0 6 3.00781 0
+ 0 3.01562 0 6 3.01562 0 0 3.02344 0 6 3.02344 0
+ 0 3.03125 0 6 3.03125 0 0 3.03906 0 6 3.03906 0
+ 0 3.04688 0 6 3.04688 0 0 3.05469 0 6 3.05469 0
+ 0 3.0625 0 6 3.0625 0 0 3.07031 0 6 3.07031 0
+ 0 3.07812 0 6 3.07812 0 0 3.08594 0 6 3.08594 0
+ 0 3.09375 0 6 3.09375 0 0 3.10156 0 6 3.10156 0
+ 0 3.10938 0 6 3.10938 0 0 3.11719 0 6 3.11719 0
+ 0 3.125 0 6 3.125 0 0 3.13281 0 6 3.13281 0
+ 0 3.14062 0 6 3.14062 0 0 3.14844 0 6 3.14844 0
+ 0 3.15625 0 6 3.15625 0 0 3.16406 0 6 3.16406 0
+ 0 3.17188 0 6 3.17188 0 0 3.17969 0 6 3.17969 0
+ 0 3.1875 0 6 3.1875 0 0 3.19531 0 6 3.19531 0
+ 0 3.20312 0 6 3.20312 0 0 3.21094 0 6 3.21094 0
+ 0 3.21875 0 6 3.21875 0 0 3.22656 0 6 3.22656 0
+ 0 3.23438 0 6 3.23438 0 0 3.24219 0 6 3.24219 0
+ 0 3.25 0 6 3.25 0 0 3.25781 0 6 3.25781 0
+ 0 3.26562 0 6 3.26562 0 0 3.27344 0 6 3.27344 0
+ 0 3.28125 0 6 3.28125 0 0 3.28906 0 6 3.28906 0
+ 0 3.29688 0 6 3.29688 0 0 3.30469 0 6 3.30469 0
+ 0 3.3125 0 6 3.3125 0 0 3.32031 0 6 3.32031 0
+ 0 3.32812 0 6 3.32812 0 0 3.33594 0 6 3.33594 0
+ 0 3.34375 0 6 3.34375 0 0 3.35156 0 6 3.35156 0
+ 0 3.35938 0 6 3.35938 0 0 3.36719 0 6 3.36719 0
+ 0 3.375 0 6 3.375 0 0 3.38281 0 6 3.38281 0
+ 0 3.39062 0 6 3.39062 0 0 3.39844 0 6 3.39844 0
+ 0 3.40625 0 6 3.40625 0 0 3.41406 0 6 3.41406 0
+ 0 3.42188 0 6 3.42188 0 0 3.42969 0 6 3.42969 0
+ 0 3.4375 0 6 3.4375 0 0 3.44531 0 6 3.44531 0
+ 0 3.45312 0 6 3.45312 0 0 3.46094 0 6 3.46094 0
+ 0 3.46875 0 6 3.46875 0 0 3.47656 0 6 3.47656 0
+ 0 3.48438 0 6 3.48438 0 0 3.49219 0 6 3.49219 0
+ 0 3.5 0 6 3.5 0 0 3.50781 0 6 3.50781 0
+ 0 3.51562 0 6 3.51562 0 0 3.52344 0 6 3.52344 0
+ 0 3.53125 0 6 3.53125 0 0 3.53906 0 6 3.53906 0
+ 0 3.54688 0 6 3.54688 0 0 3.55469 0 6 3.55469 0
+ 0 3.5625 0 6 3.5625 0 0 3.57031 0 6 3.57031 0
+ 0 3.57812 0 6 3.57812 0 0 3.58594 0 6 3.58594 0
+ 0 3.59375 0 6 3.59375 0 0 3.60156 0 6 3.60156 0
+ 0 3.60938 0 6 3.60938 0 0 3.61719 0 6 3.61719 0
+ 0 3.625 0 6 3.625 0 0 3.63281 0 6 3.63281 0
+ 0 3.64062 0 6 3.64062 0 0 3.64844 0 6 3.64844 0
+ 0 3.65625 0 6 3.65625 0 0 3.66406 0 6 3.66406 0
+ 0 3.67188 0 6 3.67188 0 0 3.67969 0 6 3.67969 0
+ 0 3.6875 0 6 3.6875 0 0 3.69531 0 6 3.69531 0
+ 0 3.70312 0 6 3.70312 0 0 3.71094 0 6 3.71094 0
+ 0 3.71875 0 6 3.71875 0 0 3.72656 0 6 3.72656 0
+ 0 3.73438 0 6 3.73438 0 0 3.74219 0 6 3.74219 0
+ 0 3.75 0 6 3.75 0 0 3.75781 0 6 3.75781 0
+ 0 3.76562 0 6 3.76562 0 0 3.77344 0 6 3.77344 0
+ 0 3.78125 0 6 3.78125 0 0 3.78906 0 6 3.78906 0
+ 0 3.79688 0 6 3.79688 0 0 3.80469 0 6 3.80469 0
+ 0 3.8125 0 6 3.8125 0 0 3.82031 0 6 3.82031 0
+ 0 3.82812 0 6 3.82812 0 0 3.83594 0 6 3.83594 0
+ 0 3.84375 0 6 3.84375 0 0 3.85156 0 6 3.85156 0
+ 0 3.85938 0 6 3.85938 0 0 3.86719 0 6 3.86719 0
+ 0 3.875 0 6 3.875 0 0 3.88281 0 6 3.88281 0
+ 0 3.89062 0 6 3.89062 0 0 3.89844 0 6 3.89844 0
+ 0 3.90625 0 6 3.90625 0 0 3.91406 0 6 3.91406 0
+ 0 3.92188 0 6 3.92188 0 0 3.92969 0 6 3.92969 0
+ 0 3.9375 0 6 3.9375 0 0 3.94531 0 6 3.94531 0
+ 0 3.95312 0 6 3.95312 0 0 3.96094 0 6 3.96094 0
+ 0 3.96875 0 6 3.96875 0 0 3.97656 0 6 3.97656 0
+ 0 3.98438 0 6 3.98438 0 0 3.99219 0 6 3.99219 0
+ 0 4 0 6 4 0
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 2 3 5 4 4 5 7 6
+ 6 7 9 8 8 9 11 10 10 11 13 12
+ 12 13 15 14 14 15 17 16 16 17 19 18
+ 18 19 21 20 20 21 23 22 22 23 25 24
+ 24 25 27 26 26 27 29 28 28 29 31 30
+ 30 31 33 32 32 33 35 34 34 35 37 36
+ 36 37 39 38 38 39 41 40 40 41 43 42
+ 42 43 45 44 44 45 47 46 46 47 49 48
+ 48 49 51 50 50 51 53 52 52 53 55 54
+ 54 55 57 56 56 57 59 58 58 59 61 60
+ 60 61 63 62 62 63 65 64 64 65 67 66
+ 66 67 69 68 68 69 71 70 70 71 73 72
+ 72 73 75 74 74 75 77 76 76 77 79 78
+ 78 79 81 80 80 81 83 82 82 83 85 84
+ 84 85 87 86 86 87 89 88 88 89 91 90
+ 90 91 93 92 92 93 95 94 94 95 97 96
+ 96 97 99 98 98 99 101 100 100 101 103 102
+ 102 103 105 104 104 105 107 106 106 107 109 108
+ 108 109 111 110 110 111 113 112 112 113 115 114
+ 114 115 117 116 116 117 119 118 118 119 121 120
+ 120 121 123 122 122 123 125 124 124 125 127 126
+ 126 127 129 128 128 129 131 130 130 131 133 132
+ 132 133 135 134 134 135 137 136 136 137 139 138
+ 138 139 141 140 140 141 143 142 142 143 145 144
+ 144 145 147 146 146 147 149 148 148 149 151 150
+ 150 151 153 152 152 153 155 154 154 155 157 156
+ 156 157 159 158 158 159 161 160 160 161 163 162
+ 162 163 165 164 164 165 167 166 166 167 169 168
+ 168 169 171 170 170 171 173 172 172 173 175 174
+ 174 175 177 176 176 177 179 178 178 179 181 180
+ 180 181 183 182 182 183 185 184 184 185 187 186
+ 186 187 189 188 188 189 191 190 190 191 193 192
+ 192 193 195 194 194 195 197 196 196 197 199 198
+ 198 199 201 200 200 201 203 202 202 203 205 204
+ 204 205 207 206 206 207 209 208 208 209 211 210
+ 210 211 213 212 212 213 215 214 214 215 217 216
+ 216 217 219 218 218 219 221 220 220 221 223 222
+ 222 223 225 224 224 225 227 226 226 227 229 228
+ 228 229 231 230 230 231 233 232 232 233 235 234
+ 234 235 237 236 236 237 239 238 238 239 241 240
+ 240 241 243 242 242 243 245 244 244 245 247 246
+ 246 247 249 248 248 249 251 250 250 251 253 252
+ 252 253 255 254 254 255 257 256 256 257 259 258
+ 258 259 261 260 260 261 263 262 262 263 265 264
+ 264 265 267 266 266 267 269 268 268 269 271 270
+ 270 271 273 272 272 273 275 274 274 275 277 276
+ 276 277 279 278 278 279 281 280 280 281 283 282
+ 282 283 285 284 284 285 287 286 286 287 289 288
+ 288 289 291 290 290 291 293 292 292 293 295 294
+ 294 295 297 296 296 297 299 298 298 299 301 300
+ 300 301 303 302 302 303 305 304 304 305 307 306
+ 306 307 309 308 308 309 311 310 310 311 313 312
+ 312 313 315 314 314 315 317 316 316 317 319 318
+ 318 319 321 320 320 321 323 322 322 323 325 324
+ 324 325 327 326 326 327 329 328 328 329 331 330
+ 330 331 333 332 332 333 335 334 334 335 337 336
+ 336 337 339 338 338 339 341 340 340 341 343 342
+ 342 343 345 344 344 345 347 346 346 347 349 348
+ 348 349 351 350 350 351 353 352 352 353 355 354
+ 354 355 357 356 356 357 359 358 358 359 361 360
+ 360 361 363 362 362 363 365 364 364 365 367 366
+ 366 367 369 368 368 369 371 370 370 371 373 372
+ 372 373 375 374 374 375 377 376 376 377 379 378
+ 378 379 381 380 380 381 383 382 382 383 385 384
+ 384 385 387 386 386 387 389 388 388 389 391 390
+ 390 391 393 392 392 393 395 394 394 395 397 396
+ 396 397 399 398 398 399 401 400 400 401 403 402
+ 402 403 405 404 404 405 407 406 406 407 409 408
+ 408 409 411 410 410 411 413 412 412 413 415 414
+ 414 415 417 416 416 417 419 418 418 419 421 420
+ 420 421 423 422 422 423 425 424 424 425 427 426
+ 426 427 429 428 428 429 431 430 430 431 433 432
+ 432 433 435 434 434 435 437 436 436 437 439 438
+ 438 439 441 440 440 441 443 442 442 443 445 444
+ 444 445 447 446 446 447 449 448 448 449 451 450
+ 450 451 453 452 452 453 455 454 454 455 457 456
+ 456 457 459 458 458 459 461 460 460 461 463 462
+ 462 463 465 464 464 465 467 466 466 467 469 468
+ 468 469 471 470 470 471 473 472 472 473 475 474
+ 474 475 477 476 476 477 479 478 478 479 481 480
+ 480 481 483 482 482 483 485 484 484 485 487 486
+ 486 487 489 488 488 489 491 490 490 491 493 492
+ 492 493 495 494 494 495 497 496 496 497 499 498
+ 498 499 501 500 500 501 503 502 502 503 505 504
+ 504 505 507 506 506 507 509 508 508 509 511 510
+ 510 511 513 512 512 513 515 514 514 515 517 516
+ 516 517 519 518 518 519 521 520 520 521 523 522
+ 522 523 525 524 524 525 527 526 526 527 529 528
+ 528 529 531 530 530 531 533 532 532 533 535 534
+ 534 535 537 536 536 537 539 538 538 539 541 540
+ 540 541 543 542 542 543 545 544 544 545 547 546
+ 546 547 549 548 548 549 551 550 550 551 553 552
+ 552 553 555 554 554 555 557 556 556 557 559 558
+ 558 559 561 560 560 561 563 562 562 563 565 564
+ 564 565 567 566 566 567 569 568 568 569 571 570
+ 570 571 573 572 572 573 575 574 574 575 577 576
+ 576 577 579 578 578 579 581 580 580 581 583 582
+ 582 583 585 584 584 585 587 586 586 587 589 588
+ 588 589 591 590 590 591 593 592 592 593 595 594
+ 594 595 597 596 596 597 599 598 598 599 601 600
+ 600 601 603 602 602 603 605 604 604 605 607 606
+ 606 607 609 608 608 609 611 610 610 611 613 612
+ 612 613 615 614 614 615 617 616 616 617 619 618
+ 618 619 621 620 620 621 623 622 622 623 625 624
+ 624 625 627 626 626 627 629 628 628 629 631 630
+ 630 631 633 632 632 633 635 634 634 635 637 636
+ 636 637 639 638 638 639 641 640 640 641 643 642
+ 642 643 645 644 644 645 647 646 646 647 649 648
+ 648 649 651 650 650 651 653 652 652 653 655 654
+ 654 655 657 656 656 657 659 658 658 659 661 660
+ 660 661 663 662 662 663 665 664 664 665 667 666
+ 666 667 669 668 668 669 671 670 670 671 673 672
+ 672 673 675 674 674 675 677 676 676 677 679 678
+ 678 679 681 680 680 681 683 682 682 683 685 684
+ 684 685 687 686 686 687 689 688 688 689 691 690
+ 690 691 693 692 692 693 695 694 694 695 697 696
+ 696 697 699 698 698 699 701 700 700 701 703 702
+ 702 703 705 704 704 705 707 706 706 707 709 708
+ 708 709 711 710 710 711 713 712 712 713 715 714
+ 714 715 717 716 716 717 719 718 718 719 721 720
+ 720 721 723 722 722 723 725 724 724 725 727 726
+ 726 727 729 728 728 729 731 730 730 731 733 732
+ 732 733 735 734 734 735 737 736 736 737 739 738
+ 738 739 741 740 740 741 743 742 742 743 745 744
+ 744 745 747 746 746 747 749 748 748 749 751 750
+ 750 751 753 752 752 753 755 754 754 755 757 756
+ 756 757 759 758 758 759 761 760 760 761 763 762
+ 762 763 765 764 764 765 767 766 766 767 769 768
+ 768 769 771 770 770 771 773 772 772 773 775 774
+ 774 775 777 776 776 777 779 778 778 779 781 780
+ 780 781 783 782 782 783 785 784 784 785 787 786
+ 786 787 789 788 788 789 791 790 790 791 793 792
+ 792 793 795 794 794 795 797 796 796 797 799 798
+ 798 799 801 800 800 801 803 802 802 803 805 804
+ 804 805 807 806 806 807 809 808 808 809 811 810
+ 810 811 813 812 812 813 815 814 814 815 817 816
+ 816 817 819 818 818 819 821 820 820 821 823 822
+ 822 823 825 824 824 825 827 826 826 827 829 828
+ 828 829 831 830 830 831 833 832 832 833 835 834
+ 834 835 837 836 836 837 839 838 838 839 841 840
+ 840 841 843 842 842 843 845 844 844 845 847 846
+ 846 847 849 848 848 849 851 850 850 851 853 852
+ 852 853 855 854 854 855 857 856 856 857 859 858
+ 858 859 861 860 860 861 863 862 862 863 865 864
+ 864 865 867 866 866 867 869 868 868 869 871 870
+ 870 871 873 872 872 873 875 874 874 875 877 876
+ 876 877 879 878 878 879 881 880 880 881 883 882
+ 882 883 885 884 884 885 887 886 886 887 889 888
+ 888 889 891 890 890 891 893 892 892 893 895 894
+ 894 895 897 896 896 897 899 898 898 899 901 900
+ 900 901 903 902 902 903 905 904 904 905 907 906
+ 906 907 909 908 908 909 911 910 910 911 913 912
+ 912 913 915 914 914 915 917 916 916 917 919 918
+ 918 919 921 920 920 921 923 922 922 923 925 924
+ 924 925 927 926 926 927 929 928 928 929 931 930
+ 930 931 933 932 932 933 935 934 934 935 937 936
+ 936 937 939 938 938 939 941 940 940 941 943 942
+ 942 943 945 944 944 945 947 946 946 947 949 948
+ 948 949 951 950 950 951 953 952 952 953 955 954
+ 954 955 957 956 956 957 959 958 958 959 961 960
+ 960 961 963 962 962 963 965 964 964 965 967 966
+ 966 967 969 968 968 969 971 970 970 971 973 972
+ 972 973 975 974 974 975 977 976 976 977 979 978
+ 978 979 981 980 980 981 983 982 982 983 985 984
+ 984 985 987 986 986 987 989 988 988 989 991 990
+ 990 991 993 992 992 993 995 994 994 995 997 996
+ 996 997 999 998 998 999 1001 1000 1000 1001 1003 1002
+ 1002 1003 1005 1004 1004 1005 1007 1006 1006 1007 1009 1008
+ 1008 1009 1011 1010 1010 1011 1013 1012 1012 1013 1015 1014
+ 1014 1015 1017 1016 1016 1017 1019 1018 1018 1019 1021 1020
+ 1020 1021 1023 1022 1022 1023 1025 1024
+ </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 1604 1608 1612 1616 1620 1624 1628 1632
+ 1636 1640 1644 1648 1652 1656 1660 1664 1668 1672 1676 1680
+ 1684 1688 1692 1696 1700 1704 1708 1712 1716 1720 1724 1728
+ 1732 1736 1740 1744 1748 1752 1756 1760 1764 1768 1772 1776
+ 1780 1784 1788 1792 1796 1800 1804 1808 1812 1816 1820 1824
+ 1828 1832 1836 1840 1844 1848 1852 1856 1860 1864 1868 1872
+ 1876 1880 1884 1888 1892 1896 1900 1904 1908 1912 1916 1920
+ 1924 1928 1932 1936 1940 1944 1948 1952 1956 1960 1964 1968
+ 1972 1976 1980 1984 1988 1992 1996 2000 2004 2008 2012 2016
+ 2020 2024 2028 2032 2036 2040 2044 2048
+ </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 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 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>