IDL> im_sum = im1+im2+im3
IDL> writefits,'sum_image.fits',im_sum
IDL> writefits,'sum_image.fits',im_sum,hdr