[Doc] Update units of solute

[ci-skip]

doc/introduction/features.rst

@@ -42,7 +42,7 @@ Solute transport Solver
 =======================
 
 The solute transport equation for unsaturated media describes movement of 
-solute by the *solute concentration* :math:`c_w`,
+solute by the *solute concentration* :math:`c_w \, [\mathrm{kg}/\mathrm{m}^3]`,
 
 .. math::
 
@@ -51,7 +51,7 @@ solute by the *solute concentration* :math:`c_w`,
               \right]
    = 0,
 
-where :math:`D \, []` is the effective hydrodynamic dispersion tensor, and :math:`\vec{j}_w` is the water flux. Currently, :math:`D` can only be used as a constant parameter by the keyword ``effHydrDips`` in the configuration file.
+where :math:`D \, [\mathrm{m}^2/\mathrm{s}]` is the effective hydrodynamic dispersion tensor, and :math:`\vec{j}_w` is the water flux. Currently, :math:`D` can only be used as a constant parameter by the keyword ``effHydrDips`` in the configuration file.
 
 .. TODO: Update effective hydrodynamic dispersion description once parameterization is merged in master.
 

doc/manual/bcfile.rst

@@ -37,20 +37,20 @@ Neumann
 
   .. object:: neumann <solute_flux>
 
-  Set a fixed volumetric solute flux [:math:`s^{-1}`] at the boundary. Positive values 
+  Set a fixed volumetric solute flux [:math:`kg/(s\,m^3)`] at the boundary. Positive values 
   imply fluxes out of the domain, negative values imply fluxes into the domain.
 
 
-Neumann
+Outflow
 +++++++
 
   .. object:: outflow <solute_outflux>
 
   Set a spatial rate of change of volumetric solute flux 
-  [:math:`s^{-1}m^{-1}`] at the boundary. Positive values imply higher rate of 
-  change out of the domain, negative values imply lower rate of change into 
+  [:math:`kg/(s\,m\,m^3)`] at the boundary. Positive values imply higher rate 
+  of change out of the domain, negative values imply lower rate of change into 
   the domain. Influx situations are not taken into account. For most of the 
-  cases an outflow of 0. [:math:`s^{-1}m^{-1}`] is a reasonable situation.
+  cases an outflow of 0. [:math:`kg/(s\,m\,m^3)`] is a reasonable situation.
 
 Datafile Structure
 ==================

doc/manual/initial.rst

@@ -84,7 +84,7 @@ Initial condition tranformations for the Transport solver.
     * ``quantity = soluteConcentration``
 
     The input data is directly interpreted as solute concentration,
-    :math:`f = c_w`.
+    :math:`f = c_w [\mathrm{kg}/\mathrm{m}^3]`.
 
 
 .. _muparser: http://beltoforion.de/article.php?a=muparser&p=features