MIRSI DATA REDUCTION PIPELINE STEPS

     -----------------------------------------------------------------
     Each step is performed by creating and executing a batch file.
     Each batch file requires one or more string vectors of filenames.
     This can be done by creating a file in the standard way (in UNIX 
     for example, "ls *.fits > batch_file.bat") or by using the 
     batch_sv.pro interface.  The function call for each step and 
     other non-string-vector steps should then be added to the created
     batch file.
     ----------------------------------------------------------------- 

     0. TO START PIPELINE: 
          a. Download and untar program files from 
               http://irtfweb.ifa.hawaii.edu/~elv
          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