Abstract. We introduce a pure-streamfunction formulation for the incompressible Navier-Stokes equations. The idea is to replace the vorticity in the vorticity- streamfunction evolution equation by the Laplacianof the streamfunction. The resulting formulation includes the streamfunction only, thus no inter-function relations need to invoked. A compact numerical scheme, which interpolates streamfunction values as well as its first order derivatives, is presented and analyzed.