There was a mix of issues in bilinear interpolation and pressure calculation but smoke is now being correctly advected. Will now code the advection for vx and vy!