From fc7d7b2003f529bce65659772c6fa642ea2938bd Mon Sep 17 00:00:00 2001 From: bangerth Date: Tue, 4 Oct 2011 18:58:33 +0000 Subject: [PATCH] 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 --- deal.II/examples/step-32/step-32.cc | 4 ++-- 1 file changed, 2 insertions(+), 2 deletions(-) diff --git a/deal.II/examples/step-32/step-32.cc b/deal.II/examples/step-32/step-32.cc index f1ad6a2718..b689c959d6 100644 --- 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