Add density-dependent gravity

jedbrown · May 27, 2011 · 9a6d48f · 9a6d48f
1 parent f6dbee1
commit 9a6d48f
Show file tree

Hide file tree

Showing 2 changed files with 1 addition and 2 deletions.
diff --git a/src/fs/tests/vht.c b/src/fs/tests/vht.c
@@ -181,7 +181,7 @@ static dErr VHTCaseSetFromOptions(VHTCase scase)
     err = dOptionsRealUnits("-rheo_T0","Reference temperature (corresponds to enthalpy=0)","",u->Temperature,rheo->T0,&rheo->T0,NULL);dCHK(err);
     err = dOptionsRealUnits("-rheo_T3","Triple point temperature","",u->Temperature,rheo->T3,&rheo->T3,NULL);dCHK(err);
     err = dOptionsRealUnits("-rheo_splice_delta","Characteristic width of split","",u->Energy,rheo->splice_delta,&rheo->splice_delta,NULL);dCHK(err);
-    err = dOptionsRealUnits("-gravity","Gravity","",u->Acceleration,scase->gravity,&scase->gravity,NULL);dCHK(err);
+    err = dOptionsRealUnits("-rheo_gravity","Gravity in z-direction (probably negative)","",u->Acceleration,rheo->gravity,&rheo->gravity,NULL);dCHK(err);
     if (scase->setfromoptions) {err = (*scase->setfromoptions)(scase);dCHK(err);}
   } err = PetscOptionsEnd();dCHK(err);
   dFunctionReturn(0);

diff --git a/src/fs/tests/vhtimpl.h b/src/fs/tests/vhtimpl.h
@@ -59,7 +59,6 @@ struct _n_VHTCase {
   dErr (*setfromoptions)(VHTCase);
   dErr (*destroy)(VHTCase);
   struct VHTRheology rheo;
-  dReal gravity;
   dReal bbox[3][2];
   dBool reality; // The "solution" is just a guess or boundary conditions
   dUnits units;