Overview

Tom Barlow has written an automated data reduction tool for reducing 2-D echelle data such as produced by HIRES and ESI. This tool, called ``makee,'' allows the user to input the file names for the object, flat, wavelength cal., bias, and a brighter trace star. and outputs a partially reduced set of 1-D spectra, one for each order on the detector.

The MAKEE manual can be found here.

ESI-specific notes

Install Location
Makee is installed at HQ for all ESI numbered accounts. Our install is for Solaris; however, it does not appear to run in Solaris10, only Sol6-Sol9. Your SA can give you the name of an appropriate machine.
Syntax Help
makee -help
[96] esi11@kealoha: makee -help
 Syntax: makee (object file) (star file) (flat file)
                     [arc file]  [2nd arc file]
  
    [arc1add=(FITS file)]  [arc2add=(FITS file)]  [arcmode=]
    [mode=]  [raw=]  [bias=(file)]  [-nobias]  [lpr=]
    [eperdn=]  [ronoise=]  [pf=(file)]  [order=n,n,..] [comm=s]
    [sl=] [hw=] [-atm] [-novac] [-noskyshift] [skyshift=]
    [-showmask]  [object=]  [observer=]  [log=]
    [oko=n,n,...]  [-nox]  [uop=]
  
    [-options] (show description of options)
    [-help   ] (show description of options)
  
  object file : Object FITS file to be reduced.
  star file   : Bright star FITS file to define trace.
  flat file   : Flat field FITS file.
  arc file    : Arclamp FITS file (optional).
  2nd arc file: Second Arclamp FITS file (optional).
  
  mode=      : 1 = Full extraction (default).
               2 = Mask CRs but do straight sum of object.
               3 = Mask CRs in sky only,straight sum of object.
               4 = Median of whole order.
  
  raw=dir   : Directory path for raw FITS files.
             (You can also set environment variable MAKEERAW,
                or do "ln -s (directory name) raw".)
  bias=file : FITS filename of bias (short dark) image.
  -nobias   : Do not apply bias correction.
  pf=file   : Use this parameter file instead of "makee.param".
 order=n,n,.. : Only extract orders n,n,.. where the bluest
 oko=n,n,..   : Only use these orders for trace centroids.
 uop=file     : Load User Object Positions data file.
       -nox : No Extraction (exit after profile computation).
   -showmask: Set pixels in masked (bad) regions to zero in
              reduced spectra. These regions will always appear
              as -1.0 in the variance spectra.
   comm=s   : Add the comment "s" to each output spectrum.
   sl=v     : Give estimate of slit length in pixels.
   hw=v     : Give object half width in pixels explicitly.
   lpr=     : Send output to this printer.
  -lesspr   : Send less output to printer (just Flux-*).
  -atm      : Apply atmospheric extinction correction.
  -novac    : No vacuum wavelength scale correction (air wave.)
 -noskyshift: Do not apply automated skyline shifted scale.
  skyshift=v: Override automated skyline wavelength scale 
              shift, set shift to v pixels.
   object=  : Reset "OBJECT" name to this string (~=blank).
  observer= : Reset "OBSERVER" name to this string (~=blank).
     log=   : Filename for log file for output messages.
  arc1add=  : Arclamp exposure to provide additional lines to
              wavelength calibration of first arclamp.
  arc2add=  : Arclamp exposure to provide additional lines to
              wavelength calibration of second arclamp.
  arcmode=  : Arclamp wavelength scale averaging mode:
                 0 : Average scales of 1st,2nd arclamps(def.).
                 1 : Shift scale to 1st arclamp.
                 2 : Shift scale to 2nd arclamp.
   eperdn=  : Inverse gain (electrons per digital number).
   ronoise= : Readout noise in electrons.
    { If "eperdn=" or "ronoise=" is not given, use (1) values
      in header, (2) values in "eperdn.dat" and "ronoise.dat",
      or (3) default values.
  
 The following are usually only used for testing purposes:
   -profpoly : Write out profile polynomial fit plots.
   -test     : Write out prior/post masking test images.
   -notwhole : Read ee.get file for trace information.
   -trace    : General testing...
              MAKEE VERSION= 3.4.0 : April 2001                      

Work Flow (Inputs and Outputs)
 Object Image    Star (trace) Image     Flat Field 
        |                    |                  |
        |                    |                  |
        |   ThAr Arc Image   |      Short Dark  |
        |    (optional)      |      (optional)  |
        |         |          |           |      |
        |         |          |           |      |
        |         |          |           |      |
        _____________________|___________________
                             |
                             |
                             |
                             |
                           makee
                             |
                             |
                             |
     ________________________|
     |
     |-- Flux-###.fits (2-D spectrum file of object in units
     |      of DN/(flat), vacuum,heliocentric wavelengths.)
     |
     |-- Err-###.fits  (1 sigma error spectrum)
     |-- Sky-###.fits  (night sky spectrum)
     |-- s2n-###.fits  (Signal-to-noise per pixel)
     |-- DN-###.fits  (object spectrum in digital number units)
     |-- Flat-###.fits  (non-normalized spectrum of the flat field)
     |
     |-- Flux-###.ps (PostScript plot).
     |-- s2n-###.ps  (PostScript plot).
     |-- DN-###.ps  (PostScript plot).
     |
     |-- trace-###_###.ps (plot of echelle trace fits)
     |-- profiles_###.ps  (plot of object profiles)
     |
     |-- Arc-###.fits   (arc lamp spectrum)
     |-- Arc-###-cf.ps  (arc correlation function)
     |-- Arc-###-fr.ps  (arc final residuals (plot))
     |-- Arc-###.res    (arc final residuals (table))
     |-- Arc-###.ids    (arc line identifications)
     |-- Arc-###-sky.ps (night sky line shift fit)

         where "###" is the observation number.
