From: Qiao Lei <qiaol618@gmail.com>
Date: Fri, 13 Feb 2015 20:20:57 +0000 (-0600)
Subject: Fix bug in energy source term according to Praveen C's report.
X-Git-Tag: v8.3.0-rc1~464^2
X-Git-Url: https://gitweb.dealii.org/cgi-bin/gitweb.cgi?a=commitdiff_plain;h=refs%2Fpull%2F546%2Fhead;p=dealii.git

Fix bug in energy source term according to Praveen C's report.

    The energy source term due to gravity should be rho*g*v=g*momentum.
---

diff --git a/examples/step-33/step-33.cc b/examples/step-33/step-33.cc
index c3162bed51..6ae98203de 100644
--- a/examples/step-33/step-33.cc
+++ b/examples/step-33/step-33.cc
@@ -314,9 +314,7 @@ namespace Step33
             forcing[c] = gravity * W[density_component];
             break;
           case energy_component:
-            forcing[c] = gravity *
-                         W[density_component] *
-                         W[first_momentum_component+dim-1];
+            forcing[c] = gravity * W[first_momentum_component+dim-1];
             break;
           default:
             forcing[c] = 0;