Commit ffb745d2 authored by Lukas Riedel's avatar Lukas Riedel

Add variable for solute mass to initial condition tip

Also sanitize the variable for the variance and fix the quantity of
the initial condition which depends on the spatial dimensions.
parent 8d6ea20f
......@@ -59,10 +59,13 @@ They are controlled by the ``initial.type`` key and available for every model.
.. tip::
The expression for a gaussian pulse of solute concentration centered at
:math:`\vec{p} = [0.5, 0.5]^T \, \mathrm{m}` and variance of
:math:`\sigma = 0.1^2 \, \mathrm{m}^2` is::
:math:`\vec{p} = [0.5, 0.5]^T \, \mathrm{m}` and variance of
:math:`\sigma^2 = \left( 0.1 \, \mathrm{m} \right)^2` is::
initial.equation = exp(-((x-0.5)^2+(y-0.5)^2)/(2.*0.1^2))/(2*pi*0.1^2)
initial.equation = <m> * exp(-((x-0.5)^2+(y-0.5)^2)/(2.*0.1^2)) / (2*pi*0.1^2)
where ``<m>`` is the total solute mass of the pulse
:math:`m_s \, [\text{kg}]`.
.. object:: Dataset
......@@ -109,7 +112,8 @@ Initial condition tranformations for the Transport solver.
* ``quantity = soluteConcentration``
The input data is directly interpreted as solute concentration,
:math:`f = c_w [\mathrm{kg}/\mathrm{m}^3]`.
:math:`f = c_w [\mathrm{kg}/\mathrm{m}^d]`, where :math:`d` indicates the
spatial dimensions of the grid.
.. _H5: https://www.h5py.org/
.. _muparser: http://beltoforion.de/article.php?a=muparser&p=features
Markdown is supported
0% or
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment