mkscr.cl

subtract A B C [D] : C=A-B[+D]
"subtract" command can be applied to both data taken by IRS1 and IRS2.
Since this command do nothing if the output image exists, you have to
remove the output data before if you want to recalculate it.

skusub A B C : C=Ax + B(1-x)
"skusub" command merges and optimizes two images to minimize the
residual sky emission.

!../subtract FMSA00004039.fits FMSA00004045.fits irs1_sample1a+1_lr.fits
!../subtract FMSA00004039.fits FMSA00004041.fits irs1_sample1a_lr.fits
!../subtract FMSA00004043.fits FMSA00004041.fits irs1_sample1b-1_lr.fits
!../subtract FMSA00004043.fits FMSA00004045.fits irs1_sample1b_lr.fits
!../subtract FMSA00004047.fits FMSA00004045.fits irs1_sample1c-1_lr.fits
!../subtract FMSA00004047.fits FMSA00004049.fits irs1_sample1c_lr.fits
!../subtract FMSA00004051.fits FMSA00004049.fits irs1_sample1d-1_lr.fits
!../subtract FMSA00004047.fits FMSA00004053.fits irs1_sample1c+1_lr.fits
!../subtract FMSA00004051.fits FMSA00004053.fits irs1_sample1d_lr.fits
!../skysub irs1_sample1a_lr.fits irs1_sample1a+1_lr.fits irs1_sample1a1_lr.fits
!../skysub irs1_sample1a_lr.fits irs1_sample1b-1_lr.fits irs1_sample1a2_lr.fits
!../skysub irs1_sample1b_lr.fits irs1_sample1b-1_lr.fits irs1_sample1b1_lr.fits
!../skysub irs1_sample1b_lr.fits irs1_sample1c-1_lr.fits irs1_sample1b2_lr.fits
!../skysub irs1_sample1c_lr.fits irs1_sample1c-1_lr.fits irs1_sample1c1_lr.fits
!../skysub irs1_sample1c_lr.fits irs1_sample1d-1_lr.fits irs1_sample1c2_lr.fits
!../skysub irs1_sample1d_lr.fits irs1_sample1d-1_lr.fits irs1_sample1d1_lr.fits
!../skysub irs1_sample1d_lr.fits irs1_sample1c+1_lr.fits irs1_sample1d2_lr.fits