[cig-commits] [commit] master: Using introspection() instead of hard coded in (7102f80)

[cig_noreply at geodynamics.org]

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

On branch  : master
Link       : https://github.com/geodynamics/aspect/compare/4a5e27e5690a0e0c1e13e19100af7398e37f36b8...d31f4e9435a1c12781f5b673b672cbd29c41c2c2

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

commit 7102f805377b11dbc63219276de212032052cd70
Author: Siqi Zhang <siqi.zhang at mq.edu.au>
Date:   Thu May 22 10:45:55 2014 -0500

    Using introspection() instead of hard coded in


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

7102f805377b11dbc63219276de212032052cd70
 source/postprocess/visualization/internal_heating.cc | 10 +++++-----
 1 file changed, 5 insertions(+), 5 deletions(-)

diff --git a/source/postprocess/visualization/internal_heating.cc b/source/postprocess/visualization/internal_heating.cc
index bf20ccf..b913dfc 100644
--- a/source/postprocess/visualization/internal_heating.cc
+++ b/source/postprocess/visualization/internal_heating.cc
@@ -56,17 +56,17 @@ namespace aspect
         const unsigned int n_quadrature_points = uh.size();
         Assert (computed_quantities.size() == n_quadrature_points,    ExcInternalError());
         Assert (computed_quantities[0].size() == 1,                   ExcInternalError());
-        Assert (uh[0].size() == dim+2+this->n_compositional_fields(), ExcInternalError());
-        Assert (duh[0].size() == dim+2+this->n_compositional_fields(),ExcInternalError());
+        Assert (uh[0].size() == this->introspection().n_components, ExcInternalError());
+        Assert (duh[0].size() == this->introspection().n_components,ExcInternalError());
 
         for (unsigned int q=0; q<n_quadrature_points; ++q)
           {
-            double temperature=uh[q][dim+1];
-            double pressure   =uh[q][dim];
+            double temperature=uh[q][this->introspection().component_indices.temperature];
+            double pressure   =uh[q][this->introspection().component_indices.pressure];
             std::vector<double> composition(this->n_compositional_fields());
 
             for (unsigned int c=0; c<this->n_compositional_fields(); ++c)
-              composition[c] = uh[q][dim+2+c];
+              composition[c] = uh[q][this->introspection().component_indices.compositional_fields[c]];
             computed_quantities[q](0) = heating_model.specific_heating_rate(temperature,
                                                                             pressure,
                                                                             composition,