I don't usually post in the middle of the day, but my code is finally successful!

My flux values for the binary spectra will now be fixed if there is errors on account of cosmic rays!

This will currently work if there are a max of 4 pixels in a row with mask =1. I might increase this later, but otherwise we're good.