Several changes to handle standard FITS spectra files

1.  Added DC-FLAG handling
2.  Dealt with BZERO (maybe not quite right..)
3.  Added error file to read_fits.read_fits_spectrum1d
This is not a problem of this PR in particular, but in general - we should have an `extension` option to use in place of assuming the data are in hdu 0.