From a1c07502ecb7030454c116d13ee63aac913c4bb2 Mon Sep 17 00:00:00 2001
From: Qiao Lei <qiaol618@gmail.com>
Date: Fri, 13 Feb 2015 14:20:57 -0600
Subject: [PATCH] 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.
---
 examples/step-33/step-33.cc | 4 +---
 1 file changed, 1 insertion(+), 3 deletions(-)

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;
-- 
2.39.5