Added the prepare function for 4th order pressure solver
Implemented offsets in advection and diffusion, small bug remaining