Fix/averagefield var p
This MR looks to sort out the filter Average field which was not working when the polynomial space was different for different variables. This is often the case when running the Incompressible Navier-Stokes equations at high Reynolds number.