Chebyshev both directions rectangle incompressible NSE

[Michael Jolly]

I'm trying to use V3 to solve the incompressible Navier-Stokes equations on a rectangle with velocity=0 on all 4 walls.   A test with a manufactured solution shows convergence except for pressure.  

# N timestep  err_u  rate_u  err_v  rate_v  err_p  rate_p   div   rate_div

  16 1.9e-06 3.6e-10        3.6e-10        7.9e+00        3.1e-17 

  32 4.9e-07 9.1e-11  2.0   9.1e-11  2.0   1.9e+02  -4.6  4.0e-16  -3.7

  64 1.2e-07 2.2e-11  2.0   2.2e-11  2.0   1.1e-01  10.7  9.5e-18   5.4

Has anyone been successful with Chebyshev in both directions for the NSE.   My script is attached.

Thanks,

Mike