This should be nvfrac, not nvfrac%mesh.
Included the heat transfer term in the MMS test.
fixed a load of bugs
added some notes in the schema about which mesh a SinkingVelocity field should be on
Added SinkingVelocity mesh checking to check_options routines of:
Advection_Diffusion_CG.F90
Advection_Diffusion_FV.F90
Field_Equations_CV.F90

Added check_options routine in Advection_Diffusion_DG.F90 with check for SinkingVelocity mesh.

not tested yet.
merge with trunk
fixed bugs in options checking routines - these now work as planned
updated test flml files
Adding PETSs 3.3 support and dropping support for 3.0.
We need to grab the harmonic frequencies each time we do an update now, not just once
fixed bug that made all runs with scalar fields, but no sinking velocity, fail
merge with trunk
better considered condition on vb slope limiter to avoid NaNs
fixed error in schema description
fixed various bugs. Main issue was only doing boundary integral on external face with a dirichlet boundary condition.

now also solved with inverse of mass matrix
- Added an options check to make sure the prognostic PhaseVolumeFraction is only in the particle phase.
- Added options "include_heat_flux_term" and "include_pressure_term" under the InternalEnergy equation_type.
Included the pressure*div(u) term in the MMS test.
Included a check for include_pressure_term in the CV code.
changed mms test schema so that bedload is not calculated