Add gravity() method.
Moved some calculations to computeExplicitData() method.
Fix initialization issue with gpot_.
Removed qs from the set of primary variables.