Merged Julianes changes and mine from svn porosity branch with latest master branch revision. Solution now includes porosity and porosity velocity fields. Pororisyt velocity does nothing for the time being, porosity is advected the same as any other compositional field.