debugging units issues when calculating additional uncertainties

the model flux arrays are coming in without units, which is a problem
when adding/subtracting them to the data.  There's a fix in bdfit
right now, but it needs to be fixed in the model file itself
eventually
plotting +/- uncertainty, errorbar was having issues