[cig-commits] commit: Fix a bug where gravity was not being used in BuoyancyForceTerm

[Mercurial]

changeset:   451:4853368a2fda
tag:         tip
user:        Walter Landry <wlandry at caltech.edu>
date:        Wed Jun 13 15:08:55 2012 -0700
files:       Utils/src/BuoyancyForceTerm.cxx
description:
Fix a bug where gravity was not being used in BuoyancyForceTerm


diff -r 9d6b2ba08b9d -r 4853368a2fda Utils/src/BuoyancyForceTerm.cxx
--- a/Utils/src/BuoyancyForceTerm.cxx	Mon May 21 15:46:35 2012 -0700
+++ b/Utils/src/BuoyancyForceTerm.cxx	Wed Jun 13 15:08:55 2012 -0700
@@ -537,7 +537,7 @@ void _BuoyancyForceTerm_CalcGravity(Buoy
                                     double* xi, double *gravity)
 {
   for(unsigned int d=0;d<mesh->generator->nDims;++d)
-    gravity[d]=self->gHat[d];
+    gravity[d]=self->gHat[d]*self->gravity;
 }
 
 void BuoyancyForceTerm_average_density_alpha(BuoyancyForceTerm *self,