Fix an actual bug: we were taking a completely meaningless temperature

field to evaluate the density in computing the rhs of the temperature
equation. Use the correctly extrapolated field.

git-svn-id: https://svn.dealii.org/trunk@24532 0785d39b-7218-0410-832d-ea1e28bc413d

diff --git a/deal.II/examples/step-32/step-32.cc b/deal.II/examples/step-32/step-32.cc
index f1ad6a271858cdb5d2331979a4fdde455b59bf2d..b689c959d6945a19f9bc8edafefa8cfc902180d7 100644 (file)
--- a/deal.II/examples/step-32/step-32.cc
+++ b/deal.II/examples/step-32/step-32.cc
@@ -3324,10 +3324,10 @@ namespace Step32
             scratch.old_strain_rates[q]);
 
        const double gamma
-         = ((EquationData::radiogenic_heating * EquationData::density(T_term_for_rhs)
+         = ((EquationData::radiogenic_heating * EquationData::density(ext_T)
              +
              2 * EquationData::eta * extrapolated_strain_rate * extrapolated_strain_rate) /
-            (EquationData::density(T_term_for_rhs) * EquationData::specific_heat));
+            (EquationData::density(ext_T) * EquationData::specific_heat));
 
        for (unsigned int i=0; i<dofs_per_cell; ++i)
          {