Add upper and lower bounds to wavefront fit parameters

[aringlis]

It would be useful to be able to define upper and lower bounds for each of the wavefront fitting parameters (at the moment these are the amplitude, position and width of a Gaussian). This would help prevent returning bad/unphysical fit results.

Add this functionality to fit_wavefront, within eitwaveutils.