Commit 3f994ad6 authored by Natalie Schroeder's avatar Natalie Schroeder
Browse files

[richards]

* change in  richards output - pressure head output only if gravity is defined
* update of reference solutions of all richards tests - needed due to change in the output (pressure head, water content and source values)

discussed and reviewed by Bernd



git-svn-id: svn://svn.iws.uni-stuttgart.de/DUMUX/dumux/trunk@15062 2fb0f335-1f38-0410-981e-8018bf24f1b0
parent 55473cfb
......@@ -127,8 +127,6 @@ class RichardsModel : public GET_PROP_TYPE(TypeTag, BaseModel)
enum { isBox = GET_PROP_VALUE(TypeTag, ImplicitIsBox) };
enum { dofCodim = isBox ? dim : 0 };
static const bool useHead = GET_PROP_VALUE(TypeTag, UseHead);
public:
/*!
* \brief All relevant primary and secondary of a given
......@@ -142,6 +140,7 @@ public:
{
typedef Dune::BlockVector<Dune::FieldVector<double, 1> > ScalarField;
typedef Dune::BlockVector<Dune::FieldVector<double, dimWorld> > VectorField;
bool gravity =GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, bool, Problem, EnableGravity);
// get the number of degrees of freedom
unsigned numDofs = this->numDofs();
......@@ -223,7 +222,8 @@ public:
(*mobN)[dofIdxGlobal] = elemVolVars[scvIdx].mobility(nPhaseIdx);
(*poro)[dofIdxGlobal] = elemVolVars[scvIdx].porosity();
(*Te)[dofIdxGlobal] = elemVolVars[scvIdx].temperature();
(*ph)[dofIdxGlobal] = elemVolVars[scvIdx].pressureHead(wPhaseIdx);
if (gravity)
(*ph)[dofIdxGlobal] = elemVolVars[scvIdx].pressureHead(wPhaseIdx);
(*wc)[dofIdxGlobal] = elemVolVars[scvIdx].waterContent(wPhaseIdx);
(*source)[dofIdxGlobal] = sourcevalues[0];
}
......@@ -244,7 +244,8 @@ public:
writer.attachDofData(*mobN, "mobN", isBox);
writer.attachDofData(*poro, "porosity", isBox);
writer.attachDofData(*Te, "temperature", isBox);
writer.attachDofData(*ph, "pressure head", isBox);
if (gravity)
writer.attachDofData(*ph, "pressure head", isBox);
writer.attachDofData(*wc, "water content", isBox);
writer.attachDofData(*source, "source", isBox);
......
......@@ -63,7 +63,13 @@ public:
*/
RichardsProblem(TimeManager &timeManager, const GridView &gridView)
: ParentType(timeManager, gridView)
{}
{
bool gravity_ =GET_RUNTIME_PARAM_FROM_GROUP(TypeTag, bool, Problem, EnableGravity);
static const bool useHead = GET_PROP_VALUE(TypeTag, UseHead);
if ( (!gravity_) && (useHead) )
DUNE_THROW(Dune::InvalidStateException, "Don't run the Richards model with pressure head without gravity!");
}
/*!
* \name Problem parameters
......
......@@ -7,6 +7,7 @@ File = ./grids/richardslens-24x16.dgf
[Problem]
Name = richardslensbox
EnableGravity = 1 # enable gravity
[Implicit]
EnablePartialReassemble = 1 # enable partial reassembly of the Jacobian matrix
......
......@@ -7,6 +7,7 @@ File = ./grids/richardslens-24x16.dgf
[Problem]
Name = richardslenscc
EnableGravity = 1 # enable gravity
[Implicit]
EnablePartialReassemble = 1 # enable partial reassembly of the Jacobian matrix
......
......@@ -212,56 +212,94 @@
290 290 290 290 290 290 290 290 290 290 290 290
290 290 290 290 290 290 290 290
</DataArray>
<DataArray type="Float32" Name="water content" NumberOfComponents="1" format="ascii">
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
</DataArray>
<DataArray type="Float32" Name="source" NumberOfComponents="1" format="ascii">
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
</DataArray>
<DataArray type="Float32" Name="velocity" NumberOfComponents="3" format="ascii">
-1.03797e-09 0 0 -1.02573e-09 0 0 -1.00183e-09 0 0 -9.69465e-10 0 0
-9.30241e-10 0 0 -8.83398e-10 0 0 -8.30174e-10 0 0 -7.7189e-10 0 0
-9.30241e-10 0 0 -8.83399e-10 0 0 -8.30174e-10 0 0 -7.7189e-10 0 0
-7.09893e-10 0 0 -6.50497e-10 0 0 -5.94981e-10 0 0 -5.39432e-10 0 0
-4.84687e-10 0 0 -4.31453e-10 0 0 -3.80308e-10 0 0 -3.31701e-10 0 0
-2.85956e-10 0 0 -2.4329e-10 0 0 -2.0382e-10 0 0 -1.6758e-10 0 0
-1.34542e-10 0 0 -1.06341e-10 0 0 -8.36012e-11 0 0 -6.4169e-11 0 0
-4.68528e-11 0 0 -3.14953e-11 0 0 -1.79351e-11 0 0 -6.01105e-12 0 0
4.43408e-12 0 0 1.35505e-11 0 0 2.14803e-11 0 0 2.83564e-11 0 0
3.4301e-11 0 0 3.94261e-11 0 0 4.38332e-11 0 0 4.76134e-11 0 0
-1.34542e-10 0 0 -1.06341e-10 0 0 -8.36011e-11 0 0 -6.41689e-11 0 0
-4.68528e-11 0 0 -3.14954e-11 0 0 -1.79352e-11 0 0 -6.01101e-12 0 0
4.43412e-12 0 0 1.35505e-11 0 0 2.14803e-11 0 0 2.83564e-11 0 0
3.4301e-11 0 0 3.94262e-11 0 0 4.38332e-11 0 0 4.76134e-11 0 0
5.08484e-11 0 0 5.36109e-11 0 0 5.59649e-11 0 0 5.7967e-11 0 0
5.96667e-11 0 0 6.11072e-11 0 0 6.23261e-11 0 0 6.33557e-11 0 0
6.42243e-11 0 0 6.4956e-11 0 0 6.55716e-11 0 0 6.60889e-11 0 0
6.6523e-11 0 0 6.6887e-11 0 0 6.71918e-11 0 0 6.74469e-11 0 0
6.766e-11 0 0 6.78381e-11 0 0 6.79866e-11 0 0 6.81105e-11 0 0
6.76601e-11 0 0 6.78381e-11 0 0 6.79866e-11 0 0 6.81105e-11 0 0
6.82136e-11 0 0 6.82995e-11 0 0 6.8371e-11 0 0 6.84304e-11 0 0
6.84798e-11 0 0 6.85208e-11 0 0 6.85548e-11 0 0 6.8583e-11 0 0
6.86064e-11 0 0 6.86258e-11 0 0 6.86419e-11 0 0 6.86551e-11 0 0
6.86662e-11 0 0 6.86753e-11 0 0 6.86828e-11 0 0 6.86891e-11 0 0
6.86064e-11 0 0 6.86258e-11 0 0 6.86419e-11 0 0 6.86552e-11 0 0
6.86662e-11 0 0 6.86753e-11 0 0 6.86829e-11 0 0 6.86891e-11 0 0
6.86942e-11 0 0 6.86984e-11 0 0 6.87019e-11 0 0 6.87048e-11 0 0
6.87072e-11 0 0 6.87092e-11 0 0 6.87109e-11 0 0 6.87121e-11 0 0
6.87132e-11 0 0 6.87142e-11 0 0 6.8715e-11 0 0 6.87156e-11 0 0
6.87161e-11 0 0 6.87165e-11 0 0 6.87169e-11 0 0 6.87171e-11 0 0
6.87174e-11 0 0 6.87176e-11 0 0 6.87177e-11 0 0 6.87179e-11 0 0
6.87072e-11 0 0 6.87092e-11 0 0 6.87108e-11 0 0 6.87122e-11 0 0
6.87133e-11 0 0 6.87142e-11 0 0 6.8715e-11 0 0 6.87156e-11 0 0
6.8716e-11 0 0 6.87165e-11 0 0 6.87169e-11 0 0 6.87172e-11 0 0
6.87174e-11 0 0 6.87175e-11 0 0 6.87177e-11 0 0 6.87179e-11 0 0
6.8718e-11 0 0 6.8718e-11 0 0 6.87181e-11 0 0 6.87182e-11 0 0
6.87182e-11 0 0 6.87182e-11 0 0 6.87183e-11 0 0 6.87183e-11 0 0
6.87182e-11 0 0 6.87183e-11 0 0 6.87183e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0
6.87184e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87184e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
6.87185e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0
6.87185e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87185e-11 0 0 6.87184e-11 0 0 6.87184e-11 0 0
6.87184e-11 0 0 6.87184e-11 0 0 6.87185e-11 0 0 6.87185e-11 0 0
</DataArray>
<DataArray type="Float32" Name="process rank" NumberOfComponents="1" format="ascii">
......
......@@ -117,11 +117,11 @@
</DataArray>
<DataArray type="Float32" Name="mobW" NumberOfComponents="1" format="ascii">
946.105 946.107 946.105 946.107 946.107 946.107 946.104 946.104 946.093 946.093 946.059 946.059
945.981 945.981 945.825 945.825 945.548 945.548 945.104 945.104 944.452 944.452 943.563 943.563
942.432 942.432 940.895 940.895 939.113 939.113 937.185 937.185 935.189 935.189 933.206 933.206
931.309 931.309 929.553 929.553 927.981 927.981 926.613 926.613 925.456 925.456 924.501 924.501
923.731 923.731 923.126 923.126 922.659 922.659 922.307 922.307 922.046 922.046 921.856 921.856
921.72 921.72 921.624 921.624 921.558 921.558 921.512 921.512 921.482 921.482 921.462 921.462
945.981 945.981 945.825 945.825 945.548 945.548 945.105 945.105 944.453 944.453 943.564 943.564
942.433 942.433 940.896 940.896 939.114 939.114 937.185 937.185 935.19 935.19 933.206 933.206
931.308 931.308 929.553 929.553 927.98 927.98 926.612 926.612 925.455 925.455 924.5 924.5
923.731 923.731 923.125 923.125 922.659 922.659 922.306 922.306 922.045 922.045 921.855 921.855
921.719 921.719 921.624 921.624 921.558 921.558 921.512 921.512 921.482 921.482 921.462 921.462
921.448 921.448 921.44 921.44 921.434 921.434 921.43 921.43 921.428 921.428 921.427 921.427
921.426 921.426 921.425 921.425 921.425 921.425 921.424 921.424 921.424 921.424 921.424 921.424
921.424 921.424 921.423 921.423 921.423 921.423 921.423 921.423 921.423 921.423 921.423 921.423
......@@ -167,7 +167,7 @@
291 291 291 291 291 291 291 291 291 291 290.998 290.998
290.995 290.995 290.988 290.988 290.975 290.975 290.956 290.956 290.926 290.926 290.887 290.887
290.836 290.836 290.775 290.775 290.706 290.706 290.63 290.63 290.552 290.552 290.473 290.473
290.398 290.398 290.328 290.328 290.265 290.265 290.211 290.211 290.164 290.164 290.126 290.126
290.398 290.398 290.328 290.328 290.265 290.265 290.21 290.21 290.164 290.164 290.125 290.125
290.094 290.094 290.07 290.07 290.051 290.051 290.037 290.037 290.026 290.026 290.019 290.019
290.013 290.013 290.009 290.009 290.007 290.007 290.005 290.005 290.003 290.003 290.003 290.003
290.002 290.002 290.002 290.002 290.002 290.002 290.001 290.001 290.001 290.001 290.001 290.001
......@@ -179,47 +179,79 @@
290.001 290.001 290.001 290.001 290.001 290.001 290.001 290.001 290.001 290.001 290.001 290.001
290.001 290.001 290.001 290.001 290 290
</DataArray>
<DataArray type="Float32" Name="water content" NumberOfComponents="1" format="ascii">
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4
</DataArray>
<DataArray type="Float32" Name="source" NumberOfComponents="1" format="ascii">
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0
</DataArray>
<DataArray type="Float32" Name="velocity" NumberOfComponents="3" format="ascii">
0.0001 -3.44191e-19 0 0.0001 -5.16287e-19 0 0.0001 -3.44191e-19 0 0.0001 -5.16287e-19 0
0.0001 0 0 0.0001 0 0 0.0001 -3.44191e-19 0 0.0001 -3.44191e-19 0
0.0001 3.44187e-19 0 0.0001 3.44187e-19 0 0.0001 0 0 0.0001 0 0
0.0001 -1.54866e-18 0 0.0001 -1.54866e-18 0 0.0001 0 0 0.0001 0 0
9.99999e-05 -1.71994e-19 0 9.99999e-05 -1.71994e-19 0 9.99998e-05 1.71913e-19 0 9.99998e-05 1.71913e-19 0
9.99996e-05 1.54615e-18 0 9.99996e-05 1.54615e-18 0 9.99993e-05 -1.71633e-19 0 9.99993e-05 -1.71633e-19 0
9.99989e-05 1.71427e-19 0 9.99989e-05 1.71427e-19 0 9.99984e-05 2.05377e-18 0 9.99984e-05 2.05377e-18 0
9.9998e-05 1.36659e-18 0 9.9998e-05 1.36659e-18 0 9.99975e-05 1.36378e-18 0 9.99975e-05 1.36378e-18 0
9.9997e-05 -3.4022e-19 0 9.9997e-05 -3.4022e-19 0 9.99964e-05 1.69749e-19 0 9.99964e-05 1.69749e-19 0
9.99959e-05 -1.69404e-19 0 9.99959e-05 -1.69404e-19 0 9.99954e-05 1.18359e-18 0 9.99954e-05 1.18359e-18 0
9.9995e-05 -1.01279e-18 0 9.9995e-05 -1.01279e-18 0 9.99946e-05 3.371e-19 0 9.99946e-05 3.371e-19 0
9.99942e-05 6.73357e-19 0 9.99942e-05 6.73357e-19 0 9.99939e-05 1.68166e-18 0 9.99939e-05 1.68166e-18 0
9.99937e-05 -1.51223e-18 0 9.99937e-05 -1.51223e-18 0 9.99935e-05 -2.01499e-18 0 9.99935e-05 -2.01499e-18 0
9.99934e-05 1.67831e-19 0 9.99934e-05 1.67831e-19 0 9.99933e-05 -3.35533e-19 0 9.99933e-05 -3.35533e-19 0
9.99932e-05 -3.35438e-19 0 9.99932e-05 -3.35438e-19 0 9.99932e-05 -5.03054e-19 0 9.99932e-05 -5.03054e-19 0
9.99931e-05 -3.3532e-19 0 9.99931e-05 -3.3532e-19 0 9.99931e-05 -6.7057e-19 0 9.99931e-05 -6.7057e-19 0
9.99931e-05 -8.38152e-19 0 9.99931e-05 -8.38152e-19 0 9.99931e-05 -1.67622e-19 0 9.99931e-05 -1.67622e-19 0
9.99931e-05 -1.67617e-19 0 9.99931e-05 -1.67617e-19 0 9.99931e-05 -3.35226e-19 0 9.99931e-05 -3.35226e-19 0
9.99931e-05 5.02831e-19 0 9.99931e-05 5.02831e-19 0 9.99931e-05 1.67609e-19 0 9.99931e-05 1.67609e-19 0
9.99931e-05 -3.35216e-19 0 9.99931e-05 -3.35216e-19 0 9.99932e-05 1.67607e-19 0 9.99932e-05 1.67607e-19 0
9.99932e-05 -8.38034e-19 0 9.99932e-05 -8.38034e-19 0 9.99932e-05 -5.0282e-19 0 9.99932e-05 -5.0282e-19 0
9.99932e-05 1.00564e-18 0 9.99932e-05 1.00564e-18 0 9.99932e-05 -1.00564e-18 0 9.99932e-05 -1.00564e-18 0
9.99932e-05 1.67606e-19 0 9.99932e-05 1.67606e-19 0 9.99932e-05 -8.3803e-19 0 9.99932e-05 -8.3803e-19 0
9.99932e-05 0 0 9.99932e-05 0 0 9.99933e-05 1.84367e-18 0 9.99933e-05 1.84367e-18 0
9.99933e-05 -1.17324e-18 0 9.99933e-05 -1.17324e-18 0 9.99933e-05 1.67606e-18 0 9.99933e-05 1.67606e-18 0
9.99933e-05 3.35212e-19 0 9.99933e-05 3.35212e-19 0 9.99933e-05 1.67606e-19 0 9.99933e-05 1.67606e-19 0
9.99933e-05 1.67606e-19 0 9.99933e-05 1.67606e-19 0 9.99933e-05 -1.67606e-19 0 9.99933e-05 -1.67606e-19 0
9.99933e-05 5.02817e-19 0 9.99933e-05 5.02817e-19 0 9.99934e-05 -5.02817e-19 0 9.99934e-05 -5.02817e-19 0
9.99934e-05 6.70423e-19 0 9.99934e-05 6.70423e-19 0 9.99934e-05 -6.70423e-19 0 9.99934e-05 -6.70423e-19 0
9.99934e-05 1.34085e-18 0 9.99934e-05 1.34085e-18 0 9.99934e-05 0 0 9.99934e-05 0 0
9.99934e-05 1.67606e-18 0 9.99934e-05 1.67606e-18 0 9.99934e-05 1.67606e-19 0 9.99934e-05 1.67606e-19 0
9.99934e-05 -3.35211e-19 0 9.99934e-05 -3.35211e-19 0 9.99935e-05 -1.50845e-18 0 9.99935e-05 -1.50845e-18 0
9.99935e-05 1.67606e-19 0 9.99935e-05 1.67606e-19 0 9.99935e-05 6.70422e-19 0 9.99935e-05 6.70422e-19 0
9.99935e-05 -1.67605e-19 0 9.99935e-05 -1.67605e-19 0 9.99935e-05 1.50845e-18 0 9.99935e-05 1.50845e-18 0
9.99935e-05 -1.17324e-18 0 9.99935e-05 -1.17324e-18 0 9.99935e-05 -1.67605e-19 0 9.99935e-05 -1.67605e-19 0
9.99936e-05 -1.67605e-18 0 9.99936e-05 -1.67605e-18 0 9.99936e-05 1.67605e-19 0 9.99936e-05 1.67605e-19 0
9.99936e-05 -5.02816e-19 0 9.99936e-05 -5.02816e-19 0 9.99936e-05 -1.67605e-18 0 9.99936e-05 -1.67605e-18 0
9.99936e-05 -6.70421e-19 0 9.99936e-05 -6.70421e-19 0 9.99936e-05 1.00563e-18 0 9.99936e-05 1.00563e-18 0
9.99936e-05 -5.02816e-19 0 9.99936e-05 -5.02816e-19 0 9.99936e-05 0 0 9.99936e-05 0 0
9.99937e-05 0 0 9.99937e-05 0 0 9.99937e-05 0 0 9.99937e-05 0 0
0.0001 -3.44191e-19 0 0.0001 0 0 0.0001 -3.44191e-19 0 0.0001 0 0
0.0001 -1.37677e-18 0 0.0001 -1.37677e-18 0 0.0001 -1.72095e-18 0 0.0001 -1.72095e-18 0
0.0001 -1.72093e-18 0 0.0001 -1.72093e-18 0 0.0001 3.44174e-19 0 0.0001 3.44174e-19 0
0.0001 6.88292e-19 0 0.0001 6.88292e-19 0 0.0001 5.16134e-19 0 0.0001 5.16134e-19 0
9.99999e-05 3.43988e-19 0 9.99999e-05 3.43988e-19 0 9.99998e-05 1.20339e-18 0 9.99998e-05 1.20339e-18 0
9.99996e-05 3.4359e-19 0 9.99996e-05 3.4359e-19 0 9.99993e-05 1.37307e-18 0 9.99993e-05 1.37307e-18 0
9.99989e-05 6.8571e-19 0 9.99989e-05 6.8571e-19 0 9.99984e-05 1.19804e-18 0 9.99984e-05 1.19804e-18 0
9.9998e-05 3.41648e-19 0 9.9998e-05 3.41648e-19 0 9.99975e-05 1.70473e-19 0 9.99975e-05 1.70473e-19 0
9.9997e-05 0 0 9.9997e-05 0 0 9.99964e-05 3.39499e-19 0 9.99964e-05 3.39499e-19 0
9.99959e-05 -1.69404e-19 0 9.99959e-05 -1.69404e-19 0 9.99954e-05 1.52176e-18 0 9.99954e-05 1.52176e-18 0
9.9995e-05 -1.68799e-18 0 9.9995e-05 -1.68799e-18 0 9.99946e-05 -1.85405e-18 0 9.99946e-05 -1.85405e-18 0
9.99942e-05 -8.41696e-19 0 9.99942e-05 -8.41696e-19 0 9.99939e-05 0 0 9.99939e-05 0 0
9.99937e-05 -1.34421e-18 0 9.99937e-05 -1.34421e-18 0 9.99935e-05 3.35831e-19 0 9.99935e-05 3.35831e-19 0
9.99934e-05 1.67831e-18 0 9.99934e-05 1.67831e-18 0 9.99933e-05 1.67767e-19 0 9.99933e-05 1.67767e-19 0
9.99932e-05 1.34175e-18 0 9.99932e-05 1.34175e-18 0 9.99932e-05 -1.50916e-18 0 9.99932e-05 -1.50916e-18 0
9.99931e-05 0 0 9.99931e-05 0 0 9.99931e-05 1.67642e-19 0 9.99931e-05 1.67642e-19 0
9.99931e-05 1.6763e-19 0 9.99931e-05 1.6763e-19 0 9.99931e-05 -1.67622e-18 0 9.99931e-05 -1.67622e-18 0
9.99931e-05 -2.51425e-18 0 9.99931e-05 -2.51425e-18 0 9.99931e-05 -1.17329e-18 0 9.99931e-05 -1.17329e-18 0
9.99931e-05 -1.6761e-19 0 9.99931e-05 -1.6761e-19 0 9.99931e-05 1.67609e-19 0 9.99931e-05 1.67609e-19 0
9.99931e-05 1.67608e-19 0 9.99931e-05 1.67608e-19 0 9.99932e-05 0 0 9.99932e-05 0 0
9.99932e-05 3.35214e-19 0 9.99932e-05 3.35214e-19 0 9.99932e-05 -1.67607e-19 0 9.99932e-05 -1.67607e-19 0
9.99932e-05 1.50846e-18 0 9.99932e-05 1.50846e-18 0 9.99932e-05 8.38031e-19 0 9.99932e-05 8.38031e-19 0
9.99932e-05 0 0 9.99932e-05 0 0 9.99932e-05 -1.34085e-18 0 9.99932e-05 -1.34085e-18 0
9.99932e-05 -1.67606e-19 0 9.99932e-05 -1.67606e-19 0 9.99933e-05 1.67606e-18 0 9.99933e-05 1.67606e-18 0
9.99933e-05 1.67606e-19 0 9.99933e-05 1.67606e-19 0 9.99933e-05 1.34085e-18 0 9.99933e-05 1.34085e-18 0
9.99933e-05 -5.02818e-19 0 9.99933e-05 -5.02818e-19 0 9.99933e-05 -2.01127e-18 0 9.99933e-05 -2.01127e-18 0
9.99933e-05 -1.34085e-18 0 9.99933e-05 -1.34085e-18 0 9.99933e-05 -1.67606e-19 0 9.99933e-05 -1.67606e-19 0
9.99934e-05 1.67606e-19 0 9.99934e-05 1.67606e-19 0 9.99934e-05 -8.38029e-19 0 9.99934e-05 -8.38029e-19 0
9.99934e-05 -6.70423e-19 0 9.99934e-05 -6.70423e-19 0 9.99934e-05 1.67606e-19 0 9.99934e-05 1.67606e-19 0
9.99934e-05 1.67606e-19 0 9.99934e-05 1.67606e-19 0 9.99934e-05 1.67606e-19 0 9.99934e-05 1.67606e-19 0
9.99934e-05 1.50845e-18 0 9.99934e-05 1.50845e-18 0 9.99934e-05 0 0 9.99934e-05 0 0
9.99935e-05 1.00563e-18 0 9.99935e-05 1.00563e-18 0 9.99935e-05 3.35211e-19 0 9.99935e-05 3.35211e-19 0
9.99935e-05 -1.67606e-19 0 9.99935e-05 -1.67606e-19 0 9.99935e-05 -3.0169e-18 0 9.99935e-05 -3.0169e-18 0
9.99935e-05 0 0 9.99935e-05 0 0 9.99935e-05 6.70422e-19 0 9.99935e-05 6.70422e-19 0
9.99935e-05 -1.50845e-18 0 9.99935e-05 -1.50845e-18 0 9.99936e-05 -3.35211e-19 0 9.99936e-05 -3.35211e-19 0
9.99936e-05 -1.17324e-18 0 9.99936e-05 -1.17324e-18 0 9.99936e-05 0 0 9.99936e-05 0 0
9.99936e-05 -5.02816e-19 0 9.99936e-05 -5.02816e-19 0 9.99936e-05 8.38026e-19 0 9.99936e-05 8.38026e-19 0
9.99936e-05 -1.67605e-18 0 9.99936e-05 -1.67605e-18 0 9.99936e-05 -1.67605e-19 0 9.99936e-05 -1.67605e-19 0
9.99936e-05 5.02816e-19 0 9.99936e-05 5.02816e-19 0 9.99937e-05 1.67605e-19 0 9.99937e-05 1.67605e-19 0
9.99937e-05 2.17887e-18 0 9.99937e-05 2.17887e-18 0 9.99937e-05 0 0 9.99937e-05 0 0
9.99937e-05 0 0 9.99937e-05 0 0
</DataArray>
<DataArray type="Float32" Name="temperatureExact" NumberOfComponents="1" format="ascii">
......
......@@ -102,6 +102,24 @@
290.001 290.001 290.001 290.001 290.001 290.001 290.001 290.001 290.001 290.001 290.001 290.001
290.001 290.001 290.001 290.001 290.001 290.001 290.001 290.001
</DataArray>
<DataArray type="Float32" Name="water content" NumberOfComponents="1" format="ascii">
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
0.4 0.4 0.4 0.4 0.4 0.4 0.4 0.4
</DataArray>
<DataArray type="Float32" Name="source" NumberOfComponents="1" format="ascii">
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0 0 0 0 0
0 0 0 0 0 0 0 0
</DataArray>
<DataArray type="Float32" Name="velocity" NumberOfComponents="3" format="ascii">
0.0001 0 0 0.0001 0 0 0.0001 0 0 0.0001 0 0
0.0001 0 0 0.0001 0 0 0.0001 0 0 0.0001 0 0
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment