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 c3162bed519f8a9512442fa5e98c0cbf08651a70..6ae98203def4c2bc844c91ab08ee645952c684f0 100644 (file)
--- 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;