We're writing fits files with pyfits and would like to include a few bit masks (in particular "is saturated" and "ignore this pixel"). We're a bit stumped on how to do it and would appreciate any advice. -- Russell