Wavelength Calibration
ESI wavelength calibration is typically done using three lamp sets: HgNe, XE, and CuAr. MAKEE can only accept two wavelength calibration files. You can take the HgNe and Xe spectra on the same image, since their intensities are fairly well-matched, or if you have the spectra in separate files, add the Xe image to the HgNe image.
Slit Length
The echellette slits are about 134 pixels in length. We find makee runs more reliably if you give a sl=134 argument in your makee call.
Multiple Objects
Multiple objects can be handled by specifying the extraction window using "uop=". See the uop help file for details.
Extended Objects
Use the "hw=P" argument to specify an extraction window, where "P" is in pixels
Example
makee e160213_0076 e160213_0048 FLAT ARC sl=134
Screen Output
............. MAKEE: Keck ESI Reduction .......... Sat Feb 13 06:02:19 2016

MAKEE VERSION= 3.4.0 : April 2001

 OBJECT:     0  2231     0  4095   ./e160213_0076.fits  Object   0.75
   STAR:     0  2231     0  4095   ./e160213_0048.fits  Object   0.75
   FLAT:     0  2231     0  4095   ./FLAT.fits  DmFlat   0.75
    ARC:     0  2231     0  4095   ./AARC.fits  Line     0.75
 Setting EPERDN =   1.2900
 Setting RONOISE =   2.7000
  Column binning=  1    Row binning=  1
 Frame#=   76
1443.3+0006  Expos= 1200s   Object=Star 27 1200s

Object to extract:
Star 27 1200s         1440+0019  16/02/13 15:19:26   76  1200s
Binning:   1x1        File: e160213_0076.fits

Mode=  1   Extract all orders.
Masking file: MaskESI_1x1_2.dat
   Bias file: /sol/apps/makee/makee//biasESI_1x1_2amp.fits
............
Reading bias : /sol/apps/makee/makee//biasESI_1x1_2amp.fits
Using 1 dimensional processed bias image.
Bias image observation date: 2000-10-21
............
Reading flat : ./FLAT.fits
 Setting EPERDN =   1.2900
 Setting RONOISE =   2.7000
  Column binning=  1    Row binning=  1
 Frame#=   18
