Improve robustness of nebular fitting

* fit-nebula.py:  New option --linear-components specifies which of the
  components should have their quadratic terms pegged at zero.
* fit_utils.py: Coefficients set to None get reset to a small floating
  value (zero caused problems with printing the relative uncertainties)
  and with "vary=False"