diff --git a/dumux/porousmediumflow/2pnc/implicit/model.hh b/dumux/porousmediumflow/2pnc/implicit/model.hh
index b91555dbfcc2545564b22013b2775a865e0341b0..93838e82e81a648a55e6c2d560b34cf58178c5e7 100644
--- a/dumux/porousmediumflow/2pnc/implicit/model.hh
+++ b/dumux/porousmediumflow/2pnc/implicit/model.hh
@@ -19,7 +19,7 @@
/*!
* \file
*
-* \brief Adaption of the fully implicit box scheme to the two-phase n-component flow model.
+* \brief Adaption of the fully implicit model to the two-phase n-component flow model.
*/
#ifndef DUMUX_2PNC_MODEL_HH
@@ -68,7 +68,7 @@ namespace Dumux
* = q_\lambda\f$
*
* All equations are discretized using a vertex-centered finite volume (box)
- * or cell-centered finite volume scheme (this is not done for 2pnc approach yet, however possible) as
+ * or cell-centered finite volume scheme as
* spatial and the implicit Euler method as time discretization.
*
* By using constitutive relations for the capillary pressure \f$p_c =
@@ -166,7 +166,7 @@ public:
{
ParentType::init(problem);
- unsigned numDofs = this->numDofs();
+ auto numDofs = this->numDofs();
staticDat_.resize(numDofs);
@@ -174,37 +174,14 @@ public:
for (const auto& element : elements(this->gridView_()))
{
- if (!isBox) // i.e. cell-centered discretization
- {
- int dofIdxGlobal = this->dofMapper().index(element);
- const GlobalPosition &globalPos = element.geometry().center();
-
- // initialize phase presence
- staticDat_[dofIdxGlobal].phasePresence
- = this->problem_().initialPhasePresence(*(this->gridView_().template begin<dim>()),
- dofIdxGlobal, globalPos);
- staticDat_[dofIdxGlobal].wasSwitched = false;
-
- staticDat_[dofIdxGlobal].oldPhasePresence
- = staticDat_[dofIdxGlobal].phasePresence;
- }
- }
-
- if (isBox) // i.e. vertex-centered discretization
- {
- for (const auto& vertex : vertices(this->gridView_()))
+ FVElementGeometry fvGeometry;
+ fvGeometry.update(this->gridView_(), element);
+ for (unsigned int scvIdx = 0; scvIdx < fvGeometry.numScv; ++scvIdx)
{
- int dofIdxGlobal = this->dofMapper().index(vertex);
- const GlobalPosition &globalPos = vertex.geometry().corner(0);
-
- // initialize phase presence
- staticDat_[dofIdxGlobal].phasePresence
- = this->problem_().initialPhasePresence(vertex, dofIdxGlobal,
- globalPos);
+ auto dofIdxGlobal = this->dofMapper().subIndex(element, scvIdx, dofCodim);
+ staticDat_[dofIdxGlobal].phasePresence = this->problem_().initialPhasePresence(element, fvGeometry, scvIdx);
staticDat_[dofIdxGlobal].wasSwitched = false;
-
- staticDat_[dofIdxGlobal].oldPhasePresence
- = staticDat_[dofIdxGlobal].phasePresence;
+ staticDat_[dofIdxGlobal].oldPhasePresence = staticDat_[dofIdxGlobal].phasePresence;
}
}
}
@@ -271,15 +248,15 @@ public:
}
/*!
- * \brief Returns the phase presence of the current or the old solution of a vertex.
+ * \brief Returns the phase presence of the current or the old solution
*
- * \param globalVertexIdx The global vertex index
+ * \param dofIdxGlobal The global DOF index
* \param oldSol Evaluate function with solution of current or previous time step
*/
- int phasePresence(int globalVertexIdx, bool oldSol) const
+ int phasePresence(int dofIdxGlobal, bool oldSol) const
{
- return oldSol ? staticDat_[globalVertexIdx].oldPhasePresence
- : staticDat_[globalVertexIdx].phasePresence;
+ return oldSol ? staticDat_[dofIdxGlobal].oldPhasePresence
+ : staticDat_[dofIdxGlobal].phasePresence;
}
/*!
@@ -299,7 +276,7 @@ public:
typedef Dune::BlockVector<Dune::FieldVector<Scalar, dim> > VectorField;
// get the number of degrees of freedom
- unsigned numDofs = this->numDofs();
+ auto numDofs = this->numDofs();
ScalarField *Sg = writer.allocateManagedBuffer (numDofs);
ScalarField *Sl = writer.allocateManagedBuffer (numDofs);
@@ -312,7 +289,6 @@ public:
ScalarField *mobG = writer.allocateManagedBuffer (numDofs);
ScalarField *temperature = writer.allocateManagedBuffer (numDofs);
ScalarField *poro = writer.allocateManagedBuffer (numDofs);
- ScalarField *boxVolume = writer.allocateManagedBuffer (numDofs);
VectorField *velocityN = writer.template allocateManagedBuffer<double, dim>(numDofs);
VectorField *velocityW = writer.template allocateManagedBuffer<double, dim>(numDofs);
ImplicitVelocityOutput<TypeTag> velocityOutput(this->problem_());
@@ -340,21 +316,17 @@ public:
for (int j = 0; j < dim; ++j) //Permeability only in main directions xx and yy
Perm[j] = writer.allocateManagedBuffer(numDofs);
- *boxVolume = 0;
-
- unsigned numElements = this->gridView_().size(0);
+ auto numElements = this->gridView_().size(0);
ScalarField *rank = writer.allocateManagedBuffer (numElements);
- FVElementGeometry fvGeometry;
- VolumeVariables volVars;
- ElementVolumeVariables elemVolVars;
-
for (const auto& element : elements(this->gridView_()))
{
- int eIdxGlobal = this->problem_().elementMapper().index(element);
+ auto eIdxGlobal = this->problem_().elementMapper().index(element);
(*rank)[eIdxGlobal] = this->gridView_().comm().rank();
+ FVElementGeometry fvGeometry;
fvGeometry.update(this->gridView_(), element);
+ ElementVolumeVariables elemVolVars;
elemVolVars.update(this->problem_(),
element,
fvGeometry,
@@ -362,14 +334,7 @@ public:
for (int scvIdx = 0; scvIdx < fvGeometry.numScv; ++scvIdx)
{
- int dofIdxGlobal = this->dofMapper().subIndex(element, scvIdx, dofCodim);
-
- volVars.update(sol[dofIdxGlobal],
- this->problem_(),
- element,
- fvGeometry,
- scvIdx,
- false);
+ auto dofIdxGlobal = this->dofMapper().subIndex(element, scvIdx, dofCodim);
GlobalPosition globalPos = fvGeometry.subContVol[scvIdx].global;
(*Sg)[dofIdxGlobal] = elemVolVars[scvIdx].saturation(nPhaseIdx);
@@ -381,25 +346,20 @@ public:
(*rhoG)[dofIdxGlobal] = elemVolVars[scvIdx].density(nPhaseIdx);
(*mobL)[dofIdxGlobal] = elemVolVars[scvIdx].mobility(wPhaseIdx);
(*mobG)[dofIdxGlobal] = elemVolVars[scvIdx].mobility(nPhaseIdx);
- (*boxVolume)[dofIdxGlobal] += fvGeometry.subContVol[scvIdx].volume;
(*poro)[dofIdxGlobal] = elemVolVars[scvIdx].porosity();
(*temperature)[dofIdxGlobal] = elemVolVars[scvIdx].temperature();
for (int phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)
- {
for (int compIdx = 0; compIdx < numComponents; ++compIdx)
- {
- (*moleFraction[phaseIdx][compIdx])[dofIdxGlobal]= volVars.moleFraction(phaseIdx,compIdx);
- Valgrind::CheckDefined((*moleFraction[phaseIdx][compIdx])[dofIdxGlobal]);
- }
- }
+ (*moleFraction[phaseIdx][compIdx])[dofIdxGlobal]= elemVolVars[scvIdx].moleFraction(phaseIdx,compIdx);
+
for (int compIdx = 0; compIdx < numComponents; ++compIdx)
- (*molarity[compIdx])[dofIdxGlobal] = (volVars.molarity(wPhaseIdx, compIdx));
+ (*molarity[compIdx])[dofIdxGlobal] = (elemVolVars[scvIdx].molarity(wPhaseIdx, compIdx));
Tensor K = perm_(this->problem_().spatialParams().intrinsicPermeability(element, fvGeometry, scvIdx));
for (int j = 0; j<dim; ++j)
- (*Perm[j])[dofIdxGlobal] = K[j][j] /* volVars.permFactor()*/;
+ (*Perm[j])[dofIdxGlobal] = K[j][j];
};
// velocity output
@@ -421,7 +381,6 @@ public:
writer.attachDofData(*mobG, "mobG", isBox);
writer.attachDofData(*poro, "porosity", isBox);
writer.attachDofData(*temperature, "temperature", isBox);
- writer.attachDofData(*boxVolume, "boxVolume", isBox);
writer.attachDofData(*Perm[0], "Kxx", isBox);
if (dim >= 2)
writer.attachDofData(*Perm[1], "Kyy", isBox);
@@ -436,7 +395,7 @@ public:
oss << "x"
<< FluidSystem::componentName(j)
<< FluidSystem::phaseName(i);
- writer.attachDofData(*moleFraction[i][j], oss.str().c_str(), isBox);
+ writer.attachDofData(*moleFraction[i][j], oss.str(), isBox);
}
}
@@ -469,10 +428,6 @@ public:
// write primary variables
ParentType::serializeEntity(outStream, entity);
int dofIdxGlobal = this->dofMapper().index(entity);
-
- if (!outStream.good())
- DUNE_THROW(Dune::IOError, "Could not serialize vertex " << dofIdxGlobal);
-
outStream << staticDat_[dofIdxGlobal].phasePresence << " ";
}
@@ -490,14 +445,9 @@ public:
ParentType::deserializeEntity(inStream, entity);
int dofIdxGlobal = this->dofMapper().index(entity);
- if (!inStream.good())
- DUNE_THROW(Dune::IOError,
- "Could not deserialize vertex " << dofIdxGlobal);
-
inStream >> staticDat_[dofIdxGlobal].phasePresence;
staticDat_[dofIdxGlobal].oldPhasePresence
= staticDat_[dofIdxGlobal].phasePresence;
-
}
/*!
@@ -514,19 +464,19 @@ public:
for (unsigned i = 0; i < staticDat_.size(); ++i)
staticDat_[i].visited = false;
- FVElementGeometry fvGeometry;
- static VolumeVariables volVars;
for (const auto& element : elements(this->gridView_()))
{
+ FVElementGeometry fvGeometry;
fvGeometry.update(this->gridView_(), element);
for (int scvIdx = 0; scvIdx < fvGeometry.numScv; ++scvIdx)
{
- int dofIdxGlobal = this->dofMapper().subIndex(element, scvIdx, dim);
+ auto dofIdxGlobal = this->dofMapper().subIndex(element, scvIdx, dofCodim);
if (staticDat_[dofIdxGlobal].visited)
continue;
staticDat_[dofIdxGlobal].visited = true;
+ VolumeVariables volVars;
volVars.update(curGlobalSol[dofIdxGlobal],
this->problem_(),
element,
@@ -542,7 +492,8 @@ public:
// make sure that if there was a variable switch in an
// other partition we will also set the switch flag
// for our partition.
- wasSwitched = this->gridView_().comm().max(wasSwitched);
+ if (this->gridView_().comm().size() > 1)
+ wasSwitched = this->gridView_().comm().max(wasSwitched);
setSwitched_(wasSwitched);
}
@@ -583,8 +534,7 @@ protected:
*/
void resetPhasePresence_()
{
- int numDofs = this->gridView_().size(dim);
- for (int i = 0; i < numDofs; ++i)
+ for (int i = 0; i < this->numDofs(); ++i)
{
staticDat_[i].phasePresence
= staticDat_[i].oldPhasePresence;
@@ -597,8 +547,7 @@ protected:
*/
void updateOldPhasePresence_()
{
- int numDofs = this->gridView_().size(dim);
- for (int i = 0; i < numDofs; ++i)
+ for (int i = 0; i < this->numDofs(); ++i)
{
staticDat_[i].oldPhasePresence
= staticDat_[i].phasePresence;
@@ -620,7 +569,8 @@ protected:
* variable switch was performed.
*/
bool primaryVarSwitch_(SolutionVector &globalSol,
- const VolumeVariables &volVars, int dofIdxGlobal,
+ const VolumeVariables &volVars,
+ int dofIdxGlobal,
const GlobalPosition &globalPos)
{
diff --git a/dumux/porousmediumflow/2pnc/implicit/volumevariables.hh b/dumux/porousmediumflow/2pnc/implicit/volumevariables.hh
index f9b1b4420fb732b96c5e3a4907a1ff466605e2c1..b0b1ac992ce0171949c296681a144930581f790d 100644
--- a/dumux/porousmediumflow/2pnc/implicit/volumevariables.hh
+++ b/dumux/porousmediumflow/2pnc/implicit/volumevariables.hh
@@ -131,18 +131,18 @@ public:
/////////////
const MaterialLawParams &materialParams = problem.spatialParams().materialLawParams(element, fvGeometry, scvIdx);
- // Second instance of a parameter cache.
+ // Second instance of a parameter cache.
// Could be avoided if diffusion coefficients also
// became part of the fluid state.
typename FluidSystem::ParameterCache paramCache;
paramCache.updateAll(fluidState_);
for (int phaseIdx = 0; phaseIdx < numPhases; ++phaseIdx)
- {// relative permeabilities
+ {
+ // relative permeabilities
Scalar kr;
if (phaseIdx == wPhaseIdx)
kr = MaterialLaw::krw(materialParams, saturation(wPhaseIdx));
- else // ATTENTION: krn requires the liquid saturation
- // as parameter!
+ else // ATTENTION: krn requires the liquid saturation // as parameter!
kr = MaterialLaw::krn(materialParams, saturation(wPhaseIdx));
mobility_[phaseIdx] = kr / fluidState_.viscosity(phaseIdx);
@@ -181,12 +181,12 @@ public:
* \param priVars The primary Variables
*/
static void completeFluidState(const PrimaryVariables& priVars,
- const Problem& problem,
- const Element& element,
- const FVElementGeometry& fvGeometry,
- int scvIdx,
- FluidState& fluidState,
- bool isOldSol = false)
+ const Problem& problem,
+ const Element& element,
+ const FVElementGeometry& fvGeometry,
+ int scvIdx,
+ FluidState& fluidState,
+ bool isOldSol = false)
{
Scalar t = Implementation::temperature_(priVars, problem, element,
diff --git a/test/porousmediumflow/2pnc/implicit/fuelcellproblem.hh b/test/porousmediumflow/2pnc/implicit/fuelcellproblem.hh
index 1f8501e8ede523db52219745fbf3806ad8399054..d49d7075e5ced201d5e0aabb89187b9370e3b278 100644
--- a/test/porousmediumflow/2pnc/implicit/fuelcellproblem.hh
+++ b/test/porousmediumflow/2pnc/implicit/fuelcellproblem.hh
@@ -188,14 +188,18 @@ public:
const int scvIdx,
const ElementVolumeVariables &elemVolVars) const
{
- const GlobalPosition globalPos = element.geometry().corner(scvIdx);
- const VolumeVariables &volVars = elemVolVars[scvIdx];
-
values = 0.0;
+ const auto& globalPos = isBox ? element.geometry().corner(scvIdx)
+ : element.geometry().center();
+
//reaction sources from electro chemistry
- if(onLowerBoundary_(globalPos))
- ElectroChemistry::reactionSource(values, volVars);
+ if(inReactionLayer_(globalPos))
+ {
+ const auto& volVars = elemVolVars[scvIdx];
+ auto currentDensity = ElectroChemistry::calculateCurrentDensity(volVars);
+ ElectroChemistry::reactionSource(values, currentDensity);
+ }
}
@@ -287,16 +291,16 @@ public:
}
/*!
- * \brief Return the initial phase state inside a control volume.
+ * \brief Return the initial phase state inside a sub control volume.
*
- * \param vertex The vertex
- * \param globalIdx The index of the global vertex
- * \param globalPos The global position
+ * \param element The element of the sub control volume
+ * \param fvGeometry The finite volume geometry
+ * \param scvIdx The sub control volume index
*/
- int initialPhasePresence(const Vertex &vertex,
- int &globalIdx,
- const GlobalPosition &globalPos) const
+ int initialPhasePresence(const Element &element,
+ const FVElementGeometry &fvGeometry,
+ int scvIdx) const
{
return Indices::bothPhases;
}
@@ -310,7 +314,7 @@ public:
typedef Dune::BlockVector<Dune::FieldVector<Scalar, 1> > ScalarField;
// get the number of degrees of freedom
- unsigned numDofs = this->model().numDofs();
+ auto numDofs = this->model().numDofs();
// create the required scalar fields
ScalarField *currentDensity = this->resultWriter().allocateManagedBuffer (numDofs);
@@ -330,28 +334,32 @@ public:
for (int scvIdx = 0; scvIdx < fvGeometry.numScv; ++scvIdx)
{
- int dofIdxGlobal = this->model().dofMapper().subIndex(element, scvIdx, dofCodim);
-
- //reactionSource Output
- PrimaryVariables source;
- this->solDependentSource(source, element, fvGeometry, scvIdx, elemVolVars);
-
- (*reactionSourceH2O)[dofIdxGlobal] = source[wPhaseIdx];
- (*reactionSourceO2)[dofIdxGlobal] = source[numComponents-1];
-
- //Current Output
- (*currentDensity)[dofIdxGlobal] = -1.0*source[numComponents-1]*4*Constant::F;
- //recorrection of the area for output
- Scalar gridYMin = 0.0;
- Scalar gridYMax = GET_RUNTIME_PARAM(TypeTag, Scalar, Grid.UpperRightY);
- Scalar nCellsY = GET_RUNTIME_PARAM(TypeTag, Scalar, Grid.NumberOfCellsY);
-
- Scalar lengthBox = (gridYMax - gridYMin)/nCellsY;
-
- // Only works for box
- (*currentDensity)[dofIdxGlobal] = (*currentDensity)[dofIdxGlobal]/2*lengthBox/10000; //i in [A/cm^2]
+ const auto& globalPos = isBox ? element.geometry().corner(scvIdx)
+ : element.geometry().center();
+
+ auto dofIdxGlobal = this->model().dofMapper().subIndex(element, scvIdx, dofCodim);
+
+ //reaction sources from electro chemistry
+ if(inReactionLayer_(globalPos))
+ {
+ //reactionSource Output
+ PrimaryVariables source;
+ auto i = ElectroChemistry::calculateCurrentDensity(elemVolVars[scvIdx]);
+ ElectroChemistry::reactionSource(source, i);
+
+ (*reactionSourceH2O)[dofIdxGlobal] = source[wPhaseIdx];
+ (*reactionSourceO2)[dofIdxGlobal] = source[numComponents-1];
+
+ //Current Output in A/cm^2
+ (*currentDensity)[dofIdxGlobal] = i/10000;
+ }
+ else
+ {
+ (*reactionSourceH2O)[dofIdxGlobal] = 0.0;
+ (*reactionSourceO2)[dofIdxGlobal] = 0.0;
+ (*currentDensity)[dofIdxGlobal] = 0.0;
+ }
}
-
}
this->resultWriter().attachDofData(*reactionSourceH2O, "reactionSourceH2O [mol/(sm^2)]", isBox);
@@ -382,6 +390,9 @@ private:
bool onUpperBoundary_(const GlobalPosition &globalPos) const
{ return globalPos[1] > this->bBoxMax()[1] - eps_; }
+ bool inReactionLayer_(const GlobalPosition& globalPos) const
+ { return globalPos[1] < 0.1*(this->bBoxMax()[1] - this->bBoxMin()[1]) + eps_; }
+
Scalar temperature_;
Scalar eps_;
int nTemperature_;
@@ -391,6 +402,7 @@ private:
Scalar temperatureLow_, temperatureHigh_;
Scalar pO2Inlet_;
};
+
} //end namespace
#endif
diff --git a/test/references/fuelcell2pncbox-reference.vtu b/test/references/fuelcell2pncbox-reference.vtu
index 03da97e532f70e568a89557e4c2ca893f1ef816a..80abd5597f0a53e04d8b3ba3d256a80b15c63be5 100644
--- a/test/references/fuelcell2pncbox-reference.vtu
+++ b/test/references/fuelcell2pncbox-reference.vtu
@@ -168,21 +168,6 @@
310 310 310 310 310 310 310 310 310 310 310 310
310 310 310 310 310 310 310 310 310 310
</DataArray>
- <DataArray type="Float32" Name="boxVolume" NumberOfComponents="1" format="ascii">
- 2.18254e-09 4.36508e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09
- 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09
- 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09
- 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 4.36508e-09 8.73016e-09 2.18254e-09 4.36508e-09 4.36508e-09 8.73016e-09 8.73016e-09 8.73016e-09
- 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09
- 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 4.36508e-09 4.36508e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09
- 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09
- 8.73016e-09 8.73016e-09 8.73016e-09 4.36508e-09 4.36508e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09
- 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09
- 8.73016e-09 4.36508e-09 4.36508e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09
- 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 8.73016e-09 4.36508e-09
- 2.18254e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09
- 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 4.36508e-09 2.18254e-09
- </DataArray>
<DataArray type="Float32" Name="reactionSourceH2O [mol/(sm^2)]" NumberOfComponents="1" format="ascii">
747.822 747.822 0 0 747.822 0 747.822 0 747.822 0 747.822 0
747.822 0 747.822 0 747.822 0 747.822 0 747.822 0 747.822 0
diff --git a/test/references/fuelcell2pnccc-reference.vtu b/test/references/fuelcell2pnccc-reference.vtu
new file mode 100644
index 0000000000000000000000000000000000000000..2071e1e1a87c785ac31d40349c80b134daf9c9f8
--- /dev/null
+++ b/test/references/fuelcell2pnccc-reference.vtu
@@ -0,0 +1,462 @@
+<?xml version="1.0"?>
+<VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian">
+ <UnstructuredGrid>
+ <Piece NumberOfCells="126" NumberOfPoints="154">
+ <CellData Scalars="Sg">
+ <DataArray type="Float32" Name="Sg" NumberOfComponents="1" format="ascii">
+ 0.69999 0.69999 0.69999 0.69999 0.69999 0.69999 0.69999 0.69999 0.69999 0.69999 0.69999 0.69999
+ 0.69999 0.69999 0.69999 0.69999 0.69999 0.69999 0.69999 0.69999 0.69999 0.699992 0.699992 0.699992
+ 0.699992 0.699992 0.699992 0.699992 0.699992 0.699992 0.699992 0.699992 0.699992 0.699992 0.699992 0.699992
+ 0.699992 0.699992 0.699992 0.699992 0.699992 0.699992 0.699994 0.699994 0.699994 0.699994 0.699994 0.699994
+ 0.699994 0.699994 0.699994 0.699994 0.699994 0.699994 0.699994 0.699994 0.699994 0.699994 0.699994 0.699994
+ 0.699994 0.699994 0.699994 0.699995 0.699995 0.699995 0.699995 0.699995 0.699995 0.699995 0.699995 0.699995
+ 0.699995 0.699995 0.699995 0.699995 0.699995 0.699995 0.699995 0.699995 0.699995 0.699995 0.699995 0.699995
+ 0.699997 0.699997 0.699997 0.699997 0.699997 0.699997 0.699997 0.699997 0.699997 0.699997 0.699997 0.699997
+ 0.699997 0.699997 0.699997 0.699997 0.699997 0.699997 0.699997 0.699997 0.699997 0.699999 0.699999 0.699999
+ 0.699999 0.699999 0.699999 0.699999 0.699999 0.699999 0.699999 0.699999 0.699999 0.699999 0.699999 0.699999
+ 0.699999 0.699999 0.699999 0.699999 0.699999 0.699999
+ </DataArray>
+ <DataArray type="Float32" Name="Sl" NumberOfComponents="1" format="ascii">
+ 0.30001 0.30001 0.30001 0.30001 0.30001 0.30001 0.30001 0.30001 0.30001 0.30001 0.30001 0.30001
+ 0.30001 0.30001 0.30001 0.30001 0.30001 0.30001 0.30001 0.30001 0.30001 0.300008 0.300008 0.300008
+ 0.300008 0.300008 0.300008 0.300008 0.300008 0.300008 0.300008 0.300008 0.300008 0.300008 0.300008 0.300008
+ 0.300008 0.300008 0.300008 0.300008 0.300008 0.300008 0.300006 0.300006 0.300006 0.300006 0.300006 0.300006
+ 0.300006 0.300006 0.300006 0.300006 0.300006 0.300006 0.300006 0.300006 0.300006 0.300006 0.300006 0.300006
+ 0.300006 0.300006 0.300006 0.300005 0.300005 0.300005 0.300005 0.300005 0.300005 0.300005 0.300005 0.300005
+ 0.300005 0.300005 0.300005 0.300005 0.300005 0.300005 0.300005 0.300005 0.300005 0.300005 0.300005 0.300005
+ 0.300003 0.300003 0.300003 0.300003 0.300003 0.300003 0.300003 0.300003 0.300003 0.300003 0.300003 0.300003
+ 0.300003 0.300003 0.300003 0.300003 0.300003 0.300003 0.300003 0.300003 0.300003 0.300001 0.300001 0.300001
+ 0.300001 0.300001 0.300001 0.300001 0.300001 0.300001 0.300001 0.300001 0.300001 0.300001 0.300001 0.300001
+ 0.300001 0.300001 0.300001 0.300001 0.300001 0.300001
+ </DataArray>
+ <DataArray type="Float32" Name="pg" NumberOfComponents="1" format="ascii">
+ 99999.8 99999.8 99999.8 99999.8 99999.8 99999.8 99999.8 99999.8 99999.8 99999.8 99999.8 99999.8
+ 99999.8 99999.8 99999.8 99999.8 99999.8 99999.8 99999.8 99999.8 99999.8 99999.9 99999.9 99999.9
+ 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9
+ 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9
+ 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9
+ 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9
+ 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9 99999.9
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000 100000
+ 100000 100000 100000 100000 100000 100000
+ </DataArray>
+ <DataArray type="Float32" Name="pl" NumberOfComponents="1" format="ascii">
+ 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440
+ 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440
+ 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440
+ 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440
+ 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440
+ 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440
+ 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440
+ 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440
+ 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440
+ 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440 111440
+ 111440 111440 111440 111440 111440 111440
+ </DataArray>
+ <DataArray type="Float32" Name="pc" NumberOfComponents="1" format="ascii">
+ -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.4
+ -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.4 -11440.3 -11440.3 -11440.3
+ -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3
+ -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3
+ -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3 -11440.3
+ -11440.3 -11440.3 -11440.3 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2
+ -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2
+ -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2
+ -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2
+ -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2
+ -11440.2 -11440.2 -11440.2 -11440.2 -11440.2 -11440.2
+ </DataArray>
+ <DataArray type="Float32" Name="rhoL" NumberOfComponents="1" format="ascii">
+ 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402
+ 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402
+ 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402
+ 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402
+ 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402 993.402
+ 993.402 993.402 993.402 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403
+ 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403
+ 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403
+ 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403
+ 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403 993.403
+ 993.403 993.403 993.403 993.403 993.403 993.403
+ </DataArray>
+ <DataArray type="Float32" Name="rhoG" NumberOfComponents="1" format="ascii">
+ 1.09924 1.09924 1.09924 1.09924 1.09924 1.09924 1.09924 1.09924 1.09924 1.09924 1.09924 1.09924
+ 1.09924 1.09924 1.09924 1.09924 1.09924 1.09924 1.09924 1.09924 1.09924 1.10315 1.10315 1.10315
+ 1.10315 1.10315 1.10315 1.10315 1.10315 1.10315 1.10315 1.10315 1.10315 1.10315 1.10315 1.10315
+ 1.10315 1.10315 1.10315 1.10315 1.10315 1.10315 1.10691 1.10691 1.10691 1.10691 1.10691 1.10691
+ 1.10691 1.10691 1.10691 1.10691 1.10691 1.10691 1.10691 1.10691 1.10691 1.10691 1.10691 1.10691
+ 1.10691 1.10691 1.10691 1.11054 1.11054 1.11054 1.11054 1.11054 1.11054 1.11054 1.11054 1.11054
+ 1.11054 1.11054 1.11054 1.11054 1.11054 1.11054 1.11054 1.11054 1.11054 1.11054 1.11054 1.11054
+ 1.11404 1.11404 1.11404 1.11404 1.11404 1.11404 1.11404 1.11404 1.11404 1.11404 1.11404 1.11404
+ 1.11404 1.11404 1.11404 1.11404 1.11404 1.11404 1.11404 1.11404 1.11404 1.11742 1.11742 1.11742
+ 1.11742 1.11742 1.11742 1.11742 1.11742 1.11742 1.11742 1.11742 1.11742 1.11742 1.11742 1.11742
+ 1.11742 1.11742 1.11742 1.11742 1.11742 1.11742
+ </DataArray>
+ <DataArray type="Float32" Name="mobL" NumberOfComponents="1" format="ascii">
+ 129.817 129.817 129.817 129.817 129.817 129.817 129.817 129.817 129.817 129.817 129.817 129.817
+ 129.817 129.817 129.817 129.817 129.817 129.817 129.817 129.817 129.817 129.815 129.815 129.815
+ 129.815 129.815 129.815 129.815 129.815 129.815 129.815 129.815 129.815 129.815 129.815 129.815
+ 129.815 129.815 129.815 129.815 129.815 129.815 129.814 129.814 129.814 129.814 129.814 129.814
+ 129.814 129.814 129.814 129.814 129.814 129.814 129.814 129.814 129.814 129.814 129.814 129.814
+ 129.814 129.814 129.814 129.812 129.812 129.812 129.812 129.812 129.812 129.812 129.812 129.812
+ 129.812 129.812 129.812 129.812 129.812 129.812 129.812 129.812 129.812 129.812 129.812 129.812
+ 129.81 129.81 129.81 129.81 129.81 129.81 129.81 129.81 129.81 129.81 129.81 129.81
+ 129.81 129.81 129.81 129.81 129.81 129.81 129.81 129.81 129.81 129.809 129.809 129.809
+ 129.809 129.809 129.809 129.809 129.809 129.809 129.809 129.809 129.809 129.809 129.809 129.809
+ 129.809 129.809 129.809 129.809 129.809 129.809
+ </DataArray>
+ <DataArray type="Float32" Name="mobG" NumberOfComponents="1" format="ascii">
+ 18729.3 18729.3 18729.3 18729.3 18729.3 18729.3 18729.3 18729.3 18729.3 18729.3 18729.3 18729.3
+ 18729.3 18729.3 18729.3 18729.3 18729.3 18729.3 18729.3 18729.3 18729.3 18654.1 18654.1 18654.1
+ 18654.1 18654.1 18654.1 18654.1 18654.1 18654.1 18654.1 18654.1 18654.1 18654.1 18654.1 18654.1
+ 18654.1 18654.1 18654.1 18654.1 18654.1 18654.1 18582.2 18582.2 18582.2 18582.2 18582.2 18582.2
+ 18582.2 18582.2 18582.2 18582.2 18582.2 18582.2 18582.2 18582.2 18582.2 18582.2 18582.2 18582.2
+ 18582.2 18582.2 18582.2 18513.5 18513.5 18513.5 18513.5 18513.5 18513.5 18513.5 18513.5 18513.5
+ 18513.5 18513.5 18513.5 18513.5 18513.5 18513.5 18513.5 18513.5 18513.5 18513.5 18513.5 18513.5
+ 18447.7 18447.7 18447.7 18447.7 18447.7 18447.7 18447.7 18447.7 18447.7 18447.7 18447.7 18447.7
+ 18447.7 18447.7 18447.7 18447.7 18447.7 18447.7 18447.7 18447.7 18447.7 18384.7 18384.7 18384.7
+ 18384.7 18384.7 18384.7 18384.7 18384.7 18384.7 18384.7 18384.7 18384.7 18384.7 18384.7 18384.7
+ 18384.7 18384.7 18384.7 18384.7 18384.7 18384.7
+ </DataArray>
+ <DataArray type="Float32" Name="porosity" NumberOfComponents="1" format="ascii">
+ 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
+ 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
+ 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
+ 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
+ 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
+ 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
+ 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
+ 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
+ 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
+ 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2 0.2
+ 0.2 0.2 0.2 0.2 0.2 0.2
+ </DataArray>
+ <DataArray type="Float32" Name="temperature" NumberOfComponents="1" format="ascii">
+ 310 310 310 310 310 310 310 310 310 310 310 310
+ 310 310 310 310 310 310 310 310 310 310 310 310
+ 310 310 310 310 310 310 310 310 310 310 310 310
+ 310 310 310 310 310 310 310 310 310 310 310 310
+ 310 310 310 310 310 310 310 310 310 310 310 310
+ 310 310 310 310 310 310 310 310 310 310 310 310
+ 310 310 310 310 310 310 310 310 310 310 310 310
+ 310 310 310 310 310 310 310 310 310 310 310 310
+ 310 310 310 310 310 310 310 310 310 310 310 310
+ 310 310 310 310 310 310 310 310 310 310 310 310
+ 310 310 310 310 310 310
+ </DataArray>
+ <DataArray type="Float32" Name="Kxx" NumberOfComponents="1" format="ascii">
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ </DataArray>
+ <DataArray type="Float32" Name="Kyy" NumberOfComponents="1" format="ascii">
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ 5e-11 5e-11 5e-11 5e-11 5e-11 5e-11
+ </DataArray>
+ <DataArray type="Float32" Name="xH2Ol" NumberOfComponents="1" format="ascii">
+ 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988
+ 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988
+ 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988
+ 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988
+ 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988
+ 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988
+ 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988
+ 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988
+ 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999988 0.999987 0.999987 0.999987
+ 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987
+ 0.999987 0.999987 0.999987 0.999987 0.999987 0.999987
+ </DataArray>
+ <DataArray type="Float32" Name="xN2l" NumberOfComponents="1" format="ascii">
+ 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06
+ 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 7.06157e-06 6.80754e-06 6.80754e-06 6.80754e-06
+ 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06
+ 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06 6.80754e-06 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06
+ 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06 6.56266e-06
+ 6.56266e-06 6.56266e-06 6.56266e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06
+ 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06 6.32658e-06
+ 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06
+ 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 6.099e-06 5.87961e-06 5.87961e-06 5.87961e-06
+ 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06
+ 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06 5.87961e-06
+ </DataArray>
+ <DataArray type="Float32" Name="xO2l" NumberOfComponents="1" format="ascii">
+ 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06
+ 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.50297e-06 4.98546e-06 4.98546e-06 4.98546e-06
+ 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06
+ 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06 4.98546e-06 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06
+ 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06 5.45059e-06
+ 5.45059e-06 5.45059e-06 5.45059e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06
+ 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06 5.89898e-06
+ 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06
+ 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.33125e-06 6.74795e-06 6.74795e-06 6.74795e-06
+ 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06
+ 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06 6.74795e-06
+ </DataArray>
+ <DataArray type="Float32" Name="xH2Og" NumberOfComponents="1" format="ascii">
+ 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062
+ 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623062 0.0623061 0.0623061 0.0623061
+ 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061
+ 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061
+ 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061
+ 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061
+ 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061
+ 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061
+ 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.0623061 0.062306 0.062306 0.062306
+ 0.062306 0.062306 0.062306 0.062306 0.062306 0.062306 0.062306 0.062306 0.062306 0.062306 0.062306 0.062306
+ 0.062306 0.062306 0.062306 0.062306 0.062306 0.062306
+ </DataArray>
+ <DataArray type="Float32" Name="xN2g" NumberOfComponents="1" format="ascii">
+ 0.702006 0.702006 0.702006 0.702006 0.702006 0.702006 0.702006 0.702006 0.702006 0.702006 0.702006 0.702006
+ 0.702006 0.702006 0.702006 0.702006 0.702006 0.702006 0.702006 0.702006 0.702006 0.676753 0.676753 0.676753
+ 0.676753 0.676753 0.676753 0.676753 0.676753 0.676753 0.676753 0.676753 0.676753 0.676753 0.676753 0.676753
+ 0.676753 0.676753 0.676753 0.676753 0.676753 0.676753 0.652408 0.652408 0.652408 0.652408 0.652408 0.652408
+ 0.652408 0.652408 0.652408 0.652408 0.652408 0.652408 0.652408 0.652408 0.652408 0.652408 0.652408 0.652408
+ 0.652408 0.652408 0.652408 0.628939 0.628939 0.628939 0.628939 0.628939 0.628939 0.628939 0.628939 0.628939
+ 0.628939 0.628939 0.628939 0.628939 0.628939 0.628939 0.628939 0.628939 0.628939 0.628939 0.628939 0.628939
+ 0.606315 0.606315 0.606315 0.606315 0.606315 0.606315 0.606315 0.606315 0.606315 0.606315 0.606315 0.606315
+ 0.606315 0.606315 0.606315 0.606315 0.606315 0.606315 0.606315 0.606315 0.606315 0.584504 0.584504 0.584504
+ 0.584504 0.584504 0.584504 0.584504 0.584504 0.584504 0.584504 0.584504 0.584504 0.584504 0.584504 0.584504
+ 0.584504 0.584504 0.584504 0.584504 0.584504 0.584504
+ </DataArray>
+ <DataArray type="Float32" Name="xO2g" NumberOfComponents="1" format="ascii">
+ 0.235687 0.235687 0.235687 0.235687 0.235687 0.235687 0.235687 0.235687 0.235687 0.235687 0.235687 0.235687
+ 0.235687 0.235687 0.235687 0.235687 0.235687 0.235687 0.235687 0.235687 0.235687 0.260941 0.260941 0.260941
+ 0.260941 0.260941 0.260941 0.260941 0.260941 0.260941 0.260941 0.260941 0.260941 0.260941 0.260941 0.260941
+ 0.260941 0.260941 0.260941 0.260941 0.260941 0.260941 0.285286 0.285286 0.285286 0.285286 0.285286 0.285286
+ 0.285286 0.285286 0.285286 0.285286 0.285286 0.285286 0.285286 0.285286 0.285286 0.285286 0.285286 0.285286
+ 0.285286 0.285286 0.285286 0.308755 0.308755 0.308755 0.308755 0.308755 0.308755 0.308755 0.308755 0.308755
+ 0.308755 0.308755 0.308755 0.308755 0.308755 0.308755 0.308755 0.308755 0.308755 0.308755 0.308755 0.308755
+ 0.331379 0.331379 0.331379 0.331379 0.331379 0.331379 0.331379 0.331379 0.331379 0.331379 0.331379 0.331379
+ 0.331379 0.331379 0.331379 0.331379 0.331379 0.331379 0.331379 0.331379 0.331379 0.35319 0.35319 0.35319
+ 0.35319 0.35319 0.35319 0.35319 0.35319 0.35319 0.35319 0.35319 0.35319 0.35319 0.35319 0.35319
+ 0.35319 0.35319 0.35319 0.35319 0.35319 0.35319
+ </DataArray>
+ <DataArray type="Float32" Name="m^w_H2O" NumberOfComponents="1" format="ascii">
+ 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4
+ 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4
+ 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4
+ 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4
+ 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4
+ 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4
+ 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4
+ 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4
+ 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4
+ 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4
+ 55141.4 55141.4 55141.4 55141.4 55141.4 55141.4
+ </DataArray>
+ <DataArray type="Float32" Name="m^w_N2" NumberOfComponents="1" format="ascii">
+ 0.38939 0.38939 0.38939 0.38939 0.38939 0.38939 0.38939 0.38939 0.38939 0.38939 0.38939 0.38939
+ 0.38939 0.38939 0.38939 0.38939 0.38939 0.38939 0.38939 0.38939 0.38939 0.375382 0.375382 0.375382
+ 0.375382 0.375382 0.375382 0.375382 0.375382 0.375382 0.375382 0.375382 0.375382 0.375382 0.375382 0.375382
+ 0.375382 0.375382 0.375382 0.375382 0.375382 0.375382 0.361879 0.361879 0.361879 0.361879 0.361879 0.361879
+ 0.361879 0.361879 0.361879 0.361879 0.361879 0.361879 0.361879 0.361879 0.361879 0.361879 0.361879 0.361879
+ 0.361879 0.361879 0.361879 0.348861 0.348861 0.348861 0.348861 0.348861 0.348861 0.348861 0.348861 0.348861
+ 0.348861 0.348861 0.348861 0.348861 0.348861 0.348861 0.348861 0.348861 0.348861 0.348861 0.348861 0.348861
+ 0.336312 0.336312 0.336312 0.336312 0.336312 0.336312 0.336312 0.336312 0.336312 0.336312 0.336312 0.336312
+ 0.336312 0.336312 0.336312 0.336312 0.336312 0.336312 0.336312 0.336312 0.336312 0.324214 0.324214 0.324214
+ 0.324214 0.324214 0.324214 0.324214 0.324214 0.324214 0.324214 0.324214 0.324214 0.324214 0.324214 0.324214
+ 0.324214 0.324214 0.324214 0.324214 0.324214 0.324214
+ </DataArray>
+ <DataArray type="Float32" Name="m^w_O2" NumberOfComponents="1" format="ascii">
+ 0.248303 0.248303 0.248303 0.248303 0.248303 0.248303 0.248303 0.248303 0.248303 0.248303 0.248303 0.248303
+ 0.248303 0.248303 0.248303 0.248303 0.248303 0.248303 0.248303 0.248303 0.248303 0.274909 0.274909 0.274909
+ 0.274909 0.274909 0.274909 0.274909 0.274909 0.274909 0.274909 0.274909 0.274909 0.274909 0.274909 0.274909
+ 0.274909 0.274909 0.274909 0.274909 0.274909 0.274909 0.300557 0.300557 0.300557 0.300557 0.300557 0.300557
+ 0.300557 0.300557 0.300557 0.300557 0.300557 0.300557 0.300557 0.300557 0.300557 0.300557 0.300557 0.300557
+ 0.300557 0.300557 0.300557 0.325282 0.325282 0.325282 0.325282 0.325282 0.325282 0.325282 0.325282 0.325282
+ 0.325282 0.325282 0.325282 0.325282 0.325282 0.325282 0.325282 0.325282 0.325282 0.325282 0.325282 0.325282
+ 0.349118 0.349118 0.349118 0.349118 0.349118 0.349118 0.349118 0.349118 0.349118 0.349118 0.349118 0.349118
+ 0.349118 0.349118 0.349118 0.349118 0.349118 0.349118 0.349118 0.349118 0.349118 0.372096 0.372096 0.372096
+ 0.372096 0.372096 0.372096 0.372096 0.372096 0.372096 0.372096 0.372096 0.372096 0.372096 0.372096 0.372096
+ 0.372096 0.372096 0.372096 0.372096 0.372096 0.372096
+ </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
+ </DataArray>
+ <DataArray type="Float32" Name="reactionSourceH2O [mol/(sm^2)]" NumberOfComponents="1" format="ascii">
+ 379.43 379.43 379.43 379.43 379.43 379.43 379.43 379.43 379.43 379.43 379.43 379.43
+ 379.43 379.43 379.43 379.43 379.43 379.43 379.43 379.43 379.43 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0
+ </DataArray>
+ <DataArray type="Float32" Name="reactionSourceO2 [mol/(sm^2)]" NumberOfComponents="1" format="ascii">
+ -129.463 -129.463 -129.463 -129.463 -129.463 -129.463 -129.463 -129.463 -129.463 -129.463 -129.463 -129.463
+ -129.463 -129.463 -129.463 -129.463 -129.463 -129.463 -129.463 -129.463 -129.463 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0
+ </DataArray>
+ <DataArray type="Float32" Name="currentDensity [A/cm^2]" NumberOfComponents="1" format="ascii">
+ 0.458013 0.458013 0.458013 0.458013 0.458013 0.458013 0.458013 0.458013 0.458013 0.458013 0.458013 0.458013
+ 0.458013 0.458013 0.458013 0.458013 0.458013 0.458013 0.458013 0.458013 0.458013 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 0 0 0 0 0 0
+ 0 0 0 0 0 0 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 9.52381e-05 0 0 0 9.16667e-05 0 9.52381e-05 9.16667e-05 0
+ 0.000190476 0 0 0.000190476 9.16667e-05 0 0.000285714 0 0 0.000285714 9.16667e-05 0
+ 0.000380952 0 0 0.000380952 9.16667e-05 0 0.00047619 0 0 0.00047619 9.16667e-05 0
+ 0.000571429 0 0 0.000571429 9.16667e-05 0 0.000666667 0 0 0.000666667 9.16667e-05 0
+ 0.000761905 0 0 0.000761905 9.16667e-05 0 0.000857143 0 0 0.000857143 9.16667e-05 0
+ 0.000952381 0 0 0.000952381 9.16667e-05 0 0.00104762 0 0 0.00104762 9.16667e-05 0
+ 0.00114286 0 0 0.00114286 9.16667e-05 0 0.0012381 0 0 0.0012381 9.16667e-05 0
+ 0.00133333 0 0 0.00133333 9.16667e-05 0 0.00142857 0 0 0.00142857 9.16667e-05 0
+ 0.00152381 0 0 0.00152381 9.16667e-05 0 0.00161905 0 0 0.00161905 9.16667e-05 0
+ 0.00171429 0 0 0.00171429 9.16667e-05 0 0.00180952 0 0 0.00180952 9.16667e-05 0
+ 0.00190476 0 0 0.00190476 9.16667e-05 0 0.002 0 0 0.002 9.16667e-05 0
+ 0 0.000183333 0 9.52381e-05 0.000183333 0 0.000190476 0.000183333 0 0.000285714 0.000183333 0
+ 0.000380952 0.000183333 0 0.00047619 0.000183333 0 0.000571429 0.000183333 0 0.000666667 0.000183333 0
+ 0.000761905 0.000183333 0 0.000857143 0.000183333 0 0.000952381 0.000183333 0 0.00104762 0.000183333 0
+ 0.00114286 0.000183333 0 0.0012381 0.000183333 0 0.00133333 0.000183333 0 0.00142857 0.000183333 0
+ 0.00152381 0.000183333 0 0.00161905 0.000183333 0 0.00171429 0.000183333 0 0.00180952 0.000183333 0
+ 0.00190476 0.000183333 0 0.002 0.000183333 0 0 0.000275 0 9.52381e-05 0.000275 0
+ 0.000190476 0.000275 0 0.000285714 0.000275 0 0.000380952 0.000275 0 0.00047619 0.000275 0
+ 0.000571429 0.000275 0 0.000666667 0.000275 0 0.000761905 0.000275 0 0.000857143 0.000275 0
+ 0.000952381 0.000275 0 0.00104762 0.000275 0 0.00114286 0.000275 0 0.0012381 0.000275 0
+ 0.00133333 0.000275 0 0.00142857 0.000275 0 0.00152381 0.000275 0 0.00161905 0.000275 0
+ 0.00171429 0.000275 0 0.00180952 0.000275 0 0.00190476 0.000275 0 0.002 0.000275 0
+ 0 0.000366667 0 9.52381e-05 0.000366667 0 0.000190476 0.000366667 0 0.000285714 0.000366667 0
+ 0.000380952 0.000366667 0 0.00047619 0.000366667 0 0.000571429 0.000366667 0 0.000666667 0.000366667 0
+ 0.000761905 0.000366667 0 0.000857143 0.000366667 0 0.000952381 0.000366667 0 0.00104762 0.000366667 0
+ 0.00114286 0.000366667 0 0.0012381 0.000366667 0 0.00133333 0.000366667 0 0.00142857 0.000366667 0
+ 0.00152381 0.000366667 0 0.00161905 0.000366667 0 0.00171429 0.000366667 0 0.00180952 0.000366667 0
+ 0.00190476 0.000366667 0 0.002 0.000366667 0 0 0.000458333 0 9.52381e-05 0.000458333 0
+ 0.000190476 0.000458333 0 0.000285714 0.000458333 0 0.000380952 0.000458333 0 0.00047619 0.000458333 0
+ 0.000571429 0.000458333 0 0.000666667 0.000458333 0 0.000761905 0.000458333 0 0.000857143 0.000458333 0
+ 0.000952381 0.000458333 0 0.00104762 0.000458333 0 0.00114286 0.000458333 0 0.0012381 0.000458333 0
+ 0.00133333 0.000458333 0 0.00142857 0.000458333 0 0.00152381 0.000458333 0 0.00161905 0.000458333 0
+ 0.00171429 0.000458333 0 0.00180952 0.000458333 0 0.00190476 0.000458333 0 0.002 0.000458333 0
+ 0 0.00055 0 9.52381e-05 0.00055 0 0.000190476 0.00055 0 0.000285714 0.00055 0
+ 0.000380952 0.00055 0 0.00047619 0.00055 0 0.000571429 0.00055 0 0.000666667 0.00055 0
+ 0.000761905 0.00055 0 0.000857143 0.00055 0 0.000952381 0.00055 0 0.00104762 0.00055 0
+ 0.00114286 0.00055 0 0.0012381 0.00055 0 0.00133333 0.00055 0 0.00142857 0.00055 0
+ 0.00152381 0.00055 0 0.00161905 0.00055 0 0.00171429 0.00055 0 0.00180952 0.00055 0
+ 0.00190476 0.00055 0 0.002 0.00055 0
+ </DataArray>
+ </Points>
+ <Cells>
+ <DataArray type="Int32" Name="connectivity" NumberOfComponents="1" format="ascii">
+ 0 1 3 2 1 4 5 3 4 6 7 5
+ 6 8 9 7 8 10 11 9 10 12 13 11
+ 12 14 15 13 14 16 17 15 16 18 19 17
+ 18 20 21 19 20 22 23 21 22 24 25 23
+ 24 26 27 25 26 28 29 27 28 30 31 29
+ 30 32 33 31 32 34 35 33 34 36 37 35
+ 36 38 39 37 38 40 41 39 40 42 43 41
+ 2 3 45 44 3 5 46 45 5 7 47 46
+ 7 9 48 47 9 11 49 48 11 13 50 49
+ 13 15 51 50 15 17 52 51 17 19 53 52
+ 19 21 54 53 21 23 55 54 23 25 56 55
+ 25 27 57 56 27 29 58 57 29 31 59 58
+ 31 33 60 59 33 35 61 60 35 37 62 61
+ 37 39 63 62 39 41 64 63 41 43 65 64
+ 44 45 67 66 45 46 68 67 46 47 69 68
+ 47 48 70 69 48 49 71 70 49 50 72 71
+ 50 51 73 72 51 52 74 73 52 53 75 74
+ 53 54 76 75 54 55 77 76 55 56 78 77
+ 56 57 79 78 57 58 80 79 58 59 81 80
+ 59 60 82 81 60 61 83 82 61 62 84 83
+ 62 63 85 84 63 64 86 85 64 65 87 86
+ 66 67 89 88 67 68 90 89 68 69 91 90
+ 69 70 92 91 70 71 93 92 71 72 94 93
+ 72 73 95 94 73 74 96 95 74 75 97 96
+ 75 76 98 97 76 77 99 98 77 78 100 99
+ 78 79 101 100 79 80 102 101 80 81 103 102
+ 81 82 104 103 82 83 105 104 83 84 106 105
+ 84 85 107 106 85 86 108 107 86 87 109 108
+ 88 89 111 110 89 90 112 111 90 91 113 112
+ 91 92 114 113 92 93 115 114 93 94 116 115
+ 94 95 117 116 95 96 118 117 96 97 119 118
+ 97 98 120 119 98 99 121 120 99 100 122 121
+ 100 101 123 122 101 102 124 123 102 103 125 124
+ 103 104 126 125 104 105 127 126 105 106 128 127
+ 106 107 129 128 107 108 130 129 108 109 131 130
+ 110 111 133 132 111 112 134 133 112 113 135 134
+ 113 114 136 135 114 115 137 136 115 116 138 137
+ 116 117 139 138 117 118 140 139 118 119 141 140
+ 119 120 142 141 120 121 143 142 121 122 144 143
+ 122 123 145 144 123 124 146 145 124 125 147 146
+ 125 126 148 147 126 127 149 148 127 128 150 149
+ 128 129 151 150 129 130 152 151 130 131 153 152
+ </DataArray>
+ <DataArray type="Int32" Name="offsets" NumberOfComponents="1" format="ascii">
+ 4 8 12 16 20 24 28 32 36 40 44 48
+ 52 56 60 64 68 72 76 80 84 88 92 96
+ 100 104 108 112 116 120 124 128 132 136 140 144
+ 148 152 156 160 164 168 172 176 180 184 188 192
+ 196 200 204 208 212 216 220 224 228 232 236 240
+ 244 248 252 256 260 264 268 272 276 280 284 288
+ 292 296 300 304 308 312 316 320 324 328 332 336
+ 340 344 348 352 356 360 364 368 372 376 380 384
+ 388 392 396 400 404 408 412 416 420 424 428 432
+ 436 440 444 448 452 456 460 464 468 472 476 480
+ 484 488 492 496 500 504
+ </DataArray>
+ <DataArray type="UInt8" Name="types" NumberOfComponents="1" format="ascii">
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9 9 9 9 9 9 9
+ 9 9 9 9 9 9
+ </DataArray>
+ </Cells>
+ </Piece>
+ </UnstructuredGrid>
+</VTKFile>