0352.0+4500  Expos=  200s   Object=Dome Flat
NOTE: Using special interpolation for ESI.
NOTE: BIW:  X binning:   1   Y binning:   1
Subtracting baseline ( 2072 --> 2151 ) from data (   24 --> 1047 ).
Subtracting baseline ( 2152 --> 2231 ) from data ( 1048 --> 2071 ).
Rotating ESI image by 90 degrees...
Flipping ESI image in columns...
Prepix:   12   Postpix:   80
Original Window:  cols=   0  2048  rows=   0  4096
New Clipped Window:  cols=   0  4094  rows=   0  2045
Subtracting bias from image...
............
Using slit length from command line.
Estimate for slit length =  134.000
Saving flat image... junk-flat_1455379346000076001.fits
............
Reading star : ./e160213_0048.fits
 Setting EPERDN =   1.2900
 Setting RONOISE =   2.7000
  Column binning=  1    Row binning=  1
 Frame#=   48
0801.4+1627  Expos=    3s   Object=HR3134 Hot Rotator
NOTE: Using special interpolation for ESI.
NOTE: BIW:  X binning:   1   Y binning:   1
Subtracting baseline ( 2072 --> 2151 ) from data (   24 --> 1047 ).
Subtracting baseline ( 2152 --> 2231 ) from data ( 1048 --> 2071 ).
Rotating ESI image by 90 degrees...
Flipping ESI image in columns...
Prepix:   12   Postpix:   80
Original Window:  cols=   0  2048  rows=   0  4096
New Clipped Window:  cols=   0  4094  rows=   0  2045
Subtracting bias from image...
Finding and fitting order traces...
Writing PostScript plot file: trace-048.ps
NOTE: EETO: EON=  1  High=  4.414  WRMS=  0.664  WPPD=  0.512  NREJ=   0
NOTE: EETO: EON=  2  High=  1.070  WRMS=  0.154  WPPD=  0.288  NREJ=  72
NOTE: EETO: EON=  3  High=  0.702  WRMS=  0.101  WPPD=  0.247  NREJ=  53
NOTE: EETO: EON=  4  High=  0.638  WRMS=  0.091  WPPD=  0.227  NREJ=  71
NOTE: EETO: EON=  5  High=  0.399  WRMS=  0.057  WPPD=  0.185  NREJ=  94
NOTE: EETO: EON=  6  High=  0.219  WRMS=  0.032  WPPD=  0.157  NREJ=  68
NOTE: EETO: EON=  7  High=  0.144  WRMS=  0.021  WPPD=  0.141  NREJ=  12
NOTE: EETO: EON=  8  High=  0.115  WRMS=  0.019  WPPD=  0.130  NREJ=   1
NOTE: EETO: EON=  9  High=  0.124  WRMS=  0.028  WPPD=  0.131  NREJ=   0
NOTE: EETO: EON= 10  High=  4.165  WRMS=  0.665  WPPD=  0.625  NREJ=   0
Number of echelle orders =  10
Recovering flat image... junk-flat_1455379346000076001.fits
Estimating background limits...
NOTE: Background between widely separated orders:  1.65382E+02
NOTE: EEFB: EON=  1  Slit Edge Lower=  86.4   Upper= -89.5   Width= 175.9
NOTE: EEFB: EON=  2  Slit Edge Lower=  81.2   Upper= -83.3   Width= 164.5
NOTE: EEFB: EON=  3  Slit Edge Lower=  76.8   Upper= -81.2   Width= 158.0
NOTE: EEFB: EON=  4  Slit Edge Lower=  73.0   Upper= -78.3   Width= 151.3
NOTE: EEFB: EON=  5  Slit Edge Lower=  69.8   Upper= -75.7   Width= 145.5
NOTE: EEFB: EON=  6  Slit Edge Lower=  66.9   Upper= -73.7   Width= 140.6
NOTE: EEFB: EON=  7  Slit Edge Lower=  64.4   Upper= -71.9   Width= 136.3
NOTE: EEFB: EON=  8  Slit Edge Lower=   0.0   Upper=   0.0   Width=   0.0
NOTE: EEFB: EON=  9  Slit Edge Lower=   0.0   Upper=   0.0   Width=   0.0
NOTE: EEFB: EON= 10  Slit Edge Lower=  57.7   Upper= -66.0   Width= 123.7
EEFB: Fitting background lower offsets.
EEFB: Fitting background upper offsets.
 Flat field scale factor=   3.778533E+03
