2.2.5 February 5, 2002 Ignore pixels with values less than MIN_PIXEL_VALUE at ends of spectra Replace pixels with values less than MIN_PIXEL_VALUE elsewhere in spectra In compbcv.x, allow use of new keywords OBS-LONG, OBS-LAT, and OBS-ELEV 2.2.4 September 19, 2001 Add 10th result reporting mode to EMSAO with line info only. Drop SPECDIR from file name in headings Fix units check in getimage() 2.2.3 July 31, 2001 Change cross-correlation velocity information in header when template changed interactively in EMSAO and XCSAO using the 't' cursor command. Fix pixel cross-correlation so it works. Add option to ignore zeroed pixels in continuum fit and removal Add option in SUMSPEC to remove continuum after rebinning (contout) Add option in SUMSPEC to write each input file to a separate output file (composite file name is "") Separate LINESPEC and SUMSPEC source directories from Makespec into Linespec and Sumspec, with shared tasks into Util 2.2.2 April 3, 2001 Fix EMSAO and plotspec.x in several places to deal with right-to-left spectra. Plot text info after spectrum in EMSAO so it all displays on xgterm. Match arguments in SUMSPEC so renormalization works. 2.2.1 February 9, 2001 Fix bug in EMSAO so emission line structures are always initialized correctly Fix bug in XCSAO to deal with right-to-left dispersion Print aperture separately in XCSAO mode 12 report and EMSAO one-line reports 2.2.0 September 13, 2000 Add options to XCSAO to cross-correlate in wavelength and pixel space by changing the "correlate" parameter into a string. Allow pixel or wavelength limits to be specified for templates in XCSAO so several pieces of a single template spectrum, such as emission lines, may be correlated against an object spectrum. Add two new reporting modes to EMSAO, 8 to print wavelength information about fit emission lines on one line and 9 to print only wavelength shifts on one line. 2.1.27 August 10, 2000 Fix bug in EMSAO which caused problems with sorting during line search 2.1.26 August 2, 2000 Fix echelle mode spectrum number tracking. Print file name and aperture as title for object and template graphs 2.1.25 July 21, 2000 In SUMSPEC, shift to first spectrum's barycentric velocity correction if svel_corr is not none. 2.1.24 July 6, 2000 In LINESPEC and SUMSPEC, specify wavelength dispersion by 3 of start and ending wavelength, wavelength/pixel, and number of pixels Automatically renormalize if maximum spectrum value is less than one. Clean up debugging code for renormalized cross-correlations. 2.1.23 June 8, 2000 Add equivalent width error to new report modes 0 and 7 in EMSAO 2.1.22 April 28, 2000 Change years from 0-120, not just 0-99, to 1900+year. Some instrumentation write 100 instead of 00 for 2000 in the old FITS format. 2.1.21 March 6, 2000 Fix bug in LINESPEC which caused wavelength to be multiplied 2+Z instead of 1+Z if redshift was entered as Z instead of velocity. 2.1.20 January 25, 2000 Read SUBCONT flag in spectrum header to find out whether to subtract continuum. If T, remove continuum; if F, do not. If it is not found, remove continuum. 2.1.19 December 8, 1999 Free spectrum WCS when freeing image stops runaway memory use when looping through large numbers of spectra in XCSAO, EMSAO, or SUMSPEC. Fix bug in BCVCORR which failed to get image observatory name correctly. 2.1.18 September 24, 1999 Fix qplot mode of XCSAO to update best template information, including R when a different best template is selected. 2.1b17 September 15, 1999 Fix tab table output in XCSAO (report_mode= -8, -10, -11, -12, or -14) 2.1b16 September 7, 1999 Write CD as well as CDELT in SUMSPEC. Add gio.h dependencies in mkpkg. 2.1b15 August 24, 1999 Add vel_plot keyword to XCSAO and velplot keyword to QPLOT, so combination velocity-redshifted lines can be displayed by default in dispmode > 0. When continuum is divided out in SUMSPEC, the new continuum is 1. A new option for cont_remove, "divzero", divides by the continuum and subtracts 1 so the spectrum can be cross-correlated. This is what "divide" did in version 2.0 and all previous 2.1 betas, but it is not what you want if you plan to compute equivalent widths. Read templates from header if NTEMP is set as well as if CZPARAMS is set. 2.1b14 July 29, 1999 Add tab table output option to XCSAO. In SUMSPEC, copy entire header if single file in and out, even if there are multiple spectra from that file. In SUMSPEC, normalize by dividing out the mean times an input constant. Add display with continuum subtracted to XCSAO. 2.1b13 July 23, 1999 Add tab table output option to EMSAO. Replace blank fields with underscore. 2.1b12 July 16, 1999 Fix bug so QPLOT can update spectrum headers in both XCSAO and EMSAO. If save_vel=yes and correlate or linefit =no the programs run in QPLOT mode, where nothing will be wriitten unless new fits or correlations are done. 2.1b11 June 30, 1999 Fix bug in XCSAO which caused template heliocentric velocity correction and header velocity to be ignored. In debug mode, note which needed parameters are missing when the heliocentric velocity correction cannot be computed. 2.1b10 June 25, 1999 Add display modes 4 and 5 in EMSAO and XCSAO to display spectrum with continuum subtracted. In EMSAO, add cursor commands to interactively modify some continuum fit parameters: C sets continuum fit parameters, S sets the VELOCITY keyword. 2.1b9 June 15, 1999 Add option in sumspec to zero or fill parts of spectra known to be bad. If they are zeroed, continuum is fit to an interpolated spectrum, and the bad pixels are zeroed after the continuum is removed. Fix fatal bug in linespec 2.1b8 May 21, 1999 Fix minor bug in filllist and change UT from end to start of observation in juldate.x, adding UTEND to UTMID and UTOPEN as possibilities for computing the Julian Date in conjunction with EXPOSURE or EXPTIME if the Julian Date is not already ain the header as GJDN or MJD-OBS. Try DATE keyword if DATE-OBS is not present 2.1b7 May 11, 1999 Fix bugs so that SUMSPEC can stack rebinned spectra 2.1b5 March 30, 1999 Search IRAF database for observatory position using OBSERVAT header keyword if observatory position is not in the spectrum header, for computation of the geocentric part of the barycentric velocity correction. 2.1b4 March 19, 1999 Fix logic in sumspec so multiple spectra are summed correctly. 2.1b3 March 18, 1999 Multiply by 1000 when renormalizing so that continuum is not less than 1 so S/N for line search in EMSAO and continuum removal in XCSAO and SUMSPEC is a usable number. 2.1b2 March 12, 1999 Fix bug so that no error message is returned if an object spectrum file cannot be opened READ_WRITE if it can be opened READ_ONLY and no writing is to be done. 2.1b1 December 30, 1998 Log normalization factor in SUMSPEC for first spectrum. Fix bug in XCSAO to properly deal with template spectrum lists with pathnames. 2.0 July 31, 1998 Add option to toggle header off in dispmode 3 summary display to create figures with a line-labelled spectrum. 2.0b17 June 12, 1998 Limit analysis to pixels within limits specified by spectral WCS in header 2.0b16 May 22, 1998 2.0b15 May 15, 1998 Allow INDEF in minvel and maxvel. Deal better with bad correlations. 2.0b14 April 22, 1998 Fix minor data type clashes. Recode code which reads velocities from headers. This version is probably the final beta before RVSAO 2.0, which is described in a paper which has been accepted for publication in the August 1998 PASP. 2.0b12 April 7, 1998 Fix SUMSPEC bug which crashed the program when automatically finding the wavelength overlap region. 2.0b11 March 4, 1998 Fix BCV calculation when no image is present 2.0b10 February 23, 1998 Fixed handling of velocities in object spectrum files. Fixed handling of velocities in template spectrum files. Fixed handling of barycentric velocity correction from files. Allow up to 16384 points in correlations. Set template velocity to 0.0 km/sec if there is no VELOCITY or APVELn keyword in the header. 2.0b January 9, 1998 Precision problems have been fixed in the barycentric velocity correction. The BCVCORR task has been added to deal with any observatory parameters which might throw off the inline computation of the barycentric velocity correction. ISO date specification (yyyy-mm-dd) is now supported in header date keywords as well as 1900-1999 standard FITS (dd/mm/yy). The new SUMSPEC task creates spectra from multiple spectra taken at different times. The new LINESPEC task creates artificial spectra from lists of emission and/or absorption lines. The unsupported Fourier quotient task FQUOT has been dropped. The option to divide by the continuum in XCSAO instead of subtracting it, which is set by a template flag, DIVCONT, has been fixed. There are several new parameters in EMSAO: MINCONT to determine whether area or equivalent width is computed, LWMIN and LWMAX for the minimum and maximum variation from the mean line width to be accepted, and LSMIN for the minimum ratio of the equivalent width (or area) to its error for a line to be accepted. All intermediate plots, except those produced by IRAF's interactive curve fitting routines during continuum fitting, have an added cursor command, h, to eliminate headings and produce clean graphs for publication. The parameters FIXBAD and BADLINES have been added to EMSAO, SUMSPEC, and XCSAO to remove specified regions of a spectrum before processing. A line list is provided to remove the regions around night sky emission lines. The parameters SPECBAND and TEMPBAND have been added to XCSAO, and SPECBAND and SKYBAND to EMSAO, to make use of multiple bands in multispec data possible. Because of all of these new parameters, users should run RELEARN on any tasks they have used in the past. *************************************************************************** 1.9 December 5, 1995 Bugs fixed, but not a major release. *************************************************************************** 1.8 September 28, 1995 The major change is the ability to use either EMSAO and XCSAO to display previous results saved in the image header. XCSAO has two final displays, the old one with both spectrum and correlation peak and complete per- template correlation data, and the EMSAO display with emission and absorption lines labelled. A quality flag can be set in the header from either final display in either XCSAO or EMSAO. If no correlation is done (correlation=no), only the final velocity information is rewritten by XCSAO. If no new line fit is done (linefit=no), EMSAO writes out only the final velocity information, not overwriting the per-line results. A new task, QPLOT, was added to check redshift results after batch runs. There were also many minor bug fixes. *************************************************************************** 1.7 June 1, 1995 Add continuum removal, explicit emission and rejection limits, and zero-padding as optional template header-parameter-directed operations. The boolean zero_pad parameter was changed to the character zeropad to allow zero-padding to be optionally controlled by the template spectrum header parameter ZEROPAD. czinit, which was a yes or no boolean in previous versions of RVSAO has turned into the character string vel_init in RVSAO 1.7 to allow an initial velocity to be set from the czguess parameter or from a previously-obtained correlation, emission, or combined velocity in the header. Parameters s_abs_reject, s_em_reject, t_abs_reject, and t_em_reject have been moved to XCSAO from CONTPARS as they are used to prepare spectra only for cross-correlation. *************************************************************************** 1.6 February 27, 1995 Optionally display the same summary in XCSAO as EMSAO and allow template- directed high-pass filtering in XCSAO so that narrow line templates, such as those containing emission lines can be used. To enable this option, square-root cosine-bell filtering is now done to the object and template spectrum transforms rather than cosine-bell filtering the correlation. *************************************************************************** 1.5 September 9, 1994 This version writes more complete results from EMSAO and XCSAO to the image header. EMSAO gets larger emission line names and a parameter to specify the dispersion error, as well as a new summary page. XCSAO adds parameters to allow emission line chopping in XCSAO to be controlled by image header parameters, allowing emission line templates to be used, and adds optional zero-padding of transforms before cross-correlation to improve errors. *************************************************************************** 1.4 June 23, 1994 Fixing many small bugs, this version was never released ouside of SAO. *************************************************************************** 1.3 April 12, 1994 Improve handling of missing or unwriteable files, add a cursor option to set peak search width in XCSAO, and improve cursor line selection in EMSAO. Add and improve cursor commands in EMSAO and add spectrum editing in XCSAO. Greatly improve handling of multispec files as both object and template spectra. Add ZVEL CL script to run both XCSAO and EMSAO on a file or list of files. Fix EMSAO line-finding bug. *************************************************************************** 1.2 December 3, 1993 Fix some bugs involving multispec header information and adds the multspec spectrum number to all logs and graphs. *************************************************************************** 1.1 July 1, 1993 Implement the IRAF MWCS spectrum file format, so multispec files can be processed. *************************************************************************** 1.02 February 16, 1993 Clean up the lesser-used XCSAO reporting modes. *************************************************************************** 1.01 February 10, 1993 Add an improved help file for EMSAO, including some hints on using cursor mode and fix some bugs in cursor line identification. The heliocentric velocity correction is now correctly applied. (Barycentric corrections from the image file header value have always worked.) *************************************************************************** 1.00 January 13, 1993 Clean up a few bugs and change the emlines.dat file to include an additional entry for each line containing the half-width in angstroms of the region to included in the fit for that line. Some emission lines are not good Gaussians, so a narrow fit near the peak gets a better center than an attempt to fit the entire line to a Gaussian profile. FQUOT no longer crashes. *************************************************************************** 0.99 August 12, 1992 Clean up multispec handling and fix some bugs in the cursor interaction section of EMSAO, especially with regards to line ID. A variety of bugs have been fixed in the way velocity results are saved and combined velocities computed for archive files and image headers. *************************************************************************** 0.98 Improve multispec capabilities, adding three new parameters to XCSAO: specnum, a list of spectra in a multispec file to be correlated, tempnum, a list of spectra in a multispec file to be used as templates echelle, which is set to "yes" if the spectrum and template are both multispec echelle files--the same list, specnum, is used. The maxshift and minshift parameters are now in km/sec, instead of pixels in Fourier space, which should make appropriate values more obvious. The following parameters have been added to EMSAO: specnum, a list of spectra to process from a multispec file, tempnum, a list of corresponding sky spectra in the multispec file, emlines, a file containing a list of emission lines, emcombine, a file containing a list of combined emission lines, npfit, the number of pixels to fit around each line peak, npcont, the number of pixels to fit around each line continuum, wspan, the wavelength to search around each redshifted line center, ablines, a file containing a list of absorption lines to label, emsearch, a file containing a list of the strongest emission lines. Local modifications of the line lists is now possible without changing the program, and changes to the parameters of the line fits allow use of EMSAO on a wide range of types of spectra. *************************************************************************** 0.97 Add multispec reading capability. *************************************************************************** 0.96 Add correct parameters for FQUOT and fix a bug in pkwidth which would cause XCSAO to bomb if a number of points rather than a peak fraction to fit was specified. *************************************************************************** 0.95 Properly remove emission lines if IRAF continuum fitting is chosen and correct several argument mismatches in subroutine calls. EMSAO no longer gives emission line centers that are one pixel off. ***************************************************************************