Add infrastructure for Gaussian fits

* plot_utils.py: Gaussian velocity components where the
  parameters (intensity, u0, sigma) are polynomial functions of position