Saving flat image... junk-flat_1455379346000076001.fits
............
Reading object : ./e160213_0076.fits
 Setting EPERDN =   1.2900
 Setting RONOISE =   2.7000
  Column binning=  1    Row binning=  1
 Frame#=   76
1443.3+0006  Expos= 1200s   Object=Star 27 1200s
NOTE: Using special interpolation for ESI.
NOTE: BIW:  X binning:   1   Y binning:   1
Subtracting baseline ( 2072 --> 2151 ) from data (   24 --> 1047 ).
Subtracting baseline ( 2152 --> 2231 ) from data ( 1048 --> 2071 ).
Rotating ESI image by 90 degrees...
Flipping ESI image in columns...
Prepix:   12   Postpix:   80
Original Window:  cols=   0  2048  rows=   0  4096
New Clipped Window:  cols=   0  4094  rows=   0  2045
Subtracting bias from image...
 Observation number=   76   UT date = 13Feb16 
Masking bad regions...
NOTE: Object x,y binning is   1 x  1 .
Recovering flat image... junk-flat_1455379346000076001.fits
Dividing object by flat field...
 Determine object profile and boundaries...
Writing FITS profile file: profiles-076.fits
Writing PostScript plot file: profiles-076.ps
 slit width, neo, eperdn, ronoise, pse_mode = 134.00  10  1.290  2.700  1
****************************** Extracting echelle order number  1 of 10
Cosmic Ray Fraction: expected= 0.00093  actual= 0.00106
Sky Pixels Masked (primary,adjacent,additional):   289   186  1827
 Sky areas ratio:   1.078.   Sky polynomial order:  1.
                Object Pixels Masked (high,low):   203    26
PSE_ObjSum: Had to deal with   256 masked pixels.
****************************** Extracting echelle order number  2 of 10
Cosmic Ray Fraction: expected= 0.00093  actual= 0.00095
Sky Pixels Masked (primary,adjacent,additional):   326   210  2240
 Sky areas ratio:   1.100.   Sky polynomial order:  1.
                Object Pixels Masked (high,low):   277    17
PSE_ObjSum: Had to deal with   328 masked pixels.
****************************** Extracting echelle order number  3 of 10
Cosmic Ray Fraction: expected= 0.00093  actual= 0.00120
Sky Pixels Masked (primary,adjacent,additional):   343   285  2713
 Sky areas ratio:   1.122.   Sky polynomial order:  1.
                Object Pixels Masked (high,low):   149    10
PSE_ObjSum: Had to deal with   213 masked pixels.
****************************** Extracting echelle order number  4 of 10
Cosmic Ray Fraction: expected= 0.00093  actual= 0.00111
Sky Pixels Masked (primary,adjacent,additional):   287   264  3480
 Sky areas ratio:   1.144.   Sky polynomial order:  1.
                Object Pixels Masked (high,low):   213     7
PSE_ObjSum: Had to deal with   279 masked pixels.
****************************** Extracting echelle order number  5 of 10
Cosmic Ray Fraction: expected= 0.00093  actual= 0.00093
Sky Pixels Masked (primary,adjacent,additional):   298   190  5140
    38 columns affected by nearby night sky lines during SkyFit.
 Sky areas ratio:   1.166.   Sky polynomial order:  1.
                Object Pixels Masked (high,low):   216     4
PSE_ObjSum: Had to deal with   285 masked pixels.
****************************** Extracting echelle order number  6 of 10
Cosmic Ray Fraction: expected= 0.00093  actual= 0.00110
Sky Pixels Masked (primary,adjacent,additional):   282   196  5601
   109 columns affected by nearby night sky lines during SkyFit.
 Sky areas ratio:   1.186.   Sky polynomial order:  1.
                Object Pixels Masked (high,low):   392    11
PSE_ObjSum: Had to deal with   475 masked pixels.
****************************** Extracting echelle order number  7 of 10
Cosmic Ray Fraction: expected= 0.00093  actual= 0.00092
Sky Pixels Masked (primary,adjacent,additional):   229   173  5302
   187 columns affected by nearby night sky lines during SkyFit.
 Sky areas ratio:   1.205.   Sky polynomial order:  1.
                Object Pixels Masked (high,low):   182     7
