From: bangerth Date: Tue, 4 Oct 2011 18:58:33 +0000 (+0000) Subject: Fix an actual bug: we were taking a completely meaningless temperature X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=fc7d7b2003f529bce65659772c6fa642ea2938bd;p=dealii-svn.git 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 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