;ngc1097_preparations.pro

;Starting from Pipeline 1 files:
   image='ngc1097_v7.phot.1.fits'
   mask='mask_ngc1097.fits'
   wt='ngc1097_v7.phot.1_wt.fits'

;Make sigma-image, new data image, and new mask  with NaNs removed
;These are used in the GALFIT inputfile ngc1097_galfit30_pasp_v3.ingal

   sigma_fits='ngc1097_v7.phot.1_sigma.fits'
   image_new='ngc1097_v7.phot.1_nonan.fits'
   mask_new='mask_ngc1097_nonan.fits'

;check the sigma-image, comparing calculated sigma image with 
;the stdev estimated from image itself

   check=[751.978-1,739.894-1,10]
   fudge=2

make_sigma_image_new, image, mask, wt, sigma_fits,$
                      image_new=image_new,mask_new=mask_new,$
                      check=check,fudge=fudge

end