PSE_ObjSum: Had to deal with   277 masked pixels.
****************************** Extracting echelle order number  8 of 10
Cosmic Ray Fraction: expected= 0.00093  actual= 0.00110
Sky Pixels Masked (primary,adjacent,additional):   252   207  4044
   664 columns affected by nearby night sky lines during SkyFit.
 Sky areas ratio:   1.221.   Sky polynomial order:  1.
                Object Pixels Masked (high,low):   145     5
PSE_ObjSum: Had to deal with   206 masked pixels.
****************************** Extracting echelle order number  9 of 10
Cosmic Ray Fraction: expected= 0.00093  actual= 0.00150
Sky Pixels Masked (primary,adjacent,additional):   342   247  3434
   817 columns affected by nearby night sky lines during SkyFit.
 Sky areas ratio:   1.233.   Sky polynomial order:  1.
                Object Pixels Masked (high,low):   213     4
PSE_ObjSum: Had to deal with   278 masked pixels.
****************************** Extracting echelle order number 10 of 10
Cosmic Ray Fraction: expected= 0.00093  actual= 0.00102
Sky Pixels Masked (primary,adjacent,additional):   227   157  2680
   649 columns affected by nearby night sky lines during SkyFit.
 Sky areas ratio:   1.239.   Sky polynomial order:  1.
                Object Pixels Masked (high,low):   210     2
PSE_ObjSum: Had to deal with   261 masked pixels.
Writing : Flux-076.fits
Writing : Var-076.fits
Writing : Sky-076.fits
Writing : Sum-076.fits
Sat Feb 13 06:06:39 2016

.....  Extracting Flat Field Spectrum  ...........

............
Reading object : ./FLAT.fits
 Setting EPERDN =   1.2900
 Setting RONOISE =   2.7000
  Column binning=  1    Row binning=  1
 Frame#=   18
0352.0+4500  Expos=  200s   Object=Dome Flat
NOTE: Using special interpolation for ESI.
NOTE: BIW:  X binning:   1   Y binning:   1
Subtracting baseline ( 2072 --> 2151 ) from data (   24 --> 1047 ).
Subtracting baseline ( 2152 --> 2231 ) from data ( 1048 --> 2071 ).
Rotating ESI image by 90 degrees...
Flipping ESI image in columns...
Prepix:   12   Postpix:   80
Original Window:  cols=   0  2048  rows=   0  4096
New Clipped Window:  cols=   0  4094  rows=   0  2045
Subtracting bias from image...
 Observation number=   18   UT date = 13Feb16 
Masking bad regions...
NOTE: Object x,y binning is   1 x  1 .
Recovering flat image... junk-flat_1455379346000076001.fits
Dividing object by flat field...
 slit width, neo, eperdn, ronoise, pse_mode = 134.00  10  1.290  2.700  4
Extracting whole slit spectrum...
****************************** Extracting echelle order number  1 of 10
****************************** Extracting echelle order number  2 of 10
****************************** Extracting echelle order number  3 of 10
****************************** Extracting echelle order number  4 of 10
****************************** Extracting echelle order number  5 of 10
****************************** Extracting echelle order number  6 of 10
****************************** Extracting echelle order number  7 of 10
****************************** Extracting echelle order number  8 of 10
****************************** Extracting echelle order number  9 of 10
****************************** Extracting echelle order number 10 of 10
Writing : Flat-018.fits
Sat Feb 13 06:06:52 2016

.....  Extracting Arc Lamp Spectrum  ...........

............
Reading object : ./AARC.fits
 Setting EPERDN =   1.2900
 Setting RONOISE =   2.7000
  Column binning=  1    Row binning=  1
 Frame#=   13
