function velmap

Go through a datacube and fit a gaussian to a line. The arguments center and fwhm are used as initial estimate for the gauss-fit.

Syntax
result = velmap(cube, center, fwhm [, threshold])

Arguments
cube:   A datacube
center:   The approximate center (wavelength) of the line to fit
fwhm:   The approximate FWHM of the line to fit
threshold:   A lower limit under which no fit is attempted (given as percentage of peak flux)


Returns
A datacube with the following slices: result[*,*,1]: fit background result[*,*,2]: fit continuum result[*,*,3]: fit line center result[*,*,4]: fit line FWHM result[*,*,5]: error in background fit result[*,*,6]: error in continuum fit result[*,*,7]: error in line center fit result[*,*,8]: error in line FWHM fit