[cig-commits] [commit] master: fix typo in formula (429e886)

[cig_noreply at geodynamics.org]

Repository : https://github.com/geodynamics/aspect

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/492d395e19a72325f7828c78924713be8ff5aa40...de71435b7d188bdffde1ace3918c3c23a76c7138

>---------------------------------------------------------------

commit 429e886f920911ee99fc4ec238b28fd4654455f7
Author: Timo Heister <timo.heister at gmail.com>
Date:   Sun May 25 20:11:15 2014 -0400

    fix typo in formula


>---------------------------------------------------------------

429e886f920911ee99fc4ec238b28fd4654455f7
 doc/manual/parameters.tex       | 2 +-
 source/material_model/simple.cc | 2 +-
 2 files changed, 2 insertions(+), 2 deletions(-)

diff --git a/doc/manual/parameters.tex b/doc/manual/parameters.tex
index 69a09f1..3ef6f73 100644
--- a/doc/manual/parameters.tex
+++ b/doc/manual/parameters.tex
@@ -2616,7 +2616,7 @@ In order to facilitate placing input files in locations relative to the ASPECT s
 
 `simple': A material model that has constant values for all coefficients but the density and viscosity. The defaults for all coefficients are chosen to be similar to what is believed to be correct for Earth's mantle. All of the values that define this model are read from a section ``Material model/Simple model'' in the input file, see Section~\ref{parameters:Material_20model/Simple_20model}.
 
-This model uses the following set of equations for the two coefficients that are non-constant: \begin{align}  \eta(p,T,\mathfrak c) &= \tau(T) \zeta(\mathfrak c) \eta_0, \\  \rho(p,T,\mathfrak c) &= \left(1-\alpha (T-T_0)\right)\rho_0 + \Deltaho \; c_0,\end{align}where $c_0$ is the first component of the compositional vector $\mathfrak c$ if the model uses compositional fields, or zero otherwise. 
+This model uses the following set of equations for the two coefficients that are non-constant: \begin{align}  \eta(p,T,\mathfrak c) &= \tau(T) \zeta(\mathfrak c) \eta_0, \\  \rho(p,T,\mathfrak c) &= \left(1-\alpha (T-T_0)\right)\rho_0 + \Delta\rho \; c_0,\end{align}where $c_0$ is the first component of the compositional vector $\mathfrak c$ if the model uses compositional fields, or zero otherwise. 
 
 The temperature pre-factor for the viscosity formula above is defined as \begin{align}  \tau(T) &= H\left(e^{\beta (T-T_0)/T_0}\right),  \qquad\qquad H(x) = \left\{                         \begin{cases}                            10^{-2} & \text{if}\; x<10^{-2}, \\                            x & \text{if}\; 10^{-2}\le x \le 10^2, \\                            10^{2} & \text{if}\; x>10^{2}, \\                         \end{cases}                       \right.\end{align} where $\beta$ corresponds to the input parameter ``Thermal viscosity exponent'' and $T_0$ to the parameter ``Reference temperature''. If you set $T_0=0$ in the input file, the thermal pre-factor $\tau(T)=1$.
 
diff --git a/source/material_model/simple.cc b/source/material_model/simple.cc
index ae78915..b654495 100644
--- a/source/material_model/simple.cc
+++ b/source/material_model/simple.cc
@@ -344,7 +344,7 @@ namespace aspect
                                    "are non-constant: "
                                    "\\begin{align}"
                                    "  \\eta(p,T,\\mathfrak c) &= \\tau(T) \\zeta(\\mathfrak c) \\eta_0, \\\\"
-                                   "  \\rho(p,T,\\mathfrak c) &= \\left(1-\\alpha (T-T_0)\\right)\\rho_0 + \\Delta\rho \\; c_0,"
+                                   "  \\rho(p,T,\\mathfrak c) &= \\left(1-\\alpha (T-T_0)\\right)\\rho_0 + \\Delta\\rho \\; c_0,"
                                    "\\end{align}"
                                    "where $c_0$ is the first component of the compositional vector "
                                    "$\\mathfrak c$ if the model uses compositional fields, or zero otherwise. "