This merge adds a pressure/density splitting technique to help to prevent spurious oscillations in the pressure field when using unbalanced finite element pairs.

When the subtract_out_reference_profile option (under the compressible equation_of_state) is enabled, the Density and Pressure fields are split up into a hydrostatic component (') and a perturbed component (''). The hydrostatic components, denoted p' and rho', satisfy the balance: grad(p') = rho'*g. The hydrostatic components (defined by the user) are then subtracted from the pressure and density used in the pressure gradient and buoyancy terms in the momentum equation. This works with both CG and DG Velocity discretisations.

This merge also includes: single-phase and multiphase tests for this functionality, a description in the manual and schema, and some extra sanity checks regarding the application of pressure boundary conditions and inlet velocity boundary conditions for the ctjacobs-multiphase flow model.