Figures improved.
Understand the spectral approx of 'div(u,w)=0' problem with no (interior) constraints pretty well.  Should look at: (1) add u_z=0 at top, (2) write up test function argument that shows -u_z+w_x=0 so u,w harmonic in interior, (3) build my own svd-based underdetermined solution.