CFL condition (centroid-midpt distances versus inscribed circle or proper error estimator)

Gradient limiters: (Reconcile Matt's optimisation with quantity based limiters. E.g. reuse same C-code)

New forcing terms:
 * Friction of Mud
 * Viscosity