0352.0+4500  Expos=  600s   Object=HgNeXeCuAr Lamp
NOTE: Using special interpolation for ESI.
NOTE: BIW:  X binning:   1   Y binning:   1
Subtracting baseline ( 2072 --> 2151 ) from data (   24 --> 1047 ).
Subtracting baseline ( 2152 --> 2231 ) from data ( 1048 --> 2071 ).
Rotating ESI image by 90 degrees...
Flipping ESI image in columns...
Prepix:   12   Postpix:   80
Original Window:  cols=   0  2048  rows=   0  4096
New Clipped Window:  cols=   0  4094  rows=   0  2045
Subtracting bias from image...
 Observation number=   13   UT date = 13Feb16 
Masking bad regions...
NOTE: Object x,y binning is   1 x  1 .
Recovering flat image... junk-flat_1455379346000076001.fits
Dividing object by flat field...
 slit width, neo, eperdn, ronoise, pse_mode = 134.00  10  1.290  2.700  4
Extracting whole slit spectrum...
****************************** Extracting echelle order number  1 of 10
****************************** Extracting echelle order number  2 of 10
****************************** Extracting echelle order number  3 of 10
****************************** Extracting echelle order number  4 of 10
****************************** Extracting echelle order number  5 of 10
****************************** Extracting echelle order number  6 of 10
****************************** Extracting echelle order number  7 of 10
****************************** Extracting echelle order number  8 of 10
****************************** Extracting echelle order number  9 of 10
****************************** Extracting echelle order number 10 of 10
Writing : Arc-013.fits
Sat Feb 13 06:07:04 2016
Compute wavelength scale for Arc-013.fits
arc2obj  Arc-013.fits  Flux-076.fits  -sky  -var  ps=Flux-076-skyshift.ps  -sum

Reading: Arc-013.fits
Arc file= Arc-013.fits                            
NOTE: Skylines from: /sol/apps/makee/makee//SkyLineDataBase/sky.results.ESI
Reading: Flux-076.fits
Second file= Flux-076.fits                           
Reading: Var-076.fits
Reading: Sum-076.fits
Reading: Sky-076.fits
*****Night sky line shift results:
#   row  wavelength  peak   pix.shft weight
    1    5   5589.128    27.6    0.37   0.28
    2    5   5889.953   482.4    0.37   1.00
    3    5   5895.923   301.1    0.53   1.00
    4    5   5915.301    60.8    0.65   0.61
    5    5   5924.723    34.5    0.75   0.34
    6    5   5932.862    47.2    0.57   0.47
    7    5   5945.932    27.4    0.89   0.27
    8    5   5953.420    39.7    1.12   0.40
    9    6   5889.953   482.4    0.17   1.00
   10    6   5895.923   301.1    0.13   1.00
   11    6   5915.301    60.8   -0.29   0.61
   12    6   5924.723    34.5   -0.17   0.34
   13    6   5932.862    47.2    0.34   0.47
   14    6   5945.932    27.4   -0.84   0.27
   15    6   5953.420    39.7    0.26   0.40
   16    6   5970.278    32.6    0.95   0.33
   17    6   6257.961    69.7    0.00   0.70
   18    6   6287.434    62.0    0.21   0.62
   19    6   6297.903    31.0    0.23   0.31
   20    6   6300.304   673.7    0.10   1.00
   21    6   6477.910    38.0    0.07   0.38
   22    6   6498.729   107.4    0.46   1.00
   23    6   6504.995    27.6    0.53   0.28
   24    6   6522.425    62.2    1.46   0.00
   25    6   6544.022    59.8    1.06   0.60
   26    6   6553.617    63.2    0.89   0.63
   27    7   6257.961    69.7   -1.55   0.00
   28    7   6287.434    62.0   -0.22   0.62
   29    7   6297.903    31.0    8.98   0.00
   30    7   6300.304   673.7   -0.01   1.00
   31    7   6477.910    38.0   -0.70   0.38
   32    7   6498.729   107.4   -0.21   1.00
   33    7   6504.995    27.6    0.26   0.28
   34    7   6522.425    62.2    0.97   0.62
   35    7   6544.022    59.8   -0.30   0.60
   36    7   6553.617    63.2   -0.07   0.63
   37    7   6568.779    62.0    0.05   0.62
   38    7   6596.643    37.6   -0.46   0.38
   39    7   6841.945    40.7   -0.35   0.41
   40    7   6863.955   255.9   -0.08   1.00
   41    7   6870.994    29.4    0.09   0.29
   42    7   6900.833    63.8    0.03   0.64
   43    7   6969.930   101.8   -0.03   1.00
   44    7   7253.127   109.6    0.60   1.00
   45    7   7276.405   492.0    0.27   1.00
   46    7   7284.439   159.0    0.83   1.00
   47    8   7253.127   109.6    0.18   1.00
   48    8   7276.405   492.0   -0.30   1.00
   49    8   7284.439   159.0    0.38   1.00
   50    8   7316.282   379.4    0.07   1.00
   51    8   7329.148   252.3    0.10   1.00
   52    8   7340.885   424.4    0.10   1.00
   53    8   7358.659   290.6    0.00   1.00
   54    8   7392.198   308.4   -0.01   1.00
   55    8   7429.781   134.5   -0.17   1.00
   56    8   7530.713   336.9   -0.01   1.00
   57    8   7557.492   240.3    0.20   1.00
   58    8   7571.746   443.1    0.11   1.00
   59    8   7586.093   238.1    0.04   1.00
   60    8   7759.996   410.1    0.16   1.00
   61    8   7808.467   526.5    0.13   1.00
   62    8   7821.503  1119.7    0.09   1.00
   63    8   7841.266   844.4    0.02   1.00
   64    8   7860.658  1044.8    0.08   1.00
   65    8   7870.730  1196.6    0.17   1.00
   66    8   8014.059  1149.3    0.26   1.00
   67    8   8052.020  1072.3    0.02   1.00
   68    9   8014.059  1149.3    0.19   1.00
   69    9   8052.020  1072.3   -0.16   1.00
   70    9   8288.603  1786.7   -0.23   1.00
   71    9   8298.896  1956.7   -0.20   1.00
   72    9   8382.392  1366.9   -0.10   1.00
   73    9   8493.389   932.7    0.27   1.00
   74    9   8903.114  1558.9    0.18   1.00
   75    9   8919.634  2112.7    0.14   1.00
