fit_2d_vf function

Estimate a 2D vector field