Commit b0149d6f authored by Timo Koch's avatar Timo Koch
Browse files

Merge branch 'fix/revert-boxdfm-flux-computation' into 'master'

[facet][box] do not use support points for coupling fluxes

See merge request !1688
parents 80afc0b7 758cb98d
......@@ -87,80 +87,80 @@ public:
const auto& insideScv = fvGeometry.scv(scvf.insideScvIdx());
const auto& insideVolVars = elemVolVars[insideScv];
// Evaluate the flux in a tpfa manner using a support point inside the element.
// Therefore, we create a segment between the integration point and another point
// in -(K*normal) direction of the face. This segment must be long enough
// to intersect the opposite face of the element. Then, we compute the
// part of the segment intersecting with the element and compute the
// point inside the element by moving half its length into the interior.
const auto& geometry = element.geometry();
auto distVec = mv(insideVolVars.permeability(), scvf.unitOuterNormal());
distVec /= -1.0*distVec.two_norm();
distVec *= diameter(geometry)*5.0; // make sure segment will be long enough
const auto p1 = scvf.ipGlobal();
auto p2 = scvf.ipGlobal();
p2 += distVec;
static constexpr int dimWorld = GridView::dimensionworld;
using Segment = Dune::MultiLinearGeometry<Scalar, 1, dimWorld>;
using Policy = IntersectionPolicy::SegmentPolicy<typename GridView::ctype, dimWorld>;
using IntersectionAlgorithm = GeometryIntersection<typename Element::Geometry, Segment, Policy>;
typename IntersectionAlgorithm::Intersection intersection;
Segment segment(Dune::GeometryTypes::line, std::vector<GlobalPosition>({p1, p2}));
if (!IntersectionAlgorithm::intersection(geometry, segment, intersection))
DUNE_THROW(Dune::InvalidStateException, "Could not compute support point for flux computation");
// use center of intersection as integration point
auto supportPoint = intersection[0];
supportPoint += intersection[1];
supportPoint /= 2.0;
std::vector< Dune::FieldVector<Scalar, 1> > supportPointShapeValues;
fvGeometry.feLocalBasis().evaluateFunction(geometry.local(supportPoint), supportPointShapeValues);
Scalar rho = 0.0;
Scalar supportPressure = 0.0;
for (const auto& scv : scvs(fvGeometry))
{
const auto& volVars = elemVolVars[scv];
rho += volVars.density(phaseIdx)*shapeValues[scv.indexInElement()][0];
supportPressure += volVars.pressure(phaseIdx)*supportPointShapeValues[scv.indexInElement()][0];
}
// the transmissibility on the matrix side
const auto d = scvf.ipGlobal() - supportPoint;
const auto tm = 1.0/d.two_norm2()*vtmv(d, insideVolVars.permeability(), scvf.unitOuterNormal());
// compute flux depending on the user's choice of boundary types
// evaluate user-defined interior boundary types
const auto bcTypes = problem.interiorBoundaryTypes(element, scvf);
const auto& facetVolVars = problem.couplingManager().getLowDimVolVars(element, scvf);
// compute tpfa flux such that flux continuity holds with flux in fracture
Scalar flux;
static const bool enableGravity = getParamFromGroup<bool>(problem.paramGroup(), "Problem.EnableGravity");
// on interior Neumann boundaries, evaluate the flux using the facet permeability
if (bcTypes.hasOnlyNeumann())
{
// On surface grids, use sqrt of aperture as distance measur
using std::sqrt;
const auto df = dim == dimWorld ? 0.5*facetVolVars.extrusionFactor() : 0.5*sqrt(facetVolVars.extrusionFactor());
const auto tf = 1.0/df*vtmv(scvf.unitOuterNormal(), facetVolVars.permeability(), scvf.unitOuterNormal());
// interpolate pressure/density to scvf integration point
Scalar p = 0.0;
Scalar rho = 0.0;
for (const auto& scv : scvs(fvGeometry))
{
const auto& volVars = elemVolVars[scv];
p += volVars.pressure(phaseIdx)*shapeValues[scv.indexInElement()][0];
rho += volVars.density(phaseIdx)*shapeValues[scv.indexInElement()][0];
}
// compute tpfa flux from integration point to facet centerline
const auto& facetVolVars = problem.couplingManager().getLowDimVolVars(element, scvf);
flux = tm*tf/(tm+tf)*(supportPressure - facetVolVars.pressure(phaseIdx))*scvf.area()*insideVolVars.extrusionFactor();
using std::sqrt;
// If this is a surface grid, use the square root of the facet extrusion factor
// as an approximate average distance from scvf ip to facet center
using std::sqrt;
const auto a = facetVolVars.extrusionFactor();
auto gradP = scvf.unitOuterNormal();
gradP *= dim == dimWorld ? 0.5*a : 0.5*sqrt(a);
gradP /= gradP.two_norm2();
gradP *= (facetVolVars.pressure(phaseIdx) - p);
if (enableGravity)
gradP.axpy(-rho, problem.spatialParams().gravity(scvf.center()));
// apply facet permeability and return the flux
return -1.0*scvf.area()
*insideVolVars.extrusionFactor()
*vtmv(scvf.unitOuterNormal(), facetVolVars.permeability(), gradP);
}
// on interior Dirichlet boundaries use the facet pressure and evaluate flux
else if (bcTypes.hasOnlyDirichlet())
flux = tm*(supportPressure - facetVolVars.pressure(phaseIdx))*scvf.area()*insideVolVars.extrusionFactor();
{
// create vector with nodal pressures
std::vector<Scalar> pressures(element.subEntities(dim));
for (const auto& scv : scvs(fvGeometry))
pressures[scv.localDofIndex()] = elemVolVars[scv].pressure(phaseIdx);
// substitute with facet pressures for those scvs touching this facet
for (const auto& scvfJ : scvfs(fvGeometry))
if (scvfJ.interiorBoundary() && scvfJ.facetIndexInElement() == scvf.facetIndexInElement())
pressures[ fvGeometry.scv(scvfJ.insideScvIdx()).localDofIndex() ]
= problem.couplingManager().getLowDimVolVars(element, scvfJ).pressure(phaseIdx);
// evaluate gradP - rho*g at integration point
Scalar rho(0.0);
Dune::FieldVector<Scalar, dimWorld> gradP(0.0);
for (const auto& scv : scvs(fvGeometry))
{
rho += elemVolVars[scv].density(phaseIdx)*shapeValues[scv.indexInElement()][0];
gradP.axpy(pressures[scv.localDofIndex()], fluxVarCache.gradN(scv.indexInElement()));
}
if (enableGravity)
gradP.axpy(-rho, problem.spatialParams().gravity(scvf.center()));
// apply matrix permeability and return the flux
return -1.0*scvf.area()
*insideVolVars.extrusionFactor()
*vtmv(scvf.unitOuterNormal(), insideVolVars.permeability(), gradP);
}
// mixed boundary types are not supported
else
DUNE_THROW(Dune::NotImplemented, "Mixed boundary types are not supported");
static const bool enableGravity = getParamFromGroup<bool>(problem.paramGroup(), "Problem.EnableGravity");
if (enableGravity)
flux -= rho*scvf.area()*insideVolVars.extrusionFactor()
*vtmv(scvf.unitOuterNormal(), insideVolVars.permeability(), problem.spatialParams().gravity(scvf.center()));
return flux;
}
// compute transmissibilities ti for analytical jacobians
......
......@@ -4,50 +4,50 @@
<Piece NumberOfCells="1927" NumberOfPoints="526">
<PointData Scalars="p">
<DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
1.88262 2 1.88288 2 2 1.84732 2 1.79699 2 1.88469 2 1.73117
1.70166 1.7318 1.80681 1.70898 1.75841 1.88424 2 2 2 2 1.80438 1.69959
1.72962 1.83826 2 2 2 1.88683 1.81704 1.72358 1.66915 1.75564 2 2
1.61471 1.65243 1.78008 1.62324 1.62079 1.56533 1.55935 1.55709 1.49838 1.50015 1.49818 1.593
1.50133 1.54324 1.39982 1.44788 1.54122 1.61656 1.65086 1.50393 1.44824 1.53603 1.50957 2
1.83213 2 2 1.64736 1.61328 1.88556 1.81227 2 2 1.75419 1.66803 2
2 1.78492 1.57067 2 2 1.57396 1.51414 1.52757 1.78362 1.55025 1.53704 1.43829
1.50902 1.46352 1.50599 1.55576 1.51955 1.51269 1.51163 1.48769 1.51323 1.51169 1.50774 1.50823
2 2 1.83296 2 1.88431 2 1.73308 1.70045 1.80442 1.79652 2 1.77358
2 1.65104 2 2 1.8832 2 1.88249 1.84749 2 1.73154 1.70205 1.73183
2 2 1.88549 2 1.80227 1.70154 1.55983 1.55143 1.44142 1.59358 1.50147 1.55767
1.49932 1.62363 1.6239 1.56648 1.50053 1.50058 1.43521 1.5935 1.50118 1.55609 2 1.76706
1.62081 2 2 1.579 2 1.83142 2 1.65357 1.81427 2 1.88593 1.76021
1.55439 1.53444 1.51241 1.51437 1.52018 1.50483 1.50972 1.55021 1.43918 1.62471 1.50512 1.43495
1.3795 1.37981 1.36782 1.41135 1.29877 1.26831 1.2687 1.30166 1.19793 1.2494 1.14189 1.28718
1.19575 1.35657 1.46474 1.38221 1.36587 1.16907 1.1777 1.32772 1.38238 1.1532 1.11834 1.11806
1.11508 1 1 1 1 1 1 1.11511 1 1 1 1
1.24519 1.11168 1 1 1 1.22245 1 1 1 1.35196 1.39448 1.3734
1.26663 1.19812 1.35053 1.46821 1.37301 1.23791 1.50231 1.48431 1.49714 1.42745 1.48753 1.50265
1.49876 1.50549 1.42397 1.48407 1.1622 1 1 1 1.11084 1.17534 1 1
1 1.21735 1 1.21842 1 1.44121 1.38054 1.27723 1.29183 1.35028 1.24402 1.38085
1.26913 1.19799 1.29934 1.26937 1.29243 1.16222 1.11577 1.19355 1 1 1 1
1.21828 1 1.11778 1.15374 1.11874 1 1 1 1 1.19695 1.11605 1
1 1 1.41312 1.27875 1.49623 1.46718 1.42273 1.5011 1.37941 1.44226 1.37201 1.27256
1.24444 1.3579 1 1 1 1 1.16157 1 1 1.18859 1.1158 1
1 2 1.83726 2 1.75787 1.72287 2 2 1.88403 1.79393 2 1.70814
1.8056 1.7745 1.64885 2 1.57224 2 2 1.62198 1.51638 1.61735 1.54054 1.55658
1.4462 1.62274 1.49823 1.55816 1.51164 1.51093 1.49554 1.53548 1.54022 1.44982 1.50187 2
1.88158 2 2 2 1.73212 2 1.8817 2 1.84648 1.8837 2 1.73153
1.70232 1.80288 1.70697 1.71877 2 1.81419 2 1.65351 1.83812 2 1.88487 2
1.72715 1.7602 1.5139 1.62248 1.49922 1.56619 1.49907 1.4994 1.62591 1.46433 1.51207 1.55785
1.44548 2 1.78049 2 2 1.76755 1.79653 1.65156 2 1.71993 2 1.83831
1.72747 2 2 1.88515 1.70846 1.55509 1.4874 1.54011 1.54089 1.44801 1.5619 1.44858
1.55879 1.49836 2 2 1.70175 1.83236 1.88571 2 2 2 2 1.88283
1.80388 1.73025 1.8823 2 2 1.8469 1.80218 1.73121 1.701 1.55053 1.59357 1.50157
1.43882 1.62077 1.4984 1.62299 1.56514 1.40024 1.50019 1.4352 1.36308 1.4099 1.36565 1.29994
1.16704 1.24792 1.42314 1.50371 1.21544 1.50008 1.43262 1.36353 1.4708 1.49056 1.37081 1
1 1 1.11442 1 1 1 1 1 1.19465 1 1.37856 1.37722
1.39239 1.19823 1.26727 1.26647 1.43381 1.40877 1.2973 1.29888 1.1976 1.34593 1.16836 1.24313
1.37324 1.14181 1.1181 1 1 1 1 1.15253 1.118 1.11453 1 1
1 1.1986 1 1 1.11011 1 1.39568 1.36852 1.22395 1.35692 1.24884 1.19116
1.36878 1.41205 1.30227 1.16491 1 1 1 1.11413 1.14134 1 1 1.27146
1.28977 1.38004 1.26851 1.19891 1.26933 1.38012 1.19706 1.29914 1 1 1.16329 1.11513
1 1 1.11886 1 1 1 1.11779 1.15345 1 1
1.87805 2 1.87832 2 2 1.84227 2 1.78511 2 1.87808 2 1.72406
1.69542 1.72491 1.80032 1.69659 1.73365 1.87729 2 2 2 2 1.79791 1.68755
1.70429 1.82501 2 2 2 1.87603 1.79257 1.70596 1.64114 1.73704 2 2
1.52532 1.63943 1.74367 1.61832 1.61566 1.56293 1.55502 1.55112 1.49855 1.50085 1.49856 1.58733
1.50216 1.53644 1.40995 1.455 1.53277 1.60146 1.63534 1.50341 1.45572 1.51054 1.50104 2
1.81908 2 2 1.63237 1.52533 1.87431 1.78542 2 2 1.73524 1.6401 2
2 1.75652 1.51384 2 2 1.51419 1.51219 1.50562 1.75748 1.54467 1.51067 1.44296
1.51013 1.47005 1.50063 1.54824 1.51209 1.50496 1.50807 1.48896 1.50816 1.5049 1.50403 1.50065
2 2 1.81974 2 1.87701 2 1.70777 1.68777 1.79713 1.7834 2 1.73577
2 1.63698 2 2 1.87843 2 1.87724 1.84183 2 1.72348 1.69479 1.72378
2 2 1.87851 2 1.79344 1.68908 1.55028 1.54525 1.45005 1.58698 1.50109 1.55244
1.49816 1.61763 1.6179 1.56312 1.49999 1.49993 1.43773 1.58677 1.50017 1.54719 2 1.73576
1.52813 2 2 1.51468 2 1.81565 2 1.63645 1.78496 2 1.8729 1.74061
1.54649 1.51061 1.50844 1.50519 1.51284 1.49669 1.50115 1.54363 1.44628 1.61732 1.50412 1.439
1.38729 1.38695 1.38118 1.42079 1.30623 1.276 1.27665 1.31462 1.20639 1.26779 1.1507 1.2986
1.20466 1.37496 1.48949 1.3976 1.39437 1.18397 1.20387 1.36002 1.47646 1.159 1.12376 1.12347
1.12199 1 1 1 1 1 1 1.12288 1 1 1 1
1.2653 1.12194 1 1 1 1.25194 1 1 1 1.37451 1.41557 1.3783
1.28058 1.22712 1.36461 1.49029 1.47527 1.25512 1.49668 1.49007 1.49304 1.48712 1.48999 1.49661
1.49293 1.49728 1.48683 1.49572 1.17902 1 1 1 1.12183 1.20401 1 1
1 1.24825 1 1.24284 1 1.44736 1.38644 1.2971 1.30599 1.36203 1.27479 1.38636
1.27604 1.21009 1.30578 1.27645 1.3058 1.17758 1.12322 1.20003 1 1 1 1
1.24669 1 1.12271 1.15928 1.12424 1 1 1 1 1.20579 1.12345 1
1 1 1.4455 1.30969 1.49286 1.48973 1.48654 1.48986 1.47528 1.44542 1.37408 1.28884
1.26554 1.37473 1 1 1 1 1.17599 1 1 1.22122 1.13039 1
1 2 1.82395 2 1.73413 1.70545 2 2 1.87739 1.78139 2 1.69602
1.79893 1.73678 1.63569 2 1.51435 2 2 1.52829 1.51268 1.60281 1.53372 1.55108
1.45153 1.61866 1.49906 1.55099 1.50829 1.5049 1.49703 1.51042 1.53443 1.45976 1.50033 2
1.87663 2 2 2 1.72594 2 1.87654 2 1.8412 1.87607 2 1.7248
1.69669 1.79393 1.69245 1.68528 2 1.78453 2 1.63724 1.82333 2 1.87105 2
1.70969 1.74108 1.51339 1.61786 1.50055 1.56466 1.50032 1.50089 1.61988 1.47302 1.5163 1.55389
1.45032 2 1.74114 2 2 1.73515 1.78274 1.636 2 1.68564 2 1.82302
1.7093 2 2 1.87797 1.69359 1.54725 1.48865 1.53215 1.53599 1.45656 1.55651 1.45628
1.55464 1.49927 2 2 1.68985 1.81701 1.87884 2 2 2 2 1.8776
1.79437 1.72178 1.87696 2 2 1.84097 1.79338 1.7233 1.69382 1.54625 1.58845 1.50327
1.44481 1.61507 1.49897 1.61791 1.56279 1.41212 1.50148 1.44008 1.38997 1.41869 1.37943 1.31284
1.18245 1.26836 1.44447 1.49665 1.24442 1.49305 1.48712 1.39152 1.49049 1.49009 1.47484 1
1 1 1.12179 1 1 1 1 1 1.23094 1 1.38676 1.38476
1.40591 1.208 1.27527 1.27395 1.43824 1.41723 1.3048 1.31021 1.20654 1.35507 1.18778 1.26718
1.37879 1.1524 1.12422 1 1 1 1 1.15874 1.124 1.12269 1 1
1 1.22937 1 1 1.12166 1 1.41562 1.40009 1.25892 1.3861 1.26849 1.21961
1.38382 1.4223 1.3158 1.17897 1 1 1 1.12129 1.15109 1 1 1.29057
1.30431 1.3889 1.27734 1.20853 1.27843 1.38853 1.2079 1.30769 1 1 1.18121 1.12308
1 1 1.12517 1 1 1 1.12381 1.16012 1 1
</DataArray>
</PointData>
<CellData Scalars="process rank">
......
......@@ -4,7 +4,7 @@
<Piece NumberOfLines="4" NumberOfPoints="5">
<PointData Scalars="p">
<DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
1.51027 1.50907 1.50652 1.50395 1.50273
1.50561 1.50403 1.50065 1.49728 1.49572
</DataArray>
</PointData>
<CellData Scalars="process rank">
......
......@@ -4,11 +4,11 @@
<Piece NumberOfCells="88" NumberOfPoints="57">
<PointData Scalars="p">
<DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
1.51056 1.51218 1.51578 1.50903 1.50266 1.49761 1.50094 1.50401 1.52013 1.51485 1.51379 1.49425
1.49944 1.49967 1.50658 1.50368 1.5097 1.5065 1.5121 1.51563 1.50101 1.49769 1.52046 1.51484
1.51372 1.49434 1.49952 1.49977 1.50963 1.50371 1.50647 1.51994 1.51541 1.51464 1.51363 1.49384
1.4991 1.49908 1.49769 1.51201 1.50087 1.50962 1.50353 1.50651 1.51217 1.51571 1.49724 1.50068
1.52053 1.51493 1.51384 1.49355 1.49904 1.49889 1.50973 1.50656 1.5035
1.50562 1.50816 1.51414 1.50403 1.49572 1.48687 1.49292 1.49728 1.52518 1.51223 1.51066 1.4754
1.48997 1.4903 1.50063 1.49661 1.5049 1.50065 1.50829 1.51431 1.49305 1.48716 1.52811 1.51273
1.51043 1.47498 1.49007 1.49049 1.50489 1.49666 1.50033 1.52517 1.51379 1.51211 1.51054 1.47658
1.48951 1.49004 1.48716 1.50807 1.49304 1.50496 1.49668 1.50104 1.50844 1.51464 1.48659 1.49286
1.52795 1.51287 1.51061 1.47541 1.48974 1.48983 1.50519 1.50114 1.4967
</DataArray>
</PointData>
<CellData Scalars="process rank">
......
......@@ -4,15 +4,15 @@
<Piece NumberOfLines="100" NumberOfPoints="101">
<PointData Scalars="p">
<DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii">
106062 106054 106037 106014 105986 105956 105924 105890 105855 105819 105782 105744
105706 105668 105622 105576 105531 105485 105440 105394 105349 105304 105258 105213
105168 105123 105079 105034 104990 104974 104959 104945 104931 104917 104905 104895
104887 104758 104758 104758 104758 104757 104757 104757 104758 104758 104759 104760
104761 104762 104781 104801 104820 104838 104857 104876 104895 104914 104933 104952
104971 104981 104973 104964 104957 104949 104943 104937 104934 105704 105704 105704
105703 105703 105702 105700 105698 105696 105694 105692 105689 105686 105682 105679
105675 105671 105670 105674 105678 104423 104487 104552 104610 104659 104700 104734
104773 104788 104806 104828 104850
106060 106053 106035 106012 105984 105954 105922 105888 105853 105817 105779 105742
105704 105665 105620 105574 105528 105483 105437 105392 105346 105301 105256 105211
105166 105121 105076 105031 104987 104972 104957 104942 104928 104915 104903 104892
104885 104755 104755 104755 104754 104754 104754 104754 104754 104755 104756 104757
104758 104759 104778 104797 104816 104835 104854 104873 104892 104911 104930 104949
104968 104978 104970 104962 104954 104946 104940 104935 104931 105702 105702 105702
105701 105700 105699 105698 105696 105694 105692 105689 105686 105683 105680 105676
105673 105669 105668 105672 105676 104419 104483 104548 104606 104655 104697 104731
104770 104784 104802 104825 104847
</DataArray>
</PointData>
<CellData Scalars="process rank">
......
This source diff could not be displayed because it is too large. You can view the blob instead.
......@@ -4,26 +4,26 @@
<Piece NumberOfLines="100" NumberOfPoints="101">
<PointData Scalars="x^tracer_0">
<DataArray type="Float32" Name="x^tracer_0" NumberOfComponents="1" format="ascii">
6.43776e-06 1.5335e-05 0.000974219 0.00570012 0.0133654 0.0193448 0.021786 0.0217631 0.0210783 0.0202995 0.0195915 0.0189528
0.0183562 0.0148532 0.0143645 0.0138346 0.013268 0.0126687 0.0120403 0.0113856 0.010708 0.0100099 0.00929513 0.00856952
0.00783964 0.00711267 0.00639621 0.0056983 0.00445852 0.00308973 0.00202863 0.0012533 0.00072174 0.000381818 0.00018018 7.20845e-05
2.48713e-05 3.30022e-31 -1.4585e-30 2.22679e-27 -2.14416e-24 1.85454e-21 -1.45043e-18 1.14701e-15 1.94754e-13 1.30445e-11 5.08534e-10 1.34303e-08
2.62169e-07 4.06127e-06 1.51536e-05 3.0443e-05 5.91849e-05 0.000111213 0.000201747 0.000352859 0.000594306 0.000962901 0.00149957 0.00224369
0.00322572 0.00245257 0.00123266 0.000561476 0.000228518 8.14487e-05 2.44043e-05 5.65291e-06 1.11135e-06 -0 0 -1.97583e-43
1.12404e-40 -6.40839e-38 3.61333e-35 -2.02264e-32 1.11411e-29 -6.22004e-27 3.41548e-24 -1.86144e-21 1.02066e-18 -5.44547e-16 2.91471e-13 -1.48958e-10
7.66076e-08 -3.47458e-05 -2.92147e-05 2.45203e-07 1.65244e-07 8.83885e-10 5.23159e-09 2.74789e-08 1.09287e-07 3.48211e-07 9.2039e-07 2.07136e-06
-2.70303e-10 1.44834e-14 1.38753e-17 6.73615e-17 2.05065e-16
0 0 0.000928025 0.00553288 0.0131449 0.0192773 0.0218448 0.0219353 0.0213199 0.0205958 0.0199429 0.0193602
0.0188177 0.0152471 0.0148085 0.0143297 0.0138137 0.0132628 0.0126796 0.012066 0.0114236 0.0107544 0.0100628 0.00935362
0.00863288 0.00790713 0.00718356 0.00646993 0.0051725 0.00365401 0.00244534 0.00154104 0.000906851 0.00049253 0.000240196 0.000103001
4.40332e-05 1.3702e-30 8.34454e-30 -7.75442e-27 7.8612e-24 -7.08269e-21 5.67853e-18 5.7979e-15 6.84508e-13 3.87419e-11 1.35923e-09 3.32544e-08
6.09866e-07 8.88015e-06 3.09885e-05 5.85639e-05 0.000107195 0.000189876 0.000325195 0.000538009 0.000859113 0.00132325 0.00196506 0.00281296
0.00388305 0.00290195 0.00148333 0.000687085 0.000284913 0.000104067 3.21961e-05 8.11928e-06 2.0242e-06 -0 0 -2.01787e-43
1.15524e-40 -6.64289e-38 3.76181e-35 -2.12179e-32 1.18104e-29 -6.54402e-27 3.61308e-24 -1.97203e-21 1.06546e-18 -5.68899e-16 2.97935e-13 -1.5088e-10
7.34868e-08 -3.36357e-05 -1.87361e-05 2.94289e-07 2.87356e-07 4.57796e-09 2.0286e-08 9.2968e-08 3.31808e-07 9.58402e-07 2.32323e-06 4.84344e-06
-4.57513e-10 2.58822e-14 5.64127e-17 2.42799e-16 1.12718e-15
</DataArray>
<DataArray type="Float32" Name="X^tracer_0" NumberOfComponents="1" format="ascii">
1.07296e-07 2.55583e-07 1.6237e-05 9.5002e-05 0.000222757 0.000322413 0.0003631 0.000362719 0.000351306 0.000338326 0.000326525 0.00031588
0.000305937 0.000247553 0.000239408 0.000230577 0.000221133 0.000211146 0.000200672 0.00018976 0.000178467 0.000166832 0.000154919 0.000142825
0.000130661 0.000118545 0.000106604 9.49716e-05 7.43087e-05 5.14955e-05 3.38105e-05 2.08884e-05 1.2029e-05 6.36364e-06 3.00301e-06 1.20141e-06
4.14522e-07 5.50036e-33 -2.43083e-32 3.71131e-29 -3.5736e-26 3.09089e-23 -2.41738e-20 1.91169e-17 3.2459e-15 2.17408e-13 8.47556e-12 2.23839e-10
4.36948e-09 6.76878e-08 2.5256e-07 5.07384e-07 9.86415e-07 1.85355e-06 3.36245e-06 5.88099e-06 9.9051e-06 1.60484e-05 2.49928e-05 3.73948e-05
5.37621e-05 4.08761e-05 2.05444e-05 9.35794e-06 3.80863e-06 1.35748e-06 4.06738e-07 9.42152e-08 1.85225e-08 -0 0 -2.8026e-45
1.87354e-42 -1.06807e-39 6.02221e-37 -3.37107e-34 1.85685e-31 -1.03667e-28 5.69247e-26 -3.1024e-23 1.7011e-20 -9.07579e-18 4.85785e-15 -2.48264e-12
1.27679e-09 -5.79097e-07 -4.86911e-07 4.08672e-09 2.75407e-09 1.47314e-11 8.71931e-11 4.57982e-10 1.82145e-09 5.80351e-09 1.53398e-08 3.45227e-08
-4.50505e-12 2.4139e-16 2.31256e-19 1.12269e-18 3.41776e-18
0 0 1.54671e-05 9.22148e-05 0.000219081 0.000321289 0.000364081 0.000365589 0.000355332 0.000343263 0.000332382 0.00032267
0.000313629 0.000254119 0.000246809 0.000238829 0.000230229 0.000221046 0.000211326 0.0002011 0.000190393 0.00017924 0.000167713 0.000155894
0.000143881 0.000131785 0.000119726 0.000107832 8.62083e-05 6.09002e-05 4.07557e-05 2.56839e-05 1.51142e-05 8.20883e-06 4.00327e-06 1.71668e-06
7.33886e-07 2.28367e-32 1.39076e-31 -1.2924e-28 1.3102e-25 -1.18045e-22 9.46421e-20 9.66317e-17 1.14085e-14 6.45699e-13 2.26538e-11 5.5424e-10
1.01644e-08 1.48003e-07 5.16474e-07 9.76066e-07 1.78658e-06 3.16459e-06 5.41991e-06 8.96682e-06 1.43186e-05 2.20541e-05 3.2751e-05 4.68826e-05
6.47175e-05 4.83658e-05 2.47222e-05 1.14514e-05 4.74855e-06 1.73446e-06 5.36602e-07 1.35321e-07 3.37366e-08 -0 0 -2.8026e-45
1.92538e-42 -1.10715e-39 6.26969e-37 -3.53631e-34 1.96841e-31 -1.09067e-28 6.0218e-26 -3.28672e-23 1.77576e-20 -9.48165e-18 4.96559e-15 -2.51466e-12
1.22478e-09 -5.60595e-07 -3.12268e-07 4.90482e-09 4.78926e-09 7.62993e-11 3.38101e-10 1.54947e-09 5.53013e-09 1.59734e-08 3.87206e-08 8.0724e-08
-7.62522e-12 4.31371e-16 9.40212e-19 4.04665e-18 1.87863e-17
</DataArray>
<DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii">
1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000 1000
......
Markdown is supported
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