1. Identify and download calibration data (biases and flux standard) from 
   Gemini Science Archive. (GSA Complete Catalog Query)
   http://www2.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/gsa/

- PI DATA PACKAGE RETRIEVAL to get science data, flats, procecessed bias, and logs.
- COMPLETE CATALOG: BIASES for flux std and DARKS for targets

   to DOWNLOAD DATA
  - save html file as file.html
  - use wget to retrieve
    wget --http-user=<user> --http-password=<pswd> -B http://www1.cadc-ccda.hia-iha.nrc-cnrc.gc.ca -F -i file1.html

to PRINT LOGS
a2ps -1r -f7 -b 

For calibration data:
- BIASES
   - NEW: archive provides already processed BIAS
   - OLD:
     - Find appropriate biases with same parameters as science data 
	CCDSEC == [1:2048,1:4608] 
	CCDSUM = 2 2
	- download all biases in date range of science data
 
- FLUX CALIBRATOR
   - find biases (need to have same binning and CCDSEC as flux obs)
       CCDSUM = 2 2
       CCDSEC == [1:2048,1792:2815]

- DARKS (For Nod and Shuffle data)
   - darks should match dark time and nodcount for objects
       NODCOUNT = 8
       NODPIX = 1536
       EXPTIME = 900..1000 
       CCDSUM = 2 2

2. Once everything is downloaded, move image frames to images/
  - mkdir images/
  - hselect > images.lis
    images = "*.fits[0]"
    fields = "$I"
    expr = grating ?= 'MIRROR'  && obstype != 'BIAS'
  - emacs images.lis to get of [0] (Cntl-x r k to use a rectangular kill)
  - imrename @images.lis images/

4. Make lists with different data types
  - OBJECTS
    hselect > objects.lis
    images = "*.fits[0]"
    fields = "$I"
    expr = obstype ?= 'OBJECT'

  - FLATS
    hselect > flats.lis
    images = "*.fits[0]"
    fields = "$I"
    expr = obstype ?= 'FLAT'

  - ARCS
    hselect > arcs.lis
    images = "*.fits[0]"
    fields = "$I"       
    expr = obstype ?= 'ARC'

  - BIASES (if one was not already processed for you)
    hselect > bias.lis 
    images = "*.fits[0]"
    fields = "$I"      
    expr =  obstype ?= 'BIAS' && nccds = 3

  - BIASES (for Nod and Shuffle data)
    hselect > dark.lis 
    images = "*.fits[0]"
    fields = "$I"      
    expr =  obstype ?= 'DARK' 

  --> don't forget to emacs *.lis to get rid of [0] (Cntl-X R K)

5. GPREPARE the data
   - OBJECTS
     inimages = @objects.lis
     rawpath = /data/hillary/4/kelle/optical_spectra_raw/Geminispec/GN_2004B_Q_10_2/RAW
     outpref = /data/hillary/2/kelle/optical_spectra/runs/gemini/GN_2004B/run2/gt
     fl_vard = no
     fl_addm = yes
     (Should get no warnings or errors)

   - FLATS
     inimages = @flats.lis
     rawpath = /data/hillary/4/kelle/optical_spectra_raw/Geminispec/GN_2004B_Q_10_2/RAW
     outpref = /data/hillary/2/kelle/optical_spectra/runs/gemini/GN_2004B/run2/gf
     fl_vard = no
     fl_addm = yes
     (Should get no warnings or errors)

   - ARCS
     inimages = @arcs.lis
     rawpath = /data/hillary/4/kelle/optical_spectra_raw/Geminispec/GN_2004B_Q_10_2/RAW
     outpref = /data/hillary/2/kelle/optical_spectra/runs/gemini/GN_2004B/run2/ga
     fl_vard = no
     fl_addm = yes
     (Should get no warnings or errors)

   - BIASES (for flats, arcs, and not Nod & Shuffle observations)
     inimages = @bais.lis
     outpref = /data/hillary/2/kelle/optical_spectra/runs/gemini/GN_2004B/run2/gb
     fl_addm = no (if set to yes, doesnt really matter, but might give error)
     (OK: WARNING - GPREPARE: ... use fl_addmdf=yes)
     (OK: WARNING - GPREPARE: ... OBSMODE=IMAGE)

  - DARKS
     inimages = @darks.lis
     rawpath = /data/hillary/4/kelle/optical_spectra_raw/Geminispec/GN_2004B_Q_10_2/RAW
     outpref = /data/hillary/2/kelle/optical_spectra/runs/gemini/GN_2004B/run2/gd
     fl_vard = no
     fl_addm = yes
     (Should get no warnings or errors)