Version with vorticity equation. Boundary condition not yet implemented.