MIRSI DATA REDUCTION STEPS

     -------------------------------------------------------------------------------
     Each step is performed by creating and executing a batch file.  Each batch file 
     requires one or more string vectors of filenames, as well as various function
     calls and keywords.
     -------------------------------------------------------------------------------

     0. TO START PIPELINE: 
          a. Download, untar, and unzip program files from 
               http://irtfweb.ifa.hawaii.edu/research/dr_resources/
          b. Download Astro Users Library from http://idlastro.gsfc.nasa.gov/homepage.html
          c. Add library paths to those searched by IDL
          d. Execute astrolib.pro in IDL (or put in startup file)
          e. Create subdirectories for each step
                    - Originals
                    - Coadded (step 1)
                    - Registration1 (step 2)
                    - Registration2 (step 3)
                    - Results (steps 4 and 5)
                    - Batches

     1. TO COADD CHOP/NOD QUADS:
          a. Create batch file with following form:                  
                 infiles = ['file1.fits', 'file2.fits', ...]
                 outfiles = ['output1.fits', 'output2.fits', ...]
                 paths = ['/input_path/', '/output_path/']
                 n = n_elements(infiles)
                 for i = 0, (n-1) do $ 
                      mirsi_coad, infiles[i], outfiles[i], paths
          b. If reducing extended objects, set ext keyword
          c. Execute batch file 

     2.  TO REGISTER COADDED QUADS OF SAME OBJECT: 
          a. Create batch file with following form:                 
                 infiles = ['coadded_file1.fits',  'coadded_file2.fits', ...]
                        (these should be files that you want to register)
                 outfile = 'output_name'
                 paths = ['/input_path/', '/output_path/']
                 mirsi_register, infiles, outfile, paths, step=1
          b. If a dither pattern was used, set the dither keyword
                 dither = [x1, y1, x2, y2, ...]
          c. Set plot keyword to display results, if desired
          d. Repeat a-c for each different object
          e. Execute batch file 
          f. If registration is not satisfactory
                i. Try adjusting maximum shift amount by setting mshft
                   keyword equal to a larger or smaller value
               ii. Enter approximate peak positions with peaks keyword
          g. Steps 2 and 3 can be accomplished in one step by not setting
             step keyword in a

     3.  TO REGISTER AN OBJECT'S CHOP/NOD IMAGES:
          a. Create batch file with following form:  
                infiles = ['file1', 'file2', ...]
                outfile = ['outfile1', 'outfile2', ...]
                paths = ['/input_path/', '/output_path/']
                n = n_elements(infiles)
                for i = 0, (n-1) do $ 
                   mirsi_register, infiles, outfile, paths, step=2
          b. Set plot keyword to display results, if desired
          c. Execute batch file 

     4.  TO CREATE A FILE WITH FLUXES AT A RANGE OF APERTURES:
            a. Create batch file with following form:
                img = ['reg_image1', 'reg_image2', ...]
                outfile = ['data_i1', 'data_i2', ...]
                path = [/img_path/', '/outfile_path/']
                aper = [minimum, maximum, # of divisions]
                skyrad = [inner sky annulus, outer sky annulus]
                n = n_elements(img)
                for i = 0, (n-1) do $ 
                     flux, img[i], outfile[i], path, aper, skyrad
            b. If multiple objects are in the field then specify 
               coordinates of desired object with peak keyword
            c. Execute batch file

     5. TO SELECT AN APERTURE AND OBTAIN RESULTS
            a. Create batch file with following form:
                img = ['data_i1', 'data_i2', ...]
                outfile = 'output_file'
                path = ['/img_path/', '/outfile_path/']
                aper = desired aperture (in pixels)
                phot_results, img, outfile, aper, path=path
            b. Execute batch file for standard stars
            c. Calculate zero point correction and extinction 
               coefficient from standard star results
            d. Create batch file as in (a) for objects
            e. Set amc keyword to [zero pt, extinction]
            f. Execute batch file for objects