It seems like we can do without this `_bandflux_errsnakesq` function, as the equivalent of 
`source._bandflux_errsnakesq(band, phase)` is `source._restframe_errsnakesq(band.wave_eff, phase)`