1. Approximate spectrophotometry using absolute fluxes of Vega

a) Fit continuum to a 2D spectrum of Vega, select to save "fit" (type=fit)

input   = raw/df_vega_005-004.ec.fits  Input images
output  =        vega-fit.fits  Output images
(lines  =                    *) Image lines to be fit
(bands  =                    *) Image bands to be fit
(type   =                  fit) Type of output
(replace=                   no) Replace rejected points by fit?
(wavesca=                  yes) Scale the X axis with wavelength?
(logscal=                   no) Take the log (base 10) of both axes?
(overrid=                   no) Override previously fit lines?
(listonl=                   no) List fit but don't modify any images?
(logfile=              logfile) List of log files
(interac=                  yes) Set fitting parameters interactively?
(sample =                    *) Sample points to use in fit
(naverag=                    1) Number of points in sample averaging
(functio=              spline3) Fitting function
(order  =                    7) Order of fitting function
(low_rej=                   3.) Low rejection in sigma of fit
(high_re=                   0.) High rejection in sigma of fit
(niterat=                   10) Number of rejection iterations
(grow   =                   1.) Rejection growing radius
(markrej=                  yes) Mark rejected points?
(graphic=             stdgraph) Graphics output device
(cursor =                     ) Graphics cursor input
ask     =                  YES                                       
(mode   =                   ql)

b) Combine continuum fit of Vega using task scombine

input   =             vega-fit  List of input spectra
output  =          vega-fit-1D  List of output spectra
(noutput=                     ) List of output number combined spectra
(logfile=               STDOUT) Log file
(apertur=                     ) Apertures to combine
(group  =               images) Grouping option
(combine=                  sum) Type of combine operation
(reject =                 none) Type of rejection
(first  =                   no) Use first spectrum for dispersion?
(w1     =                4320.) Starting wavelength of output spectra
(w2     =                7100.) Ending wavelength of output spectra
(dw     =                  0.1) Wavelength increment of output spectra
(nw     =                INDEF) Length of output spectra
(log    =                   no) Logarithmic increments?
(scale  =                 none) Image scaling
(zero   =                 none) Image zero point offset
(weight =                 none) Image weights
(sample =                     ) Wavelength sample regions for statistics
(lthresh=                INDEF) Lower threshold
(hthresh=                INDEF) Upper threshold
(nlow   =                    1) minmax: Number of low pixels to reject
(nhigh  =                    1) minmax: Number of high pixels to reject
(nkeep  =                    1) Minimum to keep (pos) or maximum to reject
(neg)
(mclip  =                  yes) Use median in sigma clipping algorithms?
(lsigma =                   3.) Lower sigma clipping factor
(hsigma =                   3.) Upper sigma clipping factor
(rdnoise=                   0.) ccdclip: CCD readout noise (electrons)
(gain   =                   1.) ccdclip: CCD gain (electrons/DN)
(snoise =                   0.) ccdclip: Sensitivity noise (fraction)
(sigscal=                  0.1) Tolerance for sigma clipping scaling
corrections
(pclip  =                 -0.5) pclip: Percentile clipping parameter
(grow   =                    0) Radius (pixels) for 1D neighbor rejection
(blank  =                   0.) Value if there are no pixels

c) Combine 2D spectra to 1D spectra of all objects using combine=sum

d) Divide 1D spectra of objects by 1D fit to the Vega spectrum using sarith

input1  =          HD128167-1D  List of input spectra
op      =                    /  Operation
input2  =          vega-fit-1D  List of input spectra or constants
output  =        HD128167-temp  List of output spectra
(w1     =                INDEF) Starting wavelength
(w2     =                INDEF) Ending wavelength
(apertur=                     ) List of input apertures or columns/lines
(bands  =                     ) List of input bands or lines/bands
(beams  =                     ) List of input beams or echelle orders
(apmodul=                    0) Input aperture modulus (0=none)
(reverse=                   no) Reverse order of operands in binary operation?
(ignorea=                   no) Ignore second operand aperture numbers?
(format =            multispec) Output spectral format
(renumbe=                   no) Renumber output apertures?
(offset =                    0) Output aperture number offset
(clobber=                   no) Modify existing output images?
(merge  =                   no) Merge with existing output images?
(rebin  =                  yes) Rebin to exact wavelength region?
(errval =                   0.) Arithmetic error replacement value
(verbose=                   no) Print operations?
(mode   =                   ql)

e) Multiple the resulting spectrum by Vega abslute flux spectrum

input1  =        HD128167-temp  List of input spectra
op      =                    *  Operation
input2  =          vega-fluxes  List of input spectra or constants
output  =      HD128167-fluxes  List of output spectra
(w1     =                INDEF) Starting wavelength
(w2     =                INDEF) Ending wavelength
(apertur=                     ) List of input apertures or columns/lines
(bands  =                     ) List of input bands or lines/bands
(beams  =                     ) List of input beams or echelle orders
(apmodul=                    0) Input aperture modulus (0=none)

(reverse=                   no) Reverse order of operands in binary
operation?
(ignorea=                   no) Ignore second operand aperture numbers?

(format =            multispec) Output spectral format
(renumbe=                   no) Renumber output apertures?
(offset =                    0) Output aperture number offset
(clobber=                   no) Modify existing output images?
(merge  =                   no) Merge with existing output images?
(rebin  =                  yes) Rebin to exact wavelength region?
(errval =                   0.) Arithmetic error replacement value
(verbose=                   no) Print operations?
(mode   =                   ql)

f) Multiply the resulting spectrum by the ratio of exposure times:
   exposure(Vega)/exposure(target)