diff --git a/test/freeflow/navierstokes/CMakeLists.txt b/test/freeflow/navierstokes/CMakeLists.txt index 172941835429ee8a0e11d215435ceb169dd38625..35794615eb04968d3c31526c95e1e05729dfb58e 100644 --- a/test/freeflow/navierstokes/CMakeLists.txt +++ b/test/freeflow/navierstokes/CMakeLists.txt @@ -97,3 +97,24 @@ dune_add_test(NAME test_navierstokes_angeli --files ${CMAKE_SOURCE_DIR}/test/references/test_angeli-reference.vtu ${CMAKE_CURRENT_BINARY_DIR}/test_angeli-00045.vtu --command "${CMAKE_CURRENT_BINARY_DIR}/test_navierstokes_angeli") + + +dune_add_test(NAME test_stokes_channel_3d + SOURCES test_stokes_channel_3d.cc + COMPILE_DEFINITIONS DIM_3D=1 + CMAKE_GUARD HAVE_UMFPACK + COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py + CMD_ARGS --script fuzzy + --files ${CMAKE_SOURCE_DIR}/test/references/stokes_channel_3d-reference.vtu + ${CMAKE_CURRENT_BINARY_DIR}/test_stokes_channel_3d-00001.vtu + --command "${CMAKE_CURRENT_BINARY_DIR}/test_stokes_channel_3d") + +dune_add_test(NAME test_stokes_channel_pseudo3d + SOURCES test_stokes_channel_3d.cc + COMPILE_DEFINITIONS DIM_3D=0 + CMAKE_GUARD HAVE_UMFPACK + COMMAND ${CMAKE_SOURCE_DIR}/bin/testing/runtest.py + CMD_ARGS --script fuzzy + --files ${CMAKE_SOURCE_DIR}/test/references/stokes_channel_pseudo3d-reference.vtu + ${CMAKE_CURRENT_BINARY_DIR}/test_stokes_channel_pseudo3d-00001.vtu + --command "${CMAKE_CURRENT_BINARY_DIR}/test_stokes_channel_pseudo3d") diff --git a/test/freeflow/navierstokes/stokeschannel3dproblem.hh b/test/freeflow/navierstokes/stokeschannel3dproblem.hh new file mode 100644 index 0000000000000000000000000000000000000000..1890daef5aaaea9806a6e869304bb7ee90617c85 --- /dev/null +++ b/test/freeflow/navierstokes/stokeschannel3dproblem.hh @@ -0,0 +1,266 @@ +// -*- 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 Channel flow test for the staggered grid (Navier-)Stokes model. + * + * The channel is either modeled in 3D or in 2D, using an additional wall friction term + * to mimic the 3D behavior of the flow. + * + */ +#ifndef DUMUX_3D_CHANNEL_PROBLEM_HH +#define DUMUX_3D_CHANNEL_PROBLEM_HH + +#include <dumux/material/fluidsystems/1pliquid.hh> +#include <dumux/material/components/constant.hh> + +#include <dumux/freeflow/navierstokes/problem.hh> +#include <dumux/discretization/staggered/freeflow/properties.hh> +#include <dumux/freeflow/navierstokes/model.hh> +#include <dune/common/float_cmp.hh> + +namespace Dumux +{ + + +template <class TypeTag> +class ThreeDChannelTestProblem; + +namespace Properties +{ +NEW_TYPE_TAG(ThreeDChannelTestTypeTag, INHERITS_FROM(StaggeredFreeFlowModel, NavierStokes)); + +// the fluid system +SET_PROP(ThreeDChannelTestTypeTag, FluidSystem) +{ + using Scalar = typename GET_PROP_TYPE(TypeTag, Scalar); + using type = FluidSystems::OnePLiquid<Scalar, Components::Constant<1, Scalar> >; +}; + +// Set the grid type +#if DIM_3D +SET_TYPE_PROP(ThreeDChannelTestTypeTag, Grid, Dune::YaspGrid<3>); +#else +SET_TYPE_PROP(ThreeDChannelTestTypeTag, Grid, Dune::YaspGrid<2>); +#endif + +// Set the problem property +SET_TYPE_PROP(ThreeDChannelTestTypeTag, Problem, ThreeDChannelTestProblem<TypeTag> ); + +SET_BOOL_PROP(ThreeDChannelTestTypeTag, EnableFVGridGeometryCache, true); +SET_BOOL_PROP(ThreeDChannelTestTypeTag, EnableGridFluxVariablesCache, true); +SET_BOOL_PROP(ThreeDChannelTestTypeTag, EnableGridVolumeVariablesCache, true); +} + +/*! + * \brief Test problem for the one-phase model: + \todo doc me! + */ +template <class TypeTag> +class ThreeDChannelTestProblem : public NavierStokesProblem<TypeTag> +{ + using ParentType = NavierStokesProblem<TypeTag>; + + using GridView = typename GET_PROP_TYPE(TypeTag, GridView); + using Scalar = typename GET_PROP_TYPE(TypeTag, Scalar); + + using Indices = typename GET_PROP_TYPE(TypeTag, ModelTraits)::Indices; + using BoundaryTypes = typename GET_PROP_TYPE(TypeTag, BoundaryTypes); + using Element = typename GridView::template Codim<0>::Entity; + + using FVGridGeometry = typename GET_PROP_TYPE(TypeTag, FVGridGeometry); + using FVElementGeometry = typename FVGridGeometry::LocalView; + using SubControlVolume = typename FVElementGeometry::SubControlVolume; + using SubControlVolumeFace = typename FVElementGeometry::SubControlVolumeFace; + + static constexpr int dim = GridView::dimension; + static constexpr int dimWorld = GridView::dimensionworld; + using GlobalPosition = Dune::FieldVector<Scalar, dimWorld>; + + using CellCenterPrimaryVariables = typename GET_PROP_TYPE(TypeTag, CellCenterPrimaryVariables); + using FacePrimaryVariables = typename GET_PROP_TYPE(TypeTag, FacePrimaryVariables); + + using PrimaryVariables = typename GET_PROP_TYPE(TypeTag, PrimaryVariables); + using NumEqVector = typename GET_PROP_TYPE(TypeTag, NumEqVector); + + static constexpr bool enablePseudoThreeDWallFriction = !DIM_3D; + +public: + ThreeDChannelTestProblem(std::shared_ptr<const FVGridGeometry> fvGridGeometry) + : ParentType(fvGridGeometry), eps_(1e-6) + { + deltaP_ = getParam<Scalar>("Problem.DeltaP"); + height_ = getParam<Scalar>("Problem.Height"); + rho_ = getParam<Scalar>("Component.LiquidDensity"); + nu_ = getParam<Scalar>("Component.LiquidKinematicViscosity"); + + if(dim == 3 && !Dune::FloatCmp::eq(height_, this->fvGridGeometry().bBoxMax()[2])) + DUNE_THROW(Dune::InvalidStateException, "z-dimension must equal height"); + + if(enablePseudoThreeDWallFriction) + extrusionFactor_ = 2.0/3.0 * height_; + else + extrusionFactor_ = 1.0; + } + + /*! + * \name Problem parameters + */ + // \{ + + /*! + * \brief Return the temperature within the domain in [K]. + * + * This problem assumes a temperature of 10 degrees Celsius. + */ + Scalar temperature() const + { return 273.15 + 10; } // 10C + + + /*! + * \brief Evaluate the source term for all phases within a given + * sub-control-volume face. + */ + using ParentType::source; + template<class ElementVolumeVariables, class ElementFaceVariables> + NumEqVector source(const Element &element, + const FVElementGeometry& fvGeometry, + const ElementVolumeVariables& elemVolVars, + const ElementFaceVariables& elemFaceVars, + const SubControlVolumeFace &scvf) const + { + auto source = NumEqVector(0.0); + +#if !DIM_3D + static const Scalar height = getParam<Scalar>("Problem.Height"); + static const Scalar factor = getParam<Scalar>("Problem.PseudoWallFractionFactor", 8.0); + source[scvf.directionIndex()] = this->pseudo3DWallFriction(scvf, elemVolVars, elemFaceVars, height, factor); +#endif + + return source; + } + + Scalar extrusionFactorAtPos(const GlobalPosition& pos) const + { return extrusionFactor_; } + + + // \} + /*! + * \name Boundary conditions + */ + // \{ + + /*! + * \brief Specifies which kind of boundary condition should be + * used for which equation on a given boundary control volume. + * + * \param globalPos The position of the center of the finite volume + */ + BoundaryTypes boundaryTypesAtPos(const GlobalPosition &globalPos) const + { + BoundaryTypes values; + + // set Dirichlet values for the velocity everywhere + values.setDirichlet(Indices::momentumXBalanceIdx); + values.setDirichlet(Indices::momentumYBalanceIdx); + if(dim == 3) + values.setDirichlet(Indices::momentumZBalanceIdx); + + // set a fixed pressure in one cell + if (isOutlet_(globalPos) || isInlet_(globalPos)) + { + values.setDirichlet(Indices::conti0EqIdx); + values.setOutflow(Indices::momentumXBalanceIdx); + values.setOutflow(Indices::momentumYBalanceIdx); + if(dim == 3) + values.setOutflow(Indices::momentumZBalanceIdx); + } + else + values.setOutflow(Indices::conti0EqIdx); + + return values; + } + + /*! + * \brief Evaluate the boundary conditions for a dirichlet + * control volume. + * + * \param globalPos The center of the finite volume which ought to be set. + */ + PrimaryVariables dirichletAtPos(const GlobalPosition &globalPos) const + { + PrimaryVariables values(0.0); + + if(isInlet_(globalPos)) + values[Indices::pressureIdx] = 1e5 + deltaP_; + + if(isOutlet_(globalPos)) + values[Indices::pressureIdx] = 1e5; + + return values; + } + + // \} + + /*! + * \brief Evaluate the initial value for a control volume. + * + * \param globalPos The global position + */ + PrimaryVariables initialAtPos(const GlobalPosition &globalPos) const + { + PrimaryVariables values(0.0); + return values; + } + + //! Returns the analytical solution for the flux through the rectangular channel + Scalar analyticalFlux() const + { + const Scalar h = height_; + const Scalar w = this->fvGridGeometry().bBoxMax()[1]; + const Scalar L = this->fvGridGeometry().bBoxMax()[0]; + + const Scalar mu = nu_*rho_; + + return h*h*h * w * deltaP_ / (12*mu*L) * (1.0 - 0.630 * h/w); + } + +private: + + bool isInlet_(const GlobalPosition& globalPos) const + { + return globalPos[0] < eps_; + } + + bool isOutlet_(const GlobalPosition& globalPos) const + { + return globalPos[0] > this->fvGridGeometry().bBoxMax()[0] - eps_; + } + + Scalar eps_; + Scalar deltaP_; + Scalar extrusionFactor_; + Scalar height_; + Scalar rho_; + Scalar nu_; +}; +} //end namespace + +#endif diff --git a/test/freeflow/navierstokes/test_stokes_channel_3d.cc b/test/freeflow/navierstokes/test_stokes_channel_3d.cc new file mode 100644 index 0000000000000000000000000000000000000000..bffaab16dac9f7384178ae88e1c7e042fc37540b --- /dev/null +++ b/test/freeflow/navierstokes/test_stokes_channel_3d.cc @@ -0,0 +1,288 @@ +// -*- 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/>. * + *****************************************************************************/ + #include <config.h> + + #include <ctime> + #include <iostream> + + #include <dune/common/parallel/mpihelper.hh> + #include <dune/common/timer.hh> + #include <dune/grid/io/file/dgfparser/dgfexception.hh> + #include <dune/grid/io/file/vtk.hh> + #include <dune/istl/io.hh> + + +#include "stokeschannel3dproblem.hh" + +#include <dumux/common/properties.hh> +#include <dumux/common/parameters.hh> +#include <dumux/common/valgrind.hh> +#include <dumux/common/dumuxmessage.hh> +#include <dumux/common/defaultusagemessage.hh> + +#include <dumux/linear/seqsolverbackend.hh> +#include <dumux/nonlinear/newtonsolver.hh> + +#include <dumux/assembly/staggeredfvassembler.hh> +#include <dumux/assembly/diffmethod.hh> + +#include <dumux/io/staggeredvtkoutputmodule.hh> +#include <dumux/freeflow/navierstokes/staggered/fluxoversurface.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 arguments 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" + "\t-SpatialParams.LensLowerLeftX x-coordinate of the lower left corner of the lens [m] \n" + "\t-SpatialParams.LensLowerLeftY y-coordinate of the lower left corner of the lens [m] \n" + "\t-SpatialParams.LensUpperRightX x-coordinate of the upper right corner of the lens [m] \n" + "\t-SpatialParams.LensUpperRightY y-coordinate of the upper right corner of the lens [m] \n" + "\t-SpatialParams.Permeability Permeability of the domain [m^2] \n" + "\t-SpatialParams.PermeabilityLens Permeability of the lens [m^2] \n"; + + std::cout << errorMessageOut + << "\n"; + } +} + +int main(int argc, char** argv) try +{ + using namespace Dumux; + + // define the type tag for this problem + using TypeTag = TTAG(ThreeDChannelTestTypeTag); + + // initialize MPI, finalize is done automatically on exit + const auto& mpiHelper = Dune::MPIHelper::instance(argc, argv); + + // print dumux start message + if (mpiHelper.rank() == 0) + DumuxMessage::print(/*firstCall=*/true); + + // parse command line arguments and input file + Parameters::init(argc, argv, usage); + + // try to create a grid (from the given grid file or the input file) + using GridCreator = typename GET_PROP_TYPE(TypeTag, GridCreator); + GridCreator::makeGrid(); + GridCreator::loadBalance(); + + //////////////////////////////////////////////////////////// + // run instationary non-linear problem on this grid + //////////////////////////////////////////////////////////// + + // we compute on the leaf grid view + const auto& leafGridView = GridCreator::grid().leafGridView(); + + // create the finite volume grid geometry + using FVGridGeometry = typename GET_PROP_TYPE(TypeTag, FVGridGeometry); + auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); + fvGridGeometry->update(); + + // the problem (boundary conditions) + using Problem = typename GET_PROP_TYPE(TypeTag, Problem); + auto problem = std::make_shared<Problem>(fvGridGeometry); + + // the solution vector + using SolutionVector = typename GET_PROP_TYPE(TypeTag, SolutionVector); + static constexpr auto cellCenterIdx = FVGridGeometry::cellCenterIdx(); + static constexpr auto faceIdx = FVGridGeometry::faceIdx(); + const auto numDofsCellCenter = leafGridView.size(0); + const auto numDofsFace = leafGridView.size(1); + SolutionVector x; + x[cellCenterIdx].resize(numDofsCellCenter); + x[faceIdx].resize(numDofsFace); + + // the grid variables + using GridVariables = typename GET_PROP_TYPE(TypeTag, GridVariables); + auto gridVariables = std::make_shared<GridVariables>(problem, fvGridGeometry); + gridVariables->init(x); + + // intialize the vtk output module + using VtkOutputFields = typename GET_PROP_TYPE(TypeTag, VtkOutputFields); + StaggeredVtkOutputModule<TypeTag> vtkWriter(*problem, *fvGridGeometry, *gridVariables, x, problem->name()); + VtkOutputFields::init(vtkWriter); //!< Add model specific output fields + vtkWriter.write(0.0); + + // the assembler with time loop for instationary problem + using Assembler = StaggeredFVAssembler<TypeTag, DiffMethod::numeric>; + auto assembler = std::make_shared<Assembler>(problem, fvGridGeometry, gridVariables); + + // the linear solver + using LinearSolver = Dumux::UMFPackBackend; + auto linearSolver = std::make_shared<LinearSolver>(); + + // the non-linear solver + using NewtonSolver = Dumux::NewtonSolver<Assembler, LinearSolver>; + NewtonSolver nonLinearSolver(assembler, linearSolver); + + // set up two planes over which fluxes are calculated + FluxOverSurface<TypeTag> flux(*assembler, x); + using GridView = typename GET_PROP_TYPE(TypeTag, GridView); + using Scalar = typename GET_PROP_TYPE(TypeTag, Scalar); + using GlobalPosition = Dune::FieldVector<Scalar, GridView::dimensionworld>; + + const Scalar xMin = fvGridGeometry->bBoxMin()[0]; + const Scalar xMax = fvGridGeometry->bBoxMax()[0]; + const Scalar yMin = fvGridGeometry->bBoxMin()[1]; + const Scalar yMax = fvGridGeometry->bBoxMax()[1]; +#if DIM_3D + const Scalar zMin = fvGridGeometry->bBoxMin()[2]; + const Scalar zMax = fvGridGeometry->bBoxMax()[2]; +#endif + + // The first plane shall be placed at the middle of the channel. + // If we have an odd number of cells in x-direction, there would not be any cell faces + // at the postion of the plane (which is required for the flux calculation). + // In this case, we add half a cell-width to the x-position in order to make sure that + // the cell faces lie on the plane. This assumes a regular cartesian grid. + // The second plane is placed at the outlet of the channel. +#if DIM_3D + const auto p0inlet = GlobalPosition{xMin, yMin, zMin}; + const auto p1inlet = GlobalPosition{xMin, yMax, zMin}; + const auto p2inlet = GlobalPosition{xMin, yMin, zMax}; + const auto p3inlet = GlobalPosition{xMin, yMax, zMax}; + flux.addSurface("inlet", p0inlet, p1inlet, p2inlet, p3inlet); +#else + const auto p0inlet = GlobalPosition{xMin, yMin}; + const auto p1inlet = GlobalPosition{xMin, yMax}; + flux.addSurface("inlet", p0inlet, p1inlet); +#endif + + const Scalar planePosMiddleX = xMin + 0.5*(xMax - xMin); + const int numCellsX = getParam<std::vector<int>>("Grid.Cells")[0]; + const Scalar offsetX = (numCellsX % 2 == 0) ? 0.0 : 0.5*((xMax - xMin) / numCellsX); + +#if DIM_3D + const auto p0middle = GlobalPosition{planePosMiddleX + offsetX, yMin, zMin}; + const auto p1middle = GlobalPosition{planePosMiddleX + offsetX, yMax, zMin}; + const auto p2middle = GlobalPosition{planePosMiddleX + offsetX, yMin, zMax}; + const auto p3middle = GlobalPosition{planePosMiddleX + offsetX, yMax, zMax}; + flux.addSurface("middle", p0middle, p1middle, p2middle, p3middle); +#else + const auto p0middle = GlobalPosition{planePosMiddleX + offsetX, yMin}; + const auto p1middle = GlobalPosition{planePosMiddleX + offsetX, yMax}; +flux.addSurface("middle", p0middle, p1middle); +#endif + + // The second plane is placed at the outlet of the channel. +#if DIM_3D + const auto p0outlet = GlobalPosition{xMax, yMin, zMin}; + const auto p1outlet = GlobalPosition{xMax, yMax, zMin}; + const auto p2outlet = GlobalPosition{xMax, yMin, zMax}; + const auto p3outlet = GlobalPosition{xMax, yMax, zMax}; + flux.addSurface("outlet", p0outlet, p1outlet, p2outlet, p3outlet); +#else + const auto p0outlet = GlobalPosition{xMax, yMin}; + const auto p1outlet = GlobalPosition{xMax, yMax}; + flux.addSurface("outlet", p0outlet, p1outlet); +#endif + + // linearize & solve + Dune::Timer timer; + nonLinearSolver.solve(x); + + // write vtk output + vtkWriter.write(1.0); + + // calculate and print mass fluxes over the planes + flux.calculateMassOrMoleFluxes(); + if(GET_PROP_TYPE(TypeTag, ModelTraits)::enableEnergyBalance()) + { + std::cout << "mass / energy flux at inlet is: " << flux.netFlux("inlet") << std::endl; + std::cout << "mass / energy flux at middle is: " << flux.netFlux("middle") << std::endl; + std::cout << "mass / energy flux at outlet is: " << flux.netFlux("outlet") << std::endl; + } + else + { + std::cout << "mass flux at inlet is: " << flux.netFlux("inlet") << std::endl; + std::cout << "mass flux at middle is: " << flux.netFlux("middle") << std::endl; + std::cout << "mass flux at outlet is: " << flux.netFlux("outlet") << std::endl; + } + + // calculate and print volume fluxes over the planes + flux.calculateVolumeFluxes(); + std::cout << "volume flux at inlet is: " << flux.netFlux("inlet")[0] << std::endl; + std::cout << "volume flux at middle is: " << flux.netFlux("middle")[0] << std::endl; + std::cout << "volume flux at outlet is: " << flux.netFlux("outlet")[0] << std::endl; + + + timer.stop(); + + std::cout << "analyticalFlux: " << problem->analyticalFlux() << std::endl; + + const auto& comm = Dune::MPIHelper::getCollectiveCommunication(); + std::cout << "Simulation took " << timer.elapsed() << " seconds on " + << comm.size() << " processes.\n" + << "The cumulative CPU time was " << timer.elapsed()*comm.size() << " seconds.\n"; + + + //////////////////////////////////////////////////////////// + // finalize, print dumux message to say goodbye + //////////////////////////////////////////////////////////// + + // print dumux end message + if (mpiHelper.rank() == 0) + { + Parameters::print(); + DumuxMessage::print(/*firstCall=*/false); + } + + return 0; +} // end main +catch (Dumux::ParameterException &e) +{ + std::cerr << std::endl << e << " ---> Abort!" << std::endl; + return 1; +} +catch (Dune::DGFException & e) +{ + std::cerr << "DGF exception thrown (" << e << + "). Most likely, the DGF file name is wrong " + "or the DGF file is corrupted, " + "e.g. missing hash at end of file or wrong number (dimensions) of entries." + << " ---> Abort!" << std::endl; + return 2; +} +catch (Dune::Exception &e) +{ + std::cerr << "Dune reported error: " << e << " ---> Abort!" << std::endl; + return 3; +} +catch (...) +{ + std::cerr << "Unknown exception thrown! ---> Abort!" << std::endl; + return 4; +} diff --git a/test/freeflow/navierstokes/test_stokes_channel_3d.input b/test/freeflow/navierstokes/test_stokes_channel_3d.input new file mode 100644 index 0000000000000000000000000000000000000000..ce28a56783b4b492a37b7bad3a5a4aa6c648addb --- /dev/null +++ b/test/freeflow/navierstokes/test_stokes_channel_3d.input @@ -0,0 +1,27 @@ +[Grid] +UpperRight = 0.005 0.002 0.00025 +Cells = 3 20 10 + +[Problem] +Name = test_stokes_channel_3d # name passed to the output routines +DeltaP = 10 +EnableGravity = false +Height = 0.00025 + +[Component] +LiquidDensity = 1e3 +LiquidKinematicViscosity = 1e-6 + +[Assembly] +NumericDifference.BaseEpsilon = 1e-8 + +[Newton] +MaxSteps = 10 +MaxRelativeShift = 1e-8 + +[Vtk] +AddVelocity = true +WriteFaceData = false + +[FluxOverPlane] +Verbose = true diff --git a/test/freeflow/navierstokes/test_stokes_channel_pseudo3d.input b/test/freeflow/navierstokes/test_stokes_channel_pseudo3d.input new file mode 100644 index 0000000000000000000000000000000000000000..462959a1c6ea09ca25b417ac58fc865659cf79c1 --- /dev/null +++ b/test/freeflow/navierstokes/test_stokes_channel_pseudo3d.input @@ -0,0 +1,25 @@ +[Grid] +UpperRight = 0.005 0.002 +Cells = 3 20 + +[Problem] +Name = test_stokes_channel_pseudo3d # name passed to the output routines +DeltaP = 10 +EnableGravity = false +Height = 0.00025 + +[Component] +LiquidDensity = 1e3 +LiquidKinematicViscosity = 1e-6 + +[Newton] +MaxSteps = 10 +MaxRelativeShift = 1e-8 + +[Assembly] +NumericDifference.BaseEpsilon = 1e-8 + + +[Vtk] +WriteFaceData = false +AddVelocity = true diff --git a/test/references/stokes_channel_3d-reference.vtu b/test/references/stokes_channel_3d-reference.vtu new file mode 100644 index 0000000000000000000000000000000000000000..6fb35f1d50050e938899bf0dbb36b78d1badba5d --- /dev/null +++ b/test/references/stokes_channel_3d-reference.vtu @@ -0,0 +1,1112 @@ +<?xml version="1.0"?> +<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian"> + <UnstructuredGrid> + <Piece NumberOfCells="600" NumberOfPoints="924"> + <CellData Scalars="p" Vectors="velocity_component (m/s)"> + <DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii"> + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + </DataArray> + <DataArray type="Float32" Name="rhoMolar" 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 1000 1000 1000 1000 + 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 + 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 + 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 + 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 + 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 + 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="rho" 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 1000 1000 1000 1000 + 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 + 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 + 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 + 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 + 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 + 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="velocity_component (m/s)" NumberOfComponents="3" format="ascii"> + 0.00187614 7.12886e-14 -1.76194e-14 0.00187614 2.46472e-14 -6.49934e-15 0.00187614 -1.70929e-14 3.30583e-15 0.00275429 3.74927e-14 2.60983e-14 + 0.00275429 -3.08739e-14 1.84864e-14 0.00275429 -6.96312e-14 7.04705e-15 0.00301195 -1.33332e-13 1.68702e-14 0.00301195 -4.75997e-14 -1.59778e-14 + 0.00301195 -5.02019e-14 -1.46059e-14 0.00309036 -2.9833e-15 -4.96805e-14 0.00309036 1.09858e-14 1.30187e-15 0.00309036 5.94208e-14 -1.34826e-14 + 0.00311437 6.34775e-14 3.31227e-14 0.00311437 3.48945e-14 -1.22321e-14 0.00311437 6.55749e-14 6.34636e-15 0.00312174 -4.03639e-14 -6.37102e-15 + 0.00312174 -8.24459e-15 1.65231e-14 0.00312174 2.05866e-14 -2.79133e-15 0.003124 -3.9773e-14 6.42325e-15 0.003124 -4.53623e-14 -9.12857e-15 + 0.003124 7.92291e-14 -1.41206e-14 0.00312469 -2.51161e-14 -1.01559e-14 0.00312469 1.2083e-13 -3.37338e-14 0.00312469 1.94482e-13 -1.54068e-14 + 0.0031249 8.99881e-14 -1.90756e-14 0.0031249 1.79763e-13 1.66617e-14 0.0031249 1.29638e-13 3.03661e-14 0.00312496 3.52859e-14 3.30418e-14 + 0.00312496 1.00381e-13 1.53798e-15 0.00312496 2.56286e-16 5.54243e-16 0.00312496 -7.64985e-14 -4.83255e-15 0.00312496 -1.01824e-13 4.60898e-14 + 0.00312496 -2.17596e-14 1.95749e-15 0.0031249 1.34939e-14 -1.79971e-14 0.0031249 -1.33717e-13 -4.19128e-14 0.0031249 1.69334e-14 -1.47088e-14 + 0.00312469 6.00216e-14 5.94536e-15 0.00312469 8.24641e-14 -1.43715e-14 0.00312469 5.24151e-14 4.99997e-15 0.003124 -8.34062e-14 3.00756e-14 + 0.003124 9.16937e-14 9.99309e-15 0.003124 5.68019e-14 -7.26045e-15 0.00312174 -1.06168e-13 -2.32739e-14 0.00312174 7.30968e-14 -7.41365e-15 + 0.00312174 4.8427e-14 6.95168e-15 0.00311437 -1.17348e-14 -4.47842e-16 0.00311437 4.90494e-14 1.13187e-14 0.00311437 8.25968e-14 -1.7257e-14 + 0.00309036 -2.94283e-14 4.08191e-15 0.00309036 -2.27261e-14 4.93497e-15 0.00309036 3.17448e-14 2.90541e-14 0.00301195 4.33373e-14 -2.29203e-14 + 0.00301195 3.73597e-14 -2.128e-14 0.00301195 -1.89696e-14 -1.71954e-14 0.00275429 1.05524e-13 6.89909e-15 0.00275429 1.10065e-13 1.36917e-15 + 0.00275429 -4.24738e-14 2.15064e-14 0.00187614 3.9095e-14 1.01013e-14 0.00187614 4.99535e-14 1.21e-14 0.00187614 -6.69758e-14 -1.76363e-14 + 0.00455804 1.77338e-14 -3.87629e-14 0.00455804 2.43883e-14 -1.93291e-14 0.00455804 -6.04549e-14 1.95251e-14 0.00705165 -3.96093e-14 7.00631e-14 + 0.00705165 -2.16851e-14 5.15753e-14 0.00705165 -1.02609e-13 6.06751e-15 0.00779704 -2.49778e-13 6.81288e-14 0.00779704 -3.99118e-15 -5.33063e-14 + 0.00779704 -5.17943e-14 -3.81483e-14 0.00802447 -2.01868e-14 -1.90804e-13 0.00802447 -2.78347e-14 3.16648e-14 0.00802447 3.41744e-14 -4.03413e-14 + 0.00809416 1.44428e-13 1.17512e-13 0.00809416 1.46376e-14 -6.80563e-14 0.00809416 4.28753e-14 1.73556e-14 0.00811554 2.63327e-14 -3.29761e-14 + 0.00811554 5.43873e-14 6.27692e-14 0.00811554 3.6842e-14 -1.52151e-14 0.0081221 5.11993e-15 3.88324e-14 0.0081221 -7.00333e-16 -3.36799e-14 + 0.0081221 1.70645e-13 -5.25927e-14 0.00812411 -8.16139e-14 -2.50439e-14 0.00812411 2.32326e-13 -1.12655e-13 0.00812411 3.50688e-13 -5.29881e-14 + 0.00812472 1.29021e-13 -8.73848e-14 0.00812472 3.09345e-13 5.51525e-14 0.00812472 2.57009e-13 1.04167e-13 0.00812489 6.76934e-14 1.30366e-13 + 0.00812489 1.46896e-13 1.82649e-14 0.00812489 4.02078e-14 9.19223e-15 0.00812489 -2.11977e-13 -4.19606e-15 0.00812489 -1.14286e-13 1.36523e-13 + 0.00812489 -2.73167e-14 7.44189e-15 0.00812472 -2.26667e-15 -9.45022e-14 0.00812472 -1.53665e-13 -1.25193e-13 0.00812472 8.83803e-14 -6.71579e-14 + 0.00812411 1.46385e-13 3.28813e-14 0.00812411 2.19507e-13 -8.59421e-14 0.00812411 2.17911e-13 1.24956e-14 0.0081221 -8.11523e-14 9.63644e-14 + 0.0081221 2.48733e-13 6.51375e-14 0.0081221 1.52938e-13 -3.68768e-15 0.00811554 -1.9816e-13 -5.25104e-14 0.00811554 1.37609e-13 -3.68362e-14 + 0.00811554 4.9101e-14 2.52331e-14 0.00809416 -1.17988e-13 -1.57045e-14 0.00809416 1.23576e-13 4.2965e-14 0.00809416 1.18081e-13 -6.68252e-14 + 0.00802447 -7.76178e-14 1.17102e-14 0.00802447 3.2678e-14 7.98572e-15 0.00802447 8.31106e-14 9.58897e-14 0.00779704 1.36251e-13 -1.0361e-13 + 0.00779704 1.09589e-13 -6.30694e-14 0.00779704 -5.4033e-15 -5.29122e-14 0.00705165 2.16253e-13 5.11512e-14 0.00705165 1.33611e-13 1.39448e-14 + 0.00705165 -4.29653e-14 6.7372e-14 0.00455804 3.39487e-14 2.93928e-14 0.00455804 4.04765e-14 3.30178e-14 0.00455804 -7.51779e-14 -5.6059e-14 + 0.00640385 -2.50306e-14 -3.43243e-14 0.00640385 2.83184e-14 -3.2196e-14 0.00640385 -1.13474e-13 5.84437e-14 0.0102083 -7.02534e-14 9.13174e-14 + 0.0102083 -6.20161e-14 9.33775e-14 0.0102083 -1.8696e-13 -1.47402e-14 0.0113645 -2.15854e-13 1.34028e-13 0.0113645 -8.12639e-14 -9.9154e-14 + 0.0113645 -9.52634e-14 -6.08745e-14 0.0117184 2.59701e-15 -3.66026e-13 0.0117184 6.17742e-15 6.67406e-14 0.0117184 -3.95425e-14 -5.54656e-14 + 0.011827 1.60385e-13 2.16417e-13 0.011827 4.18243e-14 -1.13749e-13 0.011827 -1.85366e-14 2.56392e-14 0.0118603 3.76634e-14 -6.81783e-14 + 0.0118603 1.00874e-13 9.83414e-14 0.0118603 1.17318e-13 -4.76015e-14 0.0118705 1.169e-14 8.72659e-14 0.0118705 1.20879e-13 -5.12269e-14 + 0.0118705 2.59761e-13 -8.32205e-14 0.0118736 -1.00674e-13 -2.46362e-14 0.0118736 2.13267e-13 -1.80949e-13 0.0118736 3.0159e-13 -8.13425e-14 + 0.0118746 1.1872e-13 -1.97804e-13 0.0118746 2.95068e-13 7.60078e-14 0.0118746 2.06319e-13 1.75405e-13 0.0118748 4.28208e-14 2.7255e-13 + 0.0118748 1.83415e-13 5.71773e-14 0.0118748 3.21285e-14 3.06554e-14 0.0118748 -2.48068e-13 -6.1439e-15 0.0118748 2.45648e-14 1.92349e-13 + 0.0118748 2.21155e-15 3.38071e-15 0.0118746 -2.25863e-15 -2.07561e-13 0.0118746 5.29212e-15 -1.77118e-13 0.0118746 1.59325e-13 -1.38455e-13 + 0.0118736 1.22652e-13 7.75678e-14 0.0118736 2.31298e-13 -1.96221e-13 0.0118736 2.28941e-13 2.63447e-14 0.0118705 -1.57666e-13 1.78311e-13 + 0.0118705 2.28481e-13 1.56533e-13 0.0118705 1.88709e-13 1.54803e-15 0.0118603 -2.57981e-13 -8.1474e-14 0.0118603 1.20243e-13 -8.48966e-14 + 0.0118603 1.28182e-13 5.49037e-14 0.011827 -1.35968e-13 -3.94335e-14 0.011827 1.08725e-13 8.59321e-14 0.011827 1.30929e-13 -1.21744e-13 + 0.0117184 -7.57147e-14 9.09604e-15 0.0117184 3.83531e-14 -4.32721e-15 0.0117184 6.81167e-14 1.66505e-13 0.0113645 1.02254e-13 -1.99671e-13 + 0.0113645 4.49644e-14 -7.92669e-14 0.0113645 1.23605e-15 -9.27548e-14 0.0102083 1.74554e-13 1.0923e-13 0.0102083 5.61488e-14 1.50895e-14 + 0.0102083 2.97328e-14 1.01751e-13 0.00640385 4.68377e-14 5.09889e-14 0.00640385 4.62533e-14 5.08435e-14 0.00640385 -8.74199e-15 -8.17532e-14 + 0.00756236 -1.56424e-14 -2.02269e-14 0.00756236 1.04939e-14 -3.99312e-14 0.00756236 -8.4579e-14 1.03544e-13 0.0122804 -6.74773e-14 1.0186e-13 + 0.0122804 -9.7484e-14 1.49799e-13 0.0122804 -1.37212e-13 -3.62739e-14 0.0137321 -1.06208e-13 1.65519e-13 0.0137321 -1.04337e-13 -1.50294e-13 + 0.0137321 -8.52579e-14 -7.91681e-14 0.0141777 7.08593e-14 -4.89341e-13 0.0141777 5.97154e-14 6.51825e-14 0.0141777 -5.73311e-14 -6.0403e-14 + 0.0143145 1.75302e-13 2.84073e-13 0.0143145 3.76834e-14 -6.54768e-14 0.0143145 5.02885e-14 3.35677e-14 0.0143564 6.58766e-14 -7.2917e-14 + 0.0143564 8.7557e-14 8.89017e-14 0.0143564 1.95221e-13 -6.84873e-14 0.0143693 -3.07679e-14 1.25052e-13 0.0143693 2.32897e-13 -5.80986e-14 + 0.0143693 2.60911e-13 -9.47315e-14 0.0143732 -3.39177e-14 -3.41956e-14 0.0143732 2.32945e-13 -2.01001e-13 0.0143732 2.61451e-13 -8.53046e-14 + 0.0143744 1.24e-13 -2.84293e-13 0.0143744 2.76935e-13 5.90227e-14 0.0143744 1.5294e-13 2.26157e-13 0.0143748 -2.07602e-15 4.02536e-13 + 0.0143748 2.15114e-13 1.06884e-13 0.0143748 5.03695e-14 4.57242e-14 0.0143748 -1.47832e-13 -3.41356e-14 0.0143748 2.04471e-13 1.78076e-13 + 0.0143748 8.46394e-14 -1.19077e-14 0.0143744 1.15973e-14 -2.81921e-13 0.0143744 2.50987e-13 -1.73684e-13 0.0143744 1.92872e-13 -1.89057e-13 + 0.0143732 4.5817e-15 1.23825e-13 0.0143732 2.12166e-13 -2.61865e-13 0.0143732 2.80424e-13 2.69972e-14 0.0143693 -1.68588e-13 2.42524e-13 + 0.0143693 2.20765e-13 2.03173e-13 0.0143693 2.54203e-13 8.03325e-15 0.0143564 -1.68457e-13 -1.29608e-13 0.0143564 2.0452e-13 -1.16148e-13 + 0.0143564 1.07719e-13 9.76606e-14 0.0143145 -2.41838e-14 -6.13775e-14 0.0143145 7.23042e-14 1.40737e-13 0.0143145 1.24705e-13 -1.75015e-13 + 0.0141777 -3.99148e-14 2.18632e-14 0.0141777 7.9061e-14 -5.39286e-14 0.0141777 7.7702e-14 2.35104e-13 0.0137321 -3.31388e-14 -2.57296e-13 + 0.0137321 8.57568e-14 -4.29392e-14 0.0137321 1.26877e-14 -1.40696e-13 0.0122804 4.6421e-14 1.31349e-13 0.0122804 6.1949e-14 -2.26797e-14 + 0.0122804 8.07374e-14 1.17688e-13 0.00756236 2.25298e-14 6.93795e-14 0.00756236 8.81775e-14 7.76671e-14 0.00756236 2.85604e-14 -8.28464e-14 + 0.00812128 7.43356e-15 -1.37273e-14 0.00812128 1.45548e-14 -4.31399e-14 0.00812128 -3.92156e-14 1.30752e-13 0.0133067 -3.63713e-14 1.20237e-13 + 0.0133067 -4.5809e-14 1.99219e-13 0.0133067 -7.38087e-14 -4.5258e-14 0.0149126 -1.1008e-14 1.4522e-13 0.0149126 -1.06524e-13 -1.75919e-13 + 0.0149126 -5.16997e-14 -8.39511e-14 0.0154064 1.00587e-13 -5.30305e-13 0.0154064 8.01272e-15 4.02261e-14 0.0154064 -6.0834e-14 -5.79314e-14 + 0.0155579 1.6513e-13 2.98364e-13 0.0155579 7.59611e-14 3.8158e-14 0.0155579 7.03421e-14 4.04805e-14 0.0156044 1.01149e-13 -3.70279e-14 + 0.0156044 2.28606e-13 4.72387e-14 0.0156044 2.865e-13 -6.47199e-14 0.0156187 6.2546e-15 1.3951e-13 0.0156187 3.75322e-13 -5.30269e-14 + 0.0156187 2.8966e-13 -8.60596e-14 0.0156231 7.39472e-14 -6.53881e-14 0.0156231 2.74316e-13 -1.84258e-13 0.0156231 2.31722e-13 -8.62208e-14 + 0.0156244 9.16463e-14 -2.97921e-13 0.0156244 2.67256e-13 3.15895e-14 0.0156244 1.89161e-13 2.61967e-13 0.0156248 -3.95663e-14 4.68972e-13 + 0.0156248 2.6772e-13 1.39695e-13 0.0156248 1.50469e-13 4.5513e-14 0.0156248 -6.26951e-14 -7.00955e-14 0.0156248 2.84921e-13 1.38895e-13 + 0.0156248 1.30412e-13 -8.04565e-15 0.0156244 -2.56234e-14 -2.96357e-13 0.0156244 2.94805e-13 -1.47635e-13 0.0156244 1.80537e-13 -2.26038e-13 + 0.0156231 -5.08925e-14 1.48512e-13 0.0156231 2.23525e-13 -2.77335e-13 0.0156231 3.35782e-13 -9.05548e-15 0.0156187 -5.09034e-14 2.5592e-13 + 0.0156187 2.11425e-13 1.98603e-13 0.0156187 2.58091e-13 5.89581e-14 0.0156044 -9.04661e-14 -1.4899e-13 0.0156044 1.87978e-13 -1.20459e-13 + 0.0156044 1.07535e-13 1.22054e-13 0.0155579 -1.69452e-14 -1.01496e-13 0.0155579 9.78502e-14 1.89087e-13 0.0155579 9.29791e-14 -2.01945e-13 + 0.0154064 -1.15421e-14 6.98304e-14 0.0154064 1.13428e-13 -1.17707e-13 0.0154064 6.29757e-15 2.82404e-13 0.0149126 -1.46657e-13 -2.69596e-13 + 0.0149126 8.02121e-14 1.58519e-14 0.0149126 -2.2778e-14 -1.79778e-13 0.0133067 -5.13135e-14 1.04988e-13 0.0133067 4.09643e-14 -6.87376e-14 + 0.0133067 9.47907e-14 1.02948e-13 0.00812128 3.29297e-14 8.34145e-14 0.00812128 6.02486e-14 1.05216e-13 0.00812128 9.01871e-14 -5.97045e-14 + 0.00812128 -2.30082e-15 -1.11728e-14 0.00812128 9.83696e-15 -4.69776e-14 0.00812128 4.0902e-14 1.26691e-13 0.0133067 4.29758e-14 1.18569e-13 + 0.0133067 2.2298e-14 2.20298e-13 0.0133067 -3.07939e-14 -1.68194e-14 0.0149126 1.33844e-13 1.13338e-13 0.0149126 -6.24395e-14 -1.45629e-13 + 0.0149126 -8.53351e-14 -9.21601e-14 0.0154064 5.74395e-14 -4.94639e-13 0.0154064 -1.38185e-15 -5.73415e-15 0.0154064 2.05233e-14 -6.40866e-14 + 0.0155579 4.03313e-14 2.68844e-13 0.0155579 2.30461e-13 1.12626e-13 0.0155579 1.46191e-13 6.24441e-14 0.0156044 6.28267e-14 1.07925e-14 + 0.0156044 3.45668e-13 3.29556e-14 0.0156044 2.81002e-13 -6.04222e-14 0.0156187 5.94537e-14 1.17727e-13 0.0156187 4.13698e-13 -5.23279e-14 + 0.0156187 2.4068e-13 -4.76697e-14 0.0156231 1.32852e-13 -8.03948e-14 0.0156231 3.23483e-13 -1.4085e-13 0.0156231 1.75021e-13 -1.01072e-13 + 0.0156244 1.45281e-15 -2.557e-13 0.0156244 2.54758e-13 9.02867e-15 0.0156244 2.169e-13 2.7591e-13 0.0156248 -5.36323e-14 4.5932e-13 + 0.0156248 2.91261e-13 1.44965e-13 0.0156248 2.18685e-13 4.41175e-14 0.0156248 7.81503e-14 -1.01466e-13 0.0156248 2.93326e-13 1.22872e-13 + 0.0156248 1.52601e-13 2.28943e-14 0.0156244 -4.01256e-15 -2.55167e-13 0.0156244 3.06493e-13 -1.36722e-13 0.0156244 1.22651e-13 -2.54953e-13 + 0.0156231 -7.51508e-14 1.33803e-13 0.0156231 2.02788e-13 -2.61859e-13 0.0156231 2.29354e-13 -5.84023e-14 0.0156187 4.79856e-14 2.34248e-13 + 0.0156187 1.41771e-13 1.80393e-13 0.0156187 2.14492e-13 1.19682e-13 0.0156044 1.94246e-14 -1.28408e-13 0.0156044 1.31352e-13 -1.13252e-13 + 0.0156044 1.06246e-13 1.28185e-13 0.0155579 8.18986e-14 -1.62093e-13 0.0155579 8.34543e-14 2.05542e-13 0.0155579 4.37887e-14 -1.87784e-13 + 0.0154064 2.76012e-14 1.49549e-13 0.0154064 7.34398e-14 -1.44504e-13 0.0154064 -8.93911e-15 2.91448e-13 0.0149126 -2.23799e-13 -2.47584e-13 + 0.0149126 2.34199e-14 5.07751e-14 0.0149126 -3.39893e-14 -1.9147e-13 0.0133067 -1.81788e-13 5.56597e-14 0.0133067 3.97408e-14 -1.01355e-13 + 0.0133067 8.33965e-14 4.79839e-14 0.00812128 -4.35158e-14 7.97847e-14 0.00812128 5.04956e-14 1.21982e-13 0.00812128 1.10616e-13 -1.55499e-14 + 0.00756236 -1.38598e-14 -4.81336e-15 0.00756236 -2.07736e-15 -4.9839e-14 0.00756236 7.04422e-14 9.40806e-14 0.0122804 -1.09093e-14 1.0103e-13 + 0.0122804 6.82766e-14 2.03139e-13 0.0122804 4.70626e-14 3.40675e-14 0.0137321 1.0918e-13 7.61844e-14 0.0137321 2.82506e-14 -8.14613e-14 + 0.0137321 -4.76545e-15 -1.05576e-13 0.0141777 8.38457e-15 -4.01374e-13 0.0141777 -9.32351e-16 -4.73454e-14 0.0141777 8.00428e-14 -8.48955e-14 + 0.0143145 6.445e-15 1.96084e-13 0.0143145 2.25131e-13 1.24878e-13 0.0143145 9.95269e-14 1.11655e-13 0.0143564 1.07924e-13 5.92374e-14 + 0.0143564 4.13581e-13 4.80922e-14 0.0143564 1.55656e-13 -6.56844e-14 0.0143693 6.29243e-14 8.14854e-14 0.0143693 4.229e-13 -5.24336e-14 + 0.0143693 1.21197e-13 5.34347e-15 0.0143732 1.12642e-13 -7.76065e-14 0.0143732 2.92386e-13 -8.94433e-14 0.0143732 3.61047e-14 -1.2287e-13 + 0.0143744 -3.69094e-14 -1.9084e-13 0.0143744 2.44427e-13 -1.1492e-14 0.0143744 1.61552e-13 2.55028e-13 0.0143748 -3.82341e-14 3.9533e-13 + 0.0143748 2.93375e-13 1.3739e-13 0.0143748 2.49133e-13 4.536e-14 0.0143748 1.61649e-13 -1.32928e-13 0.0143748 3.44771e-13 1.08514e-13 + 0.0143748 1.91068e-13 5.67566e-14 0.0143744 2.5422e-14 -1.70296e-13 0.0143744 3.32718e-13 -1.26162e-13 0.0143744 1.35723e-13 -2.64439e-13 + 0.0143732 -9.89592e-14 9.98331e-14 0.0143732 1.57021e-13 -2.19197e-13 0.0143732 8.61607e-14 -7.51174e-14 0.0143693 1.08308e-13 1.8083e-13 + 0.0143693 9.49249e-14 1.49908e-13 0.0143693 1.0891e-14 1.48238e-13 0.0143564 1.2064e-13 -8.58699e-14 0.0143564 9.96385e-14 -9.91337e-14 + 0.0143564 3.88395e-14 1.19984e-13 0.0143145 2.32222e-14 -2.00884e-13 0.0143145 6.40557e-14 2.00802e-13 0.0143145 3.35193e-14 -1.62271e-13 + 0.0141777 -1.42178e-14 2.17584e-13 0.0141777 4.8154e-14 -1.43051e-13 0.0141777 5.89744e-14 2.63028e-13 0.0137321 -1.04872e-13 -2.24963e-13 + 0.0137321 -3.86766e-15 6.16778e-14 0.0137321 9.41776e-14 -1.8033e-13 0.0122804 -1.39089e-13 3.87351e-14 0.0122804 9.33901e-16 -1.22611e-13 + 0.0122804 4.0731e-17 2.26882e-14 0.00756236 -7.40836e-14 4.86027e-14 0.00756236 1.46314e-14 1.27291e-13 0.00756236 -1.92729e-14 1.95366e-15 + 0.00640385 -5.49518e-15 7.8071e-16 0.00640385 -1.58899e-15 -5.31345e-14 0.00640385 8.40704e-14 4.94955e-14 0.0102083 3.9448e-14 8.29061e-14 + 0.0102083 6.45228e-14 1.6537e-13 0.0102083 1.25001e-13 6.40004e-14 0.0113645 1.38468e-13 3.96649e-14 0.0113645 5.8878e-14 -2.50071e-14 + 0.0113645 5.2214e-14 -1.02687e-13 0.0117184 -9.22502e-14 -2.73083e-13 0.0117184 -2.49196e-16 -5.82961e-14 0.0117184 -1.58167e-14 -8.16047e-14 + 0.011827 -1.69153e-13 9.75828e-14 0.011827 1.68098e-13 8.80043e-14 0.011827 1.5137e-15 1.36156e-13 0.0118603 9.40287e-15 9.14204e-14 + 0.0118603 3.60768e-13 4.76444e-14 0.0118603 5.06661e-14 -6.65054e-14 0.0118705 9.47074e-14 3.83194e-14 0.0118705 3.68296e-13 -3.52804e-14 + 0.0118705 4.33022e-14 3.37823e-14 0.0118736 1.0591e-13 -5.27393e-14 0.0118736 3.01905e-13 -6.06465e-14 0.0118736 -4.6002e-14 -1.1264e-13 + 0.0118746 -9.4447e-14 -1.32359e-13 0.0118746 2.17936e-13 -8.94871e-15 0.0118746 1.73991e-14 1.94764e-13 0.0118748 2.33094e-14 2.97853e-13 + 0.0118748 2.88561e-13 9.80936e-14 0.0118748 2.09425e-13 3.43646e-14 0.0118748 1.98795e-13 -1.3804e-13 0.0118748 3.85472e-13 1.0031e-13 + 0.0118748 2.51547e-13 7.03045e-14 0.0118746 7.96938e-15 -8.43082e-14 0.0118746 3.18918e-13 -1.06415e-13 0.0118746 1.21164e-13 -2.33546e-13 + 0.0118736 -5.65347e-14 6.2126e-14 0.0118736 1.40656e-13 -1.64966e-13 0.0118736 -1.08776e-13 -4.61156e-14 0.0118705 1.40604e-13 1.13573e-13 + 0.0118705 4.34074e-14 1.20215e-13 0.0118705 -1.37127e-13 1.36256e-13 0.0118603 1.68921e-13 -3.73067e-14 0.0118603 3.95922e-14 -8.4683e-14 + 0.0118603 3.74498e-14 7.77563e-14 0.011827 -8.67915e-15 -1.82872e-13 0.011827 4.71379e-14 1.80314e-13 0.011827 2.65486e-14 -1.19777e-13 + 0.0117184 1.57645e-14 2.06996e-13 0.0117184 6.80109e-14 -1.35745e-13 0.0117184 7.43914e-14 1.93467e-13 0.0113645 -2.51033e-14 -1.77473e-13 + 0.0113645 5.61573e-14 5.75775e-14 0.0113645 1.32201e-13 -1.45772e-13 0.0102083 -1.49433e-13 3.74736e-14 0.0102083 -6.33557e-14 -9.66112e-14 + 0.0102083 -1.29714e-14 4.01783e-14 0.00640385 -5.90166e-14 1.3644e-14 0.00640385 -9.91356e-14 9.62963e-14 0.00640385 -2.59198e-14 -1.38536e-14 + 0.00455804 -1.11364e-14 4.92061e-15 0.00455804 -6.06081e-14 -4.26763e-14 0.00455804 4.4909e-14 1.14222e-14 0.00705165 7.90152e-14 4.29541e-14 + 0.00705165 1.2653e-14 1.07899e-13 0.00705165 1.28868e-13 5.8301e-14 0.00779704 2.00748e-13 2.45361e-14 0.00779704 1.45677e-13 -3.29791e-15 + 0.00779704 5.12019e-14 -6.66959e-14 0.00802447 -3.36283e-14 -1.37268e-13 0.00802447 9.52499e-14 -4.07472e-14 0.00802447 -1.20676e-13 -5.28909e-14 + 0.00809416 -2.11892e-13 2.94264e-14 0.00809416 1.14104e-13 4.22377e-14 0.00809416 -2.68438e-14 9.14975e-14 0.00811554 -1.89349e-14 6.77487e-14 + 0.00811554 2.33574e-13 3.15598e-14 0.00811554 4.10201e-14 -3.93907e-14 0.0081221 1.47375e-13 -1.73829e-15 0.0081221 2.64027e-13 -2.13391e-14 + 0.0081221 -1.93521e-14 2.79067e-14 0.00812411 1.12954e-13 -8.99572e-15 0.00812411 1.85273e-13 -4.03612e-14 0.00812411 -1.18018e-13 -6.3687e-14 + 0.00812472 -8.00052e-14 -8.20985e-14 0.00812472 8.9538e-14 9.53826e-15 0.00812472 -4.70705e-14 1.05898e-13 0.00812489 8.12641e-14 1.72407e-13 + 0.00812489 1.70173e-13 3.42035e-14 0.00812489 1.62831e-13 1.64773e-14 0.00812489 1.91555e-13 -8.82921e-14 0.00812489 3.57367e-13 8.36596e-14 + 0.00812489 2.73526e-13 4.47823e-14 0.00812472 -2.70813e-14 -3.25579e-14 0.00812472 3.08754e-13 -7.11125e-14 0.00812472 1.04665e-13 -1.38083e-13 + 0.00812411 -3.49352e-14 2.84993e-14 0.00812411 3.89017e-14 -9.93698e-14 0.00812411 -2.06555e-13 -1.343e-14 0.0081221 1.48821e-13 4.9146e-14 + 0.0081221 -3.43757e-14 8.52822e-14 0.0081221 -8.72073e-14 7.31487e-14 0.00811554 2.43176e-13 -6.14978e-15 0.00811554 -1.89003e-14 -6.43412e-14 + 0.00811554 1.16221e-13 3.93252e-14 0.00809416 2.67495e-14 -1.14903e-13 0.00809416 6.523e-14 1.2413e-13 0.00809416 2.30712e-14 -6.26584e-14 + 0.00802447 3.43954e-14 1.3313e-13 0.00802447 6.56934e-14 -9.75439e-14 0.00802447 8.43599e-14 1.0117e-13 0.00779704 3.03481e-14 -8.95574e-14 + 0.00779704 -3.90542e-14 3.78263e-14 0.00779704 9.5015e-14 -7.8882e-14 0.00705165 -5.80149e-14 7.37755e-15 0.00705165 -8.37593e-14 -4.27035e-14 + 0.00705165 -1.91311e-16 2.73422e-14 0.00455804 2.17582e-14 3.48814e-15 0.00455804 -8.54373e-14 4.28471e-14 0.00455804 4.06358e-14 -1.3424e-14 + 0.00187614 1.50396e-14 3.81816e-15 0.00187614 -6.42425e-14 -1.4978e-14 0.00187614 -9.65387e-15 -1.25142e-15 0.00275429 6.04947e-14 8.23813e-15 + 0.00275429 2.68966e-15 3.50792e-14 0.00275429 6.34455e-14 2.24426e-14 0.00301195 1.32601e-13 1.0188e-14 0.00301195 1.2196e-13 -1.57464e-15 + 0.00301195 6.0947e-14 -1.95395e-14 0.00309036 2.48171e-14 -3.67858e-14 0.00309036 6.26504e-14 -1.26411e-14 0.00309036 -9.18672e-14 -1.82303e-14 + 0.00311437 -1.02328e-13 5.40348e-15 0.00311437 5.35003e-14 1.21955e-14 0.00311437 -6.72171e-14 2.60954e-14 0.00312174 -4.65235e-15 1.9089e-14 + 0.00312174 1.24499e-13 9.48428e-15 0.00312174 -1.86907e-14 -9.55771e-15 0.003124 4.37655e-14 -6.98462e-15 0.003124 1.26624e-13 -7.7681e-15 + 0.003124 -2.06602e-14 1.10611e-14 0.00312469 3.77741e-14 5.59117e-15 0.00312469 3.91756e-14 -1.35933e-14 0.00312469 -5.14837e-14 -1.75225e-14 + 0.0031249 -5.48191e-14 -2.80026e-14 0.0031249 5.24657e-15 8.05006e-15 0.0031249 -2.29756e-14 2.76671e-14 0.00312496 4.08777e-14 5.22257e-14 + 0.00312496 3.53083e-14 1.79597e-15 0.00312496 9.64928e-14 5.1614e-15 0.00312496 1.5123e-13 -2.44551e-14 0.00312496 1.63936e-13 3.17264e-14 + 0.00312496 1.5197e-13 1.00986e-14 0.0031249 1.66187e-14 -8.87071e-15 0.0031249 1.9396e-13 -2.19263e-14 0.0031249 3.95132e-14 -3.70088e-14 + 0.00312469 1.12404e-14 7.72238e-15 0.00312469 -2.86002e-14 -3.18616e-14 0.00312469 -1.17898e-13 -1.23313e-15 0.003124 7.13734e-14 7.57177e-15 + 0.003124 -3.91903e-14 3.14008e-14 0.003124 -7.28066e-14 1.43484e-14 0.00312174 1.1176e-13 2.07584e-15 0.00312174 -2.35946e-14 -2.53256e-14 + 0.00312174 3.22859e-14 1.43842e-14 0.00311437 -1.80553e-14 -3.51184e-14 0.00311437 3.10906e-14 4.17118e-14 0.00311437 1.49647e-14 -1.64739e-14 + 0.00309036 1.63837e-14 4.36035e-14 0.00309036 6.12373e-14 -3.15022e-14 0.00309036 4.75562e-14 2.62508e-14 0.00301195 1.09323e-13 -2.07452e-14 + 0.00301195 -2.21877e-14 1.24653e-14 0.00301195 7.50045e-14 -1.80893e-14 0.00275429 -4.41988e-15 -8.40357e-15 0.00275429 -2.13876e-14 -1.08917e-14 + 0.00275429 1.11723e-14 3.01349e-15 0.00187614 -1.81177e-14 4.3947e-15 0.00187614 -3.33164e-14 9.33429e-15 0.00187614 1.09539e-14 -2.25547e-15 + </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 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + </DataArray> + </CellData> + <Points> + <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii"> + 0 0 0 0.00166667 0 0 0 0.0001 0 0.00166667 0.0001 0 + 0 0 2.5e-05 0.00166667 0 2.5e-05 0 0.0001 2.5e-05 0.00166667 0.0001 2.5e-05 + 0.00333333 0 0 0.00333333 0.0001 0 0.00333333 0 2.5e-05 0.00333333 0.0001 2.5e-05 + 0.005 0 0 0.005 0.0001 0 0.005 0 2.5e-05 0.005 0.0001 2.5e-05 + 0 0.0002 0 0.00166667 0.0002 0 0 0.0002 2.5e-05 0.00166667 0.0002 2.5e-05 + 0.00333333 0.0002 0 0.00333333 0.0002 2.5e-05 0.005 0.0002 0 0.005 0.0002 2.5e-05 + 0 0.0003 0 0.00166667 0.0003 0 0 0.0003 2.5e-05 0.00166667 0.0003 2.5e-05 + 0.00333333 0.0003 0 0.00333333 0.0003 2.5e-05 0.005 0.0003 0 0.005 0.0003 2.5e-05 + 0 0.0004 0 0.00166667 0.0004 0 0 0.0004 2.5e-05 0.00166667 0.0004 2.5e-05 + 0.00333333 0.0004 0 0.00333333 0.0004 2.5e-05 0.005 0.0004 0 0.005 0.0004 2.5e-05 + 0 0.0005 0 0.00166667 0.0005 0 0 0.0005 2.5e-05 0.00166667 0.0005 2.5e-05 + 0.00333333 0.0005 0 0.00333333 0.0005 2.5e-05 0.005 0.0005 0 0.005 0.0005 2.5e-05 + 0 0.0006 0 0.00166667 0.0006 0 0 0.0006 2.5e-05 0.00166667 0.0006 2.5e-05 + 0.00333333 0.0006 0 0.00333333 0.0006 2.5e-05 0.005 0.0006 0 0.005 0.0006 2.5e-05 + 0 0.0007 0 0.00166667 0.0007 0 0 0.0007 2.5e-05 0.00166667 0.0007 2.5e-05 + 0.00333333 0.0007 0 0.00333333 0.0007 2.5e-05 0.005 0.0007 0 0.005 0.0007 2.5e-05 + 0 0.0008 0 0.00166667 0.0008 0 0 0.0008 2.5e-05 0.00166667 0.0008 2.5e-05 + 0.00333333 0.0008 0 0.00333333 0.0008 2.5e-05 0.005 0.0008 0 0.005 0.0008 2.5e-05 + 0 0.0009 0 0.00166667 0.0009 0 0 0.0009 2.5e-05 0.00166667 0.0009 2.5e-05 + 0.00333333 0.0009 0 0.00333333 0.0009 2.5e-05 0.005 0.0009 0 0.005 0.0009 2.5e-05 + 0 0.001 0 0.00166667 0.001 0 0 0.001 2.5e-05 0.00166667 0.001 2.5e-05 + 0.00333333 0.001 0 0.00333333 0.001 2.5e-05 0.005 0.001 0 0.005 0.001 2.5e-05 + 0 0.0011 0 0.00166667 0.0011 0 0 0.0011 2.5e-05 0.00166667 0.0011 2.5e-05 + 0.00333333 0.0011 0 0.00333333 0.0011 2.5e-05 0.005 0.0011 0 0.005 0.0011 2.5e-05 + 0 0.0012 0 0.00166667 0.0012 0 0 0.0012 2.5e-05 0.00166667 0.0012 2.5e-05 + 0.00333333 0.0012 0 0.00333333 0.0012 2.5e-05 0.005 0.0012 0 0.005 0.0012 2.5e-05 + 0 0.0013 0 0.00166667 0.0013 0 0 0.0013 2.5e-05 0.00166667 0.0013 2.5e-05 + 0.00333333 0.0013 0 0.00333333 0.0013 2.5e-05 0.005 0.0013 0 0.005 0.0013 2.5e-05 + 0 0.0014 0 0.00166667 0.0014 0 0 0.0014 2.5e-05 0.00166667 0.0014 2.5e-05 + 0.00333333 0.0014 0 0.00333333 0.0014 2.5e-05 0.005 0.0014 0 0.005 0.0014 2.5e-05 + 0 0.0015 0 0.00166667 0.0015 0 0 0.0015 2.5e-05 0.00166667 0.0015 2.5e-05 + 0.00333333 0.0015 0 0.00333333 0.0015 2.5e-05 0.005 0.0015 0 0.005 0.0015 2.5e-05 + 0 0.0016 0 0.00166667 0.0016 0 0 0.0016 2.5e-05 0.00166667 0.0016 2.5e-05 + 0.00333333 0.0016 0 0.00333333 0.0016 2.5e-05 0.005 0.0016 0 0.005 0.0016 2.5e-05 + 0 0.0017 0 0.00166667 0.0017 0 0 0.0017 2.5e-05 0.00166667 0.0017 2.5e-05 + 0.00333333 0.0017 0 0.00333333 0.0017 2.5e-05 0.005 0.0017 0 0.005 0.0017 2.5e-05 + 0 0.0018 0 0.00166667 0.0018 0 0 0.0018 2.5e-05 0.00166667 0.0018 2.5e-05 + 0.00333333 0.0018 0 0.00333333 0.0018 2.5e-05 0.005 0.0018 0 0.005 0.0018 2.5e-05 + 0 0.0019 0 0.00166667 0.0019 0 0 0.0019 2.5e-05 0.00166667 0.0019 2.5e-05 + 0.00333333 0.0019 0 0.00333333 0.0019 2.5e-05 0.005 0.0019 0 0.005 0.0019 2.5e-05 + 0 0.002 0 0.00166667 0.002 0 0 0.002 2.5e-05 0.00166667 0.002 2.5e-05 + 0.00333333 0.002 0 0.00333333 0.002 2.5e-05 0.005 0.002 0 0.005 0.002 2.5e-05 + 0 0 5e-05 0.00166667 0 5e-05 0 0.0001 5e-05 0.00166667 0.0001 5e-05 + 0.00333333 0 5e-05 0.00333333 0.0001 5e-05 0.005 0 5e-05 0.005 0.0001 5e-05 + 0 0.0002 5e-05 0.00166667 0.0002 5e-05 0.00333333 0.0002 5e-05 0.005 0.0002 5e-05 + 0 0.0003 5e-05 0.00166667 0.0003 5e-05 0.00333333 0.0003 5e-05 0.005 0.0003 5e-05 + 0 0.0004 5e-05 0.00166667 0.0004 5e-05 0.00333333 0.0004 5e-05 0.005 0.0004 5e-05 + 0 0.0005 5e-05 0.00166667 0.0005 5e-05 0.00333333 0.0005 5e-05 0.005 0.0005 5e-05 + 0 0.0006 5e-05 0.00166667 0.0006 5e-05 0.00333333 0.0006 5e-05 0.005 0.0006 5e-05 + 0 0.0007 5e-05 0.00166667 0.0007 5e-05 0.00333333 0.0007 5e-05 0.005 0.0007 5e-05 + 0 0.0008 5e-05 0.00166667 0.0008 5e-05 0.00333333 0.0008 5e-05 0.005 0.0008 5e-05 + 0 0.0009 5e-05 0.00166667 0.0009 5e-05 0.00333333 0.0009 5e-05 0.005 0.0009 5e-05 + 0 0.001 5e-05 0.00166667 0.001 5e-05 0.00333333 0.001 5e-05 0.005 0.001 5e-05 + 0 0.0011 5e-05 0.00166667 0.0011 5e-05 0.00333333 0.0011 5e-05 0.005 0.0011 5e-05 + 0 0.0012 5e-05 0.00166667 0.0012 5e-05 0.00333333 0.0012 5e-05 0.005 0.0012 5e-05 + 0 0.0013 5e-05 0.00166667 0.0013 5e-05 0.00333333 0.0013 5e-05 0.005 0.0013 5e-05 + 0 0.0014 5e-05 0.00166667 0.0014 5e-05 0.00333333 0.0014 5e-05 0.005 0.0014 5e-05 + 0 0.0015 5e-05 0.00166667 0.0015 5e-05 0.00333333 0.0015 5e-05 0.005 0.0015 5e-05 + 0 0.0016 5e-05 0.00166667 0.0016 5e-05 0.00333333 0.0016 5e-05 0.005 0.0016 5e-05 + 0 0.0017 5e-05 0.00166667 0.0017 5e-05 0.00333333 0.0017 5e-05 0.005 0.0017 5e-05 + 0 0.0018 5e-05 0.00166667 0.0018 5e-05 0.00333333 0.0018 5e-05 0.005 0.0018 5e-05 + 0 0.0019 5e-05 0.00166667 0.0019 5e-05 0.00333333 0.0019 5e-05 0.005 0.0019 5e-05 + 0 0.002 5e-05 0.00166667 0.002 5e-05 0.00333333 0.002 5e-05 0.005 0.002 5e-05 + 0 0 7.5e-05 0.00166667 0 7.5e-05 0 0.0001 7.5e-05 0.00166667 0.0001 7.5e-05 + 0.00333333 0 7.5e-05 0.00333333 0.0001 7.5e-05 0.005 0 7.5e-05 0.005 0.0001 7.5e-05 + 0 0.0002 7.5e-05 0.00166667 0.0002 7.5e-05 0.00333333 0.0002 7.5e-05 0.005 0.0002 7.5e-05 + 0 0.0003 7.5e-05 0.00166667 0.0003 7.5e-05 0.00333333 0.0003 7.5e-05 0.005 0.0003 7.5e-05 + 0 0.0004 7.5e-05 0.00166667 0.0004 7.5e-05 0.00333333 0.0004 7.5e-05 0.005 0.0004 7.5e-05 + 0 0.0005 7.5e-05 0.00166667 0.0005 7.5e-05 0.00333333 0.0005 7.5e-05 0.005 0.0005 7.5e-05 + 0 0.0006 7.5e-05 0.00166667 0.0006 7.5e-05 0.00333333 0.0006 7.5e-05 0.005 0.0006 7.5e-05 + 0 0.0007 7.5e-05 0.00166667 0.0007 7.5e-05 0.00333333 0.0007 7.5e-05 0.005 0.0007 7.5e-05 + 0 0.0008 7.5e-05 0.00166667 0.0008 7.5e-05 0.00333333 0.0008 7.5e-05 0.005 0.0008 7.5e-05 + 0 0.0009 7.5e-05 0.00166667 0.0009 7.5e-05 0.00333333 0.0009 7.5e-05 0.005 0.0009 7.5e-05 + 0 0.001 7.5e-05 0.00166667 0.001 7.5e-05 0.00333333 0.001 7.5e-05 0.005 0.001 7.5e-05 + 0 0.0011 7.5e-05 0.00166667 0.0011 7.5e-05 0.00333333 0.0011 7.5e-05 0.005 0.0011 7.5e-05 + 0 0.0012 7.5e-05 0.00166667 0.0012 7.5e-05 0.00333333 0.0012 7.5e-05 0.005 0.0012 7.5e-05 + 0 0.0013 7.5e-05 0.00166667 0.0013 7.5e-05 0.00333333 0.0013 7.5e-05 0.005 0.0013 7.5e-05 + 0 0.0014 7.5e-05 0.00166667 0.0014 7.5e-05 0.00333333 0.0014 7.5e-05 0.005 0.0014 7.5e-05 + 0 0.0015 7.5e-05 0.00166667 0.0015 7.5e-05 0.00333333 0.0015 7.5e-05 0.005 0.0015 7.5e-05 + 0 0.0016 7.5e-05 0.00166667 0.0016 7.5e-05 0.00333333 0.0016 7.5e-05 0.005 0.0016 7.5e-05 + 0 0.0017 7.5e-05 0.00166667 0.0017 7.5e-05 0.00333333 0.0017 7.5e-05 0.005 0.0017 7.5e-05 + 0 0.0018 7.5e-05 0.00166667 0.0018 7.5e-05 0.00333333 0.0018 7.5e-05 0.005 0.0018 7.5e-05 + 0 0.0019 7.5e-05 0.00166667 0.0019 7.5e-05 0.00333333 0.0019 7.5e-05 0.005 0.0019 7.5e-05 + 0 0.002 7.5e-05 0.00166667 0.002 7.5e-05 0.00333333 0.002 7.5e-05 0.005 0.002 7.5e-05 + 0 0 0.0001 0.00166667 0 0.0001 0 0.0001 0.0001 0.00166667 0.0001 0.0001 + 0.00333333 0 0.0001 0.00333333 0.0001 0.0001 0.005 0 0.0001 0.005 0.0001 0.0001 + 0 0.0002 0.0001 0.00166667 0.0002 0.0001 0.00333333 0.0002 0.0001 0.005 0.0002 0.0001 + 0 0.0003 0.0001 0.00166667 0.0003 0.0001 0.00333333 0.0003 0.0001 0.005 0.0003 0.0001 + 0 0.0004 0.0001 0.00166667 0.0004 0.0001 0.00333333 0.0004 0.0001 0.005 0.0004 0.0001 + 0 0.0005 0.0001 0.00166667 0.0005 0.0001 0.00333333 0.0005 0.0001 0.005 0.0005 0.0001 + 0 0.0006 0.0001 0.00166667 0.0006 0.0001 0.00333333 0.0006 0.0001 0.005 0.0006 0.0001 + 0 0.0007 0.0001 0.00166667 0.0007 0.0001 0.00333333 0.0007 0.0001 0.005 0.0007 0.0001 + 0 0.0008 0.0001 0.00166667 0.0008 0.0001 0.00333333 0.0008 0.0001 0.005 0.0008 0.0001 + 0 0.0009 0.0001 0.00166667 0.0009 0.0001 0.00333333 0.0009 0.0001 0.005 0.0009 0.0001 + 0 0.001 0.0001 0.00166667 0.001 0.0001 0.00333333 0.001 0.0001 0.005 0.001 0.0001 + 0 0.0011 0.0001 0.00166667 0.0011 0.0001 0.00333333 0.0011 0.0001 0.005 0.0011 0.0001 + 0 0.0012 0.0001 0.00166667 0.0012 0.0001 0.00333333 0.0012 0.0001 0.005 0.0012 0.0001 + 0 0.0013 0.0001 0.00166667 0.0013 0.0001 0.00333333 0.0013 0.0001 0.005 0.0013 0.0001 + 0 0.0014 0.0001 0.00166667 0.0014 0.0001 0.00333333 0.0014 0.0001 0.005 0.0014 0.0001 + 0 0.0015 0.0001 0.00166667 0.0015 0.0001 0.00333333 0.0015 0.0001 0.005 0.0015 0.0001 + 0 0.0016 0.0001 0.00166667 0.0016 0.0001 0.00333333 0.0016 0.0001 0.005 0.0016 0.0001 + 0 0.0017 0.0001 0.00166667 0.0017 0.0001 0.00333333 0.0017 0.0001 0.005 0.0017 0.0001 + 0 0.0018 0.0001 0.00166667 0.0018 0.0001 0.00333333 0.0018 0.0001 0.005 0.0018 0.0001 + 0 0.0019 0.0001 0.00166667 0.0019 0.0001 0.00333333 0.0019 0.0001 0.005 0.0019 0.0001 + 0 0.002 0.0001 0.00166667 0.002 0.0001 0.00333333 0.002 0.0001 0.005 0.002 0.0001 + 0 0 0.000125 0.00166667 0 0.000125 0 0.0001 0.000125 0.00166667 0.0001 0.000125 + 0.00333333 0 0.000125 0.00333333 0.0001 0.000125 0.005 0 0.000125 0.005 0.0001 0.000125 + 0 0.0002 0.000125 0.00166667 0.0002 0.000125 0.00333333 0.0002 0.000125 0.005 0.0002 0.000125 + 0 0.0003 0.000125 0.00166667 0.0003 0.000125 0.00333333 0.0003 0.000125 0.005 0.0003 0.000125 + 0 0.0004 0.000125 0.00166667 0.0004 0.000125 0.00333333 0.0004 0.000125 0.005 0.0004 0.000125 + 0 0.0005 0.000125 0.00166667 0.0005 0.000125 0.00333333 0.0005 0.000125 0.005 0.0005 0.000125 + 0 0.0006 0.000125 0.00166667 0.0006 0.000125 0.00333333 0.0006 0.000125 0.005 0.0006 0.000125 + 0 0.0007 0.000125 0.00166667 0.0007 0.000125 0.00333333 0.0007 0.000125 0.005 0.0007 0.000125 + 0 0.0008 0.000125 0.00166667 0.0008 0.000125 0.00333333 0.0008 0.000125 0.005 0.0008 0.000125 + 0 0.0009 0.000125 0.00166667 0.0009 0.000125 0.00333333 0.0009 0.000125 0.005 0.0009 0.000125 + 0 0.001 0.000125 0.00166667 0.001 0.000125 0.00333333 0.001 0.000125 0.005 0.001 0.000125 + 0 0.0011 0.000125 0.00166667 0.0011 0.000125 0.00333333 0.0011 0.000125 0.005 0.0011 0.000125 + 0 0.0012 0.000125 0.00166667 0.0012 0.000125 0.00333333 0.0012 0.000125 0.005 0.0012 0.000125 + 0 0.0013 0.000125 0.00166667 0.0013 0.000125 0.00333333 0.0013 0.000125 0.005 0.0013 0.000125 + 0 0.0014 0.000125 0.00166667 0.0014 0.000125 0.00333333 0.0014 0.000125 0.005 0.0014 0.000125 + 0 0.0015 0.000125 0.00166667 0.0015 0.000125 0.00333333 0.0015 0.000125 0.005 0.0015 0.000125 + 0 0.0016 0.000125 0.00166667 0.0016 0.000125 0.00333333 0.0016 0.000125 0.005 0.0016 0.000125 + 0 0.0017 0.000125 0.00166667 0.0017 0.000125 0.00333333 0.0017 0.000125 0.005 0.0017 0.000125 + 0 0.0018 0.000125 0.00166667 0.0018 0.000125 0.00333333 0.0018 0.000125 0.005 0.0018 0.000125 + 0 0.0019 0.000125 0.00166667 0.0019 0.000125 0.00333333 0.0019 0.000125 0.005 0.0019 0.000125 + 0 0.002 0.000125 0.00166667 0.002 0.000125 0.00333333 0.002 0.000125 0.005 0.002 0.000125 + 0 0 0.00015 0.00166667 0 0.00015 0 0.0001 0.00015 0.00166667 0.0001 0.00015 + 0.00333333 0 0.00015 0.00333333 0.0001 0.00015 0.005 0 0.00015 0.005 0.0001 0.00015 + 0 0.0002 0.00015 0.00166667 0.0002 0.00015 0.00333333 0.0002 0.00015 0.005 0.0002 0.00015 + 0 0.0003 0.00015 0.00166667 0.0003 0.00015 0.00333333 0.0003 0.00015 0.005 0.0003 0.00015 + 0 0.0004 0.00015 0.00166667 0.0004 0.00015 0.00333333 0.0004 0.00015 0.005 0.0004 0.00015 + 0 0.0005 0.00015 0.00166667 0.0005 0.00015 0.00333333 0.0005 0.00015 0.005 0.0005 0.00015 + 0 0.0006 0.00015 0.00166667 0.0006 0.00015 0.00333333 0.0006 0.00015 0.005 0.0006 0.00015 + 0 0.0007 0.00015 0.00166667 0.0007 0.00015 0.00333333 0.0007 0.00015 0.005 0.0007 0.00015 + 0 0.0008 0.00015 0.00166667 0.0008 0.00015 0.00333333 0.0008 0.00015 0.005 0.0008 0.00015 + 0 0.0009 0.00015 0.00166667 0.0009 0.00015 0.00333333 0.0009 0.00015 0.005 0.0009 0.00015 + 0 0.001 0.00015 0.00166667 0.001 0.00015 0.00333333 0.001 0.00015 0.005 0.001 0.00015 + 0 0.0011 0.00015 0.00166667 0.0011 0.00015 0.00333333 0.0011 0.00015 0.005 0.0011 0.00015 + 0 0.0012 0.00015 0.00166667 0.0012 0.00015 0.00333333 0.0012 0.00015 0.005 0.0012 0.00015 + 0 0.0013 0.00015 0.00166667 0.0013 0.00015 0.00333333 0.0013 0.00015 0.005 0.0013 0.00015 + 0 0.0014 0.00015 0.00166667 0.0014 0.00015 0.00333333 0.0014 0.00015 0.005 0.0014 0.00015 + 0 0.0015 0.00015 0.00166667 0.0015 0.00015 0.00333333 0.0015 0.00015 0.005 0.0015 0.00015 + 0 0.0016 0.00015 0.00166667 0.0016 0.00015 0.00333333 0.0016 0.00015 0.005 0.0016 0.00015 + 0 0.0017 0.00015 0.00166667 0.0017 0.00015 0.00333333 0.0017 0.00015 0.005 0.0017 0.00015 + 0 0.0018 0.00015 0.00166667 0.0018 0.00015 0.00333333 0.0018 0.00015 0.005 0.0018 0.00015 + 0 0.0019 0.00015 0.00166667 0.0019 0.00015 0.00333333 0.0019 0.00015 0.005 0.0019 0.00015 + 0 0.002 0.00015 0.00166667 0.002 0.00015 0.00333333 0.002 0.00015 0.005 0.002 0.00015 + 0 0 0.000175 0.00166667 0 0.000175 0 0.0001 0.000175 0.00166667 0.0001 0.000175 + 0.00333333 0 0.000175 0.00333333 0.0001 0.000175 0.005 0 0.000175 0.005 0.0001 0.000175 + 0 0.0002 0.000175 0.00166667 0.0002 0.000175 0.00333333 0.0002 0.000175 0.005 0.0002 0.000175 + 0 0.0003 0.000175 0.00166667 0.0003 0.000175 0.00333333 0.0003 0.000175 0.005 0.0003 0.000175 + 0 0.0004 0.000175 0.00166667 0.0004 0.000175 0.00333333 0.0004 0.000175 0.005 0.0004 0.000175 + 0 0.0005 0.000175 0.00166667 0.0005 0.000175 0.00333333 0.0005 0.000175 0.005 0.0005 0.000175 + 0 0.0006 0.000175 0.00166667 0.0006 0.000175 0.00333333 0.0006 0.000175 0.005 0.0006 0.000175 + 0 0.0007 0.000175 0.00166667 0.0007 0.000175 0.00333333 0.0007 0.000175 0.005 0.0007 0.000175 + 0 0.0008 0.000175 0.00166667 0.0008 0.000175 0.00333333 0.0008 0.000175 0.005 0.0008 0.000175 + 0 0.0009 0.000175 0.00166667 0.0009 0.000175 0.00333333 0.0009 0.000175 0.005 0.0009 0.000175 + 0 0.001 0.000175 0.00166667 0.001 0.000175 0.00333333 0.001 0.000175 0.005 0.001 0.000175 + 0 0.0011 0.000175 0.00166667 0.0011 0.000175 0.00333333 0.0011 0.000175 0.005 0.0011 0.000175 + 0 0.0012 0.000175 0.00166667 0.0012 0.000175 0.00333333 0.0012 0.000175 0.005 0.0012 0.000175 + 0 0.0013 0.000175 0.00166667 0.0013 0.000175 0.00333333 0.0013 0.000175 0.005 0.0013 0.000175 + 0 0.0014 0.000175 0.00166667 0.0014 0.000175 0.00333333 0.0014 0.000175 0.005 0.0014 0.000175 + 0 0.0015 0.000175 0.00166667 0.0015 0.000175 0.00333333 0.0015 0.000175 0.005 0.0015 0.000175 + 0 0.0016 0.000175 0.00166667 0.0016 0.000175 0.00333333 0.0016 0.000175 0.005 0.0016 0.000175 + 0 0.0017 0.000175 0.00166667 0.0017 0.000175 0.00333333 0.0017 0.000175 0.005 0.0017 0.000175 + 0 0.0018 0.000175 0.00166667 0.0018 0.000175 0.00333333 0.0018 0.000175 0.005 0.0018 0.000175 + 0 0.0019 0.000175 0.00166667 0.0019 0.000175 0.00333333 0.0019 0.000175 0.005 0.0019 0.000175 + 0 0.002 0.000175 0.00166667 0.002 0.000175 0.00333333 0.002 0.000175 0.005 0.002 0.000175 + 0 0 0.0002 0.00166667 0 0.0002 0 0.0001 0.0002 0.00166667 0.0001 0.0002 + 0.00333333 0 0.0002 0.00333333 0.0001 0.0002 0.005 0 0.0002 0.005 0.0001 0.0002 + 0 0.0002 0.0002 0.00166667 0.0002 0.0002 0.00333333 0.0002 0.0002 0.005 0.0002 0.0002 + 0 0.0003 0.0002 0.00166667 0.0003 0.0002 0.00333333 0.0003 0.0002 0.005 0.0003 0.0002 + 0 0.0004 0.0002 0.00166667 0.0004 0.0002 0.00333333 0.0004 0.0002 0.005 0.0004 0.0002 + 0 0.0005 0.0002 0.00166667 0.0005 0.0002 0.00333333 0.0005 0.0002 0.005 0.0005 0.0002 + 0 0.0006 0.0002 0.00166667 0.0006 0.0002 0.00333333 0.0006 0.0002 0.005 0.0006 0.0002 + 0 0.0007 0.0002 0.00166667 0.0007 0.0002 0.00333333 0.0007 0.0002 0.005 0.0007 0.0002 + 0 0.0008 0.0002 0.00166667 0.0008 0.0002 0.00333333 0.0008 0.0002 0.005 0.0008 0.0002 + 0 0.0009 0.0002 0.00166667 0.0009 0.0002 0.00333333 0.0009 0.0002 0.005 0.0009 0.0002 + 0 0.001 0.0002 0.00166667 0.001 0.0002 0.00333333 0.001 0.0002 0.005 0.001 0.0002 + 0 0.0011 0.0002 0.00166667 0.0011 0.0002 0.00333333 0.0011 0.0002 0.005 0.0011 0.0002 + 0 0.0012 0.0002 0.00166667 0.0012 0.0002 0.00333333 0.0012 0.0002 0.005 0.0012 0.0002 + 0 0.0013 0.0002 0.00166667 0.0013 0.0002 0.00333333 0.0013 0.0002 0.005 0.0013 0.0002 + 0 0.0014 0.0002 0.00166667 0.0014 0.0002 0.00333333 0.0014 0.0002 0.005 0.0014 0.0002 + 0 0.0015 0.0002 0.00166667 0.0015 0.0002 0.00333333 0.0015 0.0002 0.005 0.0015 0.0002 + 0 0.0016 0.0002 0.00166667 0.0016 0.0002 0.00333333 0.0016 0.0002 0.005 0.0016 0.0002 + 0 0.0017 0.0002 0.00166667 0.0017 0.0002 0.00333333 0.0017 0.0002 0.005 0.0017 0.0002 + 0 0.0018 0.0002 0.00166667 0.0018 0.0002 0.00333333 0.0018 0.0002 0.005 0.0018 0.0002 + 0 0.0019 0.0002 0.00166667 0.0019 0.0002 0.00333333 0.0019 0.0002 0.005 0.0019 0.0002 + 0 0.002 0.0002 0.00166667 0.002 0.0002 0.00333333 0.002 0.0002 0.005 0.002 0.0002 + 0 0 0.000225 0.00166667 0 0.000225 0 0.0001 0.000225 0.00166667 0.0001 0.000225 + 0.00333333 0 0.000225 0.00333333 0.0001 0.000225 0.005 0 0.000225 0.005 0.0001 0.000225 + 0 0.0002 0.000225 0.00166667 0.0002 0.000225 0.00333333 0.0002 0.000225 0.005 0.0002 0.000225 + 0 0.0003 0.000225 0.00166667 0.0003 0.000225 0.00333333 0.0003 0.000225 0.005 0.0003 0.000225 + 0 0.0004 0.000225 0.00166667 0.0004 0.000225 0.00333333 0.0004 0.000225 0.005 0.0004 0.000225 + 0 0.0005 0.000225 0.00166667 0.0005 0.000225 0.00333333 0.0005 0.000225 0.005 0.0005 0.000225 + 0 0.0006 0.000225 0.00166667 0.0006 0.000225 0.00333333 0.0006 0.000225 0.005 0.0006 0.000225 + 0 0.0007 0.000225 0.00166667 0.0007 0.000225 0.00333333 0.0007 0.000225 0.005 0.0007 0.000225 + 0 0.0008 0.000225 0.00166667 0.0008 0.000225 0.00333333 0.0008 0.000225 0.005 0.0008 0.000225 + 0 0.0009 0.000225 0.00166667 0.0009 0.000225 0.00333333 0.0009 0.000225 0.005 0.0009 0.000225 + 0 0.001 0.000225 0.00166667 0.001 0.000225 0.00333333 0.001 0.000225 0.005 0.001 0.000225 + 0 0.0011 0.000225 0.00166667 0.0011 0.000225 0.00333333 0.0011 0.000225 0.005 0.0011 0.000225 + 0 0.0012 0.000225 0.00166667 0.0012 0.000225 0.00333333 0.0012 0.000225 0.005 0.0012 0.000225 + 0 0.0013 0.000225 0.00166667 0.0013 0.000225 0.00333333 0.0013 0.000225 0.005 0.0013 0.000225 + 0 0.0014 0.000225 0.00166667 0.0014 0.000225 0.00333333 0.0014 0.000225 0.005 0.0014 0.000225 + 0 0.0015 0.000225 0.00166667 0.0015 0.000225 0.00333333 0.0015 0.000225 0.005 0.0015 0.000225 + 0 0.0016 0.000225 0.00166667 0.0016 0.000225 0.00333333 0.0016 0.000225 0.005 0.0016 0.000225 + 0 0.0017 0.000225 0.00166667 0.0017 0.000225 0.00333333 0.0017 0.000225 0.005 0.0017 0.000225 + 0 0.0018 0.000225 0.00166667 0.0018 0.000225 0.00333333 0.0018 0.000225 0.005 0.0018 0.000225 + 0 0.0019 0.000225 0.00166667 0.0019 0.000225 0.00333333 0.0019 0.000225 0.005 0.0019 0.000225 + 0 0.002 0.000225 0.00166667 0.002 0.000225 0.00333333 0.002 0.000225 0.005 0.002 0.000225 + 0 0 0.00025 0.00166667 0 0.00025 0 0.0001 0.00025 0.00166667 0.0001 0.00025 + 0.00333333 0 0.00025 0.00333333 0.0001 0.00025 0.005 0 0.00025 0.005 0.0001 0.00025 + 0 0.0002 0.00025 0.00166667 0.0002 0.00025 0.00333333 0.0002 0.00025 0.005 0.0002 0.00025 + 0 0.0003 0.00025 0.00166667 0.0003 0.00025 0.00333333 0.0003 0.00025 0.005 0.0003 0.00025 + 0 0.0004 0.00025 0.00166667 0.0004 0.00025 0.00333333 0.0004 0.00025 0.005 0.0004 0.00025 + 0 0.0005 0.00025 0.00166667 0.0005 0.00025 0.00333333 0.0005 0.00025 0.005 0.0005 0.00025 + 0 0.0006 0.00025 0.00166667 0.0006 0.00025 0.00333333 0.0006 0.00025 0.005 0.0006 0.00025 + 0 0.0007 0.00025 0.00166667 0.0007 0.00025 0.00333333 0.0007 0.00025 0.005 0.0007 0.00025 + 0 0.0008 0.00025 0.00166667 0.0008 0.00025 0.00333333 0.0008 0.00025 0.005 0.0008 0.00025 + 0 0.0009 0.00025 0.00166667 0.0009 0.00025 0.00333333 0.0009 0.00025 0.005 0.0009 0.00025 + 0 0.001 0.00025 0.00166667 0.001 0.00025 0.00333333 0.001 0.00025 0.005 0.001 0.00025 + 0 0.0011 0.00025 0.00166667 0.0011 0.00025 0.00333333 0.0011 0.00025 0.005 0.0011 0.00025 + 0 0.0012 0.00025 0.00166667 0.0012 0.00025 0.00333333 0.0012 0.00025 0.005 0.0012 0.00025 + 0 0.0013 0.00025 0.00166667 0.0013 0.00025 0.00333333 0.0013 0.00025 0.005 0.0013 0.00025 + 0 0.0014 0.00025 0.00166667 0.0014 0.00025 0.00333333 0.0014 0.00025 0.005 0.0014 0.00025 + 0 0.0015 0.00025 0.00166667 0.0015 0.00025 0.00333333 0.0015 0.00025 0.005 0.0015 0.00025 + 0 0.0016 0.00025 0.00166667 0.0016 0.00025 0.00333333 0.0016 0.00025 0.005 0.0016 0.00025 + 0 0.0017 0.00025 0.00166667 0.0017 0.00025 0.00333333 0.0017 0.00025 0.005 0.0017 0.00025 + 0 0.0018 0.00025 0.00166667 0.0018 0.00025 0.00333333 0.0018 0.00025 0.005 0.0018 0.00025 + 0 0.0019 0.00025 0.00166667 0.0019 0.00025 0.00333333 0.0019 0.00025 0.005 0.0019 0.00025 + 0 0.002 0.00025 0.00166667 0.002 0.00025 0.00333333 0.002 0.00025 0.005 0.002 0.00025 + </DataArray> + </Points> + <Cells> + <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii"> + 0 1 3 2 4 5 7 6 1 8 9 3 + 5 10 11 7 8 12 13 9 10 14 15 11 + 2 3 17 16 6 7 19 18 3 9 20 17 + 7 11 21 19 9 13 22 20 11 15 23 21 + 16 17 25 24 18 19 27 26 17 20 28 25 + 19 21 29 27 20 22 30 28 21 23 31 29 + 24 25 33 32 26 27 35 34 25 28 36 33 + 27 29 37 35 28 30 38 36 29 31 39 37 + 32 33 41 40 34 35 43 42 33 36 44 41 + 35 37 45 43 36 38 46 44 37 39 47 45 + 40 41 49 48 42 43 51 50 41 44 52 49 + 43 45 53 51 44 46 54 52 45 47 55 53 + 48 49 57 56 50 51 59 58 49 52 60 57 + 51 53 61 59 52 54 62 60 53 55 63 61 + 56 57 65 64 58 59 67 66 57 60 68 65 + 59 61 69 67 60 62 70 68 61 63 71 69 + 64 65 73 72 66 67 75 74 65 68 76 73 + 67 69 77 75 68 70 78 76 69 71 79 77 + 72 73 81 80 74 75 83 82 73 76 84 81 + 75 77 85 83 76 78 86 84 77 79 87 85 + 80 81 89 88 82 83 91 90 81 84 92 89 + 83 85 93 91 84 86 94 92 85 87 95 93 + 88 89 97 96 90 91 99 98 89 92 100 97 + 91 93 101 99 92 94 102 100 93 95 103 101 + 96 97 105 104 98 99 107 106 97 100 108 105 + 99 101 109 107 100 102 110 108 101 103 111 109 + 104 105 113 112 106 107 115 114 105 108 116 113 + 107 109 117 115 108 110 118 116 109 111 119 117 + 112 113 121 120 114 115 123 122 113 116 124 121 + 115 117 125 123 116 118 126 124 117 119 127 125 + 120 121 129 128 122 123 131 130 121 124 132 129 + 123 125 133 131 124 126 134 132 125 127 135 133 + 128 129 137 136 130 131 139 138 129 132 140 137 + 131 133 141 139 132 134 142 140 133 135 143 141 + 136 137 145 144 138 139 147 146 137 140 148 145 + 139 141 149 147 140 142 150 148 141 143 151 149 + 144 145 153 152 146 147 155 154 145 148 156 153 + 147 149 157 155 148 150 158 156 149 151 159 157 + 152 153 161 160 154 155 163 162 153 156 164 161 + 155 157 165 163 156 158 166 164 157 159 167 165 + 4 5 7 6 168 169 171 170 5 10 11 7 + 169 172 173 171 10 14 15 11 172 174 175 173 + 6 7 19 18 170 171 177 176 7 11 21 19 + 171 173 178 177 11 15 23 21 173 175 179 178 + 18 19 27 26 176 177 181 180 19 21 29 27 + 177 178 182 181 21 23 31 29 178 179 183 182 + 26 27 35 34 180 181 185 184 27 29 37 35 + 181 182 186 185 29 31 39 37 182 183 187 186 + 34 35 43 42 184 185 189 188 35 37 45 43 + 185 186 190 189 37 39 47 45 186 187 191 190 + 42 43 51 50 188 189 193 192 43 45 53 51 + 189 190 194 193 45 47 55 53 190 191 195 194 + 50 51 59 58 192 193 197 196 51 53 61 59 + 193 194 198 197 53 55 63 61 194 195 199 198 + 58 59 67 66 196 197 201 200 59 61 69 67 + 197 198 202 201 61 63 71 69 198 199 203 202 + 66 67 75 74 200 201 205 204 67 69 77 75 + 201 202 206 205 69 71 79 77 202 203 207 206 + 74 75 83 82 204 205 209 208 75 77 85 83 + 205 206 210 209 77 79 87 85 206 207 211 210 + 82 83 91 90 208 209 213 212 83 85 93 91 + 209 210 214 213 85 87 95 93 210 211 215 214 + 90 91 99 98 212 213 217 216 91 93 101 99 + 213 214 218 217 93 95 103 101 214 215 219 218 + 98 99 107 106 216 217 221 220 99 101 109 107 + 217 218 222 221 101 103 111 109 218 219 223 222 + 106 107 115 114 220 221 225 224 107 109 117 115 + 221 222 226 225 109 111 119 117 222 223 227 226 + 114 115 123 122 224 225 229 228 115 117 125 123 + 225 226 230 229 117 119 127 125 226 227 231 230 + 122 123 131 130 228 229 233 232 123 125 133 131 + 229 230 234 233 125 127 135 133 230 231 235 234 + 130 131 139 138 232 233 237 236 131 133 141 139 + 233 234 238 237 133 135 143 141 234 235 239 238 + 138 139 147 146 236 237 241 240 139 141 149 147 + 237 238 242 241 141 143 151 149 238 239 243 242 + 146 147 155 154 240 241 245 244 147 149 157 155 + 241 242 246 245 149 151 159 157 242 243 247 246 + 154 155 163 162 244 245 249 248 155 157 165 163 + 245 246 250 249 157 159 167 165 246 247 251 250 + 168 169 171 170 252 253 255 254 169 172 173 171 + 253 256 257 255 172 174 175 173 256 258 259 257 + 170 171 177 176 254 255 261 260 171 173 178 177 + 255 257 262 261 173 175 179 178 257 259 263 262 + 176 177 181 180 260 261 265 264 177 178 182 181 + 261 262 266 265 178 179 183 182 262 263 267 266 + 180 181 185 184 264 265 269 268 181 182 186 185 + 265 266 270 269 182 183 187 186 266 267 271 270 + 184 185 189 188 268 269 273 272 185 186 190 189 + 269 270 274 273 186 187 191 190 270 271 275 274 + 188 189 193 192 272 273 277 276 189 190 194 193 + 273 274 278 277 190 191 195 194 274 275 279 278 + 192 193 197 196 276 277 281 280 193 194 198 197 + 277 278 282 281 194 195 199 198 278 279 283 282 + 196 197 201 200 280 281 285 284 197 198 202 201 + 281 282 286 285 198 199 203 202 282 283 287 286 + 200 201 205 204 284 285 289 288 201 202 206 205 + 285 286 290 289 202 203 207 206 286 287 291 290 + 204 205 209 208 288 289 293 292 205 206 210 209 + 289 290 294 293 206 207 211 210 290 291 295 294 + 208 209 213 212 292 293 297 296 209 210 214 213 + 293 294 298 297 210 211 215 214 294 295 299 298 + 212 213 217 216 296 297 301 300 213 214 218 217 + 297 298 302 301 214 215 219 218 298 299 303 302 + 216 217 221 220 300 301 305 304 217 218 222 221 + 301 302 306 305 218 219 223 222 302 303 307 306 + 220 221 225 224 304 305 309 308 221 222 226 225 + 305 306 310 309 222 223 227 226 306 307 311 310 + 224 225 229 228 308 309 313 312 225 226 230 229 + 309 310 314 313 226 227 231 230 310 311 315 314 + 228 229 233 232 312 313 317 316 229 230 234 233 + 313 314 318 317 230 231 235 234 314 315 319 318 + 232 233 237 236 316 317 321 320 233 234 238 237 + 317 318 322 321 234 235 239 238 318 319 323 322 + 236 237 241 240 320 321 325 324 237 238 242 241 + 321 322 326 325 238 239 243 242 322 323 327 326 + 240 241 245 244 324 325 329 328 241 242 246 245 + 325 326 330 329 242 243 247 246 326 327 331 330 + 244 245 249 248 328 329 333 332 245 246 250 249 + 329 330 334 333 246 247 251 250 330 331 335 334 + 252 253 255 254 336 337 339 338 253 256 257 255 + 337 340 341 339 256 258 259 257 340 342 343 341 + 254 255 261 260 338 339 345 344 255 257 262 261 + 339 341 346 345 257 259 263 262 341 343 347 346 + 260 261 265 264 344 345 349 348 261 262 266 265 + 345 346 350 349 262 263 267 266 346 347 351 350 + 264 265 269 268 348 349 353 352 265 266 270 269 + 349 350 354 353 266 267 271 270 350 351 355 354 + 268 269 273 272 352 353 357 356 269 270 274 273 + 353 354 358 357 270 271 275 274 354 355 359 358 + 272 273 277 276 356 357 361 360 273 274 278 277 + 357 358 362 361 274 275 279 278 358 359 363 362 + 276 277 281 280 360 361 365 364 277 278 282 281 + 361 362 366 365 278 279 283 282 362 363 367 366 + 280 281 285 284 364 365 369 368 281 282 286 285 + 365 366 370 369 282 283 287 286 366 367 371 370 + 284 285 289 288 368 369 373 372 285 286 290 289 + 369 370 374 373 286 287 291 290 370 371 375 374 + 288 289 293 292 372 373 377 376 289 290 294 293 + 373 374 378 377 290 291 295 294 374 375 379 378 + 292 293 297 296 376 377 381 380 293 294 298 297 + 377 378 382 381 294 295 299 298 378 379 383 382 + 296 297 301 300 380 381 385 384 297 298 302 301 + 381 382 386 385 298 299 303 302 382 383 387 386 + 300 301 305 304 384 385 389 388 301 302 306 305 + 385 386 390 389 302 303 307 306 386 387 391 390 + 304 305 309 308 388 389 393 392 305 306 310 309 + 389 390 394 393 306 307 311 310 390 391 395 394 + 308 309 313 312 392 393 397 396 309 310 314 313 + 393 394 398 397 310 311 315 314 394 395 399 398 + 312 313 317 316 396 397 401 400 313 314 318 317 + 397 398 402 401 314 315 319 318 398 399 403 402 + 316 317 321 320 400 401 405 404 317 318 322 321 + 401 402 406 405 318 319 323 322 402 403 407 406 + 320 321 325 324 404 405 409 408 321 322 326 325 + 405 406 410 409 322 323 327 326 406 407 411 410 + 324 325 329 328 408 409 413 412 325 326 330 329 + 409 410 414 413 326 327 331 330 410 411 415 414 + 328 329 333 332 412 413 417 416 329 330 334 333 + 413 414 418 417 330 331 335 334 414 415 419 418 + 336 337 339 338 420 421 423 422 337 340 341 339 + 421 424 425 423 340 342 343 341 424 426 427 425 + 338 339 345 344 422 423 429 428 339 341 346 345 + 423 425 430 429 341 343 347 346 425 427 431 430 + 344 345 349 348 428 429 433 432 345 346 350 349 + 429 430 434 433 346 347 351 350 430 431 435 434 + 348 349 353 352 432 433 437 436 349 350 354 353 + 433 434 438 437 350 351 355 354 434 435 439 438 + 352 353 357 356 436 437 441 440 353 354 358 357 + 437 438 442 441 354 355 359 358 438 439 443 442 + 356 357 361 360 440 441 445 444 357 358 362 361 + 441 442 446 445 358 359 363 362 442 443 447 446 + 360 361 365 364 444 445 449 448 361 362 366 365 + 445 446 450 449 362 363 367 366 446 447 451 450 + 364 365 369 368 448 449 453 452 365 366 370 369 + 449 450 454 453 366 367 371 370 450 451 455 454 + 368 369 373 372 452 453 457 456 369 370 374 373 + 453 454 458 457 370 371 375 374 454 455 459 458 + 372 373 377 376 456 457 461 460 373 374 378 377 + 457 458 462 461 374 375 379 378 458 459 463 462 + 376 377 381 380 460 461 465 464 377 378 382 381 + 461 462 466 465 378 379 383 382 462 463 467 466 + 380 381 385 384 464 465 469 468 381 382 386 385 + 465 466 470 469 382 383 387 386 466 467 471 470 + 384 385 389 388 468 469 473 472 385 386 390 389 + 469 470 474 473 386 387 391 390 470 471 475 474 + 388 389 393 392 472 473 477 476 389 390 394 393 + 473 474 478 477 390 391 395 394 474 475 479 478 + 392 393 397 396 476 477 481 480 393 394 398 397 + 477 478 482 481 394 395 399 398 478 479 483 482 + 396 397 401 400 480 481 485 484 397 398 402 401 + 481 482 486 485 398 399 403 402 482 483 487 486 + 400 401 405 404 484 485 489 488 401 402 406 405 + 485 486 490 489 402 403 407 406 486 487 491 490 + 404 405 409 408 488 489 493 492 405 406 410 409 + 489 490 494 493 406 407 411 410 490 491 495 494 + 408 409 413 412 492 493 497 496 409 410 414 413 + 493 494 498 497 410 411 415 414 494 495 499 498 + 412 413 417 416 496 497 501 500 413 414 418 417 + 497 498 502 501 414 415 419 418 498 499 503 502 + 420 421 423 422 504 505 507 506 421 424 425 423 + 505 508 509 507 424 426 427 425 508 510 511 509 + 422 423 429 428 506 507 513 512 423 425 430 429 + 507 509 514 513 425 427 431 430 509 511 515 514 + 428 429 433 432 512 513 517 516 429 430 434 433 + 513 514 518 517 430 431 435 434 514 515 519 518 + 432 433 437 436 516 517 521 520 433 434 438 437 + 517 518 522 521 434 435 439 438 518 519 523 522 + 436 437 441 440 520 521 525 524 437 438 442 441 + 521 522 526 525 438 439 443 442 522 523 527 526 + 440 441 445 444 524 525 529 528 441 442 446 445 + 525 526 530 529 442 443 447 446 526 527 531 530 + 444 445 449 448 528 529 533 532 445 446 450 449 + 529 530 534 533 446 447 451 450 530 531 535 534 + 448 449 453 452 532 533 537 536 449 450 454 453 + 533 534 538 537 450 451 455 454 534 535 539 538 + 452 453 457 456 536 537 541 540 453 454 458 457 + 537 538 542 541 454 455 459 458 538 539 543 542 + 456 457 461 460 540 541 545 544 457 458 462 461 + 541 542 546 545 458 459 463 462 542 543 547 546 + 460 461 465 464 544 545 549 548 461 462 466 465 + 545 546 550 549 462 463 467 466 546 547 551 550 + 464 465 469 468 548 549 553 552 465 466 470 469 + 549 550 554 553 466 467 471 470 550 551 555 554 + 468 469 473 472 552 553 557 556 469 470 474 473 + 553 554 558 557 470 471 475 474 554 555 559 558 + 472 473 477 476 556 557 561 560 473 474 478 477 + 557 558 562 561 474 475 479 478 558 559 563 562 + 476 477 481 480 560 561 565 564 477 478 482 481 + 561 562 566 565 478 479 483 482 562 563 567 566 + 480 481 485 484 564 565 569 568 481 482 486 485 + 565 566 570 569 482 483 487 486 566 567 571 570 + 484 485 489 488 568 569 573 572 485 486 490 489 + 569 570 574 573 486 487 491 490 570 571 575 574 + 488 489 493 492 572 573 577 576 489 490 494 493 + 573 574 578 577 490 491 495 494 574 575 579 578 + 492 493 497 496 576 577 581 580 493 494 498 497 + 577 578 582 581 494 495 499 498 578 579 583 582 + 496 497 501 500 580 581 585 584 497 498 502 501 + 581 582 586 585 498 499 503 502 582 583 587 586 + 504 505 507 506 588 589 591 590 505 508 509 507 + 589 592 593 591 508 510 511 509 592 594 595 593 + 506 507 513 512 590 591 597 596 507 509 514 513 + 591 593 598 597 509 511 515 514 593 595 599 598 + 512 513 517 516 596 597 601 600 513 514 518 517 + 597 598 602 601 514 515 519 518 598 599 603 602 + 516 517 521 520 600 601 605 604 517 518 522 521 + 601 602 606 605 518 519 523 522 602 603 607 606 + 520 521 525 524 604 605 609 608 521 522 526 525 + 605 606 610 609 522 523 527 526 606 607 611 610 + 524 525 529 528 608 609 613 612 525 526 530 529 + 609 610 614 613 526 527 531 530 610 611 615 614 + 528 529 533 532 612 613 617 616 529 530 534 533 + 613 614 618 617 530 531 535 534 614 615 619 618 + 532 533 537 536 616 617 621 620 533 534 538 537 + 617 618 622 621 534 535 539 538 618 619 623 622 + 536 537 541 540 620 621 625 624 537 538 542 541 + 621 622 626 625 538 539 543 542 622 623 627 626 + 540 541 545 544 624 625 629 628 541 542 546 545 + 625 626 630 629 542 543 547 546 626 627 631 630 + 544 545 549 548 628 629 633 632 545 546 550 549 + 629 630 634 633 546 547 551 550 630 631 635 634 + 548 549 553 552 632 633 637 636 549 550 554 553 + 633 634 638 637 550 551 555 554 634 635 639 638 + 552 553 557 556 636 637 641 640 553 554 558 557 + 637 638 642 641 554 555 559 558 638 639 643 642 + 556 557 561 560 640 641 645 644 557 558 562 561 + 641 642 646 645 558 559 563 562 642 643 647 646 + 560 561 565 564 644 645 649 648 561 562 566 565 + 645 646 650 649 562 563 567 566 646 647 651 650 + 564 565 569 568 648 649 653 652 565 566 570 569 + 649 650 654 653 566 567 571 570 650 651 655 654 + 568 569 573 572 652 653 657 656 569 570 574 573 + 653 654 658 657 570 571 575 574 654 655 659 658 + 572 573 577 576 656 657 661 660 573 574 578 577 + 657 658 662 661 574 575 579 578 658 659 663 662 + 576 577 581 580 660 661 665 664 577 578 582 581 + 661 662 666 665 578 579 583 582 662 663 667 666 + 580 581 585 584 664 665 669 668 581 582 586 585 + 665 666 670 669 582 583 587 586 666 667 671 670 + 588 589 591 590 672 673 675 674 589 592 593 591 + 673 676 677 675 592 594 595 593 676 678 679 677 + 590 591 597 596 674 675 681 680 591 593 598 597 + 675 677 682 681 593 595 599 598 677 679 683 682 + 596 597 601 600 680 681 685 684 597 598 602 601 + 681 682 686 685 598 599 603 602 682 683 687 686 + 600 601 605 604 684 685 689 688 601 602 606 605 + 685 686 690 689 602 603 607 606 686 687 691 690 + 604 605 609 608 688 689 693 692 605 606 610 609 + 689 690 694 693 606 607 611 610 690 691 695 694 + 608 609 613 612 692 693 697 696 609 610 614 613 + 693 694 698 697 610 611 615 614 694 695 699 698 + 612 613 617 616 696 697 701 700 613 614 618 617 + 697 698 702 701 614 615 619 618 698 699 703 702 + 616 617 621 620 700 701 705 704 617 618 622 621 + 701 702 706 705 618 619 623 622 702 703 707 706 + 620 621 625 624 704 705 709 708 621 622 626 625 + 705 706 710 709 622 623 627 626 706 707 711 710 + 624 625 629 628 708 709 713 712 625 626 630 629 + 709 710 714 713 626 627 631 630 710 711 715 714 + 628 629 633 632 712 713 717 716 629 630 634 633 + 713 714 718 717 630 631 635 634 714 715 719 718 + 632 633 637 636 716 717 721 720 633 634 638 637 + 717 718 722 721 634 635 639 638 718 719 723 722 + 636 637 641 640 720 721 725 724 637 638 642 641 + 721 722 726 725 638 639 643 642 722 723 727 726 + 640 641 645 644 724 725 729 728 641 642 646 645 + 725 726 730 729 642 643 647 646 726 727 731 730 + 644 645 649 648 728 729 733 732 645 646 650 649 + 729 730 734 733 646 647 651 650 730 731 735 734 + 648 649 653 652 732 733 737 736 649 650 654 653 + 733 734 738 737 650 651 655 654 734 735 739 738 + 652 653 657 656 736 737 741 740 653 654 658 657 + 737 738 742 741 654 655 659 658 738 739 743 742 + 656 657 661 660 740 741 745 744 657 658 662 661 + 741 742 746 745 658 659 663 662 742 743 747 746 + 660 661 665 664 744 745 749 748 661 662 666 665 + 745 746 750 749 662 663 667 666 746 747 751 750 + 664 665 669 668 748 749 753 752 665 666 670 669 + 749 750 754 753 666 667 671 670 750 751 755 754 + 672 673 675 674 756 757 759 758 673 676 677 675 + 757 760 761 759 676 678 679 677 760 762 763 761 + 674 675 681 680 758 759 765 764 675 677 682 681 + 759 761 766 765 677 679 683 682 761 763 767 766 + 680 681 685 684 764 765 769 768 681 682 686 685 + 765 766 770 769 682 683 687 686 766 767 771 770 + 684 685 689 688 768 769 773 772 685 686 690 689 + 769 770 774 773 686 687 691 690 770 771 775 774 + 688 689 693 692 772 773 777 776 689 690 694 693 + 773 774 778 777 690 691 695 694 774 775 779 778 + 692 693 697 696 776 777 781 780 693 694 698 697 + 777 778 782 781 694 695 699 698 778 779 783 782 + 696 697 701 700 780 781 785 784 697 698 702 701 + 781 782 786 785 698 699 703 702 782 783 787 786 + 700 701 705 704 784 785 789 788 701 702 706 705 + 785 786 790 789 702 703 707 706 786 787 791 790 + 704 705 709 708 788 789 793 792 705 706 710 709 + 789 790 794 793 706 707 711 710 790 791 795 794 + 708 709 713 712 792 793 797 796 709 710 714 713 + 793 794 798 797 710 711 715 714 794 795 799 798 + 712 713 717 716 796 797 801 800 713 714 718 717 + 797 798 802 801 714 715 719 718 798 799 803 802 + 716 717 721 720 800 801 805 804 717 718 722 721 + 801 802 806 805 718 719 723 722 802 803 807 806 + 720 721 725 724 804 805 809 808 721 722 726 725 + 805 806 810 809 722 723 727 726 806 807 811 810 + 724 725 729 728 808 809 813 812 725 726 730 729 + 809 810 814 813 726 727 731 730 810 811 815 814 + 728 729 733 732 812 813 817 816 729 730 734 733 + 813 814 818 817 730 731 735 734 814 815 819 818 + 732 733 737 736 816 817 821 820 733 734 738 737 + 817 818 822 821 734 735 739 738 818 819 823 822 + 736 737 741 740 820 821 825 824 737 738 742 741 + 821 822 826 825 738 739 743 742 822 823 827 826 + 740 741 745 744 824 825 829 828 741 742 746 745 + 825 826 830 829 742 743 747 746 826 827 831 830 + 744 745 749 748 828 829 833 832 745 746 750 749 + 829 830 834 833 746 747 751 750 830 831 835 834 + 748 749 753 752 832 833 837 836 749 750 754 753 + 833 834 838 837 750 751 755 754 834 835 839 838 + 756 757 759 758 840 841 843 842 757 760 761 759 + 841 844 845 843 760 762 763 761 844 846 847 845 + 758 759 765 764 842 843 849 848 759 761 766 765 + 843 845 850 849 761 763 767 766 845 847 851 850 + 764 765 769 768 848 849 853 852 765 766 770 769 + 849 850 854 853 766 767 771 770 850 851 855 854 + 768 769 773 772 852 853 857 856 769 770 774 773 + 853 854 858 857 770 771 775 774 854 855 859 858 + 772 773 777 776 856 857 861 860 773 774 778 777 + 857 858 862 861 774 775 779 778 858 859 863 862 + 776 777 781 780 860 861 865 864 777 778 782 781 + 861 862 866 865 778 779 783 782 862 863 867 866 + 780 781 785 784 864 865 869 868 781 782 786 785 + 865 866 870 869 782 783 787 786 866 867 871 870 + 784 785 789 788 868 869 873 872 785 786 790 789 + 869 870 874 873 786 787 791 790 870 871 875 874 + 788 789 793 792 872 873 877 876 789 790 794 793 + 873 874 878 877 790 791 795 794 874 875 879 878 + 792 793 797 796 876 877 881 880 793 794 798 797 + 877 878 882 881 794 795 799 798 878 879 883 882 + 796 797 801 800 880 881 885 884 797 798 802 801 + 881 882 886 885 798 799 803 802 882 883 887 886 + 800 801 805 804 884 885 889 888 801 802 806 805 + 885 886 890 889 802 803 807 806 886 887 891 890 + 804 805 809 808 888 889 893 892 805 806 810 809 + 889 890 894 893 806 807 811 810 890 891 895 894 + 808 809 813 812 892 893 897 896 809 810 814 813 + 893 894 898 897 810 811 815 814 894 895 899 898 + 812 813 817 816 896 897 901 900 813 814 818 817 + 897 898 902 901 814 815 819 818 898 899 903 902 + 816 817 821 820 900 901 905 904 817 818 822 821 + 901 902 906 905 818 819 823 822 902 903 907 906 + 820 821 825 824 904 905 909 908 821 822 826 825 + 905 906 910 909 822 823 827 826 906 907 911 910 + 824 825 829 828 908 909 913 912 825 826 830 829 + 909 910 914 913 826 827 831 830 910 911 915 914 + 828 829 833 832 912 913 917 916 829 830 834 833 + 913 914 918 917 830 831 835 834 914 915 919 918 + 832 833 837 836 916 917 921 920 833 834 838 837 + 917 918 922 921 834 835 839 838 918 919 923 922 + </DataArray> + <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii"> + 8 16 24 32 40 48 56 64 72 80 88 96 + 104 112 120 128 136 144 152 160 168 176 184 192 + 200 208 216 224 232 240 248 256 264 272 280 288 + 296 304 312 320 328 336 344 352 360 368 376 384 + 392 400 408 416 424 432 440 448 456 464 472 480 + 488 496 504 512 520 528 536 544 552 560 568 576 + 584 592 600 608 616 624 632 640 648 656 664 672 + 680 688 696 704 712 720 728 736 744 752 760 768 + 776 784 792 800 808 816 824 832 840 848 856 864 + 872 880 888 896 904 912 920 928 936 944 952 960 + 968 976 984 992 1000 1008 1016 1024 1032 1040 1048 1056 + 1064 1072 1080 1088 1096 1104 1112 1120 1128 1136 1144 1152 + 1160 1168 1176 1184 1192 1200 1208 1216 1224 1232 1240 1248 + 1256 1264 1272 1280 1288 1296 1304 1312 1320 1328 1336 1344 + 1352 1360 1368 1376 1384 1392 1400 1408 1416 1424 1432 1440 + 1448 1456 1464 1472 1480 1488 1496 1504 1512 1520 1528 1536 + 1544 1552 1560 1568 1576 1584 1592 1600 1608 1616 1624 1632 + 1640 1648 1656 1664 1672 1680 1688 1696 1704 1712 1720 1728 + 1736 1744 1752 1760 1768 1776 1784 1792 1800 1808 1816 1824 + 1832 1840 1848 1856 1864 1872 1880 1888 1896 1904 1912 1920 + 1928 1936 1944 1952 1960 1968 1976 1984 1992 2000 2008 2016 + 2024 2032 2040 2048 2056 2064 2072 2080 2088 2096 2104 2112 + 2120 2128 2136 2144 2152 2160 2168 2176 2184 2192 2200 2208 + 2216 2224 2232 2240 2248 2256 2264 2272 2280 2288 2296 2304 + 2312 2320 2328 2336 2344 2352 2360 2368 2376 2384 2392 2400 + 2408 2416 2424 2432 2440 2448 2456 2464 2472 2480 2488 2496 + 2504 2512 2520 2528 2536 2544 2552 2560 2568 2576 2584 2592 + 2600 2608 2616 2624 2632 2640 2648 2656 2664 2672 2680 2688 + 2696 2704 2712 2720 2728 2736 2744 2752 2760 2768 2776 2784 + 2792 2800 2808 2816 2824 2832 2840 2848 2856 2864 2872 2880 + 2888 2896 2904 2912 2920 2928 2936 2944 2952 2960 2968 2976 + 2984 2992 3000 3008 3016 3024 3032 3040 3048 3056 3064 3072 + 3080 3088 3096 3104 3112 3120 3128 3136 3144 3152 3160 3168 + 3176 3184 3192 3200 3208 3216 3224 3232 3240 3248 3256 3264 + 3272 3280 3288 3296 3304 3312 3320 3328 3336 3344 3352 3360 + 3368 3376 3384 3392 3400 3408 3416 3424 3432 3440 3448 3456 + 3464 3472 3480 3488 3496 3504 3512 3520 3528 3536 3544 3552 + 3560 3568 3576 3584 3592 3600 3608 3616 3624 3632 3640 3648 + 3656 3664 3672 3680 3688 3696 3704 3712 3720 3728 3736 3744 + 3752 3760 3768 3776 3784 3792 3800 3808 3816 3824 3832 3840 + 3848 3856 3864 3872 3880 3888 3896 3904 3912 3920 3928 3936 + 3944 3952 3960 3968 3976 3984 3992 4000 4008 4016 4024 4032 + 4040 4048 4056 4064 4072 4080 4088 4096 4104 4112 4120 4128 + 4136 4144 4152 4160 4168 4176 4184 4192 4200 4208 4216 4224 + 4232 4240 4248 4256 4264 4272 4280 4288 4296 4304 4312 4320 + 4328 4336 4344 4352 4360 4368 4376 4384 4392 4400 4408 4416 + 4424 4432 4440 4448 4456 4464 4472 4480 4488 4496 4504 4512 + 4520 4528 4536 4544 4552 4560 4568 4576 4584 4592 4600 4608 + 4616 4624 4632 4640 4648 4656 4664 4672 4680 4688 4696 4704 + 4712 4720 4728 4736 4744 4752 4760 4768 4776 4784 4792 4800 + </DataArray> + <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii"> + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + 12 12 12 12 12 12 12 12 12 12 12 12 + </DataArray> + </Cells> + </Piece> + </UnstructuredGrid> +</VTKFile> diff --git a/test/references/stokes_channel_pseudo3d-reference.vtu b/test/references/stokes_channel_pseudo3d-reference.vtu new file mode 100644 index 0000000000000000000000000000000000000000..e46406a57191353ed1dbeaf7f1949b562e1e15c2 --- /dev/null +++ b/test/references/stokes_channel_pseudo3d-reference.vtu @@ -0,0 +1,117 @@ +<?xml version="1.0"?> +<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian"> + <UnstructuredGrid> + <Piece NumberOfCells="60" NumberOfPoints="84"> + <CellData Scalars="p" Vectors="velocity_component (m/s)"> + <DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii"> + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + 100008 100005 100002 100008 100005 100002 100008 100005 100002 100008 100005 100002 + </DataArray> + <DataArray type="Float32" Name="rhoMolar" 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 + </DataArray> + <DataArray type="Float32" Name="rho" 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 + </DataArray> + <DataArray type="Float32" Name="velocity_component (m/s)" NumberOfComponents="3" format="ascii"> + 0.00769322 -4.68964e-15 0 0.00769322 -2.4791e-15 0 0.00769322 2.90074e-15 0 0.012927 -5.09704e-15 0 + 0.012927 5.86407e-15 0 0.012927 9.61387e-15 0 0.0147073 -4.75273e-15 0 0.0147073 1.87467e-14 0 + 0.0147073 1.81825e-14 0 0.0153128 -1.67094e-14 0 0.0153128 2.37587e-14 0 0.0153128 2.78985e-14 0 + 0.0155188 -3.37813e-14 0 0.0155188 2.9366e-14 0 0.0155188 3.92226e-14 0 0.0155889 -5.15662e-14 0 + 0.0155889 3.19764e-14 0 0.0155889 5.53428e-14 0 0.0156127 -7.18326e-14 0 0.0156127 3.26904e-14 0 + 0.0156127 7.44032e-14 0 0.0156208 -9.78366e-14 0 0.0156208 3.4182e-14 0 0.0156208 9.18211e-14 0 + 0.0156235 -1.23035e-13 0 0.0156235 3.47165e-14 0 0.0156235 1.05798e-13 0 0.0156244 -1.36701e-13 0 + 0.0156244 3.28922e-14 0 0.0156244 1.14573e-13 0 0.0156244 -1.38015e-13 0 0.0156244 2.81906e-14 0 + 0.0156244 1.19985e-13 0 0.0156235 -1.32658e-13 0 0.0156235 2.33148e-14 0 0.0156235 1.19823e-13 0 + 0.0156208 -1.19409e-13 0 0.0156208 1.82874e-14 0 0.0156208 1.12109e-13 0 0.0156127 -9.86279e-14 0 + 0.0156127 1.21082e-14 0 0.0156127 9.82026e-14 0 0.0155889 -7.49714e-14 0 0.0155889 6.91552e-15 0 + 0.0155889 7.69545e-14 0 0.0155188 -5.09466e-14 0 0.0155188 3.43333e-15 0 0.0155188 5.34615e-14 0 + 0.0153128 -3.09926e-14 0 0.0153128 1.52581e-15 0 0.0153128 3.16113e-14 0 0.0147073 -1.7591e-14 0 + 0.0147073 2.762e-15 0 0.0147073 1.38826e-14 0 0.012927 -9.33011e-15 0 0.012927 8.27631e-15 0 + 0.012927 7.28086e-15 0 0.00769322 -3.07526e-15 0 0.00769322 5.94463e-15 0 0.00769322 3.82805e-15 0 + </DataArray> + <DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii"> + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + 0 0 0 0 0 0 0 0 0 0 0 0 + </DataArray> + </CellData> + <Points> + <DataArray type="Float32" Name="Coordinates" NumberOfComponents="3" format="ascii"> + 0 0 0 0.00166667 0 0 0 0.0001 0 0.00166667 0.0001 0 + 0.00333333 0 0 0.00333333 0.0001 0 0.005 0 0 0.005 0.0001 0 + 0 0.0002 0 0.00166667 0.0002 0 0.00333333 0.0002 0 0.005 0.0002 0 + 0 0.0003 0 0.00166667 0.0003 0 0.00333333 0.0003 0 0.005 0.0003 0 + 0 0.0004 0 0.00166667 0.0004 0 0.00333333 0.0004 0 0.005 0.0004 0 + 0 0.0005 0 0.00166667 0.0005 0 0.00333333 0.0005 0 0.005 0.0005 0 + 0 0.0006 0 0.00166667 0.0006 0 0.00333333 0.0006 0 0.005 0.0006 0 + 0 0.0007 0 0.00166667 0.0007 0 0.00333333 0.0007 0 0.005 0.0007 0 + 0 0.0008 0 0.00166667 0.0008 0 0.00333333 0.0008 0 0.005 0.0008 0 + 0 0.0009 0 0.00166667 0.0009 0 0.00333333 0.0009 0 0.005 0.0009 0 + 0 0.001 0 0.00166667 0.001 0 0.00333333 0.001 0 0.005 0.001 0 + 0 0.0011 0 0.00166667 0.0011 0 0.00333333 0.0011 0 0.005 0.0011 0 + 0 0.0012 0 0.00166667 0.0012 0 0.00333333 0.0012 0 0.005 0.0012 0 + 0 0.0013 0 0.00166667 0.0013 0 0.00333333 0.0013 0 0.005 0.0013 0 + 0 0.0014 0 0.00166667 0.0014 0 0.00333333 0.0014 0 0.005 0.0014 0 + 0 0.0015 0 0.00166667 0.0015 0 0.00333333 0.0015 0 0.005 0.0015 0 + 0 0.0016 0 0.00166667 0.0016 0 0.00333333 0.0016 0 0.005 0.0016 0 + 0 0.0017 0 0.00166667 0.0017 0 0.00333333 0.0017 0 0.005 0.0017 0 + 0 0.0018 0 0.00166667 0.0018 0 0.00333333 0.0018 0 0.005 0.0018 0 + 0 0.0019 0 0.00166667 0.0019 0 0.00333333 0.0019 0 0.005 0.0019 0 + 0 0.002 0 0.00166667 0.002 0 0.00333333 0.002 0 0.005 0.002 0 + </DataArray> + </Points> + <Cells> + <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii"> + 0 1 3 2 1 4 5 3 4 6 7 5 + 2 3 9 8 3 5 10 9 5 7 11 10 + 8 9 13 12 9 10 14 13 10 11 15 14 + 12 13 17 16 13 14 18 17 14 15 19 18 + 16 17 21 20 17 18 22 21 18 19 23 22 + 20 21 25 24 21 22 26 25 22 23 27 26 + 24 25 29 28 25 26 30 29 26 27 31 30 + 28 29 33 32 29 30 34 33 30 31 35 34 + 32 33 37 36 33 34 38 37 34 35 39 38 + 36 37 41 40 37 38 42 41 38 39 43 42 + 40 41 45 44 41 42 46 45 42 43 47 46 + 44 45 49 48 45 46 50 49 46 47 51 50 + 48 49 53 52 49 50 54 53 50 51 55 54 + 52 53 57 56 53 54 58 57 54 55 59 58 + 56 57 61 60 57 58 62 61 58 59 63 62 + 60 61 65 64 61 62 66 65 62 63 67 66 + 64 65 69 68 65 66 70 69 66 67 71 70 + 68 69 73 72 69 70 74 73 70 71 75 74 + 72 73 77 76 73 74 78 77 74 75 79 78 + 76 77 81 80 77 78 82 81 78 79 83 82 + </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 + </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 + </DataArray> + </Cells> + </Piece> + </UnstructuredGrid> +</VTKFile>