Archive Access Software

ESPEC

espec computes a spectrum's radial velocity by fitting individual emission line shifts and fitting those velocities weighted by their errors. The results are tabulated to the display and an optional log file, and a summary display labels the emission lines on the spectrum, also marking the positions of possible absorption lines. The following help message will appear on your terminal when you type espec, and a parameter file, espec.par, will be created in your current directory.

Parameters read from "/home/mink/Defaults/espec.par"
 ESPEC fits a radial velocity to emission lines in a
 binary spectrum archive file or files .  The range of
 the reduced file numbers (or a file containing ranges)
 is specified on the command line.
 In addition, several switches may be set:

        reduced file number of archive file (22827)
     or    reduced file numbers of archive file
     or @  file of archive reduced file numbers
  -d  sets archive root pathname (/ret)
  -inst  sets instrument directory name (zmach)
   -gz      sets initial Z for emission lines (   -1.00000)
   -gv      sets initial velocity for emission lines (   -299792.)
   -ng           zeroes initial velocity guess
   -efl(nefl)    fit (do not fit) emission line profiles
   -efv(nefv)    fit (do not fit) emission line velocity
   -warc(nwarc)  turns on(off) archive update mode
   -pe(npe)      turns on(off) plotting (on)
   -vx (c e s g)   use corr (comb emiss search guess) redshift
   -pvc (e x)   plot combination (emission correlation) redshift
   -k(nk)        turns on(off) cursor mode
   -nmenu(menu)  turns off(on) menu mode
   -efs     sets sigma limit for lines (    2.00000)
   -efp     sets no. points for line center fit (  2)
   -efc     sets no. points for continuum fit (  3)
   -rmin    sets minimum R for good xcor velocity (    4.00000)
   -nhc(hc)      turns off(on) automatic hard copies
   -hcom  sets hard copy command (imtek)
   -v(nv)        turns on(off) verbose mode

espec's default parameter file looks like this:

Parameter file for ESPEC

DIAG=F           / flag for diagnostic printing
SPECROOT='/'     / root directory of spectrum archive
SPECINST='zmach' / default instrument (zmach or mmt)
SPECRFN=22827    / default reduced file number (ignored)
SPECNUM=1        / Spectrum 1 or 2
SPECORDER=5      / Order of wavelength solution to use (5 or 7)
RFNLIST=''       / Input list of reduced file numbers
RFNFILE=''       / Input file of reduced file numbers
NPFIT=2          / Number of points to fit for line center
NPCONT=3         / Number of points to fit for continuum
VSIG=2.          / Number of sigma above which to accept lines
VELFIT=y         / Compute new emission line velocity (y or n)
LINEFIT=y        / Compute new line profiles (y or n)
WARCH=y          / Write archive records
REPLACE=y        / Replace records in archive
VELTYPE='X'      / Use C>omb, X>cor, E>miss, S>earch, G>uess vel
for initial
CZGUESS=-1.      / Guess for Cz
RMIN=4.          / Limit for R below which CZXC is not used
VELPLOT='C'      / Use C>ombination, X>cor, E>mission velocity when plo
tting
WLSIZE=2         / size of emission line labels
WLFONT='I'       / label font (T>erminal, S>implex, D>uplex, C>omple
x, I>talic)
MODE=Q           / display mode (Q>plot, P>lot, or L>ist)

NPLOTS=1         / number of data sets to plot
PCODE=1          / Current data set to modify

FMIN001=0.       /minimum pixel value for this plot (0=autoscale)
FMAX001=0.       /maximum pixel value for this plot (0=autoscale)

WAVESPACE=T      / If set to Y or T, label plot in wavelength, else in pixels
WSTART=0.        / starting wavelength in Angstroms (0=from first file)
WLENGTH=0.       / width of graph in Angstroms (0=from first file)
WEND=0.          / ending wavelength in Angstroms (0=from first file)
DATAMIN=-10000   / minimum acceptable data value
DATAMAX=50000    / maximum acceptable data value
NDSKIP=0         / number of data points (pixels) to skip before reading
NDMODE=smooth    / sample or smooth or average plot
NDTYPE='nova'    / Nova, hanning, or gaussian smoothing function
ESMOOTH=y        / Smooth spectrum before fitting lines (y or n)
NDREAD=0         / maximum number of pixels to plot
NDSAMP=10        / number by which to increment pixels while reading or to smoot
h
NDHALF=3.        / half-width of Gaussian smoothing
SAMPLE=y         / If set to Y or T, sample rather than average over increment

PLOTFILE='espec.plot'  / plot logging file
CURSOR=T         / if Y, stay in cursor mode after plotting
MENU=F           / if Y, stay in menu mode before and after plotting
KEEPZERO=Y       / If N, ignore pixels with value of 0 when auto-scaling
SPLIT=Y          / separate graphics and text screens

TIMELAB=2        / Size of time label characters; none if <0
XTIME=3300       / x-coordinate of time label
YTIME=3060       / y-coordinate of time label
DISPHRES=4096    / horizontal resolution of display
DISPVRES=3128    / vertical resolution of display
IPL=1            / plot char (1=line,-1=points,2=circles,3=squares,4=triangles
CAPSIZE=0        / size of axis coordinate labels
LABSIZE=1        / size of tick mark labels
LABFONT='S'      / label font (T>erminal, S>implex, D>uplex, C>omple
x, I>talic)
IYBOT=300        / y-coordinate of bottom of first graph
IYTOP=2900       / y-coordinate of top of last graph
IXMIN=300        / x-coordinate of left side of graph
IXMAX=4000       / x-coordinate of right side of graph

NDELREC=2        / Number of archive records to delete
DELREC001='ANALYSIS_SUMMARY'     /* correlation, emission, combined velocities
DELREC002='EMISSION_LINES'       /* emission line parameters
NADDREC=2        / Number of archive records to add
ADDREC001='ansum'        /* correlation, emission, combined velocities
ADDREC002='emlines'      /* emission line parameters

ZOOMED=F         / True if zoomed
WSTARTZ=0.       / Starting wavelength of original graph if zoomed
WLENGTHZ=0.      / width in angstroms of original graph if zoomed
WENDZ=0.         / Final wavelength of original graph if zoomed

LISTFILE='espec.log' / Name of file for logging
LOGFILE='espec.cur' / Cursor log file
HCCOM='imtek'    / Hard copy command

IXCAP=300        / Location of velocities on output page
IYCAP=2950
IXHEAD=300       / Location of file name on output page
IYHEAD=3050
REFIT=F          / Used internally by ESPEC
QUIT=F           / Used internally by ESPEC
FILEOK=T         / Used internally by ESPEC
SPECARCH=T
PCODE001=1
VERBOSE=T
END