diff --git a/lecture/efm/1p2c_2p_2p2c/CMakeLists.txt b/lecture/efm/1p2c_2p_2p2c/CMakeLists.txt index 246e7ef6d1df0aed5e11f3d2c387a67f5cbab863..4d09d31e9ab84d799025da6d80a1941a10ebb36b 100644 --- a/lecture/efm/1p2c_2p_2p2c/CMakeLists.txt +++ b/lecture/efm/1p2c_2p_2p2c/CMakeLists.txt @@ -14,7 +14,7 @@ dumux_add_test(NAME lens2p2cexercise3 CMD_ARGS --script fuzzy --files ${CMAKE_SOURCE_DIR}/lecture/references/lens-2p2c-exercise3-reference.vtu ${CMAKE_CURRENT_BINARY_DIR}/lens-2p2c-00011.vtu - --command "${CMAKE_CURRENT_BINARY_DIR}/lens2p2cexercise3 -ParameterFile ${CMAKE_CURRENT_SOURCE_DIR}/exercise3.input -Problem.Name lens-2p2c") + --command "${CMAKE_CURRENT_BINARY_DIR}/lens2p2cexercise3 -ParameterFile ${CMAKE_CURRENT_SOURCE_DIR}/exercise3.input -Problem.Name lens-2p2c -SpatialParams.Temperature 313.15") dumux_add_test(NAME lens1p2cexercise3 SOURCES lens1p2cexercise3.cc diff --git a/lecture/efm/1p2c_2p_2p2c/exercise3.input b/lecture/efm/1p2c_2p_2p2c/exercise3.input index 3cf360de721ca6259b69c1cd82c254c708183f81..85f4724779202971570fd194dade97f8f3c55617 100755 --- a/lecture/efm/1p2c_2p_2p2c/exercise3.input +++ b/lecture/efm/1p2c_2p_2p2c/exercise3.input @@ -24,6 +24,9 @@ BrooksCoreyPcEntry = 200 # entry pressure for the Brooks-Corey capillary pre Swr = 0.05 # residual saturation of the wetting phase in the coarse soil [-] Snr = 0.1 # residual saturation of the nonwetting phase in the coarse soil [-] +[SpatialParams] +Temperature = 283.15 # 10 degree Celcius for 2p problem and 1p2c problem, 40 degree Celcius for 2p2c problem + [Boundary] LowerPressure = 2.0e5 # Dirichlet pressure value for the boundary condition at the lower boundary [Pa] UpperPressure = 4.0e5 # Dirichlet pressure value for the boundary condition at the upper boundary [Pa] diff --git a/lecture/efm/1p2c_2p_2p2c/lens1p2cexercise3.cc b/lecture/efm/1p2c_2p_2p2c/lens1p2cexercise3.cc index e8d3e8445d3312637ee73e9c2f2dd70106a772f5..0170fc0b725a6656c56b29dafd06f7d51db188fc 100755 --- a/lecture/efm/1p2c_2p_2p2c/lens1p2cexercise3.cc +++ b/lecture/efm/1p2c_2p_2p2c/lens1p2cexercise3.cc @@ -114,7 +114,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/efm/1p2c_2p_2p2c/lens1p2cproblem.hh b/lecture/efm/1p2c_2p_2p2c/lens1p2cproblem.hh index aeb2e3b7760cd5153705883a9d38a19b3db7bf61..3bf7f0fefcc8c4130202f657254a5f4ef4cb3be9 100644 --- a/lecture/efm/1p2c_2p_2p2c/lens1p2cproblem.hh +++ b/lecture/efm/1p2c_2p_2p2c/lens1p2cproblem.hh @@ -91,21 +91,6 @@ public: upperPressure_ = getParam<Scalar>("Boundary.UpperPressure"); } - /*! - * \name Problem parameters - */ - // \{ - - /*! - * \brief Returns the temperature within the domain. - * This problem assumes a temperature of 10 degrees Celsius. - */ - Scalar temperature() const - { - return 273.15 + 10; // -> 10 degrees Celsius - } - - // \} /*! * \name Boundary conditions diff --git a/lecture/efm/1p2c_2p_2p2c/lens1p2cspatialparams.hh b/lecture/efm/1p2c_2p_2p2c/lens1p2cspatialparams.hh index 1d9b051b9dd386f3567b4815ab7da655e44105e8..4c6341c25db0fb9113dd0b9475ae7a6e3bce4a9f 100644 --- a/lecture/efm/1p2c_2p_2p2c/lens1p2cspatialparams.hh +++ b/lecture/efm/1p2c_2p_2p2c/lens1p2cspatialparams.hh @@ -25,7 +25,7 @@ #include <dumux/common/properties.hh> -#include <dumux/material/spatialparams/fv1p.hh> +#include <dumux/porousmediumflow/fvspatialparams1p.hh> namespace Dumux { @@ -56,10 +56,10 @@ public: * \brief Class for defining spatial parameters */ template<class FVGridGeometry, class Scalar> -class Lens1p2cSpatialParams : public FVSpatialParamsOneP<FVGridGeometry, Scalar, Lens1p2cSpatialParams<FVGridGeometry, Scalar> > +class Lens1p2cSpatialParams : public FVPorousMediumFlowSpatialParamsOneP<FVGridGeometry, Scalar, Lens1p2cSpatialParams<FVGridGeometry, Scalar> > { using ThisType = Lens1p2cSpatialParams<FVGridGeometry, Scalar>; - using ParentType = FVSpatialParamsOneP<FVGridGeometry, Scalar, ThisType>; + using ParentType = FVPorousMediumFlowSpatialParamsOneP<FVGridGeometry, Scalar, ThisType>; using GridView = typename FVGridGeometry::GridView; using Element = typename GridView::template Codim<0>::Entity; using FVElementGeometry = typename FVGridGeometry::LocalView; @@ -90,6 +90,7 @@ public: lensK_ = getParam<double>("SpatialParams.Fine.Permeability"); outerK_ = getParam<double>("SpatialParams.Coarse.Permeability"); + temperature_ = getParam<double>("SpatialParams.Temperature"); } /*! @@ -141,6 +142,15 @@ public: lensUpperRight_ = lensUpperRight; } + /*! + * \brief Returns the temperature at the domain at the given position + * \param globalPos The position in global coordinates where the temperature should be specified + */ + Scalar temperatureAtPos(const GlobalPosition& globalPos) const + { + return temperature_; + } + private: bool isInLens_(const GlobalPosition &pos) const { @@ -160,6 +170,7 @@ private: Scalar outerK_; Scalar lensPorosity_; Scalar outerPorosity_; + Scalar temperature_; }; } // end namespace Dumux diff --git a/lecture/efm/1p2c_2p_2p2c/lens2p2cexercise3.cc b/lecture/efm/1p2c_2p_2p2c/lens2p2cexercise3.cc index a2cafe1699077a3fc6d6576ada85b4a40305a470..5f8d0b2b49af247a973a31b2fd2a60d203709bdf 100644 --- a/lecture/efm/1p2c_2p_2p2c/lens2p2cexercise3.cc +++ b/lecture/efm/1p2c_2p_2p2c/lens2p2cexercise3.cc @@ -130,7 +130,7 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); + // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/efm/1p2c_2p_2p2c/lens2p2cproblem.hh b/lecture/efm/1p2c_2p_2p2c/lens2p2cproblem.hh index c4718204a4bf2aff372e7c270429ae8ac9b50875..00275d2807370f044dfc1becb11176ba892b04dc 100644 --- a/lecture/efm/1p2c_2p_2p2c/lens2p2cproblem.hh +++ b/lecture/efm/1p2c_2p_2p2c/lens2p2cproblem.hh @@ -86,12 +86,11 @@ public: LensTwoPTwoCProblem(std::shared_ptr<const FVGridGeometry> fvGridGeometry ) : ParentType(fvGridGeometry) { - temperature_ = 273.15 + 40; // [K] eps_ = 1e-6; // initialize the tables of the fluid system - FluidSystem::init(/*Tmin=*/temperature_ - 1.0, - /*Tmax=*/temperature_ + 1.0, + FluidSystem::init(/*Tmin=*/ 273.15 + getParam<Scalar>("SpatialParams.Temperature") - 1.0, + /*Tmax=*/ 273.15 + getParam<Scalar>("SpatialParams.Temperature") + 1.0, /*nT=*/3, /*pmin=*/1e5, /*pmax=*/3e7, @@ -101,20 +100,6 @@ public: upperPressure_ = getParam<Scalar>("Boundary.UpperPressure"); } - /*! - * \name Problem parameters - */ - // \{ - - /*! - * \brief Returns the temperature within the domain. - */ - Scalar temperature() const - { - return temperature_; - } - - // \} /*! * \name Boundary conditions @@ -220,7 +205,6 @@ public: private: Scalar eps_; Scalar episodeLength_; // [s] - Scalar temperature_; // [K] Scalar upperPressure_; // [Pa] Scalar lowerPressure_; // [Pa] diff --git a/lecture/efm/1p2c_2p_2p2c/lens2pexercise3.cc b/lecture/efm/1p2c_2p_2p2c/lens2pexercise3.cc index 2681a024e72b97ada8971a1723198a515ba999c0..24f073484d8515fde2e869a4f4fba61191fef456 100755 --- a/lecture/efm/1p2c_2p_2p2c/lens2pexercise3.cc +++ b/lecture/efm/1p2c_2p_2p2c/lens2pexercise3.cc @@ -132,7 +132,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/efm/1p2c_2p_2p2c/lens2pproblem.hh b/lecture/efm/1p2c_2p_2p2c/lens2pproblem.hh index b070beaba60ea7c10277d24795bea37dd1860bfa..6ef1aa5846c7737144ecb962d531a771ae63d543 100644 --- a/lecture/efm/1p2c_2p_2p2c/lens2pproblem.hh +++ b/lecture/efm/1p2c_2p_2p2c/lens2pproblem.hh @@ -85,21 +85,6 @@ public: upperPressure_ = getParam<Scalar>("Boundary.UpperPressure"); } - /*! - * \name Problem parameters - */ - // \{ - - /*! - * \brief Returns the temperature within the domain. - * This problem assumes a uniform temperature of 10 degrees Celsius. - */ - Scalar temperature() const - { - return 273.15 + 10; // -> 10°C - }; - - // \} /*! * \name Boundary conditions diff --git a/lecture/efm/1p2c_2p_2p2c/lens2pspatialparams.hh b/lecture/efm/1p2c_2p_2p2c/lens2pspatialparams.hh index 0921a384cf1dcd6df9e543f9a95d3ad516f04020..0f0b6e59d88855f47b21bdc32ba659f9b127c0c6 100644 --- a/lecture/efm/1p2c_2p_2p2c/lens2pspatialparams.hh +++ b/lecture/efm/1p2c_2p_2p2c/lens2pspatialparams.hh @@ -27,7 +27,7 @@ #define DUMUX_LENS2P_SPATIALPARAMS_HH #include <dumux/material/fluidmatrixinteractions/2p/brookscorey.hh> -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> namespace Dumux { @@ -38,14 +38,14 @@ namespace Dumux { */ template<class FVGridGeometry, class Scalar> class Lens2pSpatialParams -: public FVSpatialParams<FVGridGeometry, Scalar, Lens2pSpatialParams<FVGridGeometry, Scalar>> +: public FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, Lens2pSpatialParams<FVGridGeometry, Scalar>> { using GridView = typename FVGridGeometry::GridView; using Element = typename GridView::template Codim<0>::Entity; using FVElementGeometry = typename FVGridGeometry::LocalView; using SubControlVolume = typename FVElementGeometry::SubControlVolume; using ThisType = Lens2pSpatialParams<FVGridGeometry, Scalar>; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, ThisType>; + using ParentType = FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, ThisType>; static constexpr int dimWorld = GridView::dimensionworld; using GlobalPosition = typename Element::Geometry::GlobalCoordinate; @@ -70,8 +70,9 @@ public: lensK_ = getParam<Scalar>("SpatialParams.Fine.Permeability"); outerK_ = getParam<Scalar>("SpatialParams.Coarse.Permeability"); - } + temperature_ = getParam<Scalar>("SpatialParams.Temperature"); + } /*! * \brief Function for defining the (intrinsic) permeability \f$[m^2]\f$ * @@ -137,6 +138,15 @@ public: return FluidSystem::phase0Idx; } + /*! + * \brief Returns the temperature at the domain at the given position + * \param globalPos The position in global coordinates where the temperature should be specified + */ + Scalar temperatureAtPos(const GlobalPosition& globalPos) const + { + return temperature_; + } + /*! * \brief Set the bounding box of the fine-sand lens * \param lensLowerLeft the lower left corner coordinates vector @@ -166,6 +176,7 @@ private: Scalar outerK_; Scalar lensPorosity_; Scalar outerPorosity_; + Scalar temperature_; const PcKrSw pcKrSwFine_; const PcKrSw pcKrSwCoarse_; const Scalar eps_{1e-6}; diff --git a/lecture/efm/1p2cvs2p/exercise1.input b/lecture/efm/1p2cvs2p/exercise1.input index b36f041e503aa473da20be5e00ebe8d7e41d25f6..50c9b6f2e4591bec4cc9fcb81706c77c91928f4e 100644 --- a/lecture/efm/1p2cvs2p/exercise1.input +++ b/lecture/efm/1p2cvs2p/exercise1.input @@ -12,6 +12,7 @@ UpperRight = 5 4 # upper right corner coordinates (x, Cells = 40 64 # grid resolution in (x,y) direction [-] [SpatialParams] +Temperature = 283.15 # Temperature for the entire domain [K] Fine.Permeability = 3.1e-11 # intrinsic permeability of the fine porous medium [m^2] Coarse.Permeability = 3.1e-10 # intrinsic permeability of the coarse porous medium [m^2] diff --git a/lecture/efm/1p2cvs2p/lens1p2cexercise1.cc b/lecture/efm/1p2cvs2p/lens1p2cexercise1.cc index 3f1460614825ec682005139cae69762030de71e4..16dd4e60969b88f9a7eb1adeff055b2031a0fc26 100644 --- a/lecture/efm/1p2cvs2p/lens1p2cexercise1.cc +++ b/lecture/efm/1p2cvs2p/lens1p2cexercise1.cc @@ -116,7 +116,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/efm/1p2cvs2p/lens1p2cproblem.hh b/lecture/efm/1p2cvs2p/lens1p2cproblem.hh index 601b44a02ac0fdb213fd8608d27c795414a05f03..051b7fe2cee61f50fea5bcef3072154376923e29 100644 --- a/lecture/efm/1p2cvs2p/lens1p2cproblem.hh +++ b/lecture/efm/1p2cvs2p/lens1p2cproblem.hh @@ -107,22 +107,7 @@ public: upperPressure_ = getParam<Scalar>("Boundary.UpperPressure"); } - /*! - * \name Problem parameters - */ - // \{ - - /*! - * \brief Returns the temperature within the domain. - * \param temperature DOC ME! - * This problem assumes a temperature of 10 degrees Celsius. - */ - Scalar temperature() const - { - return 273.15 + 10; // -> 10 degrees Celsius - }; - // \} /*! * \name Boundary conditions diff --git a/lecture/efm/1p2cvs2p/lens2pexercise1.cc b/lecture/efm/1p2cvs2p/lens2pexercise1.cc index 79d723b7a939234eec4683eccd8c03af1e1ff212..a160ff705a403f92a210a793e84d1af06f07a6a1 100644 --- a/lecture/efm/1p2cvs2p/lens2pexercise1.cc +++ b/lecture/efm/1p2cvs2p/lens2pexercise1.cc @@ -134,7 +134,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/efm/1p2cvs2p/lens2pproblem.hh b/lecture/efm/1p2cvs2p/lens2pproblem.hh index 63158c7c8b9ab224df1e74b6ee91fd26132adf62..8f13223f672d133e10a0eecbbc726f005d288695 100644 --- a/lecture/efm/1p2cvs2p/lens2pproblem.hh +++ b/lecture/efm/1p2cvs2p/lens2pproblem.hh @@ -40,7 +40,7 @@ namespace Dumux { * water saturated medium. * * DOC ME! The domain is sized 5m times 4m and features a rectangular lens - * with low permeablility which spans from (1 m , 2 m) to (4 m, 3 m) + * with low permeablility which spans from (0.8 m , 2 m) to (4 m, 3 m) * and is surrounded by a medium with higher permability. Note that * this problem is discretized using only two dimensions, so from the * point of view of the two-phase model, the depth of the domain @@ -116,21 +116,6 @@ public: infiltrationEndTime_= getParam<Scalar>("Boundary.InfiltrationEndTime"); } - /*! - * \name Problem parameters - */ - // \{ - - /*! - * \brief Returns the temperature within the domain. - * This problem assumes a uniform temperature of 10 degrees Celsius. - */ - Scalar temperature() const - { - return 273.15 + 10; // -> 10°C - }; - - // \} /*! * \name Boundary conditions diff --git a/lecture/efm/2p/lens2pexercise2.cc b/lecture/efm/2p/lens2pexercise2.cc index 70afb5cd17988ac973edff1b3c54bcb0f6ca466b..5e5879fec2768c8d9be4019339e066055e7297d7 100644 --- a/lecture/efm/2p/lens2pexercise2.cc +++ b/lecture/efm/2p/lens2pexercise2.cc @@ -134,7 +134,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/efm/2p/lens2pexercise2.input b/lecture/efm/2p/lens2pexercise2.input index 3c7503222c123e3c459f01703aba4cb62ae89001..0e262798531ba71235f8f4822eb3770c0274712d 100644 --- a/lecture/efm/2p/lens2pexercise2.input +++ b/lecture/efm/2p/lens2pexercise2.input @@ -25,6 +25,9 @@ BrooksCoreyPcEntry = 200 # entry pressure for the Brooks-Corey capillary pres Swr = 0.05 # residual saturation of the wetting phase in the coarse soil [-] Snr = 0.0 # residual saturation of the nonwetting phase in the coarse soil [-] +[SpatialParams] +Temperature = 10.0 # temperature in the entire domain [°C] + [Boundary] LowerPressure = 1.19612e5 # Dirichlet pressure value for the boundary condition at the lower boundary [Pa] UpperPressure = 1.0e5 # Dirichlet pressure value for the boundary condition at the upper boundary [Pa] diff --git a/lecture/efm/2p/lens2pproblem.hh b/lecture/efm/2p/lens2pproblem.hh index c2761f12283644fdeeef6113e0d261af57758358..93c796668f9651686d3dd15b2cde41c7a21b8bec 100644 --- a/lecture/efm/2p/lens2pproblem.hh +++ b/lecture/efm/2p/lens2pproblem.hh @@ -130,22 +130,6 @@ public: } } - /*! - * \name Problem parameters - */ - // \{ - - /*! - * \brief Returns the temperature within the domain. - * - * This problem assumes a uniform temperature of 10 degrees Celsius. - */ - Scalar temperature() const - { - return 273.15 + 10; // -> 10°C - }; - - // \} /*! * \name Boundary conditions @@ -184,7 +168,7 @@ public: { PrimaryVariables values(0.0); FluidState fluidState; - fluidState.setTemperature(this->temperature()); + fluidState.setTemperature(this->spatialParams().temperatureAtPos(globalPos)); fluidState.setPressure(FluidSystem::phase0Idx, /*pressure=*/1e5); fluidState.setPressure(FluidSystem::phase1Idx, /*pressure=*/1e5); @@ -260,7 +244,7 @@ public: { PrimaryVariables values; FluidState fluidState; - fluidState.setTemperature(this->temperature()); + fluidState.setTemperature(this->spatialParams().temperatureAtPos(globalPos)); fluidState.setPressure(FluidSystem::phase0Idx, /*pressure=*/1e5); fluidState.setPressure(FluidSystem::phase1Idx, /*pressure=*/1e5); diff --git a/lecture/mm/co2plume/co2plumeshapeexercise.cc b/lecture/mm/co2plume/co2plumeshapeexercise.cc index 12162f20d3e211b862530ba92e4fd52236cb60f9..7b4097a96296c438a5ba7681c164f9d5bafcdc6f 100644 --- a/lecture/mm/co2plume/co2plumeshapeexercise.cc +++ b/lecture/mm/co2plume/co2plumeshapeexercise.cc @@ -79,7 +79,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/co2plume/co2plumeshapespatialparameters.hh b/lecture/mm/co2plume/co2plumeshapespatialparameters.hh index 509bafbc7e24cc2a41ec6480aac305db8abfd2fc..ab6371da3b2845af65a3b7c4ebd17f91aeb4c06f 100644 --- a/lecture/mm/co2plume/co2plumeshapespatialparameters.hh +++ b/lecture/mm/co2plume/co2plumeshapespatialparameters.hh @@ -27,7 +27,7 @@ #define DUMUX_PLUMESHAPE_SPATIAL_PARAMS_HH #include <dumux/io/grid/griddata.hh> -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include <dumux/material/fluidmatrixinteractions/2p/brookscorey.hh> #include <dumux/porousmediumflow/co2/model.hh> @@ -41,7 +41,7 @@ namespace Dumux { */ template<class FVGridGeometry, class Scalar> class PlumeShapeSpatialParams -: public FVSpatialParams<FVGridGeometry, +: public FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, PlumeShapeSpatialParams<FVGridGeometry, Scalar>> { @@ -50,7 +50,7 @@ class PlumeShapeSpatialParams using FVElementGeometry = typename FVGridGeometry::LocalView; using SubControlVolume = typename FVElementGeometry::SubControlVolume; using Element = typename GridView::template Codim<0>::Entity; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, + using ParentType = FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, PlumeShapeSpatialParams<FVGridGeometry, Scalar>>; using GlobalPosition = typename SubControlVolume::GlobalPosition; using PcKrSwCurve = FluidMatrix::BrooksCoreyDefault<Scalar>; diff --git a/lecture/mm/columnxylene/columnxyleneexercise.cc b/lecture/mm/columnxylene/columnxyleneexercise.cc index 09cc76c28084970b0b00884036587b883d10a70d..01fcb54cc589ed69568be67d4a9d17f0c7730d74 100644 --- a/lecture/mm/columnxylene/columnxyleneexercise.cc +++ b/lecture/mm/columnxylene/columnxyleneexercise.cc @@ -102,7 +102,7 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); + // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/columnxylene/columnxylenespatialparams.hh b/lecture/mm/columnxylene/columnxylenespatialparams.hh index b4b6eb1d5793090c2c1b56a09ae13c06aa1223a4..1246233f5ad07d500635e9cced951e59d473e34b 100644 --- a/lecture/mm/columnxylene/columnxylenespatialparams.hh +++ b/lecture/mm/columnxylene/columnxylenespatialparams.hh @@ -26,7 +26,7 @@ #include <dumux/porousmediumflow/3p3c/indices.hh> #include <dumux/common/properties.hh> -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include <dumux/material/fluidmatrixinteractions/3p/parkervangenuchten.hh> #include <dumux/material/fluidmatrixinteractions/3p/napladsorption.hh> @@ -37,7 +37,7 @@ namespace Dumux { */ template<class FVGridGeometry, class Scalar> class ColumnSpatialParams -: public FVSpatialParams<FVGridGeometry, Scalar, +: public FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, ColumnSpatialParams<FVGridGeometry, Scalar>> { using GridView = typename FVGridGeometry::GridView; @@ -45,7 +45,7 @@ class ColumnSpatialParams using SubControlVolume = typename FVElementGeometry::SubControlVolume; using Element = typename GridView::template Codim<0>::Entity; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, + using ParentType = FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, ColumnSpatialParams<FVGridGeometry, Scalar>>; using GlobalPosition = typename SubControlVolume::GlobalPosition; diff --git a/lecture/mm/convectivemixing/convmixexercise.cc b/lecture/mm/convectivemixing/convmixexercise.cc index 3cfabe4bc7c9c86983eb3220fdf9a0d539551a01..77197227fac0a2044d73b2d1449a774993226530 100644 --- a/lecture/mm/convectivemixing/convmixexercise.cc +++ b/lecture/mm/convectivemixing/convmixexercise.cc @@ -84,7 +84,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/convectivemixing/convmixexercise.input b/lecture/mm/convectivemixing/convmixexercise.input index 0f005f7cdfb30f73380ca6f198140fee250fede5..3c0e4f94b84d247f6e7910802d988f8f2c23a424 100644 --- a/lecture/mm/convectivemixing/convmixexercise.input +++ b/lecture/mm/convectivemixing/convmixexercise.input @@ -11,11 +11,12 @@ Cells = 30 30 [SpatialParams] Tortuosity = 0.706 Permeability = 3e-13 +DepthBOR = 1400 # [m] depth below ground surface [Problem] Name = convmix # name passed to the output routines EnableGravity = 1 # enable gravity -DepthBOR = 1400 # [m] depth below ground surface + [Brine] Salinity = 0.05 # brine salinity diff --git a/lecture/mm/convectivemixing/convmixproblem.hh b/lecture/mm/convectivemixing/convmixproblem.hh index 07c72d8a4d3fe1dfdbf0cc11732f14e40c10febc..abbae292783513b807aed0a82b5a6f773abb9b7b 100644 --- a/lecture/mm/convectivemixing/convmixproblem.hh +++ b/lecture/mm/convectivemixing/convmixproblem.hh @@ -90,7 +90,6 @@ public: ConvmixProblem(std::shared_ptr<const FVGridGeometry> fvGridGeometry) : ParentType(fvGridGeometry) { - depthBOR_ = getParam<Scalar>("Problem.DepthBOR"); name_ = getParam<std::string>("Problem.Name"); //initialize fluid system @@ -108,17 +107,6 @@ public: */ // \{ - /*! - * \brief Returns the temperature within the domain [K]. - * - * This problem assumes a temperature of 20 degrees Celsius. - */ - Scalar temperatureAtPos(const GlobalPosition &globalPos) const - { - return (283.15 + (depthBOR_ - globalPos[1])*0.032); // geothermal temperatur in [K] - // assuming 10°C at surface and a - // geothermal gradient of 0.032 K/m - }; // \} @@ -232,10 +220,10 @@ private: { PrimaryVariables priVars; - const Scalar temp = temperatureAtPos(globalPos); + const Scalar temp = this->spatialParams().temperatureAtPos(globalPos); const Scalar salinity = getParam<Scalar>("Brine.Salinity"); - priVars[pressureIdx] = 1.013e5 + (depthBOR_ - globalPos[1]) * 1100 * 9.81; // hydrostatic pressure distribution + priVars[pressureIdx] = 1.013e5 + (this->spatialParams().depthBorAtPos(globalPos) - globalPos[1]) * 1100 * 9.81; // hydrostatic pressure distribution priVars[CO2Idx] = 0.; // initial condition for the CO2 massfraction Scalar moleFracCO2, xgH2O; @@ -254,8 +242,6 @@ private: return priVars; } static constexpr Scalar eps_ = 1e-6; - - Scalar depthBOR_; // bottom of reservoir [m] std::string name_ ; }; diff --git a/lecture/mm/convectivemixing/convmixspatialparams.hh b/lecture/mm/convectivemixing/convmixspatialparams.hh index fc6a76c6d95d300c461c311cac4095fb0171c133..114b993e268272006d19ec63475ae31818254335 100644 --- a/lecture/mm/convectivemixing/convmixspatialparams.hh +++ b/lecture/mm/convectivemixing/convmixspatialparams.hh @@ -25,7 +25,7 @@ #define DUMUX_CONVMIX_SPATIAL_PARAMS_HH #include <dumux/porousmediumflow/properties.hh> -#include <dumux/material/spatialparams/fv1p.hh> +#include <dumux/porousmediumflow/fvspatialparams1p.hh> namespace Dumux { @@ -35,14 +35,14 @@ namespace Dumux { */ template<class FVGridGeometry, class Scalar> class ConvmixSpatialParams -: public FVSpatialParamsOneP<FVGridGeometry, Scalar, +: public FVPorousMediumFlowSpatialParamsOneP<FVGridGeometry, Scalar, ConvmixSpatialParams<FVGridGeometry, Scalar>> { using GridView = typename FVGridGeometry::GridView; using FVElementGeometry = typename FVGridGeometry::LocalView; using SubControlVolume = typename FVElementGeometry::SubControlVolume; using Element = typename GridView::template Codim<0>::Entity; - using ParentType = FVSpatialParamsOneP<FVGridGeometry, Scalar, + using ParentType = FVPorousMediumFlowSpatialParamsOneP<FVGridGeometry, Scalar, ConvmixSpatialParams<FVGridGeometry, Scalar>>; static const int dimWorld = GridView::dimensionworld; @@ -57,6 +57,7 @@ public: { permeability_ = getParam<Scalar>("SpatialParams.Permeability", 3e-13); porosity_ = 0.2; + depthBOR_ = getParam<Scalar>("SpatialParams.DepthBOR"); } /*! @@ -88,9 +89,30 @@ public: const ElementSolution& elemSol) const { return 0; } + /*! + * \brief Returns the temperature at the domain at the given position + * \param globalPos The position in global coordinates where the temperature should be specified + * // geothermal temperatur in [K] + // assuming 10°C at surface and a + // geothermal gradient of 0.032 K/m + */ + Scalar temperatureAtPos(const GlobalPosition& globalPos) const + { + return 283.15 + (depthBOR_ - globalPos[1])*0.032; + } + + /*! + * \brief Define bottom of reservoir \f$\mathrm{[m]}\f$. + * + * \param globalPos The global position + */ + Scalar depthBorAtPos(const GlobalPosition& globalPos) const + { return depthBOR_; } + private: Scalar permeability_; Scalar porosity_; + Scalar depthBOR_; }; } // end namespace Dumux diff --git a/lecture/mm/fractures/fracture_exercise.input b/lecture/mm/fractures/fracture_exercise.input index 9fe6a9cc427a1526516dbecb75855f2915194443..bf1d0e3736f8c6354e5554c28070a726864e3190 100644 --- a/lecture/mm/fractures/fracture_exercise.input +++ b/lecture/mm/fractures/fracture_exercise.input @@ -22,6 +22,7 @@ Problem.Name = matrix Problem.Name = fractures [SpatialParams] +Temperature = 283.15 #[K] Aperture = 1e-1 Permeability = 1e-14 Porosity = 0.1 diff --git a/lecture/mm/fractures/fractureproblem.hh b/lecture/mm/fractures/fractureproblem.hh index 7f32eacfea1ad8d3685f88816e8808d03d50957c..34e49b651ef10ebad7a0d8463e188bdce55ffa9f 100644 --- a/lecture/mm/fractures/fractureproblem.hh +++ b/lecture/mm/fractures/fractureproblem.hh @@ -111,14 +111,6 @@ public: return source; } - //! Set the aperture as extrusion factor. - Scalar extrusionFactorAtPos(const GlobalPosition& globalPos) const - { - // We treat the fractures as lower-dimensional in the grid, - // but we have to give it the aperture as extrusion factor - // such that the dimensions are correct in the end. - return aperture_; - } //! evaluates the Dirichlet boundary condition for a given position PrimaryVariables dirichletAtPos(const GlobalPosition& globalPos) const @@ -167,16 +159,12 @@ public: if (scvf.boundary()) flux -= neumannAtPos(scvf.ipGlobal())[contiH2OEqIdx] *scvf.area() - *extrusionFactorAtPos(scvf.ipGlobal()); + *this->spatialParams().extrusionFactorAtPos(scvf.ipGlobal()); } return flux; } - //! returns the temperature in \f$\mathrm{[K]}\f$ in the domain - Scalar temperature() const - { return 283.15; /*10°*/ } - //! sets the pointer to the coupling manager. void setCouplingManager(std::shared_ptr<CouplingManager> cm) { couplingManagerPtr_ = cm; } diff --git a/lecture/mm/fractures/fractures.cc b/lecture/mm/fractures/fractures.cc index 739de2748b674dd7cb366917baf450c90fdf90af..c4f05e1dfa0984293892642ba1502bc8e6791471 100644 --- a/lecture/mm/fractures/fractures.cc +++ b/lecture/mm/fractures/fractures.cc @@ -99,8 +99,7 @@ int main(int argc, char** argv) try // create the finite volume grid geometries auto matrixFvGridGeometry = std::make_shared<MatrixFVGridGeometry>(matrixGridView); auto fractureFvGridGeometry = std::make_shared<FractureFVGridGeometry>(fractureGridView); - matrixFvGridGeometry->update(); - fractureFvGridGeometry->update(); + // the problems (boundary/initial conditions etc) using MatrixProblem = GetPropType<MatrixTypeTag, Properties::Problem>; diff --git a/lecture/mm/fractures/fracturespatialparams.hh b/lecture/mm/fractures/fracturespatialparams.hh index 8f6277bace3dabfe3a65d6ba5750f88156397f53..38ff9d3b3bdbe8d172721524842e7e2d1bcaffaa 100644 --- a/lecture/mm/fractures/fracturespatialparams.hh +++ b/lecture/mm/fractures/fracturespatialparams.hh @@ -24,7 +24,7 @@ #ifndef DUMUX_COURSE_FRACTURESEXERCISE_FRACTURE_SPATIALPARAMS_HH #define DUMUX_COURSE_FRACTURESEXERCISE_FRACTURE_SPATIALPARAMS_HH -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include <dumux/material/fluidmatrixinteractions/2p/vangenuchten.hh> namespace Dumux { @@ -34,10 +34,10 @@ namespace Dumux { */ template< class FVGridGeometry, class Scalar > class FractureSpatialParams -: public FVSpatialParams< FVGridGeometry, Scalar, FractureSpatialParams<FVGridGeometry, Scalar> > +: public FVPorousMediumFlowSpatialParamsMP< FVGridGeometry, Scalar, FractureSpatialParams<FVGridGeometry, Scalar> > { using ThisType = FractureSpatialParams< FVGridGeometry, Scalar >; - using ParentType = FVSpatialParams< FVGridGeometry, Scalar, ThisType >; + using ParentType = FVPorousMediumFlowSpatialParamsMP< FVGridGeometry, Scalar, ThisType >; using SubControlVolume = typename FVGridGeometry::SubControlVolume; using GridView = typename FVGridGeometry::GridView; @@ -64,6 +64,8 @@ public: { porosity_ = getParamFromGroup<Scalar>(paramGroup, "SpatialParams.Fracture.Porosity"); permeability_ = getParamFromGroup<Scalar>(paramGroup, "SpatialParams.Fracture.Permeability"); + temperature_ = getParamFromGroup<Scalar>(paramGroup, "SpatialParams.Temperature"); + aperture_ = getParamFromGroup<Scalar>(paramGroup, "SpatialParams.Aperture"); } //! Function for defining the (intrinsic) permeability \f$[m^2]\f$. @@ -93,8 +95,28 @@ public: return FluidSystem::phase0Idx; } + /*! + * \brief Returns the temperature at the domain at the given position + * \param globalPos The position in global coordinates where the temperature should be specified + */ + Scalar temperatureAtPos(const GlobalPosition& globalPos) const + { + return temperature_; + } + + //! Set the aperture as extrusion factor. + Scalar extrusionFactorAtPos(const GlobalPosition& globalPos) const + { + // We treat the fractures as lower-dimensional in the grid, + // but we have to give it the aperture as extrusion factor + // such that the dimensions are correct in the end. + return aperture_; + } + private: Scalar porosity_; + Scalar temperature_; + Scalar aperture_; PermeabilityType permeability_; const PcKrSw pcKrSwFracture_; }; diff --git a/lecture/mm/fractures/matrixproblem.hh b/lecture/mm/fractures/matrixproblem.hh index bb556bfb69ad625fc38c067f43d0eb70e0c76bdf..bb20ee1d2a28c8a5adadacb780a55f2719ce021f 100644 --- a/lecture/mm/fractures/matrixproblem.hh +++ b/lecture/mm/fractures/matrixproblem.hh @@ -244,10 +244,6 @@ public: return flux; } - //! returns the temperature in \f$\mathrm{[K]}\f$ in the domain - Scalar temperature() const - { return 283.15; /*10°*/ } - //! sets the pointer to the coupling manager. void setCouplingManager(std::shared_ptr<CouplingManager> cm) { couplingManagerPtr_ = cm; } diff --git a/lecture/mm/fractures/matrixspatialparams.hh b/lecture/mm/fractures/matrixspatialparams.hh index 29c37e202429baf08ce84c2e4e0a53ef957b6d9c..d1cc419114e68762793e1ee814ac05501bbdec8c 100644 --- a/lecture/mm/fractures/matrixspatialparams.hh +++ b/lecture/mm/fractures/matrixspatialparams.hh @@ -24,7 +24,8 @@ #ifndef DUMUX_COURSE_FRACTURESEXERCISE_MATRIX_SPATIALPARAMS_HH #define DUMUX_COURSE_FRACTURESEXERCISE_MATRIX_SPATIALPARAMS_HH -#include <dumux/material/spatialparams/fv.hh> + +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include <dumux/material/fluidmatrixinteractions/2p/vangenuchten.hh> namespace Dumux { @@ -34,10 +35,10 @@ namespace Dumux { */ template< class FVGridGeometry, class Scalar > class MatrixSpatialParams -: public FVSpatialParams< FVGridGeometry, Scalar, MatrixSpatialParams<FVGridGeometry, Scalar> > +: public FVPorousMediumFlowSpatialParamsMP< FVGridGeometry, Scalar, MatrixSpatialParams<FVGridGeometry, Scalar> > { using ThisType = MatrixSpatialParams< FVGridGeometry, Scalar >; - using ParentType = FVSpatialParams< FVGridGeometry, Scalar, ThisType >; + using ParentType = FVPorousMediumFlowSpatialParamsMP< FVGridGeometry, Scalar, ThisType >; using GridView = typename FVGridGeometry::GridView; using Grid = typename GridView::Grid; @@ -60,6 +61,7 @@ public: { porosity_ = getParamFromGroup<Scalar>(paramGroup, "SpatialParams.Porosity"); permeability_ = getParamFromGroup<Scalar>(paramGroup, "SpatialParams.Permeability"); + temperature_ = getParamFromGroup<Scalar>(paramGroup, "SpatialParams.Temperature"); porosityOverburden_ = getParamFromGroup<Scalar>(paramGroup, "SpatialParams.Overburden.Porosity"); permeabilityOverburden_ = getParamFromGroup<Scalar>(paramGroup, "SpatialParams.Overburden.Permeability"); } @@ -97,6 +99,15 @@ public: return FluidSystem::phase0Idx; } + /*! + * \brief Returns the temperature at the domain at the given position + * \param globalPos The position in global coordinates where the temperature should be specified + */ + Scalar temperatureAtPos(const GlobalPosition& globalPos) const + { + return temperature_; + } + private: Scalar porosity_; Scalar porosityOverburden_; @@ -104,6 +115,7 @@ private: PermeabilityType permeabilityOverburden_; const PcKrSw pcKrSw_; const PcKrSw pcKrSwOverburden_; + Scalar temperature_; }; } // end namespace Dumux diff --git a/lecture/mm/fuelcell/fuelcell.cc b/lecture/mm/fuelcell/fuelcell.cc index 2ca67fa3b3522becbd4998cb414590accfc4148a..19599b3ddb2c9a6f797768234fd97c7ea0acbb1a 100644 --- a/lecture/mm/fuelcell/fuelcell.cc +++ b/lecture/mm/fuelcell/fuelcell.cc @@ -116,7 +116,7 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); + // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/fuelcell/fuelcell.input b/lecture/mm/fuelcell/fuelcell.input index 464cb9dd4796cebd7d9552948c010e4f0ed66066..49cac53e207a2ee4d5b56012c4725c102e03798c 100644 --- a/lecture/mm/fuelcell/fuelcell.input +++ b/lecture/mm/fuelcell/fuelcell.input @@ -45,6 +45,7 @@ MaxRelativeShift = 1e-5 [OperationalConditions] GasPressureInlet1 = 2.013e5 # [Pa] pressure at the upper inlet GasPressureInlet2 = 2.013e5 # [Pa] pressure at the lower inlet # 2.013e5 +# for interdigitated conditions, GasPressureInlet1 may be set, e.g. to 2.015e5 LiquidWaterSaturationInlet = 0.055 #LiquidWaterSaturationInlet = 0.41 # drainage diff --git a/lecture/mm/fuelcell/fuelcellspatialparams.hh b/lecture/mm/fuelcell/fuelcellspatialparams.hh index 44e0f652037c41c920aed12d4248aea1bb92bafb..2ec4f0217adac15be2804fbaa46a2635caedd00b 100644 --- a/lecture/mm/fuelcell/fuelcellspatialparams.hh +++ b/lecture/mm/fuelcell/fuelcellspatialparams.hh @@ -27,7 +27,7 @@ #include <dumux/common/properties.hh> -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include "./material/acosta.hh" #include "./material/thermalconductivityconstant.hh" @@ -61,14 +61,14 @@ struct ThermalConductivityModel<TypeTag, TTag::FuelCellLectureSpatialParams> { u * problem which uses the (non-)isothermal 2pnc box model */ template<class TypeTag> -class FuelCellLectureSpatialParams: public FVSpatialParams<GetPropType<TypeTag, Properties::GridGeometry>, +class FuelCellLectureSpatialParams: public FVPorousMediumFlowSpatialParamsMP<GetPropType<TypeTag, Properties::GridGeometry>, GetPropType<TypeTag, Properties::Scalar>, FuelCellLectureSpatialParams<TypeTag>> { using ThisType = FuelCellLectureSpatialParams<TypeTag>; using Scalar = GetPropType<TypeTag, Properties::Scalar>; using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, ThisType>; + using ParentType = FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, ThisType>; using GridView = typename FVGridGeometry::GridView; using FluidSystem = GetPropType<TypeTag, Properties::FluidSystem>; using Element = typename GridView::template Codim<0>::Entity; diff --git a/lecture/mm/heatpipe/heatpipe.cc b/lecture/mm/heatpipe/heatpipe.cc index bc47cb037250fb05ebf8a0811525d049144ea9eb..d78b26b4574d0da92908b17366411acbb732cc2a 100644 --- a/lecture/mm/heatpipe/heatpipe.cc +++ b/lecture/mm/heatpipe/heatpipe.cc @@ -81,7 +81,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/heatpipe/heatpipespatialparams.hh b/lecture/mm/heatpipe/heatpipespatialparams.hh index e05d26d5078b61c7782d75494054f77e511a95a9..e6a525975634b3302a0b4bdd43a3f0b30b3aa55c 100644 --- a/lecture/mm/heatpipe/heatpipespatialparams.hh +++ b/lecture/mm/heatpipe/heatpipespatialparams.hh @@ -19,7 +19,7 @@ #ifndef DUMUX_HEATPIPE_SPATIAL_PARAMS_HH #define DUMUX_HEATPIPE_SPATIAL_PARAMS_HH -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include "krpcheatpipe.hh" @@ -28,11 +28,11 @@ namespace Dumux template<class FVGridGeometry, class Scalar> class HeatPipeSpatialParams -: public FVSpatialParams<FVGridGeometry, Scalar, HeatPipeSpatialParams<FVGridGeometry, Scalar>> +: public FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, HeatPipeSpatialParams<FVGridGeometry, Scalar>> { using ThisType = HeatPipeSpatialParams<FVGridGeometry, Scalar>; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, ThisType>; + using ParentType = FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, ThisType>; using GridView = typename FVGridGeometry::GridView; using Element = typename GridView::template Codim<0>::Entity; diff --git a/lecture/mm/heavyoil/sagd/sagd.cc b/lecture/mm/heavyoil/sagd/sagd.cc index d1b6da09a7491415e985185ca988d42f080ec65f..1260cbf8fa8bc54c75da3dd3b00f8d4e9b649ef2 100644 --- a/lecture/mm/heavyoil/sagd/sagd.cc +++ b/lecture/mm/heavyoil/sagd/sagd.cc @@ -106,7 +106,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/heavyoil/sagd/spatialparams.hh b/lecture/mm/heavyoil/sagd/spatialparams.hh index a1185077e944fefe7ca3f85cf544132fd398be60..b98c0afcfa56a48b11f31ba31e1bf215f2a9366e 100644 --- a/lecture/mm/heavyoil/sagd/spatialparams.hh +++ b/lecture/mm/heavyoil/sagd/spatialparams.hh @@ -26,7 +26,7 @@ #include <dumux/porousmediumflow/properties.hh> #include <dumux/porousmediumflow/3pwateroil/indices.hh> -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include <dumux/material/fluidmatrixinteractions/3p/parkervangenuchten.hh> namespace Dumux { @@ -37,7 +37,7 @@ namespace Dumux { */ template<class FVGridGeometry, class Scalar> class SagdSpatialParams -: public FVSpatialParams<FVGridGeometry, Scalar, +: public FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, SagdSpatialParams<FVGridGeometry, Scalar>> { using GridView = typename FVGridGeometry::GridView; @@ -49,7 +49,7 @@ class SagdSpatialParams using GlobalPosition = typename SubControlVolume::GlobalPosition; using Element = typename GridView::template Codim<0>::Entity; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, + using ParentType = FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, SagdSpatialParams<FVGridGeometry, Scalar>>; using ThreePhasePcKrSw = FluidMatrix::ParkerVanGenuchten3PDefault<Scalar>; diff --git a/lecture/mm/heavyoil/sagdcyclic/sagd_cyclic.cc b/lecture/mm/heavyoil/sagdcyclic/sagd_cyclic.cc index d36dda2daff3794143eded07180b5575fde82f1c..ff10528cb42e2f94902ab196e27ca1f3ed3e625d 100644 --- a/lecture/mm/heavyoil/sagdcyclic/sagd_cyclic.cc +++ b/lecture/mm/heavyoil/sagdcyclic/sagd_cyclic.cc @@ -106,7 +106,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/heavyoil/sagdcyclic/spatialparams.hh b/lecture/mm/heavyoil/sagdcyclic/spatialparams.hh index 0fc1bd1e6f63666bda8b3978f8a45415592a7929..fe7545686ec770f70f7b71fcac83e0bf19f2a3b7 100644 --- a/lecture/mm/heavyoil/sagdcyclic/spatialparams.hh +++ b/lecture/mm/heavyoil/sagdcyclic/spatialparams.hh @@ -26,7 +26,7 @@ #include <dumux/porousmediumflow/properties.hh> #include <dumux/porousmediumflow/3pwateroil/indices.hh> -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include <dumux/material/fluidmatrixinteractions/3p/parkervangenuchten.hh> namespace Dumux { @@ -36,7 +36,7 @@ namespace Dumux { */ template<class FVGridGeometry, class Scalar> class SagdSpatialParams -: public FVSpatialParams<FVGridGeometry, Scalar, +: public FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, SagdSpatialParams<FVGridGeometry, Scalar>> { using GridView = typename FVGridGeometry::GridView; @@ -48,7 +48,7 @@ class SagdSpatialParams using GlobalPosition = typename SubControlVolume::GlobalPosition; using Element = typename GridView::template Codim<0>::Entity; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, + using ParentType = FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, SagdSpatialParams<FVGridGeometry, Scalar>>; using ThreePhasePcKrSw = FluidMatrix::ParkerVanGenuchten3PDefault<Scalar>; diff --git a/lecture/mm/heavyoil/sagdcyclichyst/sagd_cyclic_hyst.cc b/lecture/mm/heavyoil/sagdcyclichyst/sagd_cyclic_hyst.cc index 52f9b326ee35efb55a71d1ce8cce4f0a9a7028fe..54a49221d2fedd5b403ed3a8c719fec0d0f7fd3b 100644 --- a/lecture/mm/heavyoil/sagdcyclichyst/sagd_cyclic_hyst.cc +++ b/lecture/mm/heavyoil/sagdcyclichyst/sagd_cyclic_hyst.cc @@ -106,7 +106,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/heavyoil/sagdcyclichyst/spatialparams.hh b/lecture/mm/heavyoil/sagdcyclichyst/spatialparams.hh index acae9de6ad4e8bb65d1fa51905b97d0c2df4acfc..135b718c413c347932c1a90f8d4f65944cf5cd2d 100644 --- a/lecture/mm/heavyoil/sagdcyclichyst/spatialparams.hh +++ b/lecture/mm/heavyoil/sagdcyclichyst/spatialparams.hh @@ -25,7 +25,7 @@ #define DUMUX_SAGD_SPATIAL_PARAMS_HH #include <dumux/porousmediumflow/properties.hh> -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include <dumux/porousmediumflow/3pwateroil/indices.hh> @@ -38,7 +38,7 @@ namespace Dumux { */ template<class FVGridGeometry, class Scalar> class SagdSpatialParams -: public FVSpatialParams<FVGridGeometry, Scalar, +: public FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, SagdSpatialParams<FVGridGeometry, Scalar>> { using GridView = typename FVGridGeometry::GridView; @@ -50,7 +50,7 @@ class SagdSpatialParams using GlobalPosition = typename SubControlVolume::GlobalPosition; using Element = typename GridView::template Codim<0>::Entity; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, + using ParentType = FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, SagdSpatialParams<FVGridGeometry, Scalar>>; static constexpr bool isBox = FVGridGeometry::discMethod == DiscretizationMethod::box; diff --git a/lecture/mm/henryproblem/henry1p2c/henry1p2c.cc b/lecture/mm/henryproblem/henry1p2c/henry1p2c.cc index 918dc03af54792c0beabb3ab7fa32137202a751e..8118c802c83072126693d34298ff584303c5343e 100644 --- a/lecture/mm/henryproblem/henry1p2c/henry1p2c.cc +++ b/lecture/mm/henryproblem/henry1p2c/henry1p2c.cc @@ -109,7 +109,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/henryproblem/henry1p2c/henry1p2c.input b/lecture/mm/henryproblem/henry1p2c/henry1p2c.input index 0517f7ee771f3ebc57a78bb7bf71dc32b4e2e698..d9e32d1082ae723a054d20f25f3d057e3dc14591 100644 --- a/lecture/mm/henryproblem/henry1p2c/henry1p2c.input +++ b/lecture/mm/henryproblem/henry1p2c/henry1p2c.input @@ -7,7 +7,7 @@ File = ./grids/henry.dgf [Problem] freshWaterFluxRate = 6.6E-2 # [kg/sm^2] fresh water flux rate at the left boundary -diffusionCoefficient = 6.6e-6 # 1.2e-9 #[m2/s] Diffusion coefficient +diffusionCoefficient = 0.12e-9 #[m2/s] Diffusion coefficient Name = Henry1p2c AlphaL = 0.1 # longitual dispersivity AlphaT = 0.01 # transversal dispersivity diff --git a/lecture/mm/henryproblem/henry1p2c/henry1p2cproblem.hh b/lecture/mm/henryproblem/henry1p2c/henry1p2cproblem.hh index 2c0a22ba377f7d72eba1a9cd510296fb6fda5e95..d6a449fd5ba896c3219849cac0033f6f6f9c5f28 100644 --- a/lecture/mm/henryproblem/henry1p2c/henry1p2cproblem.hh +++ b/lecture/mm/henryproblem/henry1p2c/henry1p2cproblem.hh @@ -70,23 +70,6 @@ public: freshWaterFluxRate_= getParam<Scalar>("Problem.freshWaterFluxRate"); } - /*! - * \name Problem parameters - */ - // \{ - - /*! - * \brief Returns the temperature within the domain. - * - * This problem assumes a temperature of 36 degrees Celsius. - */ - Scalar temperature() const - { - return 273.15 + 20; // in [K] - }; - - // \} - /*! * \name Boundary conditions */ diff --git a/lecture/mm/henryproblem/henry1p2c/henry1p2cspatialparameters.hh b/lecture/mm/henryproblem/henry1p2c/henry1p2cspatialparameters.hh index d4fc7f5b1dbb5851a4a6b0e003ae6eb1fe70922e..d39cf35baa88b0612b44f949060b95e5d57e7bf8 100644 --- a/lecture/mm/henryproblem/henry1p2c/henry1p2cspatialparameters.hh +++ b/lecture/mm/henryproblem/henry1p2c/henry1p2cspatialparameters.hh @@ -25,8 +25,7 @@ #ifndef DUMUX_HENRY1P2C_SPATIAL_PARAMETERS_HH #define DUMUX_HENRY1P2C_SPATIAL_PARAMETERS_HH -#include <dumux/material/spatialparams/fv.hh> -#include <dumux/material/fluidmatrixinteractions/2p/linearmaterial.hh> +#include <dumux/porousmediumflow/fvspatialparams1p.hh> namespace Dumux { @@ -37,7 +36,7 @@ namespace Dumux { */ template<class FVGridGeometry, class Scalar> class Henry1p2cSpatialParams : -public FVSpatialParams<FVGridGeometry, Scalar, +public FVPorousMediumFlowSpatialParamsOneP<FVGridGeometry, Scalar, Henry1p2cSpatialParams<FVGridGeometry, Scalar>> { using GridView = typename FVGridGeometry::GridView; @@ -45,7 +44,7 @@ public FVSpatialParams<FVGridGeometry, Scalar, using SubControlVolume = typename FVElementGeometry::SubControlVolume; using Element = typename GridView::template Codim<0>::Entity; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, + using ParentType = FVPorousMediumFlowSpatialParamsOneP<FVGridGeometry, Scalar, Henry1p2cSpatialParams<FVGridGeometry, Scalar>>; using GlobalPosition = typename SubControlVolume::GlobalPosition; @@ -59,6 +58,7 @@ public: // intrinsic permeabilities perm_ = 1.019368e-9; porosity_ = 0.35; + temperature_ = 293.15; } /*! @@ -91,9 +91,21 @@ public: return std::array<Scalar, 2> {alphaL, alphaT}; } + /*! + * \brief Returns the temperature within the domain. + * + * This problem assumes a temperature of 20 degrees Celsius. + */ + Scalar temperatureAtPos(const GlobalPosition& globalPos) const + { + return temperature_; // in [K] + }; + + private: - Scalar perm_; + PermeabilityType perm_; Scalar porosity_; + Scalar temperature_; }; } // end namespace Dumux diff --git a/lecture/mm/henryproblem/henry2p/henry2p.cc b/lecture/mm/henryproblem/henry2p/henry2p.cc index 1ee74a01885417a13d24d78973aecdfd583608f9..050d33f78fcdc56bb19387735a26c1705165df5e 100644 --- a/lecture/mm/henryproblem/henry2p/henry2p.cc +++ b/lecture/mm/henryproblem/henry2p/henry2p.cc @@ -111,7 +111,7 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); + // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/henryproblem/henry2p/henry2pproblem.hh b/lecture/mm/henryproblem/henry2p/henry2pproblem.hh index ecaebbac76b402bcb8eaca4e8c9c8df810e905cd..5a733fea73fa67afd62abc2baa861011ffa6233d 100644 --- a/lecture/mm/henryproblem/henry2p/henry2pproblem.hh +++ b/lecture/mm/henryproblem/henry2p/henry2pproblem.hh @@ -67,7 +67,6 @@ public: { FluidSystem::init(); eps_ = 3e-6; - temperature_ = 273.15 + 20; // -> 20°C freshWaterFluxRate_= getParam<Scalar>("Problem.freshWaterFluxRate"); } @@ -76,17 +75,6 @@ public: */ // \{ - /*! - * \brief Returns the temperature within the domain. - * - * This problem assumes a uniform temperature of 10 degrees Celsius. - */ - Scalar temperature() const - { - return temperature_; - } - - NumEqVector sourceAtPos( const GlobalPosition &globalPos) const { return NumEqVector(0.0); @@ -235,7 +223,6 @@ private: } static constexpr Scalar depthBOR_ = 1.0; - Scalar temperature_; Scalar eps_; Scalar freshWaterFluxRate_; }; diff --git a/lecture/mm/henryproblem/henry2p/henry2pspatialparams.hh b/lecture/mm/henryproblem/henry2p/henry2pspatialparams.hh index f55a5a522f3591d68fdbbca606702c9d9212614f..29aa4f19f8b0e188aed163af7e100378139c3a65 100644 --- a/lecture/mm/henryproblem/henry2p/henry2pspatialparams.hh +++ b/lecture/mm/henryproblem/henry2p/henry2pspatialparams.hh @@ -25,7 +25,7 @@ #ifndef DUMUX_HENRY2P_SPATIAL_PARAMS_HH #define DUMUX_HENRY2P_SPATIAL_PARAMS_HH -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include <dumux/material/fluidmatrixinteractions/2p/linearmaterial.hh> #include <dumux/porousmediumflow/2p/model.hh> @@ -38,14 +38,14 @@ namespace Dumux { template<class FVGridGeometry, class Scalar> class Henry2pSpatialParams : -public FVSpatialParams<FVGridGeometry, Scalar, Henry2pSpatialParams<FVGridGeometry, Scalar>> +public FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, Henry2pSpatialParams<FVGridGeometry, Scalar>> { using ThisType = Henry2pSpatialParams<FVGridGeometry, Scalar>; using GridView = typename FVGridGeometry::GridView; using Element = typename GridView::template Codim<0>::Entity; using FVElementGeometry = typename FVGridGeometry::LocalView; using SubControlVolume = typename FVElementGeometry::SubControlVolume; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, ThisType>; + using ParentType = FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, ThisType>; static constexpr int dimWorld = GridView::dimensionworld; using GlobalPosition = typename Element::Geometry::GlobalCoordinate; @@ -71,8 +71,9 @@ public: typename PcKrSwCurve::BasicParams params(0.0/*pcEntry*/, 0.0/*pcMax*/); pcKrSwCurve_ = std::make_unique<PcKrSwCurve>(params); - K_ = 1.019368e-9; + k_ = 1.019368e-9; porosity_=0.35; + temperature_ = 293.15; } /*! @@ -81,14 +82,14 @@ public: * \param globalPos The global position */ PermeabilityType permeabilityAtPos(const GlobalPosition& globalPos) const - { return K_; } + { return k_; } /*! * \brief Define the porosity \f$\mathrm{[-]}\f$. * * \param globalPos The global position */ - double porosityAtPos(const GlobalPosition& globalPos) const + Scalar porosityAtPos(const GlobalPosition& globalPos) const { return porosity_; } /*! @@ -111,9 +112,20 @@ public: int wettingPhaseAtPos(const GlobalPosition& globalPos) const { return FluidSystem::phase0Idx; } + /*! + * \brief Returns the temperature within the domain. + * + * This problem assumes a temperature of 20 degrees Celsius. + */ + Scalar temperatureAtPos(const GlobalPosition& globalPos) const + { + return temperature_; // in [K] + }; + private: - Scalar K_; + PermeabilityType k_; Scalar porosity_; + Scalar temperature_; std::unique_ptr<PcKrSwCurve> pcKrSwCurve_; diff --git a/lecture/mm/naplinfiltration/3p/naplinfiltration3p.cc b/lecture/mm/naplinfiltration/3p/naplinfiltration3p.cc index 09c200f2f0d9814dd95a49e1feae615798bbbb81..b6ccfb425d4687b6638070730e7525f5b9c90878 100644 --- a/lecture/mm/naplinfiltration/3p/naplinfiltration3p.cc +++ b/lecture/mm/naplinfiltration/3p/naplinfiltration3p.cc @@ -112,7 +112,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/naplinfiltration/3p/naplinfiltration3p.input b/lecture/mm/naplinfiltration/3p/naplinfiltration3p.input index 770afb7a82c4da72efc54fb081c7104eb84a672f..d4353182b7f674cdfa3ff2160cfd9d9bbbaf3b91 100644 --- a/lecture/mm/naplinfiltration/3p/naplinfiltration3p.input +++ b/lecture/mm/naplinfiltration/3p/naplinfiltration3p.input @@ -9,6 +9,7 @@ UpperRight = 500 10 # coordinates of upper right grid corner [m] Cells = 250 10 # number of cells in (x, y) direction [-] [SpatialParams] +temperature = 283.15 # [K] permeability = 1.e-11 # [m^2] porosity = 0.40 # [-] Swr = 0.12 diff --git a/lecture/mm/naplinfiltration/3p/problem.hh b/lecture/mm/naplinfiltration/3p/problem.hh index 6bc6273e2573a481fd270fe8988bc4f8707fad2b..954445562aa73afffe0adacffed86001d1f62472 100644 --- a/lecture/mm/naplinfiltration/3p/problem.hh +++ b/lecture/mm/naplinfiltration/3p/problem.hh @@ -95,7 +95,6 @@ public: InfiltrationThreePProblem(std::shared_ptr<const FVGridGeometry> fvGridGeometry) : ParentType(fvGridGeometry) { - temperature_ = 273.15 + 10.0; // -> 10 degrees Celsius FluidSystem::init(282.15, 284.15, 3, 8e4, 3e5, 200); name_ = getParam<std::string>("Problem.Name"); @@ -119,18 +118,6 @@ public: const std::string& name() const { return name_; } - /*! - * \brief Returns the temperature within the domain. - * - * \param globalPos The position - * - * This problem assumes a temperature of 10 degrees Celsius. - */ - Scalar temperatureAtPos(const GlobalPosition &globalPos) const - { - return temperature_; - } - // \} /*! @@ -217,16 +204,6 @@ public: return initial_(globalPos); } - /*! - * \brief Returns the temperature within the domain. - * - * This problem assumes a uniform temperature of 10 degrees Celsius. - */ - Scalar temperature() const - { return temperature_; } - - - private: // internal method for the initial condition (reused for the // dirichlet conditions!) @@ -289,7 +266,6 @@ private: return sw; } - Scalar temperature_; static constexpr Scalar eps_ = 1e-6; std::string name_; Scalar time_; diff --git a/lecture/mm/naplinfiltration/3p3c/naplinfiltration3p3c.cc b/lecture/mm/naplinfiltration/3p3c/naplinfiltration3p3c.cc index 0a8816fba49c2dd2b6a4225517d595019084df61..9c6722e70890030dd136b355569c10134805956d 100644 --- a/lecture/mm/naplinfiltration/3p3c/naplinfiltration3p3c.cc +++ b/lecture/mm/naplinfiltration/3p3c/naplinfiltration3p3c.cc @@ -109,7 +109,6 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/naplinfiltration/3p3c/naplinfiltration3p3c.input b/lecture/mm/naplinfiltration/3p3c/naplinfiltration3p3c.input index b5500a54546065bb1dfe228f4077ee2ec461420c..b5fede1bdb21ea708acf6c7ad17e0877d72d686a 100644 --- a/lecture/mm/naplinfiltration/3p3c/naplinfiltration3p3c.input +++ b/lecture/mm/naplinfiltration/3p3c/naplinfiltration3p3c.input @@ -9,6 +9,7 @@ UpperRight = 500 10 # coordinates of upper right grid corner [m] Cells = 250 10 # number of cells in (x, y) direction [-] [SpatialParams] +temperature = 283.15 #[K] permeability = 1.e-11 # [m^2] porosity = 0.40 # [-] Swr = 0.12 diff --git a/lecture/mm/naplinfiltration/3p3c/problem.hh b/lecture/mm/naplinfiltration/3p3c/problem.hh index 3139c392e9708478662726235db211ad8faeb065..e2eb41786b95a48e1180d6e4669f2d3fb8147aae 100644 --- a/lecture/mm/naplinfiltration/3p3c/problem.hh +++ b/lecture/mm/naplinfiltration/3p3c/problem.hh @@ -107,7 +107,6 @@ public: InfiltrationThreePThreeCProblem(std::shared_ptr<const FVGridGeometry> fvGridGeometry) : ParentType(fvGridGeometry) { - temperature_ = 273.15 + 10.0; // -> 10 degrees Celsius FluidSystem::init(282.15, 284.15, 3, 8e4, 3e5, 200); name_ = getParam<std::string>("Problem.Name"); @@ -125,18 +124,6 @@ public: const std::string& name() const { return name_; } - /*! - * \brief Returns the temperature within the domain. - * - * \param globalPos The position - * - * This problem assumes a temperature of 10 degrees Celsius. - */ - Scalar temperatureAtPos(const GlobalPosition &globalPos) const - { - return temperature_; - } - // \} /*! @@ -283,7 +270,6 @@ private: return sw; } - Scalar temperature_; static constexpr Scalar eps_ = 1e-6; std::string name_; Scalar time_; diff --git a/lecture/mm/naplinfiltration/spatialparams.hh b/lecture/mm/naplinfiltration/spatialparams.hh index 7722719d7d6fd496e8bf243e22e2b079fcef394e..44e54fb5bca5d8e3217f0df329fc42f96779a3f0 100644 --- a/lecture/mm/naplinfiltration/spatialparams.hh +++ b/lecture/mm/naplinfiltration/spatialparams.hh @@ -25,7 +25,7 @@ #define DUMUX_NAPLINFILTRATION_SPATIAL_PARAMS_HH #include <dumux/porousmediumflow/properties.hh> -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include <dumux/material/fluidmatrixinteractions/3p/parkervangenuchten.hh> #include <dumux/material/fluidmatrixinteractions/3p/napladsorption.hh> #include <dumux/io/gnuplotinterface.hh> @@ -37,14 +37,14 @@ namespace Dumux { */ template<class FVGridGeometry, class Scalar> class InfiltrationSpatialParams -: public FVSpatialParams<FVGridGeometry, Scalar, +: public FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, InfiltrationSpatialParams<FVGridGeometry, Scalar>> { using GridView = typename FVGridGeometry::GridView; using FVElementGeometry = typename FVGridGeometry::LocalView; using SubControlVolume = typename FVElementGeometry::SubControlVolume; using Element = typename GridView::template Codim<0>::Entity; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, + using ParentType = FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, InfiltrationSpatialParams<FVGridGeometry, Scalar>>; using GlobalPosition = typename SubControlVolume::GlobalPosition; using ThreePhasePcKrSw = FluidMatrix::ParkerVanGenuchten3PDefault<Scalar>; @@ -70,6 +70,9 @@ public: // porosities porosity_ = getParam<Scalar>("SpatialParams.porosity"); + // temperature + temperature_ = getParam<Scalar>("SpatialParams.temperature"); + plotFluidMatrixInteractions_ = getParam<bool>("Output.PlotFluidMatrixInteractions"); } @@ -147,6 +150,16 @@ public: */ auto fluidMatrixInteractionAtPos(const GlobalPosition& globalPos) const { return makeFluidMatrixInteraction(threePhasePcKrSw_, adsorption_); } + + /*! + * \brief Returns the temperature at the domain at the given position + * \param globalPos The position in global coordinates where the temperature should be specified + */ + Scalar temperatureAtPos(const GlobalPosition& globalPos) const + { + return temperature_; + } + private: /* bool isFineMaterial_(const GlobalPosition &globalPos) const @@ -158,6 +171,7 @@ private: Scalar permeability_; Scalar porosity_; + Scalar temperature_; const ThreePhasePcKrSw threePhasePcKrSw_; const AdsorptionModel adsorption_; diff --git a/lecture/mm/remediationscenarios/remediationscenariosexercise.cc b/lecture/mm/remediationscenarios/remediationscenariosexercise.cc index e6c91f142dbc4d9b74fedb00cc3e711204ed0eb9..c53ffd77421d202acd029d259c765a52a8045ff1 100644 --- a/lecture/mm/remediationscenarios/remediationscenariosexercise.cc +++ b/lecture/mm/remediationscenarios/remediationscenariosexercise.cc @@ -111,7 +111,7 @@ int main(int argc, char** argv) try // create the finite volume grid geometry using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; auto fvGridGeometry = std::make_shared<FVGridGeometry>(leafGridView); - fvGridGeometry->update(); + // the problem (initial and boundary conditions) using Problem = GetPropType<TypeTag, Properties::Problem>; diff --git a/lecture/mm/remediationscenarios/remediationscenariosspatialparams.hh b/lecture/mm/remediationscenarios/remediationscenariosspatialparams.hh index 1eaabaa221861fd7682884e9005a39b21a772257..bdebd20fe7ce45afd1525c022dc352d136fb6a02 100644 --- a/lecture/mm/remediationscenarios/remediationscenariosspatialparams.hh +++ b/lecture/mm/remediationscenarios/remediationscenariosspatialparams.hh @@ -26,7 +26,7 @@ #include <dune/common/float_cmp.hh> #include <dumux/porousmediumflow/properties.hh> -#include <dumux/material/spatialparams/fv.hh> +#include <dumux/porousmediumflow/fvspatialparamsmp.hh> #include <dumux/material/fluidmatrixinteractions/3p/parkervangenuchten.hh> namespace Dumux { @@ -36,7 +36,7 @@ namespace Dumux { */ template<class TypeTag> class KuevetteSpatialParams -: public FVSpatialParams<GetPropType<TypeTag, Properties::GridGeometry>, GetPropType<TypeTag, Properties::Scalar>, KuevetteSpatialParams<TypeTag>> +: public FVPorousMediumFlowSpatialParamsMP<GetPropType<TypeTag, Properties::GridGeometry>, GetPropType<TypeTag, Properties::Scalar>, KuevetteSpatialParams<TypeTag>> { using Scalar = GetPropType<TypeTag, Properties::Scalar>; using FVGridGeometry = GetPropType<TypeTag, Properties::GridGeometry>; @@ -45,7 +45,7 @@ class KuevetteSpatialParams using FVElementGeometry = typename FVGridGeometry::LocalView; using SubControlVolume = typename FVElementGeometry::SubControlVolume; using GlobalPosition = typename SubControlVolume::GlobalPosition; - using ParentType = FVSpatialParams<FVGridGeometry, Scalar, KuevetteSpatialParams<TypeTag>>; + using ParentType = FVPorousMediumFlowSpatialParamsMP<FVGridGeometry, Scalar, KuevetteSpatialParams<TypeTag>>; using ThreePhasePcKrSw = FluidMatrix::ParkerVanGenuchten3PDefault<Scalar>; diff --git a/lecture/references/Henry1p2c-reference.vtu b/lecture/references/Henry1p2c-reference.vtu index 728dc6d62858e0d61a796ee6310a627493150083..bdfebab5ae2e74832d0bbaac62b0335a86a511bc 100644 --- a/lecture/references/Henry1p2c-reference.vtu +++ b/lecture/references/Henry1p2c-reference.vtu @@ -2,524 +2,154 @@ <VTKFile type="UnstructuredGrid" version="0.1" byte_order="LittleEndian"> <UnstructuredGrid> <Piece NumberOfCells="800" NumberOfPoints="861"> - <PointData Scalars="P"> + <PointData Scalars="p"> <DataArray type="Float32" Name="p" NumberOfComponents="1" format="ascii"> - 111395 111390 110903 110900 111386 110896 111382 110893 111379 110889 111376 110886 - 111372 110883 111369 110879 111366 110876 111363 110873 111360 110870 111357 110867 - 111354 110865 111352 110862 111350 110860 111348 110857 111346 110855 111344 110854 - 111344 110853 111343 110852 111344 110851 111345 110851 111347 110852 111348 110852 - 111350 110853 111352 110855 111354 110856 111356 110858 111359 110859 111361 110861 - 111363 110863 111366 110865 111368 110867 111370 110868 111372 110870 111375 110872 - 111377 110874 111379 110876 111381 110878 111383 110880 111384 110882 110413 110410 - 110406 110403 110399 110396 110393 110390 110386 110383 110380 110378 110375 110372 - 110369 110367 110365 110363 110361 110360 110358 110358 110357 110357 110357 110357 - 110358 110359 110360 110361 110362 110364 110365 110367 110369 110370 110372 110374 - 110376 110378 110380 109923 109920 109916 109913 109910 109906 109903 109900 109897 - 109894 109890 109888 109885 109882 109879 109877 109874 109872 109870 109868 109866 - 109864 109863 109862 109861 109861 109861 109861 109861 109862 109862 109863 109864 - 109866 109867 109869 109870 109872 109873 109875 109877 109433 109430 109426 109423 - 109420 109416 109413 109410 109407 109404 109401 109397 109395 109392 109389 109386 - 109383 109381 109378 109376 109374 109371 109369 109368 109366 109365 109364 109363 - 109363 109363 109363 109363 109364 109365 109366 109367 109368 109370 109371 109373 - 109374 108943 108940 108936 108933 108930 108927 108923 108920 108917 108914 108910 - 108907 108904 108901 108898 108895 108893 108890 108887 108884 108881 108879 108876 - 108874 108872 108870 108868 108867 108865 108865 108864 108864 108864 108864 108865 - 108865 108866 108868 108869 108870 108871 108453 108450 108447 108443 108440 108437 - 108433 108430 108427 108424 108420 108417 108414 108411 108408 108405 108402 108399 - 108396 108393 108390 108387 108384 108381 108378 108375 108373 108371 108369 108367 - 108366 108365 108364 108364 108364 108364 108365 108366 108367 108368 108369 107963 - 107960 107957 107953 107950 107947 107944 107940 107937 107934 107930 107927 107924 - 107921 107917 107914 107911 107908 107904 107901 107898 107895 107891 107888 107885 - 107882 107879 107876 107873 107871 107869 107867 107866 107865 107864 107864 107864 - 107864 107864 107865 107866 107473 107470 107467 107464 107460 107457 107454 107450 - 107447 107444 107440 107437 107434 107430 107427 107424 107420 107417 107413 107410 - 107406 107403 107399 107396 107392 107388 107385 107382 107378 107375 107372 107370 - 107368 107366 107364 107363 107363 107362 107362 107363 107363 106984 106980 106977 + 111395 111390 110904 110900 111386 110896 111382 110893 111379 110889 111376 110886 + 111372 110883 111369 110880 111366 110877 111363 110874 111361 110871 111358 110868 + 111356 110866 111353 110863 111352 110861 111351 110860 111352 110858 111353 110857 + 111354 110857 111356 110858 111357 110858 111359 110859 111360 110860 111362 110861 + 111364 110862 111365 110863 111366 110864 111368 110866 111369 110867 111371 110868 + 111372 110869 111373 110871 111375 110872 111376 110873 111377 110875 111379 110876 + 111380 110877 111381 110879 111382 110880 111384 110881 111384 110882 110413 110410 + 110406 110403 110400 110396 110393 110390 110387 110384 110381 110378 110375 110373 + 110371 110368 110366 110364 110363 110361 110361 110360 110360 110361 110361 110362 + 110363 110364 110365 110366 110367 110368 110370 110371 110372 110373 110375 110376 + 110377 110378 110380 109923 109920 109916 109913 109910 109906 109903 109900 109897 + 109894 109891 109888 109885 109882 109880 109877 109875 109872 109870 109867 109866 + 109864 109863 109862 109862 109862 109862 109863 109863 109864 109865 109866 109867 + 109868 109869 109871 109872 109873 109874 109876 109877 109433 109430 109427 109423 + 109420 109417 109413 109410 109407 109404 109401 109398 109395 109392 109389 109386 + 109383 109380 109378 109375 109372 109370 109368 109366 109365 109364 109363 109363 + 109363 109363 109364 109364 109365 109366 109367 109368 109369 109371 109372 109373 + 109374 108943 108940 108937 108933 108930 108927 108923 108920 108917 108914 108911 + 108908 108904 108901 108898 108895 108892 108889 108886 108883 108880 108877 108874 + 108872 108869 108867 108866 108865 108864 108863 108863 108863 108864 108864 108865 + 108866 108867 108868 108869 108870 108871 108453 108450 108447 108443 108440 108437 + 108434 108430 108427 108424 108421 108417 108414 108411 108408 108404 108401 108398 + 108395 108391 108388 108385 108381 108378 108375 108373 108370 108368 108366 108365 + 108364 108363 108363 108363 108363 108364 108365 108366 108367 108368 108369 107963 + 107960 107957 107954 107950 107947 107944 107940 107937 107934 107930 107927 107924 + 107920 107917 107914 107910 107907 107903 107900 107896 107893 107889 107886 107882 + 107879 107876 107873 107870 107868 107866 107864 107863 107862 107862 107862 107863 + 107863 107864 107865 107866 107473 107470 107467 107464 107460 107457 107454 107450 + 107447 107444 107440 107437 107434 107430 107427 107423 107420 107416 107412 107409 + 107405 107401 107397 107393 107390 107386 107382 107379 107375 107372 107369 107367 + 107365 107363 107362 107361 107361 107361 107362 107362 107363 106984 106980 106977 106974 106970 106967 106964 106960 106957 106954 106950 106947 106943 106940 106936 - 106933 106929 106926 106922 106919 106915 106911 106907 106903 106899 106896 106892 - 106888 106884 106880 106877 106874 106871 106868 106866 106864 106862 106861 106861 - 106860 106860 106494 106490 106487 106484 106480 106477 106474 106470 106467 106464 - 106460 106457 106453 106450 106446 106442 106439 106435 106431 106428 106424 106420 - 106416 106411 106407 106403 106399 106395 106390 106386 106382 106378 106374 106371 - 106368 106365 106362 106361 106359 106358 106358 106004 106001 105997 105994 105991 - 105987 105984 105980 105977 105974 105970 105967 105963 105959 105956 105952 105948 - 105944 105941 105937 105933 105928 105924 105920 105915 105911 105906 105902 105897 - 105892 105888 105883 105879 105874 105870 105867 105863 105860 105858 105856 105855 + 106933 106929 106925 106921 106918 106914 106910 106906 106902 106897 106893 106889 + 106885 106881 106878 106874 106871 106868 106865 106863 106861 106860 106860 106860 + 106860 106860 106494 106490 106487 106484 106481 106477 106474 106470 106467 106464 + 106460 106457 106453 106449 106446 106442 106438 106435 106431 106427 106423 106418 + 106414 106410 106406 106401 106397 106392 106388 106384 106380 106375 106372 106368 + 106365 106362 106360 106359 106358 106358 106358 106004 106001 105997 105994 105991 + 105987 105984 105980 105977 105974 105970 105966 105963 105959 105955 105952 105948 + 105944 105940 105936 105932 105927 105923 105919 105914 105909 105905 105900 105895 + 105890 105886 105881 105876 105872 105868 105864 105861 105859 105857 105856 105855 105514 105511 105508 105504 105501 105497 105494 105491 105487 105484 105480 105476 - 105473 105469 105465 105462 105458 105454 105450 105446 105442 105437 105433 105428 - 105424 105419 105414 105409 105404 105399 105394 105389 105384 105379 105374 105369 - 105365 105361 105357 105354 105352 105024 105021 105018 105014 105011 105008 105004 + 105473 105469 105465 105461 105457 105453 105449 105445 105441 105437 105432 105427 + 105423 105418 105413 105408 105403 105398 105392 105387 105382 105377 105372 105367 + 105363 105359 105356 105354 105352 105024 105021 105018 105015 105011 105008 105004 105001 104997 104994 104990 104986 104983 104979 104975 104971 104967 104963 104959 - 104955 104951 104946 104942 104937 104932 104927 104922 104917 104912 104906 104901 - 104895 104889 104883 104878 104872 104866 104861 104857 104853 104849 104535 104531 + 104955 104950 104946 104941 104936 104932 104927 104921 104916 104911 104905 104900 + 104894 104888 104882 104876 104871 104865 104860 104856 104852 104849 104535 104532 104528 104525 104521 104518 104514 104511 104507 104504 104500 104496 104493 104489 - 104485 104481 104477 104473 104469 104465 104460 104456 104451 104446 104441 104436 - 104431 104425 104420 104414 104408 104402 104395 104389 104382 104375 104369 104362 - 104356 104351 104347 104045 104042 104038 104035 104032 104028 104025 104021 104018 - 104014 104010 104006 104003 103999 103995 103991 103987 103983 103979 103974 103970 - 103965 103960 103955 103950 103945 103939 103934 103928 103922 103915 103909 103902 + 104485 104481 104477 104473 104468 104464 104460 104455 104450 104446 104441 104435 + 104430 104425 104419 104413 104407 104401 104395 104388 104381 104375 104368 104362 + 104356 104351 104347 104045 104042 104039 104035 104032 104028 104025 104021 104017 + 104014 104010 104006 104003 103999 103995 103991 103987 103983 103978 103974 103969 + 103965 103960 103955 103950 103945 103939 103933 103928 103922 103915 103909 103902 103895 103887 103880 103872 103864 103856 103849 103844 103556 103552 103549 103545 - 103542 103538 103535 103531 103528 103524 103520 103517 103513 103509 103505 103501 - 103497 103493 103488 103484 103479 103475 103470 103465 103460 103454 103449 103443 - 103437 103430 103424 103417 103409 103402 103394 103385 103376 103367 103357 103348 - 103341 103066 103063 103059 103056 103052 103049 103045 103042 103038 103034 103031 - 103027 103023 103019 103015 103011 103007 103003 102998 102994 102989 102984 102979 - 102974 102969 102964 102958 102952 102946 102939 102932 102925 102917 102909 102901 - 102892 102882 102872 102860 102849 102838 102577 102573 102570 102566 102563 102559 - 102556 102552 102548 102545 102541 102537 102533 102529 102525 102521 102517 102513 - 102508 102504 102499 102494 102489 102484 102479 102473 102467 102461 102455 102448 - 102441 102434 102426 102418 102409 102399 102389 102378 102365 102351 102336 102088 - 102084 102080 102077 102073 102070 102066 102062 102059 102055 102051 102047 102043 + 103542 103538 103535 103531 103528 103524 103520 103516 103513 103509 103505 103501 + 103497 103492 103488 103484 103479 103474 103469 103464 103459 103454 103448 103443 + 103437 103430 103424 103417 103410 103402 103394 103385 103376 103367 103357 103348 + 103341 103066 103063 103059 103056 103052 103049 103045 103042 103038 103034 103030 + 103027 103023 103019 103015 103011 103007 103002 102998 102994 102989 102984 102979 + 102974 102969 102963 102958 102952 102946 102939 102932 102925 102918 102910 102901 + 102892 102882 102872 102861 102849 102838 102577 102574 102570 102566 102563 102559 + 102556 102552 102548 102544 102541 102537 102533 102529 102525 102521 102517 102513 + 102508 102504 102499 102494 102489 102484 102479 102473 102468 102462 102455 102449 + 102442 102434 102427 102418 102410 102400 102390 102378 102366 102352 102336 102088 + 102084 102081 102077 102073 102070 102066 102062 102059 102055 102051 102047 102043 102039 102035 102031 102027 102023 102018 102014 102009 102004 101999 101994 101989 - 101983 101977 101971 101965 101958 101951 101943 101935 101927 101918 101908 101897 - 101885 101872 101857 101833 101600 101595 101591 101587 101584 101580 101576 101573 - 101569 101565 101561 101558 101554 101550 101546 101542 101538 101533 101529 101524 - 101520 101515 101510 101504 101499 101493 101487 101481 101475 101468 101461 101453 - 101445 101437 101427 101417 101407 101395 101382 101367 101359 - </DataArray> - <DataArray type="Float32" Name="delp" NumberOfComponents="1" format="ascii"> - 11395.1 11389.8 10903.4 10899.9 11385.9 10896.1 11382.3 10892.6 11378.9 10889.2 11375.5 10885.8 - 11372.2 10882.6 11369 10879.4 11365.9 10876.2 11362.9 10873.2 11360 10870.2 11357.1 10867.4 - 11354.4 10864.6 11351.9 10862 11349.6 10859.6 11347.5 10857.4 11345.8 10855.5 11344.4 10853.8 - 11343.6 10852.5 11343.5 10851.7 11344.1 10851.3 11345.2 10851.3 11346.6 10851.7 11348.3 10852.4 - 11350.2 10853.4 11352.2 10854.6 11354.3 10856 11356.5 10857.5 11358.7 10859.2 11361 10860.9 - 11363.2 10862.7 11365.5 10864.6 11367.8 10866.5 11370.1 10868.5 11372.3 10870.5 11374.5 10872.4 - 11376.7 10874.4 11378.8 10876.4 11380.9 10878.3 11382.8 10880.2 11383.7 10882.5 10413 10409.7 - 10406.2 10402.8 10399.4 10396.1 10392.8 10389.6 10386.5 10383.4 10380.4 10377.5 10374.7 10372 - 10369.5 10367.1 10364.9 10362.9 10361.2 10359.7 10358.5 10357.6 10357 10356.9 10357 10357.4 - 10358 10358.8 10359.9 10361 10362.4 10363.8 10365.3 10367 10368.7 10370.4 10372.2 10374.1 - 10375.9 10377.8 10379.7 9922.94 9919.66 9916.31 9912.95 9909.6 9906.3 9903.03 9899.82 9896.65 - 9893.54 9890.5 9887.54 9884.65 9881.86 9879.18 9876.62 9874.19 9871.9 9869.75 9867.75 9865.91 - 9864.28 9862.92 9861.86 9861.12 9860.7 9860.55 9860.67 9861.02 9861.6 9862.37 9863.32 9864.43 - 9865.69 9867.07 9868.56 9870.14 9871.78 9873.47 9875.2 9876.96 9432.95 9429.68 9426.38 9423.06 - 9419.75 9416.45 9413.19 9409.96 9406.77 9403.62 9400.53 9397.49 9394.52 9391.61 9388.79 9386.04 - 9383.38 9380.81 9378.32 9375.9 9373.58 9371.38 9369.36 9367.56 9366.02 9364.78 9363.83 9363.18 - 9362.82 9362.73 9362.89 9363.29 9363.91 9364.72 9365.72 9366.88 9368.18 9369.58 9371.08 9372.62 - 9374.2 8943.01 8939.75 8936.47 8933.17 8929.87 8926.58 8923.31 8920.06 8916.84 8913.65 8910.5 - 8907.39 8904.33 8901.31 8898.34 8895.42 8892.55 8889.73 8886.94 8884.19 8881.47 8878.83 8876.28 - 8873.88 8871.67 8869.69 8867.97 8866.54 8865.41 8864.58 8864.04 8863.81 8863.84 8864.15 8864.69 - 8865.45 8866.4 8867.51 8868.75 8870.08 8871.44 8453.1 8449.84 8446.57 8443.28 8439.99 8436.7 - 8433.41 8430.14 8426.88 8423.65 8420.44 8417.26 8414.1 8410.97 8407.87 8404.79 8401.73 8398.68 - 8395.64 8392.6 8389.57 8386.56 8383.59 8380.71 8377.94 8375.34 8372.95 8370.79 8368.89 8367.29 - 8366 8365.02 8364.38 8364.05 8364.04 8364.31 8364.85 8365.6 8366.53 8367.57 8368.67 7963.21 - 7959.96 7956.69 7953.4 7950.1 7946.8 7943.5 7940.2 7936.91 7933.63 7930.36 7927.1 7923.85 - 7920.62 7917.39 7914.16 7910.93 7907.69 7904.43 7901.15 7897.84 7894.53 7891.22 7887.94 7884.72 - 7881.6 7878.62 7875.82 7873.23 7870.88 7868.82 7867.07 7865.66 7864.6 7863.9 7863.56 7863.57 - 7863.87 7864.42 7865.12 7865.91 7473.35 7470.09 7466.82 7463.53 7460.22 7456.91 7453.58 7450.25 - 7446.93 7443.6 7440.27 7436.94 7433.6 7430.27 7426.92 7423.56 7420.18 7416.76 7413.31 7409.83 - 7406.29 7402.72 7399.12 7395.51 7391.91 7388.36 7384.88 7381.52 7378.31 7375.28 7372.49 7369.97 - 7367.76 7365.9 7364.41 7363.33 7362.66 7362.39 7362.44 7362.74 7363.15 6983.5 6980.25 6976.97 - 6973.67 6970.35 6967.02 6963.67 6960.31 6956.94 6953.57 6950.18 6946.78 6943.37 6939.94 6936.48 - 6933 6929.48 6925.92 6922.3 6918.63 6914.9 6911.11 6907.26 6903.37 6899.45 6895.53 6891.63 - 6887.78 6884.02 6880.38 6876.92 6873.66 6870.66 6867.98 6865.66 6863.75 6862.28 6861.27 6860.68 - 6860.44 6860.39 6493.68 6490.43 6487.14 6483.83 6480.5 6477.15 6473.77 6470.38 6466.97 6463.55 - 6460.1 6456.64 6453.15 6449.64 6446.09 6442.49 6438.85 6435.16 6431.4 6427.57 6423.66 6419.68 - 6415.61 6411.48 6407.29 6403.05 6398.79 6394.52 6390.28 6386.09 6382.01 6378.06 6374.31 6370.81 - 6367.63 6364.83 6362.47 6360.61 6359.24 6358.3 6357.62 6003.89 6000.63 5997.34 5994.02 5990.67 - 5987.29 5983.89 5980.47 5977.02 5973.55 5970.05 5966.53 5962.97 5959.38 5955.74 5952.05 5948.3 - 5944.49 5940.6 5936.63 5932.57 5928.42 5924.17 5919.83 5915.4 5910.89 5906.31 5901.67 5897.01 - 5892.34 5887.69 5883.11 5878.63 5874.33 5870.28 5866.54 5863.23 5860.41 5858.12 5856.33 5854.86 - 5514.13 5510.87 5507.57 5504.24 5500.87 5497.47 5494.04 5490.58 5487.09 5483.58 5480.03 5476.45 - 5472.83 5469.17 5465.45 5461.68 5457.83 5453.92 5449.92 5445.82 5441.63 5437.33 5432.92 5428.4 - 5423.76 5419.01 5414.15 5409.2 5404.16 5399.06 5393.9 5388.73 5383.58 5378.5 5373.57 5368.87 - 5364.51 5360.62 5357.28 5354.5 5352.1 5024.41 5021.15 5017.84 5014.48 5011.09 5007.67 5004.21 - 5000.72 4997.2 4993.64 4990.05 4986.42 4982.74 4979.01 4975.23 4971.38 4967.45 4963.44 4959.34 - 4955.14 4950.83 4946.41 4941.86 4937.17 4932.36 4927.4 4922.31 4917.07 4911.71 4906.22 4900.61 - 4894.9 4889.12 4883.3 4877.5 4871.8 4866.32 4861.22 4856.68 4852.78 4849.34 4534.74 4531.47 - 4528.15 4524.77 4521.36 4517.9 4514.41 4510.89 4507.33 4503.74 4500.1 4496.43 4492.7 4488.92 - 4485.07 4481.15 4477.16 4473.07 4468.89 4464.6 4460.19 4455.65 4450.98 4446.16 4441.19 4436.05 - 4430.75 4425.28 4419.63 4413.8 4407.8 4401.61 4395.25 4388.74 4382.11 4375.4 4368.73 4362.26 - 4356.27 4351.08 4346.58 4045.12 4041.85 4038.5 4035.1 4031.66 4028.17 4024.65 4021.09 4017.5 - 4013.88 4010.21 4006.49 4002.72 3998.89 3994.99 3991.02 3986.96 3982.81 3978.56 3974.19 3969.69 - 3965.06 3960.28 3955.35 3950.25 3944.96 3939.49 3933.81 3927.92 3921.81 3915.46 3908.86 3902 - 3894.87 3887.47 3879.82 3871.96 3864 3856.25 3849.3 3843.81 3555.58 3552.29 3548.92 3545.48 - 3541.99 3538.47 3534.92 3531.33 3527.71 3524.06 3520.35 3516.6 3512.8 3508.93 3504.99 3500.97 - 3496.86 3492.66 3488.35 3483.91 3479.35 3474.64 3469.78 3464.75 3459.54 3454.13 3448.51 3442.67 - 3436.58 3430.23 3423.59 3416.65 3409.37 3401.72 3393.68 3385.21 3376.3 3366.96 3357.34 3348.04 - 3341.05 3066.13 3062.82 3059.4 3055.9 3052.37 3048.81 3045.22 3041.61 3037.96 3034.28 3030.55 - 3026.78 3022.94 3019.04 3015.07 3011.01 3006.87 3002.62 2998.26 2993.78 2989.16 2984.39 2979.46 - 2974.36 2969.06 2963.56 2957.83 2951.85 2945.61 2939.07 2932.21 2925 2917.39 2909.34 2900.79 - 2891.68 2881.95 2871.51 2860.35 2848.73 2838.29 2576.84 2573.47 2569.94 2566.35 2562.76 2559.16 - 2555.55 2551.92 2548.25 2544.55 2540.8 2537.01 2533.15 2529.23 2525.23 2521.15 2516.98 2512.7 - 2508.31 2503.79 2499.13 2494.32 2489.34 2484.18 2478.82 2473.24 2467.43 2461.36 2455 2448.33 - 2441.3 2433.89 2426.05 2417.71 2408.8 2399.25 2388.93 2377.71 2365.34 2351.27 2335.52 2087.85 - 2084.3 2080.5 2076.81 2073.17 2069.54 2065.9 2062.25 2058.58 2054.86 2051.1 2047.3 2043.43 - 2039.49 2035.48 2031.38 2027.19 2022.89 2018.48 2013.94 2009.25 2004.41 1999.4 1994.21 1988.81 - 1983.19 1977.32 1971.18 1964.75 1957.99 1950.86 1943.33 1935.33 1926.81 1917.69 1907.86 1897.19 - 1885.48 1872.4 1857.27 1832.76 1600.28 1594.95 1590.97 1587.22 1583.55 1579.91 1576.27 1572.62 - 1568.94 1565.22 1561.46 1557.64 1553.77 1549.83 1545.81 1541.71 1537.51 1533.21 1528.78 1524.23 - 1519.54 1514.68 1509.66 1504.45 1499.03 1493.38 1487.49 1481.32 1474.85 1468.04 1460.86 1453.27 - 1445.2 1436.6 1427.38 1417.44 1406.64 1394.77 1381.54 1366.82 1359.16 - </DataArray> - <DataArray type="Float32" Name="x^Water_liq" NumberOfComponents="1" format="ascii"> - 0.999999 0.999998 0.999999 0.999998 0.999998 0.999998 0.999996 0.999997 0.999995 0.999995 0.999993 0.999993 - 0.999989 0.99999 0.999985 0.999985 0.999978 0.999978 0.999968 0.999969 0.999953 0.999955 0.999932 0.999934 - 0.999901 0.999905 0.999856 0.999861 0.999787 0.999796 0.999685 0.9997 0.999529 0.999555 0.999286 0.999332 - 0.998893 0.99897 0.998122 0.998345 0.99706 0.997375 0.996027 0.996297 0.995037 0.995323 0.994115 0.994403 - 0.993265 0.993529 0.992487 0.992727 0.991783 0.991999 0.991151 0.991343 0.990588 0.990754 0.990089 0.990232 - 0.989649 0.989771 0.989265 0.989366 0.988931 0.989014 0.988644 0.98871 0.988398 0.98845 0.988191 0.988229 - 0.988017 0.988044 0.987875 0.987891 0.987762 0.987765 0.987676 0.987663 0.987619 0.987583 0.999999 0.999998 - 0.999998 0.999997 0.999995 0.999993 0.99999 0.999986 0.99998 0.999971 0.999958 0.99994 0.999913 0.999874 - 0.999817 0.999733 0.999609 0.999419 0.999122 0.998651 0.997948 0.997042 0.996046 0.995073 0.994156 0.993351 - 0.992591 0.991878 0.99123 0.990648 0.990128 0.989667 0.989263 0.988911 0.988607 0.988347 0.988128 0.987944 - 0.987794 0.987674 0.987583 0.999999 0.999998 0.999998 0.999997 0.999996 0.999994 0.999991 0.999987 0.999982 - 0.999974 0.999963 0.999946 0.999923 0.99989 0.999842 0.999771 0.999668 0.999513 0.999279 0.998929 0.998424 - 0.997751 0.996936 0.99604 0.995125 0.994233 0.993404 0.992623 0.991943 0.991302 0.990705 0.990166 0.989684 - 0.989258 0.988883 0.988558 0.98828 0.988047 0.987855 0.987701 0.987583 0.999999 0.999999 0.999998 0.999997 - 0.999996 0.999994 0.999992 0.999989 0.999984 0.999977 0.999968 0.999954 0.999934 0.999907 0.999867 0.999809 - 0.999725 0.999602 0.999421 0.99916 0.998793 0.998303 0.997688 0.996965 0.996164 0.995321 0.994469 0.993636 - 0.992841 0.992101 0.991414 0.990824 0.990265 0.98975 0.989287 0.988876 0.988517 0.988209 0.987952 0.987745 - 0.987583 0.999999 0.999999 0.999998 0.999997 0.999996 0.999995 0.999993 0.99999 0.999986 0.99998 0.999972 - 0.999961 0.999945 0.999922 0.999889 0.999843 0.999776 0.99968 0.999542 0.999347 0.999077 0.998719 0.998262 - 0.997708 0.997063 0.996344 0.995572 0.994767 0.99395 0.993143 0.992365 0.991628 0.990942 0.990351 0.989809 - 0.989303 0.988846 0.988443 0.988097 0.987811 0.987583 0.999999 0.999999 0.999998 0.999998 0.999997 0.999996 - 0.999994 0.999992 0.999988 0.999984 0.999977 0.999968 0.999955 0.999936 0.99991 0.999873 0.99982 0.999746 - 0.999641 0.999495 0.999296 0.999033 0.998694 0.998276 0.997775 0.997196 0.996545 0.995833 0.995074 0.994282 - 0.993475 0.992668 0.991881 0.991131 0.990427 0.989816 0.989266 0.988754 0.988296 0.987904 0.987583 0.999999 - 0.999999 0.999999 0.999998 0.999997 0.999996 0.999995 0.999993 0.99999 0.999986 0.999981 0.999973 0.999963 - 0.999948 0.999927 0.999898 0.999858 0.999801 0.999721 0.999612 0.999465 0.99927 0.99902 0.998706 0.998323 - 0.997868 0.997342 0.996746 0.996085 0.995367 0.994599 0.993796 0.992969 0.992136 0.991315 0.990523 0.989778 - 0.989134 0.988556 0.988031 0.987583 0.999999 0.999999 0.999999 0.999998 0.999998 0.999997 0.999996 0.999994 - 0.999992 0.999989 0.999985 0.999978 0.99997 0.999958 0.999942 0.99992 0.999888 0.999845 0.999785 0.999704 - 0.999595 0.999451 0.999265 0.999031 0.998741 0.99839 0.997975 0.997492 0.99694 0.99632 0.995633 0.994885 - 0.994083 0.993235 0.992354 0.991456 0.990561 0.989688 0.988872 0.988189 0.987583 1 0.999999 0.999999 - 0.999999 0.999998 0.999997 0.999997 0.999995 0.999993 0.999991 0.999987 0.999983 0.999976 0.999967 0.999954 - 0.999937 0.999913 0.99988 0.999836 0.999775 0.999695 0.999588 0.999451 0.999276 0.999058 0.998791 0.998468 - 0.998085 0.997637 0.99712 0.996531 0.995867 0.995129 0.994317 0.993436 0.992493 0.991498 0.990468 0.989435 - 0.988441 0.987583 1 0.999999 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999995 0.999993 - 0.99999 0.999986 0.999981 0.999974 0.999964 0.999951 0.999933 0.999908 0.999875 0.99983 0.999771 0.999693 - 0.999591 0.999462 0.999299 0.999097 0.99885 0.998551 0.998195 0.997775 0.997284 0.996715 0.996061 0.995318 - 0.994477 0.993536 0.992492 0.991348 0.990118 0.988843 0.987583 1 1 0.999999 0.999999 0.999999 - 0.999998 0.999998 0.999997 0.999996 0.999994 0.999992 0.999989 0.999985 0.99998 0.999972 0.999962 0.999949 - 0.99993 0.999906 0.999873 0.999829 0.999771 0.999697 0.999602 0.999481 0.99933 0.999143 0.998913 0.998635 - 0.998301 0.997901 0.997426 0.996866 0.996206 0.995432 0.994527 0.993471 0.992246 0.990843 0.989279 0.987583 - 1 1 1 0.999999 0.999999 0.999999 0.999998 0.999998 0.999997 0.999995 0.999994 0.999992 - 0.999989 0.999984 0.999979 0.999971 0.999961 0.999947 0.999929 0.999905 0.999873 0.999831 0.999777 0.999707 - 0.999618 0.999506 0.999366 0.999193 0.998979 0.998718 0.9984 0.998014 0.997545 0.996977 0.996286 0.995442 - 0.994405 0.99313 0.991574 0.989726 0.987583 1 1 1 0.999999 0.999999 0.999999 0.999999 - 0.999998 0.999997 0.999996 0.999995 0.999993 0.999991 0.999988 0.999984 0.999978 0.999971 0.999961 0.999947 - 0.999929 0.999906 0.999876 0.999836 0.999786 0.999721 0.999639 0.999535 0.999406 0.999245 0.999045 0.998798 - 0.998493 0.998114 0.997641 0.997047 0.996288 0.995305 0.994011 0.992321 0.990183 0.987583 1 1 - 1 1 0.999999 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999995 0.999993 0.999991 - 0.999988 0.999983 0.999978 0.99997 0.999961 0.999948 0.999931 0.999909 0.999881 0.999844 0.999797 0.999738 - 0.999663 0.999568 0.999449 0.9993 0.999113 0.998879 0.998585 0.99821 0.997726 0.997086 0.996213 0.994972 - 0.993172 0.990702 0.987583 1 1 1 1 1 0.999999 0.999999 0.999999 0.999998 - 0.999998 0.999997 0.999996 0.999995 0.999993 0.999991 0.999988 0.999983 0.999978 0.999971 0.999961 0.999949 - 0.999934 0.999913 0.999887 0.999854 0.999811 0.999757 0.999689 0.999603 0.999495 0.999359 0.999187 0.998969 - 0.99869 0.998328 0.997844 0.997166 0.996147 0.994462 0.991516 0.987583 1 1 1 1 - 1 0.999999 0.999999 0.999999 0.999999 0.999998 0.999998 0.999997 0.999996 0.999995 0.999993 0.999991 - 0.999987 0.999983 0.999978 0.999972 0.999963 0.999951 0.999937 0.999918 0.999894 0.999864 0.999826 0.999778 - 0.999717 0.999641 0.999545 0.999425 0.999273 0.999082 0.998838 0.998524 0.998108 0.997515 0.996535 0.994395 - 0.987583 1 1 1 1 1 1 0.999999 0.999999 0.999999 0.999999 0.999998 - 0.999998 0.999997 0.999996 0.999994 0.999993 0.99999 0.999987 0.999984 0.999979 0.999972 0.999964 0.999954 - 0.99994 0.999923 0.999902 0.999875 0.999842 0.9998 0.999747 0.999682 0.999601 0.999502 0.99938 0.999233 - 0.999059 0.998861 0.998657 0.998513 0.99862 1 1 1 1 1 1 1 - 1 0.999999 0.999999 0.999999 0.999999 0.999998 0.999997 0.999997 0.999996 0.999994 0.999992 0.99999 - 0.999987 0.999983 0.999978 0.999972 0.999965 0.999955 0.999943 0.999927 0.999908 0.999884 0.999855 0.999819 - 0.999776 0.999723 0.99966 0.999587 0.999504 0.999415 0.999326 0.999257 0.999244 0.999375 1 1 - 1 1 1 1 1 1 0.999999 0.999999 0.999999 0.999999 0.999998 0.999998 - 0.999997 0.999996 0.999995 0.999993 0.999992 0.999989 0.999986 0.999982 0.999977 0.999971 0.999963 0.999954 - 0.999942 0.999927 0.999909 0.999888 0.999862 0.999831 0.999795 0.999754 0.99971 0.999663 0.999619 0.999585 - 0.999574 0.99961 0.999734 1 1 1 1 1 1 1 1 0.999999 - 0.999999 0.999999 0.999999 0.999998 0.999998 0.999997 0.999996 0.999995 0.999994 0.999992 0.99999 0.999987 - 0.999983 0.999979 0.999973 0.999966 0.999957 0.999947 0.999933 0.999918 0.999898 0.999876 0.99985 0.99982 - 0.999786 0.999751 0.999717 0.999688 0.99967 0.999674 0.999711 0.999788 0.999879 - </DataArray> - <DataArray type="Float32" Name="x^Salt_liq" NumberOfComponents="1" format="ascii"> - 1.127e-06 1.63795e-06 1.18651e-06 1.68277e-06 2.39909e-06 2.40463e-06 3.48966e-06 3.45936e-06 5.05674e-06 4.9865e-06 7.31528e-06 7.19181e-06 - 1.05771e-05 1.03758e-05 1.52963e-05 1.49763e-05 2.21374e-05 2.16325e-05 3.20768e-05 3.12798e-05 4.6558e-05 4.52938e-05 6.77297e-05 6.57089e-05 - 9.88203e-05 9.55562e-05 0.000144741 0.000139399 0.000213098 0.000204211 0.000315976 0.000300891 0.000473379 0.000447135 0.00072078 0.000673852 - 0.001132 0.00104323 0.00199344 0.00168423 0.0030488 0.002705 0.00407528 0.00379944 0.00506893 0.00477495 0.00598889 0.00569122 - 0.00683387 0.00656279 0.0076051 0.00736056 0.00830137 0.00808126 0.00892499 0.00873055 0.00948022 0.00931148 0.00997149 0.00982687 - 0.0104033 0.010281 0.0107804 0.0106785 0.0111072 0.0110241 0.0113882 0.0113221 0.0116278 0.0115767 0.0118298 0.011792 - 0.011998 0.0119721 0.0121357 0.0121208 0.0122454 0.012242 0.0123281 0.0123404 0.0123828 0.0124171 1.16028e-06 1.64434e-06 - 2.33581e-06 3.33164e-06 4.76788e-06 6.83628e-06 9.81263e-06 1.40951e-05 2.02599e-05 2.91425e-05 4.19583e-05 6.04827e-05 8.73264e-05 0.000126365 - 0.000183436 0.00026755 0.000393133 0.000584479 0.000883837 0.00135796 0.00207013 0.00299604 0.00402043 0.00502476 0.00593628 0.00673424 - 0.00748746 0.00819479 0.00883641 0.0094107 0.00992323 0.010377 0.010775 0.0111212 0.0114197 0.0116746 0.0118892 0.0120676 - 0.0122134 0.0123293 0.0124171 1.09016e-06 1.54314e-06 2.18726e-06 3.10581e-06 4.41888e-06 6.29713e-06 8.98336e-06 1.28244e-05 1.83172e-05 - 2.6175e-05 3.74252e-05 5.35539e-05 7.67258e-05 0.000110126 0.00015851 0.000229121 0.000333236 0.000488669 0.000722933 0.00107387 0.00157938 - 0.00225431 0.00307431 0.00398457 0.00492374 0.00583635 0.00668717 0.00745698 0.00812934 0.00876102 0.00935048 0.00988293 0.0103572 - 0.0107775 0.0111459 0.011465 0.0117377 0.0119665 0.0121538 0.0123027 0.0124171 9.94037e-07 1.40564e-06 1.98794e-06 2.8127e-06 - 3.98254e-06 5.64326e-06 8.00151e-06 1.13503e-05 1.6106e-05 2.28616e-05 3.24648e-05 4.61328e-05 6.56259e-05 9.35127e-05 0.000133583 0.000191494 - 0.000275747 0.000399022 0.000579462 0.00084047 0.00120659 0.00169503 0.00230771 0.00302971 0.00383437 0.00468951 0.00556147 0.00641838 - 0.00723015 0.00797767 0.00864599 0.00923266 0.00978375 0.0102911 0.0107469 0.011151 0.011504 0.0118059 0.0120573 0.0122596 - 0.0124171 8.84735e-07 1.24985e-06 1.76356e-06 2.48675e-06 3.5056e-06 4.9419e-06 6.96743e-06 9.82448e-06 1.38552e-05 1.9544e-05 2.75788e-05 - 3.89412e-05 5.50396e-05 7.79087e-05 0.000110506 0.000157148 0.000224113 0.000320369 0.000458166 0.000652877 0.00092124 0.00127784 0.00173125 - 0.00228193 0.00292261 0.00364008 0.004417 0.0052337 0.00606874 0.00689784 0.00769471 0.00843973 0.00911044 0.00969589 0.0102316 - 0.0107293 0.011179 0.0115752 0.0119145 0.0121948 0.0124171 7.718e-07 1.08924e-06 1.53342e-06 2.15515e-06 3.02578e-06 4.24545e-06 - 5.95471e-06 8.35079e-06 1.17108e-05 1.64248e-05 2.30435e-05 3.23473e-05 4.54463e-05 6.39256e-05 9.00538e-05 0.000127073 0.000179571 0.000253909 - 0.000358536 0.000503953 0.000701964 0.000964132 0.00129987 0.00171496 0.00221099 0.00278555 0.00343286 0.00414434 0.00490885 0.00571244 - 0.00653816 0.00736555 0.00816911 0.00892581 0.00961591 0.0102193 0.010764 0.0112678 0.0117175 0.0121024 0.0124171 6.62001e-07 - 9.33338e-07 1.31099e-06 1.83674e-06 2.5689e-06 3.58883e-06 5.01006e-06 6.99116e-06 9.75376e-06 1.36082e-05 1.89897e-05 2.65106e-05 3.70338e-05 - 5.17769e-05 7.24568e-05 0.000101482 0.000142193 0.000199117 0.000278183 0.000386742 0.000533281 0.000726748 0.000975612 0.00128697 0.00166598 - 0.00211557 0.0026366 0.00322786 0.00388618 0.00460625 0.00538021 0.00619694 0.00704126 0.0078933 0.00872818 0.00952242 0.0102532 - 0.0108902 0.0114611 0.0119773 0.0124171 5.59711e-07 7.88263e-07 1.10475e-06 1.54312e-06 2.15047e-06 2.99215e-06 4.15892e-06 5.77686e-06 - 8.02129e-06 1.11362e-05 1.54618e-05 2.14724e-05 2.98307e-05 4.14612e-05 5.76498e-05 8.01732e-05 0.000111455 0.000154736 0.000214217 0.000295125 - 0.000403626 0.00054655 0.00073096 0.000963664 0.00125081 0.00159763 0.0020083 0.00248595 0.00303255 0.00364879 0.00433383 0.00508474 - 0.0058958 0.0067574 0.00765505 0.00856916 0.0094743 0.0103464 0.011146 0.0118225 0.0124171 4.67405e-07 6.5746e-07 9.19371e-07 - 1.2804e-06 1.7782e-06 2.46473e-06 3.41182e-06 4.71872e-06 6.52267e-06 9.01354e-06 1.24541e-05 1.72082e-05 2.37791e-05 3.28616e-05 4.5411e-05 - 6.27304e-05 8.65756e-05 0.00011927 0.000163814 0.000223952 0.000304181 0.000409667 0.000546061 0.000719283 0.000935289 0.0011999 0.0015187 - 0.00189697 0.00233963 0.00285121 0.00343567 0.00409616 0.00483453 0.00565062 0.00654116 0.00749844 0.00850886 0.00955143 0.010592 - 0.0115828 0.0124171 3.86156e-07 5.4241e-07 7.56745e-07 1.05084e-06 1.45449e-06 2.00864e-06 2.7696e-06 3.81476e-06 5.25051e-06 7.22311e-06 - 9.93361e-06 1.36581e-05 1.87754e-05 2.58031e-05 3.54457e-05 4.86543e-05 6.66981e-05 9.12451e-05 0.000124442 0.000168983 0.000228147 0.000305803 - 0.000406358 0.000534686 0.000696036 0.000895958 0.00114027 0.00143502 0.00178657 0.00220156 0.00268698 0.00325014 0.00389865 0.00464011 - 0.00548166 0.00642888 0.00748387 0.00864168 0.00988342 0.0111629 0.0124171 3.16083e-07 4.43255e-07 6.16914e-07 8.54141e-07 1.17831e-06 - 1.62142e-06 2.22721e-06 3.0555e-06 4.188e-06 5.73638e-06 7.85299e-06 1.07455e-05 1.46964e-05 2.00885e-05 2.74379e-05 3.74354e-05 5.09958e-05 - 6.9315e-05 9.39301e-05 0.000126778 0.000170243 0.000227194 0.000301 0.000395531 0.000515155 0.000664746 0.00084971 0.00107604 0.00135043 - 0.0016804 0.00207452 0.00254269 0.00309651 0.00374972 0.00451874 0.00542273 0.00648241 0.00771618 0.00913095 0.0107058 0.0124171 - 2.56706e-07 3.59296e-07 4.98769e-07 6.88479e-07 9.46632e-07 1.29802e-06 1.77638e-06 2.42755e-06 3.31379e-06 4.5196e-06 6.15946e-06 8.38827e-06 - 1.1415e-05 1.55208e-05 2.10812e-05 2.85955e-05 3.87206e-05 5.23116e-05 7.04665e-05 9.45732e-05 0.000126357 0.000167923 0.000221801 0.000290978 - 0.000378946 0.000489753 0.000628078 0.000799348 0.0010099 0.00126721 0.00158027 0.0019601 0.00242059 0.00297973 0.00366152 0.00449838 - 0.00553262 0.00681275 0.00838012 0.0102445 0.0124171 2.07198e-07 2.89352e-07 4.00555e-07 5.51189e-07 7.55354e-07 1.03215e-06 1.4074e-06 - 1.91601e-06 2.60505e-06 3.53795e-06 4.8001e-06 6.50616e-06 8.80965e-06 1.19154e-05 1.60954e-05 2.17083e-05 2.92236e-05 3.92503e-05 5.25693e-05 - 7.01703e-05 9.32911e-05 0.000123461 0.000162544 0.000212792 0.000276903 0.000358094 0.000460201 0.000587819 0.000746493 0.000943007 0.00118582 - 0.00148575 0.00185709 0.00231957 0.00290167 0.00364636 0.00461913 0.00591217 0.00761783 0.00977872 0.0124171 1.66553e-07 2.31998e-07 - 3.20217e-07 4.39264e-07 6.0003e-07 8.17165e-07 1.11035e-06 1.506e-06 2.03949e-06 2.75817e-06 3.7253e-06 5.02522e-06 6.77e-06 9.10809e-06 - 1.22351e-05 1.64076e-05 2.19592e-05 2.93213e-05 3.90462e-05 5.18344e-05 6.8566e-05 9.03365e-05 0.000118497 0.000154705 0.000200978 0.000259775 - 0.000334091 0.000427591 0.000544802 0.000691386 0.000874559 0.00110377 0.00139183 0.00175702 0.00222713 0.00284765 0.00369824 0.00492153 - 0.00673844 0.0092532 0.0124171 1.33696e-07 1.85677e-07 2.55591e-07 3.49707e-07 4.76334e-07 6.46727e-07 8.75868e-07 1.18369e-06 1.59671e-06 - 2.15015e-06 2.89072e-06 3.88021e-06 5.20015e-06 6.95763e-06 9.29286e-06 1.23883e-05 1.64802e-05 2.18725e-05 2.8953e-05 3.82138e-05 5.02738e-05 - 6.59069e-05 8.60743e-05 0.000111964 0.000145042 0.000187112 0.000240401 0.000307666 0.000392348 0.000498787 0.000632539 0.000800874 0.00101361 - 0.00128467 0.00163515 0.0021004 0.00274737 0.00372066 0.00536192 0.00839598 0.0124171 1.08124e-07 1.49656e-07 2.05278e-07 2.8e-07 - 3.80474e-07 5.15218e-07 6.95665e-07 9.36912e-07 1.25889e-06 1.68787e-06 2.25842e-06 3.01591e-06 4.01969e-06 5.34712e-06 7.09864e-06 9.40404e-06 - 1.24302e-05 1.63906e-05 2.15566e-05 2.82713e-05 3.69657e-05 4.81788e-05 6.25814e-05 8.10051e-05 0.000104479 0.000134274 0.000171964 0.000219496 - 0.000279294 0.000354387 0.000448588 0.00056674 0.000715085 0.000901842 0.00113826 0.00144096 0.00183825 0.00239162 0.0032847 0.00526134 - 0.0124171 8.88494e-08 1.22567e-07 1.67632e-07 2.28113e-07 3.0927e-07 4.17905e-07 5.62779e-07 7.55489e-07 1.01123e-06 1.3499e-06 1.79744e-06 - 2.38761e-06 3.16416e-06 4.18363e-06 5.51875e-06 7.26272e-06 9.53438e-06 1.24846e-05 1.63038e-05 2.12313e-05 2.75664e-05 3.56815e-05 4.60381e-05 - 5.92059e-05 7.58862e-05 9.69399e-05 0.000123421 0.000156618 0.0001981 0.000249768 0.00031391 0.000393237 0.000490881 0.000610258 0.000754632 - 0.000925864 0.00112099 0.00132247 0.00146886 0.00135855 0 7.49564e-08 1.03155e-07 1.40943e-07 1.91823e-07 2.59998e-07 3.50913e-07 - 4.71637e-07 6.31391e-07 8.422e-07 1.11968e-06 1.48402e-06 1.96127e-06 2.58484e-06 3.39752e-06 4.45378e-06 5.82275e-06 7.59173e-06 9.87045e-06 - 1.27962e-05 1.65397e-05 2.13127e-05 2.73759e-05 3.50488e-05 4.47211e-05 5.68652e-05 7.20496e-05 9.09532e-05 0.000114378 0.000143255 0.000178642 - 0.000221696 0.000273593 0.000335363 0.00040755 0.00048954 0.000578281 0.000665833 0.000734631 0.000747402 0.000617605 0 6.53992e-08 - 9.01373e-08 1.23901e-07 1.69602e-07 2.3058e-07 3.11434e-07 4.18258e-07 5.58936e-07 7.43678e-07 9.85647e-07 1.30174e-06 1.71354e-06 2.24857e-06 - 2.9417e-06 3.83697e-06 4.98974e-06 6.46924e-06 8.36166e-06 1.07737e-05 1.38368e-05 1.77117e-05 2.25942e-05 2.87203e-05 3.63727e-05 4.58861e-05 - 5.76515e-05 7.21168e-05 8.9782e-05 0.000111182 0.000136851 0.000167251 0.000202655 0.000242935 0.000287228 0.000333401 0.000377216 0.000411048 - 0.000421768 0.000386612 0.000263634 0 5.83715e-08 8.27107e-08 1.16694e-07 1.61483e-07 2.20489e-07 2.98245e-07 4.0059e-07 5.35027e-07 - 7.11196e-07 9.41464e-07 1.24164e-06 1.63185e-06 2.13766e-06 2.79133e-06 3.63345e-06 4.71481e-06 6.09867e-06 7.8634e-06 1.01055e-05 1.29433e-05 - 1.65205e-05 2.10109e-05 2.66221e-05 3.36004e-05 4.22329e-05 5.28489e-05 6.58155e-05 8.15263e-05 0.000100376 0.000122711 0.000148753 0.00017846 - 0.000211315 0.000246019 0.000280039 0.000309044 0.000326285 0.000322293 0.000285948 0.00021 0.000119861 - </DataArray> - <DataArray type="Float32" Name="X^Water_liq" NumberOfComponents="1" format="ascii"> - 0.999996 0.999995 0.999996 0.999995 0.999992 0.999992 0.999989 0.999989 0.999984 0.999984 0.999976 0.999977 - 0.999966 0.999966 0.99995 0.999951 0.999928 0.99993 0.999896 0.999898 0.999849 0.999853 0.99978 0.999787 - 0.99968 0.99969 0.999531 0.999548 0.99931 0.999339 0.998978 0.999028 0.998472 0.998558 0.997685 0.997835 - 0.996416 0.996664 0.993927 0.994648 0.990517 0.991526 0.987215 0.988077 0.984064 0.984974 0.981143 0.982053 - 0.978459 0.979291 0.976012 0.976766 0.973805 0.974483 0.971831 0.972428 0.970075 0.970594 0.968522 0.968967 - 0.967158 0.967535 0.965968 0.966282 0.964936 0.965192 0.964048 0.964253 0.96329 0.963449 0.96265 0.962768 - 0.962117 0.962198 0.961679 0.961726 0.96133 0.96134 0.961065 0.961026 0.96089 0.96078 0.999996 0.999995 - 0.999992 0.999989 0.999985 0.999978 0.999968 0.999954 0.999934 0.999905 0.999864 0.999804 0.999717 0.999591 - 0.999406 0.999135 0.99873 0.998115 0.997155 0.995635 0.993368 0.99046 0.987275 0.984178 0.981272 0.978732 - 0.97634 0.974103 0.972077 0.970261 0.968643 0.967213 0.96596 0.964872 0.963934 0.963133 0.962456 0.961892 - 0.961429 0.961061 0.96078 0.999996 0.999995 0.999993 0.99999 0.999986 0.99998 0.999971 0.999958 0.999941 - 0.999915 0.999879 0.999826 0.999751 0.999643 0.999486 0.999258 0.998922 0.99842 0.997664 0.996532 0.994902 - 0.992734 0.990121 0.987257 0.984344 0.981515 0.978897 0.976441 0.974305 0.972302 0.970439 0.968761 0.967268 - 0.965945 0.964786 0.963783 0.962926 0.962207 0.961615 0.961144 0.96078 0.999997 0.999995 0.999994 0.999991 - 0.999987 0.999982 0.999974 0.999963 0.999948 0.999926 0.999895 0.99985 0.999787 0.999697 0.999567 0.999379 - 0.999107 0.998708 0.998124 0.997278 0.996092 0.994511 0.992532 0.990212 0.987653 0.984967 0.982262 0.979629 - 0.977124 0.974803 0.972652 0.970812 0.969069 0.967469 0.966036 0.964765 0.963655 0.962706 0.961916 0.961278 - 0.96078 0.999997 0.999996 0.999994 0.999992 0.999989 0.999984 0.999977 0.999968 0.999955 0.999936 0.99991 - 0.999874 0.999821 0.999747 0.999641 0.99949 0.999273 0.998961 0.998514 0.997882 0.997011 0.995853 0.994381 - 0.992595 0.990527 0.988228 0.985764 0.983207 0.980622 0.978075 0.975631 0.97332 0.971177 0.969339 0.967655 - 0.966086 0.964672 0.963427 0.96236 0.96148 0.96078 0.999997 0.999996 0.999995 0.999993 0.99999 0.999986 - 0.999981 0.999973 0.999962 0.999947 0.999925 0.999895 0.999852 0.999792 0.999708 0.999587 0.999417 0.999176 - 0.998836 0.998363 0.997719 0.996866 0.995774 0.994423 0.992812 0.990952 0.988868 0.986597 0.984182 0.981672 - 0.979121 0.97658 0.974114 0.971767 0.969576 0.967676 0.965972 0.964387 0.962975 0.961767 0.96078 0.999998 - 0.999997 0.999996 0.999994 0.999992 0.999988 0.999984 0.999977 0.999968 0.999956 0.999938 0.999914 0.99988 - 0.999832 0.999765 0.99967 0.999538 0.999353 0.999096 0.998743 0.998265 0.997635 0.996824 0.99581 0.994575 - 0.993113 0.991422 0.989513 0.987402 0.985112 0.982676 0.980133 0.977528 0.97491 0.972341 0.969872 0.967558 - 0.965562 0.963778 0.962157 0.96078 0.999998 0.999997 0.999996 0.999995 0.999993 0.99999 0.999986 0.999981 - 0.999974 0.999964 0.99995 0.99993 0.999903 0.999865 0.999812 0.999739 0.999637 0.999497 0.999303 0.99904 - 0.998686 0.99822 0.997619 0.996859 0.995923 0.994792 0.993454 0.991902 0.990133 0.98815 0.985961 0.983583 - 0.98104 0.978364 0.975595 0.972782 0.969991 0.967277 0.964754 0.962646 0.96078 0.999998 0.999998 0.999997 - 0.999996 0.999994 0.999992 0.999989 0.999985 0.999979 0.999971 0.999959 0.999944 0.999923 0.999893 0.999852 - 0.999796 0.999718 0.999612 0.999467 0.999271 0.999009 0.998665 0.99822 0.997654 0.996949 0.996084 0.995044 - 0.99381 0.99237 0.990711 0.988825 0.986706 0.984357 0.981783 0.979 0.976032 0.972915 0.969703 0.966495 - 0.963422 0.96078 0.999999 0.999998 0.999998 0.999997 0.999995 0.999993 0.999991 0.999988 0.999983 0.999976 - 0.999968 0.999955 0.999939 0.999916 0.999885 0.999842 0.999783 0.999703 0.999594 0.999449 0.999256 0.999003 - 0.998674 0.998255 0.997727 0.997074 0.996275 0.995312 0.994164 0.992812 0.991234 0.989412 0.987325 0.984956 - 0.982289 0.979315 0.976029 0.972445 0.968614 0.964663 0.96078 0.999999 0.999999 0.999998 0.999997 0.999996 - 0.999995 0.999993 0.99999 0.999986 0.999981 0.999974 0.999965 0.999952 0.999934 0.999911 0.999878 0.999834 - 0.999774 0.999694 0.999586 0.999444 0.999258 0.999017 0.998708 0.998317 0.997827 0.997222 0.996481 0.995583 - 0.994504 0.993217 0.991692 0.989895 0.987785 0.985319 0.982446 0.979108 0.975257 0.970871 0.966011 0.96078 - 0.999999 0.999999 0.999998 0.999998 0.999997 0.999996 0.999994 0.999992 0.999989 0.999985 0.99998 0.999973 - 0.999963 0.999949 0.999931 0.999907 0.999874 0.999829 0.99977 0.999691 0.999587 0.999452 0.999275 0.999049 - 0.998761 0.998398 0.997945 0.997383 0.996693 0.995849 0.994823 0.99358 0.992074 0.990253 0.988043 0.985352 - 0.98206 0.978034 0.973153 0.967396 0.96078 0.999999 0.999999 0.999999 0.999998 0.999998 0.999997 0.999995 - 0.999994 0.999991 0.999988 0.999984 0.999979 0.999971 0.999961 0.999947 0.999929 0.999905 0.999872 0.999828 - 0.999771 0.999695 0.999596 0.999469 0.999304 0.999094 0.998828 0.998493 0.998074 0.997553 0.996907 0.996109 - 0.995123 0.993902 0.992382 0.990475 0.988049 0.984915 0.980815 0.975492 0.968816 0.96078 0.999999 0.999999 - 0.999999 0.999999 0.999998 0.999997 0.999996 0.999995 0.999993 0.999991 0.999988 0.999984 0.999978 0.99997 - 0.99996 0.999946 0.999928 0.999904 0.999872 0.999831 0.999776 0.999705 0.999612 0.999494 0.999342 0.999149 - 0.998905 0.998598 0.998213 0.99773 0.997127 0.996371 0.995419 0.994211 0.992654 0.990602 0.987807 0.983857 - 0.978165 0.970429 0.96078 1 0.999999 0.999999 0.999999 0.999998 0.999998 0.999997 0.999996 0.999995 - 0.999993 0.999991 0.999987 0.999983 0.999977 0.99997 0.999959 0.999946 0.999928 0.999905 0.999875 0.999835 - 0.999784 0.999718 0.999633 0.999525 0.999387 0.999212 0.998991 0.998712 0.998362 0.99792 0.997364 0.99666 - 0.99576 0.994591 0.993033 0.990858 0.987597 0.982242 0.972971 0.96078 1 1 0.999999 0.999999 - 0.999999 0.999998 0.999998 0.999997 0.999996 0.999994 0.999993 0.99999 0.999987 0.999982 0.999977 0.999969 - 0.999959 0.999946 0.999929 0.999907 0.999879 0.999842 0.999795 0.999735 0.999658 0.99956 0.999436 0.99928 - 0.999083 0.998836 0.998525 0.998135 0.997644 0.997025 0.996238 0.995225 0.993882 0.991977 0.988837 0.982029 - 0.96078 1 1 0.999999 0.999999 0.999999 0.999999 0.999998 0.999998 0.999997 0.999996 0.999994 - 0.999992 0.99999 0.999986 0.999982 0.999976 0.999969 0.999959 0.999947 0.99993 0.99991 0.999883 0.999849 - 0.999806 0.999751 0.999682 0.999595 0.999486 0.999349 0.999179 0.998968 0.998707 0.998385 0.997991 0.997515 - 0.996951 0.99631 0.995654 0.995187 0.995532 1 1 1 1 0.999999 0.999999 0.999999 - 0.999998 0.999998 0.999997 0.999996 0.999995 0.999994 0.999992 0.999989 0.999985 0.999981 0.999975 0.999968 - 0.999958 0.999946 0.99993 0.99991 0.999885 0.999853 0.999813 0.999764 0.999701 0.999624 0.99953 0.999413 - 0.999272 0.999101 0.998898 0.998661 0.998392 0.998102 0.997815 0.997591 0.997549 0.997973 1 1 - 1 1 0.999999 0.999999 0.999999 0.999999 0.999998 0.999998 0.999997 0.999996 0.999994 0.999993 - 0.99999 0.999987 0.999984 0.999979 0.999973 0.999965 0.999955 0.999942 0.999926 0.999906 0.999881 0.999849 - 0.999811 0.999763 0.999705 0.999635 0.999551 0.999451 0.999335 0.999203 0.999058 0.998906 0.998763 0.998653 - 0.998618 0.998733 0.999137 1 1 1 1 0.999999 0.999999 0.999999 0.999999 0.999998 - 0.999998 0.999997 0.999996 0.999995 0.999993 0.999991 0.999988 0.999985 0.99998 0.999974 0.999967 0.999958 - 0.999946 0.999931 0.999913 0.99989 0.999861 0.999826 0.999784 0.999732 0.99967 0.999597 0.999512 0.999414 - 0.999307 0.999193 0.999082 0.998987 0.998931 0.998944 0.999063 0.999312 0.999608 - </DataArray> - <DataArray type="Float32" Name="X^Salt_liq" NumberOfComponents="1" format="ascii"> - 3.659e-06 5.31786e-06 3.8522e-06 5.46336e-06 7.78902e-06 7.807e-06 1.13297e-05 1.12313e-05 1.64174e-05 1.61893e-05 2.37499e-05 2.3349e-05 - 3.43395e-05 3.3686e-05 4.96604e-05 4.86215e-05 7.18691e-05 7.02301e-05 0.000104135 0.000101548 0.000151142 0.000147039 0.000219862 0.000213303 - 0.000320765 0.000310173 0.000469774 0.000452441 0.000691526 0.000662702 0.00102514 0.000976232 0.00153527 0.00145024 0.00233635 0.00218447 - 0.00366592 0.00337909 0.00644318 0.00544751 0.0098311 0.0087292 0.013111 0.0122311 0.0162718 0.0153381 0.0191858 0.0182442 - 0.0218518 0.0209976 0.0242764 0.0235085 0.0264583 0.0257693 0.0284069 0.0277999 0.0301372 0.0296118 0.0316647 0.0312154 - 0.0330047 0.0326253 0.0341726 0.0338574 0.0351833 0.0349266 0.0360514 0.0358472 0.0367904 0.0366329 0.0374131 0.0372966 - 0.0379311 0.0378513 0.0383548 0.038309 0.0386921 0.0386819 0.0389464 0.0389844 0.0391147 0.03922 3.76702e-06 5.33862e-06 - 7.58355e-06 1.08166e-05 1.54795e-05 2.21948e-05 3.18576e-05 4.57607e-05 6.57743e-05 9.46099e-05 0.000136212 0.00019634 0.000283464 0.000410147 - 0.00059531 0.000868123 0.00127524 0.00189512 0.00286384 0.00439543 0.00668991 0.00966211 0.0129361 0.0161316 0.0190195 0.021538 - 0.0239071 0.0261248 0.0281304 0.0299208 0.0315148 0.0329232 0.034156 0.0352265 0.0361485 0.0369347 0.0375961 0.0381452 - 0.0385939 0.0389502 0.03922 3.53937e-06 5.01004e-06 7.10128e-06 1.00834e-05 1.43465e-05 2.04444e-05 2.91654e-05 4.16354e-05 5.94674e-05 - 8.49765e-05 0.000121497 0.000173851 0.00024906 0.000357454 0.000514447 0.000743496 0.0010811 0.0015848 0.00234332 0.00347812 0.00510958 - 0.00728212 0.00991278 0.0128218 0.0158108 0.0187034 0.0213897 0.0238114 0.0259199 0.027895 0.0297333 0.0313896 0.0328618 - 0.0341637 0.035303 0.0362882 0.0371293 0.0378342 0.0384106 0.0388685 0.03922 3.2273e-06 4.56362e-06 6.45416e-06 9.13185e-06 - 1.29299e-05 1.83216e-05 2.59778e-05 3.68497e-05 5.22889e-05 7.42202e-05 0.000105395 0.000149762 0.000213034 0.000303541 0.00043357 0.000621449 - 0.000894703 0.00129433 0.00187887 0.00272358 0.0039068 0.00548232 0.00745373 0.00976994 0.0123426 0.0150665 0.0178334 0.0205421 - 0.0230987 0.0254448 0.0275358 0.0293662 0.0310814 0.0326569 0.0340691 0.0353188 0.0364086 0.0373394 0.0381137 0.038736 - 0.03922 2.87243e-06 4.05782e-06 5.72567e-06 8.0736e-06 1.13814e-05 1.60445e-05 2.26206e-05 3.18961e-05 4.49818e-05 6.345e-05 8.95335e-05 - 0.000126418 0.000178673 0.000252899 0.000358687 0.000510027 0.000727255 0.00103938 0.00148598 0.00211657 0.00298478 0.00413685 0.00559901 - 0.00737087 0.00942685 0.0117223 0.0141996 0.0167946 0.0194382 0.0220532 0.0245576 0.0268911 0.0289853 0.0308082 0.0324721 - 0.0340145 0.0354052 0.0366282 0.037674 0.0385366 0.03922 2.50577e-06 3.53639e-06 4.9785e-06 6.99704e-06 9.82363e-06 1.37834e-05 - 1.93327e-05 2.71117e-05 3.802e-05 5.33237e-05 7.48107e-05 0.000105013 0.000147534 0.000207515 0.000292316 0.000412444 0.000582774 0.000823887 - 0.00116311 0.00163432 0.00227546 0.00312345 0.00420795 0.00554654 0.00714286 0.00898751 0.0110601 0.0133312 0.0157636 0.0183114 - 0.0209199 0.0235242 0.0260444 0.0284094 0.0305595 0.0324341 0.034122 0.0356794 0.037067 0.0382525 0.03922 2.14929e-06 - 3.03023e-06 4.25633e-06 5.96327e-06 8.34033e-06 1.16516e-05 1.62658e-05 2.26976e-05 3.16665e-05 4.41799e-05 6.16507e-05 8.60661e-05 0.000120226 - 0.000168083 0.000235205 0.000329404 0.000461505 0.000646179 0.000902603 0.00125453 0.00172932 0.00235566 0.00316056 0.00416633 0.0053887 - 0.00683607 0.00850975 0.0104043 0.0125079 0.0148018 0.0172591 0.0198431 0.0225046 0.0251804 0.0277925 0.0302686 0.0325391 - 0.0345125 0.0362764 0.0378673 0.03922 1.81719e-06 2.55922e-06 3.58675e-06 5.00998e-06 6.98181e-06 9.71444e-06 1.35025e-05 1.87553e-05 - 2.6042e-05 3.61547e-05 5.01974e-05 6.97104e-05 9.6844e-05 0.000134598 0.000187146 0.000260249 0.000361767 0.000502201 0.000695155 0.000957537 - 0.00130925 0.00177229 0.00236929 0.00312194 0.00404958 0.00516841 0.006491 0.00802623 0.00977904 0.0117501 0.0139348 0.016322 - 0.0188915 0.0216109 0.0244332 0.0272957 0.0301188 0.0328282 0.0353032 0.0373907 0.03922 1.51751e-06 2.13455e-06 2.98488e-06 - 4.15703e-06 5.77318e-06 8.00211e-06 1.1077e-05 1.532e-05 2.11766e-05 2.92634e-05 4.04332e-05 5.58671e-05 7.71986e-05 0.000106683 0.000147419 - 0.000203636 0.000281027 0.000387128 0.000531653 0.00072673 0.000986901 0.00132883 0.00177071 0.0023315 0.0030302 0.00388521 0.00491395 - 0.00613268 0.00755627 0.009198 0.011069 0.0131776 0.0155275 0.0181157 0.0209294 0.0239416 0.0271072 0.0303588 0.0335894 - 0.0366518 0.03922 1.25372e-06 1.76102e-06 2.45689e-06 3.41172e-06 4.72222e-06 6.52136e-06 8.99192e-06 1.23851e-05 1.70464e-05 2.34506e-05 - 3.22504e-05 4.43421e-05 6.09548e-05 8.37691e-05 0.000115071 0.000157947 0.000216514 0.000296182 0.000403909 0.000548422 0.000740339 0.000992159 - 0.00131811 0.00173386 0.00225627 0.00290303 0.0036926 0.00464407 0.00577721 0.00711254 0.00867137 0.0104756 0.0125477 0.0149095 - 0.0175806 0.0205753 0.0238959 0.0275223 0.0313911 0.0353554 0.03922 1.02622e-06 1.4391e-06 2.00291e-06 2.77311e-06 3.82558e-06 - 5.2642e-06 7.23098e-06 9.92011e-06 1.35969e-05 1.86239e-05 2.54956e-05 3.48863e-05 4.77128e-05 6.52178e-05 8.90764e-05 0.00012153 0.000165547 - 0.000225008 0.000304895 0.000411487 0.00055251 0.000737248 0.000976587 0.00128302 0.0016706 0.00215499 0.00275347 0.00348512 0.00437113 - 0.00543517 0.00670403 0.00820838 0.00998387 0.0120724 0.0145234 0.0173939 0.0207441 0.024625 0.0290492 0.0339419 0.03922 - 8.33437e-07 1.16651e-06 1.61934e-06 2.23526e-06 3.07339e-06 4.21423e-06 5.76729e-06 7.8814e-06 1.07587e-05 1.46735e-05 1.99974e-05 2.72334e-05 - 3.70599e-05 5.0389e-05 6.84404e-05 9.28341e-05 0.000125702 0.000169818 0.000228745 0.000306983 0.000410122 0.000544986 0.000719756 0.000944092 - 0.00122927 0.00158832 0.00203629 0.00259056 0.00327137 0.00410252 0.00511245 0.00633589 0.00781633 0.00960984 0.0117907 0.0144586 - 0.017742 0.0217853 0.0267047 0.0325123 0.03922 6.72702e-07 9.39427e-07 1.30047e-06 1.78953e-06 2.45238e-06 3.35103e-06 4.56936e-06 - 6.22063e-06 8.45767e-06 1.14865e-05 1.55842e-05 2.1123e-05 2.86014e-05 3.86843e-05 5.22544e-05 7.04761e-05 9.48732e-05 0.000127422 0.000170655 - 0.000227784 0.000302822 0.000400724 0.000527533 0.000690535 0.000898452 0.00116168 0.00149258 0.00190593 0.00241955 0.00305516 0.00383974 - 0.00480769 0.00600432 0.00749183 0.00935975 0.0117423 0.0148428 0.0189432 0.0243164 0.0310658 0.03922 5.40741e-07 7.5322e-07 - 1.03964e-06 1.42614e-06 1.94809e-06 2.65306e-06 3.60494e-06 4.88947e-06 6.62152e-06 8.95481e-06 1.20947e-05 1.6315e-05 2.19796e-05 2.95703e-05 - 3.97223e-05 5.32679e-05 7.12907e-05 9.51904e-05 0.000126759 0.000168269 0.000222577 0.000293233 0.000384619 0.0005021 0.000652214 0.000842912 - 0.00108387 0.00138691 0.00176663 0.00224122 0.00283383 0.0035747 0.00450471 0.00568203 0.00719476 0.0091866 0.011908 0.0158038 - 0.0215512 0.0294302 0.03922 4.34065e-07 6.02833e-07 8.29818e-07 1.13538e-06 1.54649e-06 2.0997e-06 2.84365e-06 3.84304e-06 5.18397e-06 - 6.98078e-06 9.38513e-06 1.25976e-05 1.68829e-05 2.25888e-05 3.01702e-05 4.02196e-05 5.35038e-05 7.10091e-05 9.39946e-05 0.000124057 0.000163204 - 0.000213946 0.0002794 0.000363419 0.00047075 0.000607236 0.000780081 0.000998199 0.0012727 0.00161758 0.00205073 0.0025955 0.00328339 - 0.00415888 0.00528935 0.00678726 0.00886507 0.0119796 0.0172011 0.0267543 0.03922 3.51044e-07 4.85882e-07 6.66469e-07 9.09066e-07 - 1.23527e-06 1.67274e-06 2.25859e-06 3.04183e-06 4.08718e-06 5.47992e-06 7.3323e-06 9.7916e-06 1.30505e-05 1.73601e-05 2.30465e-05 3.05311e-05 - 4.03557e-05 5.3213e-05 6.99838e-05 9.17816e-05 0.000120005 0.000156404 0.000203152 0.000262949 0.000339128 0.000435812 0.000558093 0.000712278 - 0.000906205 0.00114966 0.00145495 0.00183768 0.00231792 0.00292206 0.00368614 0.00466324 0.00594363 0.00772329 0.0105862 0.0168823 - 0.03922 2.88464e-07 3.97935e-07 5.44244e-07 7.40608e-07 1.0041e-06 1.3568e-06 1.82715e-06 2.45282e-06 3.28313e-06 4.38267e-06 5.83568e-06 - 7.75173e-06 1.02729e-05 1.35827e-05 1.79173e-05 2.35792e-05 3.09543e-05 4.05321e-05 5.2931e-05 6.89278e-05 8.94935e-05 0.000115837 0.000149455 - 0.000192196 0.000246335 0.000314663 0.000400597 0.000508308 0.000642878 0.000810458 0.00101844 0.00127558 0.00159197 0.00197859 0.00244589 - 0.00299973 0.00363034 0.0042809 0.00475322 0.00439733 0 2.43358e-07 3.3491e-07 4.57594e-07 6.22784e-07 8.44126e-07 1.1393e-06 - 1.53125e-06 2.04991e-06 2.73434e-06 3.63521e-06 4.81811e-06 6.36755e-06 8.39207e-06 1.10305e-05 1.44598e-05 1.89043e-05 2.46474e-05 3.20454e-05 - 4.15437e-05 5.3697e-05 6.9192e-05 8.88749e-05 0.000113783 0.00014518 0.000184599 0.000233883 0.000295234 0.000371251 0.000464951 0.00057976 - 0.000719416 0.00088772 0.00108799 0.00132197 0.00158763 0.00187505 0.00215851 0.00238117 0.0024225 0.00200238 0 2.12329e-07 - 2.92646e-07 4.02266e-07 5.50642e-07 7.48615e-07 1.01112e-06 1.35794e-06 1.81468e-06 2.41447e-06 3.20006e-06 4.22629e-06 5.56328e-06 7.30033e-06 - 9.55067e-06 1.24573e-05 1.61998e-05 2.10032e-05 2.7147e-05 3.49778e-05 4.4922e-05 5.75018e-05 7.33521e-05 9.32392e-05 0.00011808 0.000148962 - 0.000187151 0.000234101 0.000291434 0.000360881 0.000444172 0.000542806 0.000657655 0.0007883 0.000931933 0.00108163 0.00122366 0.0013333 - 0.00136804 0.00125411 0.000855425 0 1.89513e-07 2.68534e-07 3.78866e-07 5.24281e-07 7.15854e-07 9.68301e-07 1.30058e-06 1.73705e-06 - 2.30901e-06 3.05661e-06 4.03117e-06 5.29806e-06 6.94024e-06 9.06246e-06 1.17965e-05 1.53073e-05 1.98001e-05 2.55294e-05 3.28086e-05 4.20214e-05 - 5.36347e-05 6.8212e-05 8.6428e-05 0.000109081 0.000137103 0.000171562 0.00021365 0.00026464 0.000325812 0.000398293 0.000482791 0.000579166 - 0.000685744 0.0007983 0.000908623 0.00100267 0.00105856 0.00104562 0.000927783 0.000681477 0.000389043 + 101983 101978 101971 101965 101958 101951 101944 101936 101928 101919 101909 101898 + 101886 101873 101858 101833 101600 101595 101591 101587 101584 101580 101576 101573 + 101569 101565 101561 101558 101554 101550 101546 101542 101537 101533 101529 101524 + 101519 101515 101510 101504 101499 101494 101488 101482 101475 101469 101462 101454 + 101446 101438 101429 101419 101408 101396 101382 101367 101360 </DataArray> <DataArray type="Float32" Name="rho" NumberOfComponents="1" format="ascii"> - 998.213 998.214 998.213 998.214 998.216 998.216 998.219 998.218 998.222 998.222 998.227 998.226 - 998.234 998.233 998.244 998.243 998.259 998.258 998.281 998.279 998.312 998.309 998.358 998.353 - 998.425 998.418 998.524 998.513 998.672 998.653 998.894 998.861 999.233 999.175 999.76 999.66 - 1000.61 1000.44 1002.28 1001.8 1004.58 1003.9 1006.82 1006.23 1008.96 1008.34 1010.95 1010.33 - 1012.78 1012.21 1014.46 1013.94 1015.98 1015.51 1017.34 1016.92 1018.55 1018.19 1019.62 1019.31 - 1020.57 1020.3 1021.39 1021.17 1022.11 1021.93 1022.72 1022.58 1023.25 1023.14 1023.7 1023.61 - 1024.07 1024.01 1024.37 1024.34 1024.62 1024.61 1024.8 1024.83 1024.92 1025 998.213 998.214 - 998.216 998.218 998.221 998.225 998.232 998.241 998.254 998.274 998.301 998.342 998.4 998.484 - 998.608 998.789 999.06 999.472 1000.12 1001.14 1002.66 1004.62 1006.78 1008.88 1010.86 1012.6 - 1014.23 1015.77 1017.17 1018.42 1019.54 1020.53 1021.4 1022.15 1022.8 1023.36 1023.83 1024.23 - 1024.55 1024.8 1025 998.213 998.214 998.215 998.217 998.22 998.224 998.23 998.238 998.25 - 998.267 998.291 998.326 998.377 998.449 998.554 998.706 998.931 999.267 999.773 1000.53 1001.63 - 1003.09 1004.85 1006.79 1008.77 1010.69 1012.48 1014.16 1015.63 1017.01 1018.3 1019.46 1020.49 - 1021.41 1022.21 1022.91 1023.5 1024.01 1024.42 1024.75 1025 998.212 998.213 998.214 998.216 - 998.219 998.222 998.227 998.235 998.245 998.26 998.281 998.31 998.352 998.413 998.5 998.625 - 998.807 999.074 999.465 1000.03 1000.83 1001.89 1003.22 1004.79 1006.52 1008.34 1010.18 1011.98 - 1013.7 1015.29 1016.77 1018.04 1019.24 1020.35 1021.34 1022.23 1023 1023.66 1024.21 1024.65 - 1025 998.212 998.213 998.214 998.215 998.217 998.221 998.225 998.231 998.24 998.252 998.27 - 998.294 998.329 998.379 998.45 998.551 998.696 998.905 999.204 999.627 1000.21 1000.99 1001.98 - 1003.18 1004.57 1006.13 1007.8 1009.54 1011.3 1013.04 1014.72 1016.31 1017.79 1019.06 1020.22 - 1021.31 1022.29 1023.16 1023.9 1024.51 1025 998.211 998.212 998.213 998.214 998.216 998.219 - 998.223 998.228 998.235 998.245 998.26 998.28 998.308 998.348 998.405 998.485 998.599 998.761 - 998.988 999.305 999.736 1000.31 1001.04 1001.95 1003.03 1004.29 1005.7 1007.24 1008.88 1010.59 - 1012.33 1014.07 1015.76 1017.38 1018.89 1020.21 1021.39 1022.49 1023.47 1024.31 1025 998.211 - 998.211 998.212 998.213 998.215 998.217 998.22 998.225 998.231 998.239 998.251 998.267 998.29 - 998.322 998.367 998.43 998.518 998.642 998.814 999.05 999.37 999.792 1000.34 1001.02 1001.85 - 1002.83 1003.97 1005.26 1006.69 1008.24 1009.9 1011.64 1013.42 1015.21 1016.98 1018.69 1020.29 - 1021.67 1022.91 1024.04 1025 998.21 998.211 998.212 998.213 998.214 998.216 998.218 998.222 - 998.227 998.233 998.243 998.256 998.274 998.299 998.334 998.383 998.451 998.546 998.675 998.851 - 999.088 999.4 999.803 1000.31 1000.94 1001.7 1002.6 1003.65 1004.84 1006.18 1007.67 1009.28 - 1011.02 1012.85 1014.74 1016.68 1018.6 1020.48 1022.23 1023.7 1025 998.21 998.21 998.211 - 998.212 998.213 998.214 998.216 998.219 998.223 998.229 998.236 998.246 998.261 998.28 998.308 - 998.345 998.397 998.469 998.566 998.697 998.872 999.102 999.4 999.779 1000.25 1000.83 1001.53 - 1002.36 1003.33 1004.45 1005.73 1007.16 1008.76 1010.51 1012.41 1014.44 1016.59 1018.8 1021.02 - 1023.16 1025 998.21 998.21 998.21 998.211 998.212 998.213 998.215 998.217 998.22 998.224 - 998.23 998.238 998.25 998.265 998.286 998.315 998.354 998.407 998.48 998.577 998.706 998.876 - 999.096 999.376 999.73 1000.17 1000.7 1001.35 1002.12 1003.03 1004.1 1005.33 1006.74 1008.35 - 1010.16 1012.19 1014.45 1016.91 1019.56 1022.3 1025 998.209 998.209 998.21 998.21 998.211 - 998.212 998.213 998.215 998.218 998.221 998.226 998.232 998.24 998.252 998.268 998.29 998.32 - 998.36 998.413 998.485 998.58 998.704 998.866 999.073 999.335 999.662 1000.07 1000.57 1001.17 - 1001.89 1002.76 1003.79 1005 1006.43 1008.1 1010.06 1012.34 1014.98 1018 1021.36 1025 - 998.209 998.209 998.209 998.21 998.21 998.211 998.212 998.214 998.215 998.218 998.222 998.227 - 998.233 998.242 998.254 998.271 998.293 998.322 998.362 998.415 998.484 998.575 998.693 998.844 - 999.037 999.28 999.583 999.96 1000.42 1000.99 1001.68 1002.52 1003.53 1004.76 1006.25 1008.08 - 1010.32 1013.07 1016.42 1020.4 1025 998.208 998.209 998.209 998.209 998.21 998.21 998.211 - 998.212 998.214 998.216 998.219 998.222 998.227 998.234 998.243 998.255 998.272 998.294 998.323 - 998.361 998.412 998.478 998.563 998.673 998.814 998.992 999.216 999.497 999.846 1000.28 1000.81 - 1001.48 1002.3 1003.32 1004.61 1006.25 1008.38 1011.17 1014.81 1019.42 1025 998.208 998.208 - 998.208 998.209 998.209 998.21 998.21 998.211 998.212 998.214 998.216 998.219 998.223 998.228 - 998.234 998.244 998.256 998.272 998.293 998.321 998.358 998.405 998.467 998.546 998.648 998.777 - 998.94 999.146 999.404 999.727 1000.13 1000.64 1001.28 1002.09 1003.14 1004.52 1006.41 1009.1 - 1012.98 1018.3 1025 998.208 998.208 998.208 998.208 998.209 998.209 998.209 998.21 998.211 - 998.212 998.214 998.216 998.219 998.223 998.228 998.235 998.244 998.255 998.271 998.291 998.318 - 998.352 998.396 998.453 998.525 998.618 998.735 998.883 999.069 999.304 999.599 999.972 1000.44 - 1001.05 1001.83 1002.88 1004.35 1006.56 1010.19 1016.55 1025 998.208 998.208 998.208 998.208 - 998.208 998.208 998.209 998.209 998.21 998.211 998.212 998.214 998.216 998.219 998.223 998.228 - 998.234 998.243 998.254 998.269 998.288 998.313 998.344 998.385 998.436 998.502 998.584 998.689 - 998.821 998.986 999.194 999.455 999.784 1000.2 1000.73 1001.41 1002.31 1003.59 1005.72 1010.34 - 1025 998.207 998.207 998.207 998.208 998.208 998.208 998.208 998.209 998.209 998.21 998.211 - 998.212 998.214 998.216 998.219 998.223 998.228 998.234 998.243 998.254 998.267 998.285 998.308 - 998.337 998.373 998.42 998.478 998.551 998.642 998.756 998.897 999.072 999.287 999.551 999.87 - 1000.25 1000.68 1001.12 1001.43 1001.2 998.207 998.207 998.207 998.207 998.207 998.207 998.208 - 998.208 998.208 998.209 998.209 998.21 998.211 998.213 998.214 998.217 998.22 998.224 998.229 - 998.235 998.243 998.254 998.267 998.284 998.305 998.332 998.365 998.406 998.458 998.521 998.599 - 998.694 998.808 998.944 999.102 999.282 999.477 999.669 999.819 999.847 999.563 998.207 998.207 - 998.207 998.207 998.207 998.207 998.207 998.208 998.208 998.208 998.209 998.21 998.21 998.212 - 998.213 998.215 998.218 998.221 998.225 998.23 998.237 998.245 998.256 998.27 998.286 998.307 - 998.333 998.365 998.404 998.451 998.507 998.574 998.651 998.74 998.837 998.938 999.034 999.108 - 999.131 999.054 998.784 998.207 998.207 998.207 998.207 998.207 998.207 998.207 998.207 998.208 - 998.208 998.208 998.209 998.21 998.211 998.213 998.214 998.217 998.22 998.224 998.229 998.235 - 998.243 998.253 998.265 998.28 998.299 998.322 998.351 998.385 998.427 998.476 998.533 998.598 - 998.67 998.746 998.821 998.884 998.922 998.913 998.833 998.667 998.469 + 998.198 998.199 998.198 998.199 998.201 998.2 998.203 998.202 998.207 998.206 998.213 998.212 + 998.223 998.221 998.238 998.234 998.263 998.256 998.302 998.289 998.366 998.343 998.471 998.427 + 998.648 998.564 998.959 998.791 999.544 999.198 1000.97 1000.1 1007.94 1002.62 1011.9 1007.08 + 1015.16 1010.39 1017.67 1013.15 1019.48 1015.85 1020.75 1018.16 1021.64 1020.05 1022.33 1021.46 + 1022.91 1022.43 1023.38 1023.07 1023.76 1023.5 1024.04 1023.82 1024.26 1024.08 1024.43 1024.29 + 1024.56 1024.45 1024.67 1024.58 1024.75 1024.69 1024.81 1024.77 1024.86 1024.83 1024.9 1024.88 + 1024.93 1024.91 1024.95 1024.94 1024.96 1024.96 1024.97 1024.98 1024.98 1024.99 998.197 998.198 + 998.199 998.201 998.204 998.209 998.216 998.226 998.242 998.267 998.305 998.364 998.458 998.612 + 998.876 999.363 1000.31 1002.06 1004.78 1007.94 1010.95 1013.41 1015.2 1017.23 1018.92 1020.43 + 1021.68 1022.61 1023.26 1023.7 1024.02 1024.26 1024.45 1024.6 1024.71 1024.79 1024.86 1024.91 + 1024.94 1024.97 1024.99 998.197 998.197 998.198 998.2 998.202 998.205 998.21 998.218 998.229 + 998.247 998.274 998.315 998.381 998.486 998.66 998.953 999.447 1000.25 1001.51 1003.27 1005.45 + 1007.96 1010.59 1013.07 1015.24 1016.84 1018.53 1019.87 1021.1 1022.12 1022.9 1023.48 1023.89 + 1024.2 1024.43 1024.61 1024.74 1024.83 1024.9 1024.95 1024.99 998.196 998.197 998.197 998.198 + 998.2 998.202 998.206 998.211 998.22 998.232 998.251 998.279 998.324 998.395 998.509 998.692 + 998.983 999.434 1000.1 1001.05 1002.32 1003.92 1005.85 1008.02 1010.35 1012.72 1014.95 1016.92 + 1018.5 1019.76 1020.86 1021.91 1022.76 1023.4 1023.87 1024.23 1024.49 1024.69 1024.82 1024.92 + 1024.98 998.196 998.196 998.196 998.197 998.198 998.2 998.203 998.207 998.212 998.221 998.234 + 998.254 998.284 998.332 998.406 998.521 998.7 998.971 999.366 999.924 1000.68 1001.67 1002.91 + 1004.42 1006.19 1008.19 1010.35 1012.61 1014.84 1016.9 1018.67 1020.14 1021.06 1022.04 1022.89 + 1023.54 1024.04 1024.41 1024.68 1024.86 1024.98 998.195 998.195 998.196 998.196 998.197 998.198 + 998.2 998.203 998.207 998.213 998.222 998.235 998.256 998.287 998.336 998.41 998.522 998.689 + 998.931 999.274 999.744 1000.37 1001.17 1002.18 1003.41 1004.87 1006.56 1008.48 1010.58 1012.8 + 1015.05 1017.17 1019.06 1020.59 1021.52 1022.47 1023.31 1023.95 1024.42 1024.76 1024.98 998.195 + 998.195 998.195 998.196 998.196 998.197 998.198 998.2 998.203 998.207 998.213 998.222 998.236 + 998.257 998.288 998.336 998.407 998.512 998.664 998.879 999.176 999.576 1000.1 1000.78 1001.62 + 1002.65 1003.89 1005.35 1007.04 1008.94 1011.04 1013.29 1015.59 1017.79 1019.76 1021.3 1022.29 + 1023.24 1024.03 1024.59 1024.98 998.194 998.195 998.195 998.195 998.195 998.196 998.197 998.198 + 998.2 998.203 998.207 998.213 998.222 998.236 998.257 998.287 998.333 998.399 998.495 998.631 + 998.82 999.078 999.422 999.872 1000.45 1001.17 1002.06 1003.13 1004.42 1005.92 1007.66 1009.62 + 1011.8 1014.14 1016.54 1018.85 1020.96 1022.44 1023.46 1024.36 1024.98 998.194 998.194 998.194 + 998.195 998.195 998.195 998.196 998.197 998.198 998.2 998.203 998.207 998.213 998.222 998.235 + 998.255 998.284 998.326 998.387 998.474 998.595 998.761 998.986 999.284 999.673 1000.17 1000.8 + 1001.58 1002.53 1003.68 1005.05 1006.65 1008.51 1010.62 1012.97 1015.5 1018.11 1020.62 1023 + 1024.16 1024.98 998.194 998.194 998.194 998.194 998.194 998.195 998.195 998.196 998.196 998.198 + 998.2 998.202 998.206 998.212 998.221 998.234 998.252 998.279 998.318 998.373 998.451 998.558 + 998.704 998.9 999.16 999.5 999.937 1000.49 1001.19 1002.04 1003.09 1004.36 1005.88 1007.67 + 1009.76 1012.15 1014.83 1017.72 1020.71 1023.58 1024.98 998.194 998.194 998.194 998.194 998.194 + 998.194 998.194 998.195 998.195 998.196 998.197 998.199 998.202 998.206 998.211 998.22 998.232 + 998.249 998.274 998.309 998.358 998.427 998.521 998.65 998.822 999.051 999.35 999.738 1000.23 + 1000.86 1001.65 1002.63 1003.83 1005.31 1007.1 1009.25 1011.81 1014.78 1018.13 1021.68 1024.98 + 998.193 998.193 998.193 998.194 998.194 998.194 998.194 998.194 998.195 998.195 998.196 998.197 + 998.199 998.201 998.205 998.21 998.218 998.229 998.245 998.267 998.299 998.343 998.403 998.487 + 998.6 998.752 998.954 999.22 999.568 1000.02 1000.6 1001.33 1002.27 1003.45 1004.95 1006.84 + 1009.22 1012.19 1015.85 1020.16 1024.98 998.193 998.193 998.193 998.193 998.193 998.193 998.193 + 998.194 998.194 998.194 998.195 998.196 998.197 998.198 998.201 998.204 998.209 998.216 998.226 + 998.24 998.26 998.288 998.327 998.381 998.454 998.554 998.688 998.868 999.107 999.421 999.835 + 1000.37 1001.08 1002 1003.21 1004.81 1006.96 1009.89 1013.81 1018.85 1024.98 998.193 998.193 + 998.193 998.193 998.193 998.193 998.193 998.193 998.193 998.194 998.194 998.195 998.195 998.196 + 998.198 998.2 998.203 998.208 998.214 998.223 998.236 998.253 998.278 998.312 998.36 998.424 + 998.512 998.631 998.791 999.005 999.292 999.673 1000.18 1000.86 1001.79 1003.06 1004.88 1007.58 + 1011.71 1017.56 1024.98 998.193 998.193 998.193 998.193 998.193 998.193 998.193 998.193 998.193 + 998.193 998.193 998.194 998.194 998.195 998.196 998.197 998.199 998.202 998.206 998.212 998.22 + 998.231 998.247 998.268 998.298 998.34 998.396 998.473 998.577 998.719 998.91 999.167 999.514 + 999.985 1000.63 1001.54 1002.86 1004.95 1008.62 1015.63 1024.98 998.192 998.192 998.192 998.193 + 998.193 998.193 998.193 998.193 998.193 998.193 998.193 998.193 998.193 998.194 998.195 998.195 + 998.197 998.199 998.201 998.205 998.21 998.217 998.227 998.24 998.259 998.285 998.32 998.369 + 998.435 998.525 998.647 998.811 999.032 999.33 999.732 1000.28 1001.04 1002.17 1004.11 1008.63 + 1024.98 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.193 + 998.193 998.193 998.193 998.194 998.194 998.195 998.196 998.198 998.2 998.203 998.208 998.214 + 998.222 998.234 998.25 998.272 998.302 998.343 998.397 998.471 998.57 998.701 998.873 999.096 + 999.379 999.725 1000.11 1000.42 1000.31 998.192 998.192 998.192 998.192 998.192 998.192 998.192 + 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.193 998.193 998.193 998.194 998.195 + 998.196 998.197 998.199 998.202 998.206 998.211 998.219 998.229 998.242 998.261 998.285 998.318 + 998.361 998.417 998.489 998.581 998.695 998.83 998.98 999.124 999.209 999.094 998.192 998.192 + 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.192 + 998.192 998.192 998.193 998.193 998.194 998.194 998.196 998.197 998.199 998.202 998.205 998.21 + 998.217 998.226 998.238 998.253 998.273 998.299 998.331 998.371 998.418 998.473 998.531 998.585 + 998.621 998.61 998.494 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.192 + 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.192 998.193 998.193 998.194 998.195 + 998.196 998.198 998.2 998.203 998.207 998.213 998.22 998.23 998.242 998.258 998.278 998.302 + 998.33 998.363 998.398 998.431 998.456 998.463 998.441 998.379 998.282 </DataArray> <DataArray type="Float32" Name="mu" NumberOfComponents="1" format="ascii"> 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 @@ -595,6 +225,376 @@ 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 0.001 </DataArray> + <DataArray type="Float32" Name="delp" NumberOfComponents="1" format="ascii"> + 11395.2 11389.9 10903.5 10900 11386 10896.3 11382.4 10892.8 11379 10889.3 11375.7 10886 + 11372.5 10882.8 11369.3 10879.6 11366.3 10876.6 11363.3 10873.6 11360.5 10870.8 11357.9 10868.1 + 11355.5 10865.5 11353.4 10863.2 11351.8 10861.3 11351 10859.6 11351.6 10858.2 11352.8 10857.5 + 11354.2 10857.3 11355.7 10857.5 11357.3 10858.1 11358.9 10858.8 11360.5 10859.7 11362 10860.8 + 11363.5 10861.9 11365 10863.1 11366.5 10864.4 11367.9 10865.6 11369.3 10866.9 11370.6 10868.2 + 11372 10869.5 11373.4 10870.8 11374.7 10872.1 11376 10873.4 11377.4 10874.7 11378.7 10876 + 11380 10877.3 11381.3 10878.5 11382.5 10879.8 11383.7 10881 11384.2 10882.5 10413.2 10409.8 + 10406.4 10403 10399.6 10396.3 10393 10389.9 10386.8 10383.8 10380.9 10378.1 10375.4 10372.9 + 10370.5 10368.3 10366.1 10364.1 10362.6 10361.5 10360.8 10360.4 10360.4 10360.7 10361.2 10361.9 + 10362.8 10363.8 10364.8 10365.9 10367.1 10368.3 10369.5 10370.8 10372 10373.3 10374.6 10375.8 + 10377.1 10378.4 10379.7 9923.07 9919.79 9916.45 9913.1 9909.76 9906.47 9903.23 9900.04 9896.92 + 9893.85 9890.87 9887.96 9885.14 9882.41 9879.75 9877.14 9874.56 9872.03 9869.63 9867.46 9865.64 + 9864.2 9863.17 9862.47 9862.12 9862.07 9862.29 9862.73 9863.36 9864.14 9865.05 9866.07 9867.15 + 9868.29 9869.48 9870.69 9871.92 9873.17 9874.42 9875.69 9876.96 9433.08 9429.81 9426.51 9423.2 + 9419.9 9416.62 9413.37 9410.16 9406.99 9403.87 9400.8 9397.79 9394.83 9391.91 9389.03 9386.15 + 9383.28 9380.41 9377.57 9374.83 9372.26 9369.92 9367.87 9366.16 9364.82 9363.85 9363.22 9362.92 + 9362.93 9363.19 9363.67 9364.33 9365.13 9366.06 9367.1 9368.2 9369.35 9370.54 9371.75 9372.98 + 9374.2 8943.13 8939.87 8936.59 8933.3 8930.01 8926.73 8923.46 8920.23 8917.02 8913.84 8910.69 + 8907.58 8904.49 8901.42 8898.35 8895.28 8892.18 8889.06 8885.94 8882.84 8879.81 8876.89 8874.14 + 8871.62 8869.36 8867.42 8865.83 8864.6 8863.76 8863.29 8863.14 8863.28 8863.67 8864.26 8865.01 + 8865.91 8866.92 8868 8869.13 8870.28 8871.44 8453.21 8449.95 8446.68 8443.4 8440.11 8436.82 + 8433.54 8430.27 8427.01 8423.78 8420.56 8417.35 8414.15 8410.94 8407.73 8404.49 8401.22 8397.91 + 8394.57 8391.21 8387.86 8384.56 8381.34 8378.24 8375.3 8372.59 8370.13 8367.98 8366.19 8364.77 + 8363.75 8363.14 8362.9 8362.99 8363.36 8363.95 8364.7 8365.6 8366.59 8367.63 8368.67 7963.31 + 7960.06 7956.79 7953.5 7950.2 7946.9 7943.6 7940.3 7937 7933.71 7930.41 7927.12 7923.82 + 7920.5 7917.16 7913.78 7910.36 7906.9 7903.38 7899.83 7896.25 7892.66 7889.09 7885.58 7882.15 + 7878.85 7875.71 7872.8 7870.15 7867.82 7865.85 7864.29 7863.16 7862.48 7862.24 7862.35 7862.76 + 7863.39 7864.16 7865.02 7865.91 7473.43 7470.18 7466.9 7463.61 7460.31 7456.99 7453.66 7450.32 + 7446.98 7443.63 7440.27 7436.9 7433.51 7430.1 7426.65 7423.15 7419.61 7416 7412.34 7408.63 + 7404.86 7401.05 7397.23 7393.4 7389.59 7385.84 7382.19 7378.66 7375.31 7372.19 7369.35 7366.86 + 7364.77 7363.13 7361.98 7361.34 7361.18 7361.39 7361.88 7362.49 7363.15 6983.58 6980.32 6977.04 + 6973.74 6970.42 6967.08 6963.73 6960.36 6956.97 6953.57 6950.15 6946.71 6943.24 6939.73 6936.19 + 6932.59 6928.94 6925.22 6921.43 6917.58 6913.65 6909.67 6905.63 6901.56 6897.46 6893.36 6889.28 + 6885.27 6881.35 6877.57 6873.99 6870.67 6867.66 6865.05 6862.9 6861.28 6860.23 6859.74 6859.73 + 6860.01 6860.39 6493.75 6490.49 6487.2 6483.89 6480.55 6477.19 6473.81 6470.4 6466.98 6463.53 + 6460.05 6456.55 6453 6449.42 6445.79 6442.1 6438.35 6434.53 6430.63 6426.65 6422.6 6418.46 + 6414.25 6409.97 6405.64 6401.25 6396.84 6392.43 6388.04 6383.71 6379.5 6375.44 6371.62 6368.09 + 6364.96 6362.31 6360.23 6358.77 6357.93 6357.64 6357.62 6003.94 6000.69 5997.39 5994.07 5990.71 + 5987.33 5983.91 5980.47 5977.01 5973.51 5969.98 5966.42 5962.81 5959.16 5955.45 5951.68 5947.84 + 5943.93 5939.94 5935.86 5931.68 5927.42 5923.06 5918.61 5914.07 5909.45 5904.76 5900.01 5895.23 + 5890.44 5885.67 5880.98 5876.42 5872.05 5867.98 5864.31 5861.15 5858.62 5856.77 5855.58 5854.86 + 5514.18 5510.91 5507.61 5504.27 5500.9 5497.49 5494.04 5490.57 5487.07 5483.53 5479.95 5476.33 + 5472.67 5468.95 5465.18 5461.34 5457.42 5453.43 5449.35 5445.17 5440.9 5436.52 5432.03 5427.43 + 5422.72 5417.9 5412.97 5407.95 5402.83 5397.64 5392.4 5387.14 5381.9 5376.76 5371.79 5367.09 + 5362.83 5359.14 5356.15 5353.87 5352.1 5024.45 5021.18 5017.87 5014.51 5011.11 5007.68 5004.21 + 5000.7 4997.16 4993.58 4989.96 4986.3 4982.58 4978.81 4974.98 4971.07 4967.09 4963.02 4958.86 + 4954.61 4950.24 4945.76 4941.16 4936.43 4931.57 4926.58 4921.45 4916.18 4910.78 4905.24 4899.58 + 4893.82 4887.98 4882.11 4876.27 4870.56 4865.12 4860.15 4855.85 4852.32 4849.34 4534.77 4531.5 + 4528.17 4524.79 4521.37 4517.91 4514.4 4510.87 4507.29 4503.67 4500.02 4496.31 4492.55 4488.73 + 4484.84 4480.89 4476.84 4472.72 4468.49 4464.16 4459.71 4455.14 4450.44 4445.61 4440.62 4435.48 + 4430.17 4424.7 4419.04 4413.21 4407.19 4400.99 4394.61 4388.08 4381.42 4374.7 4368.03 4361.62 + 4355.77 4350.81 4346.58 4045.15 4041.87 4038.53 4035.12 4031.67 4028.17 4024.64 4021.07 4017.46 + 4013.81 4010.12 4006.38 4002.58 3998.72 3994.79 3990.78 3986.69 3982.51 3978.22 3973.83 3969.31 + 3964.67 3959.89 3954.95 3949.86 3944.59 3939.14 3933.49 3927.63 3921.54 3915.21 3908.63 3901.78 + 3894.66 3887.27 3879.61 3871.74 3863.78 3856.05 3849.21 3843.81 3555.61 3552.32 3548.94 3545.5 + 3542 3538.47 3534.91 3531.31 3527.67 3524 3520.27 3516.5 3512.67 3508.77 3504.81 3500.76 + 3496.63 3492.4 3488.07 3483.62 3479.05 3474.35 3469.49 3464.48 3459.3 3453.93 3448.35 3442.56 + 3436.52 3430.22 3423.64 3416.75 3409.51 3401.9 3393.88 3385.42 3376.5 3367.12 3357.43 3348.03 + 3341.05 3066.16 3062.85 3059.42 3055.92 3052.38 3048.81 3045.21 3041.58 3037.92 3034.22 3030.48 + 3026.68 3022.82 3018.9 3014.91 3010.83 3006.67 3002.4 2998.03 2993.55 2988.93 2984.18 2979.27 + 2974.19 2968.94 2963.48 2957.81 2951.9 2945.73 2939.27 2932.49 2925.35 2917.81 2909.81 2901.31 + 2892.23 2882.48 2871.99 2860.71 2848.88 2838.29 2576.87 2573.5 2569.97 2566.37 2562.78 2559.17 + 2555.54 2551.9 2548.22 2544.5 2540.73 2536.92 2533.04 2529.1 2525.09 2520.99 2516.8 2512.51 + 2508.12 2503.6 2498.95 2494.16 2489.21 2484.09 2478.78 2473.27 2467.53 2461.54 2455.27 2448.69 + 2441.76 2434.45 2426.69 2418.42 2409.57 2400.04 2389.71 2378.41 2365.91 2351.61 2335.52 2087.89 + 2084.33 2080.53 2076.83 2073.18 2069.55 2065.9 2062.24 2058.55 2054.82 2051.04 2047.22 2043.33 + 2039.38 2035.35 2031.24 2027.04 2022.73 2018.32 2013.79 2009.12 2004.3 1999.33 1994.18 1988.84 + 1983.29 1977.51 1971.46 1965.13 1958.48 1951.47 1944.04 1936.15 1927.72 1918.66 1908.85 1898.16 + 1886.35 1873.1 1857.71 1832.76 1600.32 1594.99 1591 1587.25 1583.58 1579.93 1576.28 1572.61 + 1568.91 1565.18 1561.4 1557.57 1553.68 1549.73 1545.69 1541.58 1537.37 1533.06 1528.65 1524.1 + 1519.43 1514.61 1509.63 1504.47 1499.12 1493.55 1487.75 1481.69 1475.33 1468.65 1461.6 1454.13 + 1446.18 1437.69 1428.55 1418.64 1407.82 1395.84 1382.4 1367.38 1359.51 + </DataArray> + <DataArray type="Float32" Name="x^Water_liq" NumberOfComponents="1" format="ascii"> + 0.999999 0.999999 0.999999 0.999999 0.999998 0.999998 0.999997 0.999997 0.999995 0.999995 0.999992 0.999993 + 0.999988 0.999989 0.999981 0.999982 0.999969 0.999972 0.999951 0.999957 0.999922 0.999932 0.999873 0.999893 + 0.999792 0.999831 0.999648 0.999726 0.999379 0.999538 0.998722 0.999122 0.9955 0.997959 0.993669 0.995897 + 0.992157 0.994365 0.99099 0.993087 0.990149 0.991833 0.989556 0.990761 0.989141 0.98988 0.988819 0.989227 + 0.988551 0.988773 0.98833 0.988478 0.988157 0.988277 0.988024 0.988127 0.987921 0.988007 0.987842 0.987909 + 0.98778 0.987832 0.987732 0.98777 0.987694 0.987723 0.987664 0.987685 0.987641 0.987656 0.987624 0.987634 + 0.98761 0.987617 0.987601 0.987604 0.987593 0.987595 0.987588 0.987588 0.987584 0.987583 0.999999 0.999999 + 0.999998 0.999997 0.999996 0.999994 0.999991 0.999986 0.999979 0.999967 0.99995 0.999922 0.999879 0.999808 + 0.999686 0.999462 0.999025 0.998219 0.996961 0.995499 0.994107 0.992968 0.992136 0.991193 0.990408 0.989704 + 0.989125 0.988689 0.988388 0.988182 0.988033 0.98792 0.987832 0.987764 0.987712 0.987672 0.987642 0.987619 + 0.987603 0.987591 0.987583 0.999999 0.999999 0.999999 0.999998 0.999997 0.999995 0.999993 0.99999 0.999984 + 0.999976 0.999964 0.999945 0.999915 0.999866 0.999786 0.999651 0.999423 0.999051 0.998473 0.997661 0.99665 + 0.995493 0.994274 0.993125 0.992116 0.991373 0.990589 0.989968 0.989395 0.988917 0.988553 0.988287 0.988093 + 0.987948 0.98784 0.987759 0.987698 0.987653 0.987621 0.987598 0.987583 1 0.999999 0.999999 0.999999 + 0.999998 0.999997 0.999995 0.999992 0.999989 0.999983 0.999974 0.999961 0.99994 0.999908 0.999855 0.999771 + 0.999637 0.999429 0.99912 0.998683 0.998098 0.997357 0.996466 0.995462 0.994386 0.993285 0.992254 0.991339 + 0.990604 0.990017 0.989503 0.989017 0.98862 0.988324 0.988101 0.987935 0.987811 0.987722 0.987658 0.987613 + 0.987583 1 1 0.999999 0.999999 0.999998 0.999998 0.999996 0.999995 0.999992 0.999988 0.999982 + 0.999973 0.999959 0.999937 0.999903 0.999849 0.999767 0.999642 0.99946 0.999203 0.998854 0.998398 0.997824 + 0.997128 0.996312 0.995386 0.994383 0.993339 0.992301 0.991347 0.990522 0.98984 0.98941 0.988954 0.988559 + 0.988255 0.988023 0.98785 0.987727 0.987641 0.987583 1 1 0.999999 0.999999 0.999999 0.999998 + 0.999997 0.999996 0.999994 0.999992 0.999988 0.999981 0.999972 0.999957 0.999935 0.999901 0.999849 0.999772 + 0.99966 0.999502 0.999286 0.998998 0.998627 0.998162 0.997594 0.996919 0.996137 0.995252 0.994279 0.993247 + 0.992206 0.99122 0.990344 0.98963 0.989197 0.988753 0.988365 0.988067 0.987846 0.987689 0.987583 1 + 1 1 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 0.999991 0.999987 0.999981 + 0.999971 0.999957 0.999935 0.999902 0.999854 0.999784 0.999685 0.999548 0.999363 0.999121 0.99881 0.998421 + 0.997944 0.997371 0.996697 0.995918 0.995036 0.994063 0.993023 0.991956 0.990933 0.990019 0.9893 0.988841 + 0.988396 0.988029 0.987766 0.987583 1 1 1 1 0.999999 0.999999 0.999999 0.999998 + 0.999997 0.999996 0.999994 0.999991 0.999987 0.999981 0.999971 0.999957 0.999936 0.999906 0.999861 0.999799 + 0.999711 0.999593 0.999434 0.999227 0.998962 0.998629 0.998218 0.997721 0.997129 0.996434 0.995631 0.994722 + 0.993713 0.992629 0.991511 0.99044 0.989459 0.988769 0.988294 0.987873 0.987583 1 1 1 + 1 1 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 0.999991 0.999987 0.999981 + 0.999972 0.999958 0.999939 0.999911 0.999871 0.999815 0.999739 0.999635 0.999498 0.999318 0.999089 0.998799 + 0.99844 0.998 0.997469 0.996838 0.996096 0.995237 0.99426 0.993172 0.991998 0.990786 0.989615 0.988508 + 0.987969 0.987583 1 1 1 1 1 1 0.999999 0.999999 0.999999 0.999998 + 0.999997 0.999996 0.999994 0.999991 0.999987 0.999982 0.999973 0.999961 0.999943 0.999917 0.999882 0.999832 + 0.999765 0.999674 0.999554 0.999398 0.999196 0.998941 0.998621 0.998225 0.99774 0.997155 0.996454 0.995626 + 0.994659 0.99355 0.992307 0.990963 0.989574 0.988236 0.987583 1 1 1 1 1 + 1 1 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 0.999992 0.999988 0.999982 + 0.999974 0.999963 0.999947 0.999924 0.999893 0.999849 0.99979 0.99971 0.999605 0.999467 0.999288 0.999059 + 0.99877 0.998406 0.997955 0.997398 0.996717 0.995889 0.994893 0.993709 0.992329 0.990774 0.989123 0.987583 + 1 1 1 1 1 1 1 1 0.999999 0.999999 0.999999 0.999998 + 0.999997 0.999996 0.999995 0.999992 0.999989 0.999984 0.999976 0.999966 0.999951 0.999931 0.999903 0.999865 + 0.999813 0.999743 0.99965 0.999527 0.999367 0.999159 0.998892 0.998552 0.99812 0.997574 0.996883 0.99601 + 0.994909 0.99353 0.991835 0.989832 0.987583 1 1 1 1 1 1 1 + 1 1 0.999999 0.999999 0.999999 0.999998 0.999998 0.999996 0.999995 0.999993 0.999989 0.999985 + 0.999978 0.999969 0.999956 0.999938 0.999913 0.99988 0.999834 0.999772 0.999689 0.999579 0.999434 0.999243 + 0.998994 0.998669 0.998244 0.997686 0.996947 0.995952 0.994598 0.992779 0.990441 0.987583 1 1 + 1 1 1 1 1 1 1 1 0.999999 0.999999 0.999999 0.999998 + 0.999998 0.999997 0.999995 0.999993 0.99999 0.999986 0.99998 0.999972 0.999961 0.999945 0.999923 0.999893 + 0.999853 0.999798 0.999724 0.999626 0.999494 0.999318 0.999083 0.998768 0.998343 0.997754 0.996915 0.995666 + 0.993752 0.99104 0.987583 1 1 1 1 1 1 1 1 1 + 1 1 1 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 0.999991 0.999987 + 0.999982 0.999975 0.999965 0.999951 0.999932 0.999906 0.999871 0.999823 0.999758 0.99967 0.999551 0.999391 + 0.999174 0.998876 0.998458 0.997847 0.996883 0.995185 0.991937 0.987583 1 1 1 1 + 1 1 1 1 1 1 1 1 1 0.999999 0.999999 0.999999 + 0.999998 0.999997 0.999996 0.999994 0.999992 0.999989 0.999984 0.999978 0.999969 0.999957 0.999941 0.999919 + 0.999888 0.999847 0.999791 0.999715 0.999613 0.999476 0.99929 0.999038 0.998686 0.998165 0.997269 0.99518 + 0.987583 1 1 1 1 1 1 1 1 1 1 1 + 1 1 1 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999995 0.999993 0.99999 + 0.999986 0.999981 0.999973 0.999963 0.999949 0.999931 0.999905 0.999871 0.999826 0.999766 0.999686 0.999583 + 0.999453 0.999294 0.999118 0.998973 0.999022 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 1 0.999999 0.999999 0.999999 + 0.999998 0.999998 0.999997 0.999995 0.999994 0.999991 0.999988 0.999983 0.999977 0.999968 0.999957 0.999942 + 0.999922 0.999896 0.999863 0.999821 0.999768 0.999706 0.999637 0.99957 0.999531 0.999584 1 1 + 1 1 1 1 1 1 1 1 1 1 1 1 + 1 1 1 0.999999 0.999999 0.999999 0.999998 0.999998 0.999997 0.999995 0.999994 0.999991 + 0.999988 0.999984 0.999979 0.999972 0.999963 0.999951 0.999936 0.999918 0.999896 0.999871 0.999844 0.999819 + 0.999802 0.999807 0.999861 1 1 1 1 1 1 1 1 1 + 1 1 1 1 1 1 1 1 0.999999 0.999999 0.999999 0.999998 + 0.999998 0.999997 0.999996 0.999995 0.999993 0.99999 0.999987 0.999982 0.999977 0.999969 0.99996 0.999949 + 0.999936 0.999921 0.999905 0.99989 0.999878 0.999875 0.999885 0.999913 0.999958 + </DataArray> + <DataArray type="Float32" Name="x^Salt_liq" NumberOfComponents="1" format="ascii"> + 8.52212e-07 1.30914e-06 8.41594e-07 1.29453e-06 2.05758e-06 1.9826e-06 3.22967e-06 3.05903e-06 5.05913e-06 4.73552e-06 7.92094e-06 7.34256e-06 + 1.24148e-05 1.1399e-05 1.95078e-05 1.77222e-05 3.07785e-05 2.76074e-05 4.88476e-05 4.31227e-05 7.81679e-05 6.76093e-05 0.000126549 0.000106576 + 0.000208329 0.000169485 0.000351781 0.000274158 0.000621386 0.000462155 0.00127839 0.000877606 0.00450035 0.00204069 0.00633134 0.00410303 + 0.0078434 0.00563491 0.00901013 0.00691296 0.00985139 0.00816661 0.0104445 0.0092391 0.0108591 0.0101201 0.0111814 0.0107734 + 0.0114492 0.0112273 0.0116697 0.0115221 0.0118432 0.0117226 0.0119764 0.0118728 0.0120785 0.0119933 0.0121577 0.0120907 + 0.0122197 0.0121684 0.0122683 0.0122295 0.0123063 0.0122773 0.0123359 0.0123146 0.0123587 0.0123436 0.0123762 0.0123659 + 0.0123895 0.0123829 0.0123995 0.0123957 0.0124068 0.0124053 0.012412 0.0124124 0.0124163 0.0124171 7.24416e-07 1.11691e-06 + 1.70387e-06 2.59319e-06 3.95149e-06 6.02769e-06 9.20061e-06 1.40505e-05 2.1472e-05 3.28584e-05 5.04187e-05 7.77641e-05 0.000121102 0.000191933 + 0.00031371 0.00053805 0.000974531 0.00178114 0.00303864 0.00450119 0.00589304 0.00703248 0.0078644 0.00880738 0.009592 0.0102965 + 0.0108749 0.011311 0.0116115 0.0118179 0.0119666 0.0120798 0.0121678 0.0122358 0.0122881 0.0123279 0.0123581 0.0123807 + 0.0123974 0.0124092 0.0124171 5.7553e-07 8.83337e-07 1.34248e-06 2.02785e-06 3.05536e-06 4.60045e-06 6.92758e-06 1.04384e-05 1.57494e-05 + 2.38202e-05 3.61751e-05 5.53019e-05 8.54041e-05 0.000133863 0.000213984 0.000349195 0.000577089 0.000949216 0.00152679 0.00233867 0.00334997 + 0.00450684 0.00572646 0.00687486 0.00788443 0.00862691 0.00941115 0.0100323 0.0106055 0.0110829 0.0114472 0.0117133 0.0119074 + 0.0120517 0.01216 0.0122415 0.0123022 0.0123467 0.012379 0.0124016 0.0124171 4.34189e-07 6.63892e-07 1.00475e-06 1.50985e-06 + 2.25958e-06 3.37541e-06 5.04088e-06 7.53464e-06 1.12839e-05 1.69521e-05 2.5586e-05 3.88693e-05 5.95643e-05 9.22719e-05 0.000144638 0.000228909 + 0.000363157 0.000571122 0.00088025 0.00131729 0.00190177 0.00264315 0.00353444 0.00453791 0.00561417 0.00671469 0.00774598 0.0086609 + 0.00939609 0.00998342 0.0104974 0.0109826 0.0113805 0.0116762 0.0118987 0.0120653 0.0121886 0.0122784 0.0123422 0.0123867 + 0.0124171 3.17301e-07 4.83803e-07 7.29815e-07 1.09279e-06 1.62876e-06 2.42192e-06 3.59927e-06 5.35337e-06 7.97844e-06 1.19285e-05 1.79123e-05 + 2.70491e-05 4.11237e-05 6.29895e-05 9.71547e-05 0.000150505 0.000232935 0.000357502 0.000539843 0.000797025 0.00114608 0.00160207 0.00217582 + 0.00287206 0.00368823 0.00461352 0.00561656 0.00666096 0.00769892 0.0086527 0.00947788 0.0101598 0.0105896 0.0110463 0.0114413 + 0.0117446 0.0119769 0.0121498 0.0122734 0.0123589 0.0124171 2.27399e-07 3.45981e-07 5.20612e-07 7.77523e-07 1.15569e-06 1.71337e-06 + 2.53808e-06 3.76198e-06 5.5859e-06 8.3172e-06 1.24298e-05 1.86592e-05 2.81504e-05 4.26805e-05 6.49694e-05 9.90645e-05 0.000150726 0.000227671 + 0.000339556 0.000497652 0.000714308 0.00100218 0.00137333 0.00183844 0.00240592 0.00308076 0.00386338 0.00474769 0.00572054 0.00675301 + 0.00779391 0.00877978 0.00965614 0.0103705 0.0108025 0.0112472 0.0116352 0.0119326 0.012154 0.0123109 0.0124171 1.60965e-07 + 2.44445e-07 3.67054e-07 5.47011e-07 8.11274e-07 1.19995e-06 1.773e-06 2.62045e-06 3.87823e-06 5.75245e-06 8.55718e-06 1.27723e-05 1.91307e-05 + 2.87458e-05 4.32868e-05 6.52005e-05 9.79559e-05 0.000146267 0.000216229 0.000315324 0.000452282 0.000636806 0.000879192 0.00118989 0.00157903 + 0.0020559 0.00262852 0.00330288 0.00408195 0.0049637 0.00593658 0.00697684 0.00804398 0.00906696 0.00998128 0.0107004 0.0111592 + 0.0116037 0.0119706 0.0122344 0.0124171 1.12985e-07 1.7126e-07 2.5664e-07 3.81678e-07 5.64891e-07 8.33706e-07 1.22892e-06 1.81144e-06 + 2.67253e-06 3.94939e-06 5.84874e-06 8.6823e-06 1.29191e-05 1.92604e-05 2.87428e-05 4.28716e-05 6.3778e-05 9.43885e-05 0.000138578 0.000201279 + 0.00028853 0.000407438 0.00056607 0.000773274 0.00103845 0.00137131 0.00178156 0.00227862 0.00287113 0.00356631 0.00436878 0.00527838 + 0.00628654 0.0073711 0.00848863 0.00955996 0.0105415 0.0112308 0.011706 0.0121274 0.0124171 7.88122e-08 1.19214e-07 1.78262e-07 + 2.64542e-07 3.90676e-07 5.75293e-07 8.45961e-07 1.24359e-06 1.82907e-06 2.69314e-06 3.97119e-06 5.86511e-06 8.67502e-06 1.28444e-05 1.90222e-05 + 2.81444e-05 4.15358e-05 6.10287e-05 8.90903e-05 0.000128947 0.000184689 0.000261346 0.000364922 0.000502383 0.000681612 0.000911317 0.00120092 + 0.00156042 0.00200018 0.00253071 0.00316222 0.00390393 0.00476259 0.00574009 0.00682844 0.0080019 0.00921431 0.010385 0.0114916 + 0.0120314 0.0124171 5.46962e-08 8.25363e-08 1.2312e-07 1.82276e-07 2.68552e-07 3.94512e-07 5.78664e-07 8.4832e-07 1.24384e-06 1.82495e-06 + 2.67995e-06 3.93928e-06 5.79501e-06 8.52827e-06 1.25473e-05 1.8438e-05 2.70294e-05 3.94738e-05 5.73416e-05 8.27256e-05 0.000118351 0.00016768 + 0.000235012 0.000325562 0.000445522 0.000602102 0.000803558 0.00105919 0.00137936 0.00177542 0.0022597 0.00284531 0.00354578 0.00437415 + 0.00534095 0.00645007 0.00769251 0.00903749 0.0104264 0.0117643 0.0124171 3.77909e-08 5.68645e-08 8.45916e-08 1.24904e-07 1.8355e-07 + 2.68954e-07 3.93465e-07 5.75211e-07 8.40815e-07 1.22939e-06 1.79833e-06 2.63175e-06 3.85248e-06 5.63899e-06 8.24882e-06 1.20502e-05 1.75637e-05 + 2.55154e-05 3.69032e-05 5.30754e-05 7.5822e-05 0.000107476 0.000151023 0.000210215 0.000289684 0.000395065 0.000533107 0.00071181 0.000940556 + 0.00123027 0.00159363 0.00204528 0.00260212 0.00328346 0.00411091 0.00510697 0.00629109 0.00767081 0.00922613 0.0108774 0.0124171 + 2.60047e-08 3.89982e-08 5.78302e-08 8.51363e-08 1.24763e-07 1.82324e-07 2.66018e-07 3.87815e-07 5.65192e-07 8.23662e-07 1.20042e-06 1.74957e-06 + 2.54961e-06 3.71386e-06 5.40503e-06 7.85494e-06 1.1391e-05 1.6471e-05 2.37274e-05 3.40235e-05 4.85219e-05 6.87668e-05 9.67812e-05 0.000135178 + 0.00018729 0.000257312 0.000350474 0.000473231 0.000633506 0.000840989 0.00110754 0.00144774 0.00187974 0.00242644 0.00311719 0.00398972 + 0.00509062 0.00646961 0.00816456 0.0101684 0.0124171 1.7827e-08 2.66284e-08 3.93437e-08 5.77317e-08 8.43532e-08 1.22935e-07 1.78892e-07 + 2.60095e-07 3.77969e-07 5.49095e-07 7.97496e-07 1.15791e-06 1.68043e-06 2.43699e-06 3.53037e-06 5.10661e-06 7.37166e-06 1.06138e-05 1.52327e-05 + 2.1778e-05 3.09972e-05 4.38972e-05 6.18205e-05 8.65388e-05 0.000120368 0.000166308 0.000228215 0.000311021 0.000421001 0.000566148 0.000756683 + 0.00100582 0.00133097 0.00175575 0.00231353 0.00305328 0.00404841 0.00540205 0.00722053 0.00955923 0.0124171 1.21755e-08 1.81034e-08 + 2.66393e-08 3.89535e-08 5.67484e-08 8.24883e-08 1.19739e-07 1.73654e-07 2.51676e-07 3.64545e-07 5.27738e-07 7.63507e-07 1.10377e-06 1.59412e-06 + 2.29947e-06 3.31172e-06 4.76028e-06 6.82627e-06 9.76146e-06 1.39134e-05 1.97582e-05 2.79435e-05 3.93437e-05 5.51308e-05 7.6865e-05 0.000106611 + 0.000147087 0.000201857 0.000275592 0.000374413 0.000506404 0.00068237 0.000917097 0.0012316 0.00165744 0.00224556 0.00308454 0.004334 + 0.00624803 0.00896007 0.0124171 8.26525e-09 1.22228e-08 1.79272e-08 2.6156e-08 3.80224e-08 5.51585e-08 7.99069e-08 1.15634e-07 1.67177e-07 + 2.41479e-07 3.48488e-07 5.0244e-07 7.23637e-07 1.04096e-06 1.49534e-06 2.14453e-06 3.06969e-06 4.38425e-06 6.24601e-06 8.87322e-06 1.25661e-05 + 1.77355e-05 2.49403e-05 3.49373e-05 4.8746e-05 6.77342e-05 9.37301e-05 0.000129171 0.000177305 0.000242459 0.000330437 0.00044909 0.000609264 + 0.00082647 0.00112424 0.00154171 0.00215267 0.0031168 0.00481487 0.00806333 0.0124171 5.64704e-09 8.29164e-09 1.20873e-08 1.75692e-08 + 2.55033e-08 3.69455e-08 5.34374e-08 7.71804e-08 1.1132e-07 1.6035e-07 2.30673e-07 3.314e-07 4.75462e-07 6.81155e-07 9.74289e-07 1.39114e-06 + 1.98251e-06 2.81924e-06 3.99969e-06 5.65991e-06 7.98716e-06 1.12381e-05 1.57633e-05 2.20389e-05 3.07101e-05 4.26471e-05 5.90203e-05 8.13984e-05 + 0.00011188 0.000153265 0.000209291 0.000284948 0.000386933 0.000524342 0.000709888 0.000962496 0.0013144 0.0018353 0.00273118 0.0048197 + 0.0124171 3.93131e-09 5.72716e-09 8.29825e-09 1.20088e-08 1.73759e-08 2.51273e-08 3.62764e-08 5.22776e-08 7.51987e-08 1.07977e-07 1.54774e-07 + 2.21476e-07 3.16387e-07 4.51189e-07 6.42263e-07 9.12515e-07 1.29385e-06 1.83057e-06 2.58393e-06 3.63834e-06 5.10962e-06 7.1561e-06 9.99337e-06 + 1.39138e-05 1.9312e-05 2.6719e-05 3.6845e-05 5.06359e-05 6.93434e-05 9.4612e-05 0.000128581 0.000173992 0.000234285 0.000313606 0.000416563 + 0.000547263 0.000706361 0.00088238 0.00102649 0.000977728 0 2.84383e-09 4.11099e-09 5.9293e-09 8.56226e-09 1.23638e-08 1.7838e-08 + 2.56885e-08 3.69111e-08 5.29144e-08 7.56865e-08 1.08026e-07 1.53864e-07 2.18704e-07 3.10234e-07 4.39155e-07 6.20317e-07 8.74254e-07 1.22926e-06 + 1.72416e-06 2.41204e-06 3.36515e-06 4.68138e-06 6.49277e-06 8.97641e-06 1.23685e-05 1.69818e-05 2.32275e-05 3.16406e-05 4.29091e-05 5.79028e-05 + 7.76964e-05 0.000103571 0.000136966 0.000179323 0.000231705 0.000293975 0.000363086 0.000429492 0.000468795 0.000415863 0 2.20205e-09 + 3.17964e-09 4.59463e-09 6.65422e-09 9.61884e-09 1.38645e-08 1.99257e-08 2.85492e-08 4.07859e-08 5.81091e-08 8.25787e-08 1.17066e-07 1.65562e-07 + 2.33593e-07 3.28792e-07 4.61655e-07 6.46564e-07 9.03141e-07 1.25803e-06 1.74725e-06 2.4192e-06 3.33852e-06 4.59098e-06 6.2894e-06 8.58085e-06 + 1.16547e-05 1.57515e-05 2.11707e-05 2.82763e-05 3.74945e-05 4.92981e-05 6.41636e-05 8.24824e-05 0.000104393 0.000129483 0.000156284 0.000181408 + 0.000198038 0.000192886 0.000139253 0 1.8957e-09 2.74374e-09 4.07032e-09 5.97835e-09 8.69917e-09 1.25706e-08 1.807e-08 2.58666e-08 + 3.68972e-08 5.24697e-08 7.44053e-08 1.05233e-07 1.48452e-07 2.08889e-07 2.93179e-07 4.104e-07 5.72926e-07 7.97541e-07 1.10689e-06 1.53137e-06 + 2.11149e-06 2.90083e-06 3.96968e-06 5.40927e-06 7.33647e-06 9.89865e-06 1.32778e-05 1.76923e-05 2.3394e-05 3.0656e-05 3.97437e-05 5.08602e-05 + 6.40495e-05 7.90407e-05 9.50114e-05 0.000110264 0.000121847 0.000125291 0.000114979 8.65827e-05 4.19137e-05 + </DataArray> + <DataArray type="Float32" Name="X^Water_liq" NumberOfComponents="1" format="ascii"> + 0.999997 0.999996 0.999997 0.999996 0.999993 0.999994 0.99999 0.99999 0.999984 0.999985 0.999974 0.999976 + 0.99996 0.999963 0.999937 0.999942 0.9999 0.99991 0.999841 0.99986 0.999746 0.999781 0.999589 0.999654 + 0.999324 0.99945 0.998859 0.99911 0.997985 0.998501 0.995861 0.997156 0.985535 0.993405 0.979733 0.9868 + 0.974976 0.981934 0.971328 0.977899 0.968708 0.973963 0.966868 0.970614 0.965584 0.967874 0.964587 0.965849 + 0.96376 0.964446 0.96308 0.963535 0.962546 0.962917 0.962135 0.962454 0.961821 0.962083 0.961577 0.961783 + 0.961387 0.961545 0.961237 0.961357 0.96112 0.96121 0.96103 0.961095 0.960959 0.961006 0.960905 0.960937 + 0.960865 0.960885 0.960834 0.960846 0.960812 0.960816 0.960796 0.960794 0.960782 0.96078 0.999998 0.999996 + 0.999994 0.999992 0.999987 0.99998 0.99997 0.999954 0.99993 0.999893 0.999836 0.999748 0.999607 0.999377 + 0.998982 0.998255 0.996843 0.99424 0.990201 0.985532 0.981117 0.977523 0.97491 0.97196 0.969515 0.967327 + 0.965535 0.964187 0.96326 0.962624 0.962166 0.961817 0.961547 0.961337 0.961176 0.961054 0.960961 0.960892 + 0.960841 0.960804 0.96078 0.999998 0.999997 0.999996 0.999993 0.99999 0.999985 0.999978 0.999966 0.999949 + 0.999923 0.999883 0.99982 0.999723 0.999566 0.999306 0.998867 0.998129 0.996925 0.99506 0.992447 0.989205 + 0.985514 0.981644 0.978019 0.974847 0.972524 0.970078 0.968146 0.966369 0.964892 0.963767 0.962946 0.962348 + 0.961904 0.961571 0.96132 0.961133 0.960996 0.960897 0.960828 0.96078 0.999999 0.999998 0.999997 0.999995 + 0.999993 0.999989 0.999984 0.999976 0.999963 0.999945 0.999917 0.999874 0.999807 0.9997 0.999531 0.999257 + 0.998822 0.998148 0.997148 0.995736 0.993852 0.991469 0.988615 0.985416 0.982 0.978524 0.975282 0.972418 + 0.970125 0.968298 0.966704 0.965202 0.963973 0.96306 0.962375 0.961862 0.961482 0.961206 0.96101 0.960873 + 0.96078 0.999999 0.999998 0.999998 0.999996 0.999995 0.999992 0.999988 0.999983 0.999974 0.999961 0.999942 + 0.999912 0.999866 0.999795 0.999685 0.999512 0.999244 0.99884 0.998249 0.997417 0.996289 0.994817 0.99297 + 0.990735 0.988124 0.985175 0.981992 0.978693 0.975429 0.972443 0.96987 0.967751 0.966418 0.965005 0.963785 + 0.962849 0.962134 0.961602 0.961222 0.960959 0.96078 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 + 0.999992 0.999988 0.999982 0.999973 0.99996 0.999939 0.999909 0.999861 0.999789 0.999678 0.999511 0.999261 + 0.998898 0.998386 0.997685 0.996754 0.995555 0.994056 0.992231 0.990067 0.987565 0.984748 0.981663 0.978403 + 0.975131 0.972046 0.969315 0.967097 0.965759 0.964384 0.963187 0.96227 0.961589 0.961106 0.96078 0.999999 + 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 0.999991 0.999987 0.999981 0.999972 0.999959 0.999938 + 0.999907 0.999859 0.999788 0.999682 0.999525 0.999298 0.998977 0.998533 0.997935 0.997151 0.996147 0.994892 + 0.993356 0.991516 0.989356 0.986868 0.984062 0.98098 0.977698 0.974347 0.97115 0.968305 0.966075 0.964656 + 0.963284 0.962153 0.961342 0.96078 1 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 + 0.999991 0.999987 0.999981 0.999972 0.999958 0.999937 0.999907 0.999861 0.999793 0.999694 0.99955 0.999347 + 0.999064 0.998678 0.998164 0.997494 0.996636 0.995561 0.994239 0.99264 0.990738 0.988513 0.985954 0.983064 + 0.979874 0.976458 0.972956 0.969615 0.966567 0.964435 0.962969 0.961671 0.96078 1 1 0.999999 + 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 0.999991 0.999987 0.999981 0.999972 0.999958 0.999938 + 0.999909 0.999865 0.999802 0.999711 0.999581 0.999401 0.999152 0.998816 0.998371 0.99779 0.997047 0.996112 + 0.994952 0.993535 0.99183 0.989806 0.987435 0.984701 0.981601 0.978165 0.974479 0.970691 0.967052 0.96363 + 0.961966 0.96078 1 1 1 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 + 0.999991 0.999987 0.999981 0.999972 0.999959 0.99994 0.999912 0.999872 0.999814 0.999731 0.999616 0.999456 + 0.999237 0.998944 0.998555 0.998048 0.997396 0.996569 0.995535 0.994259 0.992701 0.990821 0.988579 0.985937 + 0.982865 0.979358 0.975449 0.971242 0.966924 0.962789 0.96078 1 1 1 1 0.999999 + 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 0.999991 0.999987 0.999982 0.999973 0.999961 0.999943 + 0.999917 0.99988 0.999828 0.999754 0.999651 0.99951 0.999318 0.99906 0.998719 0.998271 0.997693 0.996953 + 0.996017 0.994844 0.99339 0.991601 0.989418 0.986775 0.983607 0.97986 0.975517 0.970654 0.965527 0.96078 + 1 1 1 1 1 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 + 0.999992 0.999988 0.999982 0.999974 0.999963 0.999947 0.999923 0.99989 0.999842 0.999777 0.999686 0.999561 + 0.999392 0.999165 0.998863 0.998465 0.997946 0.997275 0.996413 0.995315 0.993923 0.992165 0.98995 0.987162 + 0.983659 0.979296 0.97397 0.967724 0.96078 1 1 1 1 1 1 0.999999 + 0.999999 0.999999 0.999998 0.999997 0.999996 0.999995 0.999992 0.999989 0.999983 0.999976 0.999966 0.999951 + 0.999929 0.999899 0.999857 0.999799 0.999719 0.999609 0.99946 0.999259 0.998991 0.998634 0.998164 0.997547 + 0.996742 0.995692 0.994322 0.992528 0.990155 0.986975 0.982672 0.976932 0.969617 0.96078 1 1 + 1 1 1 1 1 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999995 + 0.999993 0.999989 0.999985 0.999978 0.999968 0.999955 0.999936 0.999909 0.999872 0.999821 0.99975 0.999654 + 0.999523 0.999345 0.999106 0.998785 0.998358 0.997788 0.997029 0.996012 0.994639 0.992746 0.990054 0.986065 + 0.979996 0.971484 0.96078 1 1 1 1 1 1 1 1 0.999999 + 0.999999 0.999999 0.999998 0.999998 0.999997 0.999995 0.999993 0.99999 0.999986 0.99998 0.999971 0.999959 + 0.999942 0.999919 0.999887 0.999842 0.99978 0.999696 0.999581 0.999425 0.999213 0.998928 0.998543 0.998025 + 0.997322 0.996359 0.995012 0.993045 0.989951 0.984535 0.974287 0.96078 1 1 1 1 + 1 1 1 1 1 0.999999 0.999999 0.999999 0.999998 0.999998 0.999997 0.999995 + 0.999994 0.999991 0.999987 0.999982 0.999974 0.999964 0.999949 0.999928 0.9999 0.999862 0.999808 0.999736 + 0.999637 0.999503 0.999321 0.999075 0.998745 0.9983 0.997699 0.996882 0.995745 0.994066 0.991187 0.98452 + 0.96078 1 1 1 1 1 1 1 1 1 1 1 + 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 0.999992 0.999988 0.999983 0.999977 0.999968 + 0.999955 0.999937 0.999913 0.99988 0.999836 0.999775 0.999693 0.999583 0.999435 0.99924 0.998983 0.998649 + 0.998225 0.99771 0.997141 0.996675 0.996833 1 1 1 1 1 1 1 + 1 1 1 1 1 1 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 + 0.999994 0.999992 0.999989 0.999985 0.999979 0.999971 0.99996 0.999945 0.999925 0.999897 0.999861 0.999812 + 0.999748 0.999664 0.999555 0.999418 0.999248 0.999046 0.998822 0.998607 0.99848 0.998651 1 1 + 1 1 1 1 1 1 1 1 1 1 1 0.999999 + 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999994 0.999992 0.999989 0.999985 0.99998 0.999972 + 0.999962 0.999949 0.999931 0.999908 0.999878 0.99984 0.999792 0.999732 0.999661 0.99958 0.999493 0.999411 + 0.999357 0.999374 0.999548 1 1 1 1 1 1 1 1 1 + 1 1 1 1 1 0.999999 0.999999 0.999999 0.999998 0.999997 0.999996 0.999995 + 0.999993 0.999991 0.999987 0.999982 0.999976 0.999968 0.999957 0.999943 0.999924 0.9999 0.999871 0.999835 + 0.999792 0.999743 0.999692 0.999642 0.999605 0.999593 0.999627 0.999719 0.999864 + </DataArray> + <DataArray type="Float32" Name="X^Salt_liq" NumberOfComponents="1" format="ascii"> + 2.76684e-06 4.25033e-06 2.73237e-06 4.20291e-06 6.68025e-06 6.43681e-06 1.04856e-05 9.93159e-06 1.64251e-05 1.53745e-05 2.57162e-05 2.38385e-05 + 4.03057e-05 3.70077e-05 6.33327e-05 5.75358e-05 9.99205e-05 8.96266e-05 0.000158575 0.000139992 0.000253741 0.000219471 0.000410745 0.000345935 + 0.00067606 0.000550052 0.00114121 0.000889551 0.00201462 0.00149891 0.00413863 0.00284369 0.0144649 0.0065952 0.0202675 0.0131995 + 0.025024 0.018066 0.0286725 0.0221008 0.0312916 0.0260365 0.0331323 0.0293863 0.0344161 0.032126 0.0354127 0.0341512 + 0.0362395 0.0355544 0.0369198 0.0364646 0.0374545 0.0370828 0.0378645 0.0375456 0.0381789 0.0379166 0.0384226 0.0382165 + 0.0386133 0.0384555 0.0387628 0.0386434 0.0388796 0.0387904 0.0389704 0.0389051 0.0390407 0.0389941 0.0390945 0.0390626 + 0.0391353 0.0391149 0.0391659 0.0391544 0.0391883 0.0391839 0.0392044 0.0392057 0.0392176 0.03922 2.35193e-06 3.62624e-06 + 5.53187e-06 8.41916e-06 1.2829e-05 1.95696e-05 2.98707e-05 4.56159e-05 6.97092e-05 0.000106673 0.000163674 0.00025243 0.000393069 0.000622875 + 0.00101779 0.00174476 0.00315707 0.00575973 0.00979856 0.0144675 0.0188827 0.022477 0.0250898 0.0280398 0.0304851 0.0326733 + 0.0344651 0.0358131 0.0367403 0.0373764 0.0378344 0.0381829 0.0384535 0.0386628 0.0388236 0.038946 0.0390387 0.0391082 + 0.0391594 0.0391957 0.03922 1.86855e-06 2.8679e-06 4.35857e-06 6.58372e-06 9.91968e-06 1.4936e-05 2.24912e-05 3.38893e-05 5.11314e-05 + 7.73322e-05 0.000117439 0.000179524 0.000277225 0.000434479 0.0006944 0.00113283 0.00187119 0.00307523 0.00494002 0.0075532 0.010795 + 0.0144855 0.0183558 0.0219809 0.0251526 0.0274762 0.0299222 0.0318536 0.0336311 0.0351082 0.0362334 0.0370541 0.0376521 + 0.0380962 0.0384294 0.0386802 0.0388668 0.0390038 0.039103 0.0391725 0.03922 1.40967e-06 2.15543e-06 3.26209e-06 4.90198e-06 + 7.33607e-06 1.09588e-05 1.63659e-05 2.44621e-05 3.66342e-05 5.50357e-05 8.30645e-05 0.000126185 0.000193359 0.000299514 0.000469439 0.00074281 + 0.00117809 0.00185187 0.00285224 0.00426417 0.00614815 0.00853077 0.0113847 0.0145844 0.0180003 0.0214764 0.0247185 0.0275824 + 0.0298753 0.0317018 0.0332962 0.0347982 0.0360275 0.0369397 0.0376252 0.0381382 0.0385176 0.0387936 0.0389898 0.0391265 + 0.03922 1.03017e-06 1.57074e-06 2.36946e-06 3.54793e-06 5.28802e-06 7.86312e-06 1.16855e-05 1.73804e-05 2.59029e-05 3.87269e-05 5.81528e-05 + 8.7814e-05 0.000133503 0.000204477 0.00031536 0.000488474 0.000755866 0.00115976 0.00175057 0.00258305 0.00371138 0.00518272 0.00702979 + 0.00926484 0.011876 0.0148249 0.0180079 0.0213071 0.0245708 0.0275567 0.0301299 0.0322493 0.0335821 0.0349952 0.0362152 + 0.0371505 0.0378662 0.0383982 0.0387784 0.0390412 0.03922 7.38289e-07 1.12328e-06 1.69025e-06 2.52436e-06 3.75213e-06 5.56272e-06 + 8.24026e-06 1.22138e-05 1.81353e-05 2.70027e-05 4.03543e-05 6.05776e-05 9.13893e-05 0.000138556 0.000210903 0.000321558 0.00048919 0.000738795 + 0.00110159 0.00161391 0.0023154 0.00324644 0.00444503 0.00594424 0.00776922 0.00993345 0.0124352 0.0152515 0.018337 0.0215971 + 0.0248688 0.0279536 0.0306846 0.0329028 0.0342411 0.035616 0.0368133 0.0377297 0.0384112 0.0388936 0.03922 5.22599e-07 + 7.9363e-07 1.1917e-06 1.77596e-06 2.63393e-06 3.89583e-06 5.7563e-06 8.50769e-06 1.25912e-05 1.86761e-05 2.77818e-05 4.14661e-05 6.21082e-05 + 9.33219e-05 0.000140524 0.000211653 0.00031796 0.000474724 0.000701683 0.00102303 0.00146692 0.00206454 0.00284882 0.00385288 0.00510847 + 0.00664415 0.00848382 0.0106444 0.0131323 0.0159377 0.0190204 0.0223019 0.0256525 0.0288497 0.0316951 0.0339249 0.0353441 + 0.0367162 0.0378468 0.0386583 0.03922 3.66823e-07 5.56023e-07 8.33223e-07 1.23918e-06 1.83401e-06 2.70676e-06 3.98988e-06 5.88111e-06 + 8.67676e-06 1.28222e-05 1.89887e-05 2.8188e-05 4.19427e-05 6.25293e-05 9.33124e-05 0.000139176 0.000207036 0.000306383 0.000449775 0.00065319 + 0.000936154 0.00132161 0.00183551 0.00250621 0.00336367 0.00443851 0.00576107 0.00736024 0.00926186 0.0114866 0.0140461 0.0169363 + 0.020126 0.0235416 0.027044 0.0303854 0.033433 0.0355653 0.0370314 0.0383293 0.03922 2.55877e-07 3.87048e-07 5.78758e-07 + 8.58878e-07 1.26839e-06 1.86778e-06 2.74655e-06 4.03753e-06 5.93836e-06 8.74367e-06 1.2893e-05 1.90418e-05 2.81644e-05 4.17004e-05 6.17562e-05 + 9.13697e-05 0.00013484 0.000198113 0.000289189 0.000418526 0.000599374 0.000848005 0.00118381 0.00162923 0.00220958 0.0029527 0.0038885 + 0.00504847 0.00646488 0.00816991 0.0101943 0.0125646 0.0152988 0.0183989 0.0218347 0.0255207 0.029309 0.032948 0.0363705 + 0.0380338 0.03922 1.7758e-07 2.67968e-07 3.9973e-07 5.9179e-07 8.71898e-07 1.28085e-06 1.87873e-06 2.75421e-06 4.03834e-06 5.92498e-06 + 8.70084e-06 1.27894e-05 1.88142e-05 2.76879e-05 4.07357e-05 5.98596e-05 8.775e-05 0.000128147 0.000186145 0.000268532 0.000384143 0.000544196 + 0.000762604 0.00105622 0.00144501 0.00195218 0.00260418 0.00343069 0.00446449 0.0057413 0.00729942 0.00917909 0.011421 0.0140632 + 0.0171347 0.0206421 0.0245507 0.0287578 0.0330762 0.0372111 0.03922 1.22695e-07 1.8462e-07 2.74641e-07 4.0552e-07 5.95926e-07 + 8.73204e-07 1.27745e-06 1.86751e-06 2.72984e-06 3.9914e-06 5.83856e-06 8.54437e-06 1.25076e-05 1.83077e-05 2.67807e-05 3.91219e-05 5.70211e-05 + 8.28353e-05 0.000119803 0.000172298 0.000246127 0.000348855 0.000490157 0.000682176 0.000939897 0.00128151 0.00172875 0.00230732 0.00304723 + 0.00398328 0.00515553 0.00660998 0.00839911 0.0105822 0.0132246 0.0163926 0.0201404 0.0244826 0.0293459 0.0344727 0.03922 + 8.44284e-08 1.26614e-07 1.87755e-07 2.76409e-07 4.05063e-07 5.91946e-07 8.63672e-07 1.2591e-06 1.83499e-06 2.67415e-06 3.89734e-06 5.68025e-06 + 8.27768e-06 1.20576e-05 1.75481e-05 2.55019e-05 3.69819e-05 5.34738e-05 7.70308e-05 0.000110454 0.000157517 0.000223229 0.000314148 0.000438745 + 0.000607811 0.000834924 0.00113698 0.00153479 0.00205386 0.00272526 0.00358688 0.00468509 0.00607722 0.00783512 0.0100501 0.0128382 + 0.0163407 0.0207037 0.0260301 0.0322761 0.03922 5.78783e-08 8.64534e-08 1.27736e-07 1.87435e-07 2.73867e-07 3.99128e-07 5.80803e-07 + 8.44441e-07 1.22714e-06 1.78273e-06 2.5892e-06 3.75935e-06 5.45579e-06 7.91205e-06 1.14619e-05 1.65793e-05 2.39329e-05 3.44585e-05 4.94539e-05 + 7.07024e-05 0.00010063 0.000142506 0.000200683 0.000280908 0.000390688 0.000539744 0.00074056 0.00100908 0.00136556 0.00183576 0.00245253 + 0.0032582 0.00430832 0.00567795 0.00747241 0.00984546 0.0130254 0.0173284 0.0230684 0.0303831 0.03922 3.95299e-08 5.87756e-08 + 8.6489e-08 1.26469e-07 1.84243e-07 2.67812e-07 3.88753e-07 5.63796e-07 8.17107e-07 1.18356e-06 1.71339e-06 2.47885e-06 3.58355e-06 5.17556e-06 + 7.46557e-06 1.0752e-05 1.54549e-05 2.21623e-05 3.16915e-05 4.51706e-05 6.41453e-05 9.07175e-05 0.000127725 0.000178969 0.000249512 0.000346048 + 0.000477384 0.000655066 0.0008942 0.00121457 0.00164226 0.00221204 0.00297139 0.00398756 0.00536118 0.00725397 0.00994554 0.0139354 + 0.0200045 0.0285163 0.03922 2.68345e-08 3.96834e-08 5.82035e-08 8.49197e-08 1.23446e-07 1.79081e-07 2.59431e-07 3.75426e-07 5.42769e-07 + 7.84e-07 1.13142e-06 1.63125e-06 2.34941e-06 3.37965e-06 4.85485e-06 6.96254e-06 9.96618e-06 1.42341e-05 2.02784e-05 2.88078e-05 4.07968e-05 + 5.75789e-05 8.09682e-05 0.000113421 0.000158245 0.000219877 0.000304246 0.000419255 0.00057542 0.000786756 0.00107202 0.00145658 0.00197537 + 0.0026783 0.00364083 0.00498815 0.00695537 0.0100488 0.015465 0.0257131 0.03922 1.8334e-08 2.69202e-08 3.92435e-08 5.70413e-08 + 8.28006e-08 1.1995e-07 1.73493e-07 2.50579e-07 3.6142e-07 5.20602e-07 7.48917e-07 1.07595e-06 1.54367e-06 2.21148e-06 3.16319e-06 4.51656e-06 + 6.43652e-06 9.15306e-06 1.29856e-05 1.83756e-05 2.59312e-05 3.64856e-05 5.11763e-05 7.15495e-05 9.96986e-05 0.000138448 0.000191594 0.000264225 + 0.000363145 0.000497429 0.000679178 0.000924539 0.00125515 0.00170036 0.0023011 0.00311816 0.00425487 0.00593413 0.00881315 0.0154803 + 0.03922 1.27636e-08 1.85942e-08 2.69417e-08 3.89887e-08 5.64139e-08 8.15799e-08 1.17778e-07 1.69728e-07 2.44145e-07 3.50564e-07 5.02498e-07 + 7.19058e-07 1.0272e-06 1.46486e-06 2.08521e-06 2.96262e-06 4.20069e-06 5.94323e-06 8.38912e-06 1.18124e-05 1.6589e-05 2.32331e-05 3.24444e-05 + 4.51719e-05 6.2697e-05 8.67424e-05 0.000119613 0.000164379 0.0002251 0.000307108 0.000417338 0.000564672 0.000760245 0.00101746 0.00135118 + 0.0017746 0.00228969 0.00285913 0.003325 0.0031674 0 9.23297e-09 1.3347e-08 1.92505e-08 2.77988e-08 4.01412e-08 5.7914e-08 + 8.34019e-08 1.19838e-07 1.71795e-07 2.45729e-07 3.50725e-07 4.99545e-07 7.10059e-07 1.00723e-06 1.42579e-06 2.01396e-06 2.83841e-06 3.99098e-06 + 5.59776e-06 7.83105e-06 1.09254e-05 1.51987e-05 2.10796e-05 2.91428e-05 4.01552e-05 5.51321e-05 7.54079e-05 0.000102719 0.000139298 0.000187967 + 0.00025221 0.000336182 0.000444546 0.000581967 0.000751877 0.000953808 0.00117786 0.00139307 0.00152042 0.00134891 0 7.14931e-09 + 1.03232e-08 1.49172e-08 2.1604e-08 3.12292e-08 4.50135e-08 6.4692e-08 9.26896e-08 1.32418e-07 1.88661e-07 2.68106e-07 3.80074e-07 5.37523e-07 + 7.58398e-07 1.06748e-06 1.49884e-06 2.09918e-06 2.93219e-06 4.0844e-06 5.67273e-06 7.85429e-06 1.0839e-05 1.49052e-05 2.04193e-05 2.78586e-05 + 3.7838e-05 5.11381e-05 6.87311e-05 9.17979e-05 0.000121722 0.000160037 0.000208288 0.000267743 0.00033885 0.000420267 0.000507224 0.000588732 + 0.000642677 0.000625966 0.000451966 0 6.15471e-09 8.90801e-09 1.3215e-08 1.94097e-08 2.82433e-08 4.08126e-08 5.86672e-08 8.39803e-08 + 1.19793e-07 1.70351e-07 2.41569e-07 3.41655e-07 4.81973e-07 6.78194e-07 9.51854e-07 1.33243e-06 1.8601e-06 2.58934e-06 3.59371e-06 4.97184e-06 + 6.85526e-06 9.41796e-06 1.28881e-05 1.75619e-05 2.38187e-05 3.21369e-05 4.31073e-05 5.74387e-05 7.59487e-05 9.9523e-05 0.000129023 0.000165107 + 0.000207918 0.000256573 0.000308404 0.000357902 0.00039549 0.000406664 0.000373203 0.000281051 0.000136067 + </DataArray> </PointData> <CellData Scalars="process rank"> <DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">