iraf.net

I'm trying to add two images while removing cosmic rays in the process. If I use crreject in imcombine in IMAGES the pixels in the final image will only have half the pixel value as the surrounding pixels if there was a cosmic ray in one of the images as pixels were just used from one image. I've tried using ccdclip to reject the cosmic ray pixels and the corresponding pixels in the other image and set blank to twice the mode. This works for cosmic rays that aren't against any object but it doesn't work if it overlaps an object. I've also tried using cosmicrays in NOAO.IMRED.CRUTIL but the program just runs and runs without doing anything. Ideally I would like to use the pixel values from the image without the cosmic ray twice so that the combined image doesn't have any spots with low pixel values.I'm using IRAF 2.14 on MAC OS 10.6

I discovered if I average the two images using crreject and then multiply the pixels in the combined image by 2, I get the combined image that I want.