Commit 51cfe0da authored by Timo Koch's avatar Timo Koch
Browse files

[3p3c] Removed permeability function in volumevariables and moved perm-output...

[3p3c] Removed permeability function in volumevariables and moved perm-output to problem file - so that the model gives no error when choosing anisotropic/tensor-valued permeabilities.

This changed order of output fields, so reference vtus where adjusted too

(Reviewed by bernd)



git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@13572 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent d21d8d30
......@@ -311,7 +311,6 @@ public:
moleFraction[i][j] = writer.allocateManagedBuffer (numDofs);
ScalarField *temperature = writer.allocateManagedBuffer (numDofs);
ScalarField *poro = writer.allocateManagedBuffer(numDofs);
ScalarField *perm = writer.allocateManagedBuffer(numDofs);
VectorField *velocityN = writer.template allocateManagedBuffer<double, dim>(numDofs);
VectorField *velocityW = writer.template allocateManagedBuffer<double, dim>(numDofs);
VectorField *velocityG = writer.template allocateManagedBuffer<double, dim>(numDofs);
......@@ -371,7 +370,6 @@ public:
}
(*poro)[globalIdx] = elemVolVars[scvIdx].porosity();
(*perm)[globalIdx] = elemVolVars[scvIdx].permeability();
(*temperature)[globalIdx] = elemVolVars[scvIdx].temperature();
(*phasePresence)[globalIdx] = staticDat_[globalIdx].phasePresence;
}
......@@ -406,7 +404,6 @@ public:
}
}
writer.attachDofData(*poro, "porosity", isBox);
writer.attachDofData(*perm, "permeability", isBox);
writer.attachDofData(*temperature, "temperature", isBox);
writer.attachDofData(*phasePresence, "phase presence", isBox);
......
......@@ -574,12 +574,6 @@ public:
scvIdx);
Valgrind::CheckDefined(porosity_);
// permeability
permeability_ = problem.spatialParams().intrinsicPermeability(element,
fvGeometry,
scvIdx);
Valgrind::CheckDefined(permeability_);
// energy related quantities not contained in the fluid state
asImp_().updateEnergy_(priVars, problem, element, fvGeometry, scvIdx, isOldSol);
// compute and set the enthalpy
......@@ -685,12 +679,6 @@ public:
Scalar porosity() const
{ return porosity_; }
/*!
* \brief Returns the permeability within the control volume.
*/
Scalar permeability() const
{ return permeability_; }
/*!
* \brief Returns the diffusivity coefficient matrix.
*/
......@@ -740,7 +728,6 @@ protected:
Scalar massFrac_[numPhases][numComponents];
Scalar porosity_; //!< Effective porosity within the control volume
Scalar permeability_; //!< Effective porosity within the control volume
Scalar mobility_[numPhases]; //!< Effective mobility within the control volume
Scalar bulkDensTimesAdsorpCoeff_; //!< the basis for calculating adsorbed NAPL
/* We need a tensor here !! */
......
......@@ -138,6 +138,7 @@ class InfiltrationProblem : public ImplicitPorousMediaProblem<TypeTag>
typedef typename GET_PROP_TYPE(TypeTag, TimeManager) TimeManager;
typedef typename GridView::template Codim<0>::Entity Element;
typedef typename GridView::template Codim<0>::Iterator ElementIterator;
typedef typename GridView::template Codim<dim>::Entity Vertex;
typedef typename GridView::Intersection Intersection;
......@@ -147,6 +148,7 @@ class InfiltrationProblem : public ImplicitPorousMediaProblem<TypeTag>
typedef Dune::FieldVector<Scalar, dimWorld> GlobalPosition;
enum { isBox = GET_PROP_VALUE(TypeTag, ImplicitIsBox) };
enum { dofCodim = isBox ? dim : 0 };
public:
/*!
......@@ -345,6 +347,38 @@ public:
return wgPhaseOnly;
}
/*!
* \brief Append all quantities of interest which can be derived
* from the solution of the current time step to the VTK
* writer. Adjust this in case of anisotropic permeabilities.
*/
void addOutputVtkFields()
{
// get the number of degrees of freedom
unsigned numDofs = this->model().numDofs();
// create the scalar field required for the permeabilities
typedef Dune::BlockVector<Dune::FieldVector<double, 1> > ScalarField;
ScalarField *Kxx = this->resultWriter().allocateManagedBuffer(numDofs);
FVElementGeometry fvGeometry;
ElementIterator eIt = this->gridView().template begin<0>();
ElementIterator eEndIt = this->gridView().template end<0>();
for (; eIt != eEndIt; ++eIt)
{
fvGeometry.update(this->gridView(), *eIt);
for (int scvIdx = 0; scvIdx < fvGeometry.numScv; ++scvIdx)
{
int globalIdx = this->model().dofMapper().map(*eIt, scvIdx, dofCodim);
(*Kxx)[globalIdx] = this->spatialParams().intrinsicPermeability(*eIt, fvGeometry, scvIdx);
}
}
this->resultWriter().attachDofData(*Kxx, "permeability", isBox);
}
private:
// internal method for the initial condition (reused for the
// dirichlet conditions!)
......
......@@ -139,6 +139,7 @@ class ColumnProblem : public ImplicitPorousMediaProblem<TypeTag>
typedef typename GET_PROP_TYPE(TypeTag, TimeManager) TimeManager;
typedef typename GridView::template Codim<0>::Entity Element;
typedef typename GridView::template Codim<0>::Iterator ElementIterator;
typedef typename GridView::template Codim<dim>::Entity Vertex;
typedef typename GridView::Intersection Intersection;
......@@ -148,7 +149,8 @@ class ColumnProblem : public ImplicitPorousMediaProblem<TypeTag>
typedef Dune::FieldVector<Scalar, dimWorld> GlobalPosition;
enum { isBox = GET_PROP_VALUE(TypeTag, ImplicitIsBox) };
enum { dofCodim = isBox ? dim : 0 };
public:
/*!
* \brief The constructor
......@@ -295,6 +297,38 @@ public:
return threePhases;
}
/*!
* \brief Append all quantities of interest which can be derived
* from the solution of the current time step to the VTK
* writer. Adjust this in case of anisotropic permeabilities.
*/
void addOutputVtkFields()
{
// get the number of degrees of freedom
unsigned numDofs = this->model().numDofs();
// create the scalar field required for the permeabilities
typedef Dune::BlockVector<Dune::FieldVector<double, 1> > ScalarField;
ScalarField *Kxx = this->resultWriter().allocateManagedBuffer(numDofs);
FVElementGeometry fvGeometry;
ElementIterator eIt = this->gridView().template begin<0>();
ElementIterator eEndIt = this->gridView().template end<0>();
for (; eIt != eEndIt; ++eIt)
{
fvGeometry.update(this->gridView(), *eIt);
for (int scvIdx = 0; scvIdx < fvGeometry.numScv; ++scvIdx)
{
int globalIdx = this->model().dofMapper().map(*eIt, scvIdx, dofCodim);
(*Kxx)[globalIdx] = this->spatialParams().intrinsicPermeability(*eIt, fvGeometry, scvIdx);
}
}
this->resultWriter().attachDofData(*Kxx, "permeability", isBox);
}
private:
// internal method for the initial condition (reused for the
// dirichlet conditions!)
......
......@@ -151,6 +151,7 @@ class KuevetteProblem : public ImplicitPorousMediaProblem<TypeTag>
typedef typename GET_PROP_TYPE(TypeTag, TimeManager) TimeManager;
typedef typename GridView::template Codim<0>::Entity Element;
typedef typename GridView::template Codim<0>::Iterator ElementIterator;
typedef typename GridView::template Codim<dim>::Entity Vertex;
typedef typename GridView::Intersection Intersection;
......@@ -160,6 +161,7 @@ class KuevetteProblem : public ImplicitPorousMediaProblem<TypeTag>
typedef Dune::FieldVector<Scalar, dimWorld> GlobalPosition;
enum { isBox = GET_PROP_VALUE(TypeTag, ImplicitIsBox) };
enum { dofCodim = isBox ? dim : 0 };
public:
/*!
......@@ -318,6 +320,38 @@ public:
else return wgPhaseOnly;
}
/*!
* \brief Append all quantities of interest which can be derived
* from the solution of the current time step to the VTK
* writer. Adjust this in case of anisotropic permeabilities.
*/
void addOutputVtkFields()
{
// get the number of degrees of freedom
unsigned numDofs = this->model().numDofs();
// create the scalar field required for the permeabilities
typedef Dune::BlockVector<Dune::FieldVector<double, 1> > ScalarField;
ScalarField *Kxx = this->resultWriter().allocateManagedBuffer(numDofs);
FVElementGeometry fvGeometry;
ElementIterator eIt = this->gridView().template begin<0>();
ElementIterator eEndIt = this->gridView().template end<0>();
for (; eIt != eEndIt; ++eIt)
{
fvGeometry.update(this->gridView(), *eIt);
for (int scvIdx = 0; scvIdx < fvGeometry.numScv; ++scvIdx)
{
int globalIdx = this->model().dofMapper().map(*eIt, scvIdx, dofCodim);
(*Kxx)[globalIdx] = this->spatialParams().intrinsicPermeability(*eIt, fvGeometry, scvIdx);
}
}
this->resultWriter().attachDofData(*Kxx, "permeability", isBox);
}
private:
// internal method for the initial condition (reused for the
......
......@@ -254,7 +254,7 @@
2.62566e-05 2.62566e-05 2.62566e-05 2.62566e-05 2.62566e-05 2.62566e-05 2.62566e-05 2.62566e-05 2.62566e-05 2.62566e-05 2.62567e-05 2.62567e-05
2.62567e-05 2.62567e-05 2.62567e-05 2.62567e-05 2.62567e-05 2.62567e-05 2.62567e-05 2.62567e-05 2.62574e-05 2.62574e-05 2.62655e-05 2.62655e-05
2.63671e-05 2.63671e-05 2.76023e-05 2.76023e-05 4.3287e-05 4.3287e-05 0.000374854 0.000374854 0.000635214 0.000635214 4.80611e-05 4.80611e-05
-2.3834e-31 8.8195e-30
6.24415e-31 -1.28444e-30
</DataArray>
<DataArray type="Float32" Name="x^w_Air" NumberOfComponents="1" format="ascii">
1.37483e-05 1.37483e-05 1.37482e-05 1.37482e-05 1.37482e-05 1.37482e-05 1.37482e-05 1.37482e-05 1.37482e-05 1.37482e-05 1.37482e-05 1.37482e-05
......@@ -323,7 +323,7 @@
1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 1 1
1 1 1 1 1 1 1 1 1 1 0.0593977 0.0593977
-2.89386e-28 1.07084e-26
7.58148e-28 -1.55953e-27
</DataArray>
<DataArray type="Float32" Name="x^n_Air" NumberOfComponents="1" format="ascii">
1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10 1e-10
......@@ -392,7 +392,7 @@
0.00983708 0.00983708 0.00983701 0.00983701 0.00983695 0.00983695 0.0098369 0.0098369 0.00983686 0.00983686 0.00983683 0.00983683
0.0098368 0.0098368 0.00983677 0.00983677 0.00983673 0.00983673 0.00983672 0.00983672 0.00983693 0.00983693 0.00983993 0.00983993
0.009878 0.009878 0.0103407 0.0103407 0.0162166 0.0162166 0.140428 0.140428 0.237483 0.237483 0.0179472 0.0179472
-8.88902e-29 3.28929e-27
2.32879e-28 -4.79039e-28
</DataArray>
<DataArray type="Float32" Name="x^g_Air" NumberOfComponents="1" format="ascii">
0.961908 0.961908 0.961908 0.961908 0.961908 0.961908 0.961908 0.961908 0.961908 0.961908 0.961907 0.961907
......@@ -440,29 +440,6 @@
0.46 0.46 0.46 0.46 0.46 0.46 0.46 0.46 0.46 0.46 0.46 0.46
0.46 0.46
</DataArray>
<DataArray type="Float32" Name="permeability" NumberOfComponents="1" format="ascii">
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11
1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11
1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11
1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11
1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11
1.4e-11 1.4e-11
</DataArray>
<DataArray type="Float32" Name="temperature" NumberOfComponents="1" format="ascii">
296.15 296.15 296.15 296.15 296.15 296.15 296.15 296.15 296.15 296.15 296.15 296.15
296.15 296.15 296.15 296.15 296.15 296.15 296.15 296.15 296.15 296.15 296.15 296.15
......@@ -509,6 +486,29 @@
1 1 1 1 1 1 1 1 1 1 6 6
6 6
</DataArray>
<DataArray type="Float32" Name="permeability" NumberOfComponents="1" format="ascii">
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08 1.4e-08
1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11
1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11
1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11
1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11
1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11 1.4e-11
1.4e-11 1.4e-11
</DataArray>
</PointData>
<CellData Scalars="process rank">
<DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
......
......@@ -383,26 +383,6 @@
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4
</DataArray>
<DataArray type="Float32" Name="permeability" NumberOfComponents="1" format="ascii">
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11
</DataArray>
<DataArray type="Float32" Name="temperature" NumberOfComponents="1" format="ascii">
283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15 283.15
......@@ -443,6 +423,26 @@
6 6 6 6 6 6 6 6 6 6 6 6
6
</DataArray>
<DataArray type="Float32" Name="permeability" NumberOfComponents="1" format="ascii">
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11 1e-11
1e-11
</DataArray>
</PointData>
<CellData Scalars="process rank">
<DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment