Use the flat-fielded orders to extract the stamps
Transfer bad pixels (CR removal) into extracted orders
Make sure there is only on sigma image in FITS file

Write a separate file (with extension "n") from find-pixel-noise.py,
rather than modifying in place.