Skyline shift=   0.1441  +/-  0.3817  Npt=  72  NRej=  3
Writing to device: Flux-076-skyshift.ps/PS
 Applying skyline shift of    0.1441 pixels to arclamp scale.
Sky Shift used=    0.1441   measured=    0.1441   RMS=  0.3817  #Pt=  72
 Heliocentric velocity correction (km/s) =  28.1344
Air to vacuum correction applied...
Heliocentric velocity correction applied...
Over-writing file: Flux-076.fits
Over-writing file: Sum-076.fits
Over-writing file: Sky-076.fits
Writing file: Err-076.fits
arc2obj: All done.
.....................................................

arc2arc  Flux-076.fits  Sum-076.fits  pxsh=0
First file= Flux-076.fits                           
Second file= Sum-076.fits                            
 Shift in pixels=   0.000000
 Zero shift, no refitting necessary.
arc2arc: All done.

opim  mul  Flux-076.fits  Flat-018.fits  DN-076.fits  -clobber
  Writing file : DN-076.fits
opim  div  DN-076.fits  const=+3.7785327E+03  -clobber
  Writing file : DN-076.fits
spim0  DN-076.fits  min=0.  bin=2  ps=DN-076.ps
 Output is in : DN-076.ps

opim  div  Flux-076.fits  Err-076.fits  s2n-076.fits  -clobber
  Writing file : s2n-076.fits
opim  min  s2n-076.fits  const=0.  -clobber
  Writing file : s2n-076.fits
spim0  s2n-076.fits  min=0.  bin=2  ps=s2n-076.ps
 Output is in : s2n-076.ps

spim0  Flux-076.fits  minp=-0.08  bin=2  ps=Flux-076.ps
 Output is in : Flux-076.ps
NOTE: spectrum plot binning factor is 2 pixels.

rhh  Flux-076.fits  -wave
  76  Flux-076 13Feb16 15:19  Ec:Obje       0.75  1x1  1200s  Star 27 1200s (Flu

Sat Feb 13 06:07:05 2016
makee: All done.