Amazing: Requiring the vector field to have zero divergence, known surface w, zero basal velocity already gives great result with standard use of least squares.  To demo: run noconstraint.py.  TODO: add isochrone-type info, improve (=speed up) matrix assembly.