iraf.net

One of my users is having trouble when using splot on some WIYN Hydra
spectra. Simply writing the spectrum out with "i" in splot is changing
the CRVAL1 and CDELT1 values in the header, as well as the EPOCH, HJD,
and VLSR keywords(!) and adding an WAT1_002 keyword.I also notice that the DC-FLAG keyword has a slightly odd appearance
(the multiple "/" characters).We are running V2.12.2-a EXPORT Wed Jul 14 15:14:35 MST 2004
under Solaris 2.8I've included the pre and post splot headers below, and can upload the
image(s) somewhere if that would help.Thanks,
RickyHere is the "diff" between the pre and post-splot headers (followed by
the headers in their entirety) :1c1
<
/home/harvest/jc4qn/monoceros/decfeb_final_spectra/HD107328a.0085.fits[1884][real]:
HD107328
---
> test2.fits[1884][real]: HD107328
3,5c3,5
< Line storage mode, physdim [1884], length of user area 8546 s.u.
< Created Wed 20:10:07 26-Jul-2006, Last modified Fri 13:24:36 24-Feb-2006
< Pixel file
"/home/harvest/jc4qn/monoceros/decfeb_final_spectra/HD107328a.0085.fits"
[ok]
---
> Line storage mode, physdim [1884], length of user area 4172 s.u.
> Created Mon 08:31:57 27-Nov-2006, Last modified Mon 08:31:57 27-Nov-2006
> Pixel file "test2.fits" [ok]
8,9c8,9
< DATE = '2006-02-22T20:29:47' / Date FITS file was generated
< IRAF-TLM= '13:24:36 (24/02/2006)' / Time of last modification
---
> DATE = '2006-11-27T13:31:57' / Date FITS file was generated
> IRAF-TLM= '08:31:57 (27/11/2006)' / Time of last modification
83c83
< WAT1_001= 'wtype=linear label=Wavelength units=angstroms'
---
> WAT1_001= 'wtype=linear label=Wavelength units=angstroms
units_display=angstrom'
91c91
< CRVAL1 = 3.64443859422733
---
> CRVAL1 = 3.64044143551894
93,94c93,94
< CDELT1 = 1.05884953326994E-4
< CD1_1 = 1.05884953326994E-4
---
> CDELT1 = 1.07872600621072E-4
> CD1_1 = 1.07872600621072E-4
104a105
> WAT1_002= 's '
Pre-splot spectral header:/home/harvest/jc4qn/monoceros/decfeb_final_spectra/HD107328a.0085.fits[1884][real]:
HD107328
No bad pixels, min=0., max=0. (old)
Line storage mode, physdim [1884], length of user area 8546 s.u.
Created Wed 20:10:07 26-Jul-2006, Last modified Fri 13:24:36 24-Feb-2006
Pixel file
"/home/harvest/jc4qn/monoceros/decfeb_final_spectra/HD107328a.0085.fits"
[ok]
EXTEND = F / File may contain extensions
ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator
DATE = '2006-02-22T20:29:47' / Date FITS file was generated
IRAF-TLM= '13:24:36 (24/02/2006)' / Time of last modification
OBJECT = 'HD107328' / Name of the
OBSERVAT= 'KPNO ' / observatory
OBSERVER= 'R.Munoz ' / observers
PROPID = '06A-0347 ' / observing proposal ID
EXPTIME = 56. / actual integration time, seconds
DARKTIME= 740. / total elapsed time, seconds
IMAGETYP= 'object ' / object, flat, bias, etc.
TELESCOP= 'kpno ' / telescope name
COMMENT
DATE-OBS= '2006-02-08 ' / Y2K compliant (yyyy-mm-ddThh:mm:ss)
UT = '10:30:47.802 ' / universal time (start of exposure)
JD = '2453774.93837729' / julian date of observation
INITTIME= 0.182 / software/hardware initialization (sec)
PREPTIME= 8.144 / CCD preparation time (sec)
OBSTIME = 739.384 / active observing time (sec)
READTIME= 121.392 / CCD readout time (sec)
TOTTIME = 869.106 / total time spent in task (sec)
COMMENT
DETECTOR= 't2ka ' / detector name
PREFLASH= 0 / preflash time, seconds
GAIN = 2.1 / gain, electrons per adu
DWELL = 5 / sample integration time, microseconds
RDNOISE = 4. / read noise, electrons
CAMTEMP = -87 / camera temperature, Celsius
DEWTEMP = -155 / dewar temperature, Celsius
CCDSEC = '[600:1528,1:2048]' / orientation to full frame
ORIGSEC = '[1:2048,1:2048] ' / original size full frame
CCDSUM = '1 1 ' / on chip summation
INSTRUME= 'test ' / instrument name
MJDHDR = 2.4537749481550E+06 / MJD
LSTHDR = '12:32:33.719 ' / local sidereal time
ELMAP = '61:10:58.848 ' / mapper elevation
AZMAP = '186:07:36.399 ' / mapper azimuth
TRACK = 'Tracking ' / tracking
COORDS = ' Helio. mean FK5 ' / coordinate system
EPOCH = 2006.100 / same as EQUINOX (for back compat.)
TARGRA = '12:20:40.798 ' / Telescope target RA
TARGDEC = '03:16:35.047 ' / Telescope target DEC
EQUINOX = 2006.100 / equinox of position
RA = '12:20:40.8' / Telescope RA
DEC = '+03:16:35' / Telescope DEC
RAOFFST = '-23:59:59.939 ' / Telescope RA offset
DECOFFST= '-00:00:00.23 ' / Telescope DEC offset
ZD = '28:48:27.472 ' / zenith distance
AIRMASS = 1.141 / airmass at end of exposure
ROTANGLE= 0.000 / MNIR rotator angle, degrees CW from N
ROTOFF = 0.000 / MNIR rotator angle offset
FOCUS = 8143.062 / Telescope focus (microns)
ROTPORT = 'mnir ' / Tertiary fold position
FOLDPOS = 'On ' / Tertiary port selection
GRATNAME= '600@10.1 ' / grating name
GRATANGL= 25.288 / grating angle (degrees)
DEWARANG= 0.200 / dewar azimuth (degrees)
CAMANGLE= 30.000 / camera-collimator angle (degrees)
FIBFOCUS= -0.212 / fiber focus (inches)
FIBRNAME= 'Red ' / fiber name
CAMFOCUS= -0.010 / camera focus (inches)
CAMNAME = 'Bench Spect (Red) Camera' / camera name
FILTNAMA= 'none ' / bench filter a (in beam)
FILTNAMB= 'none ' / bench filter b (in beam)
FILTNAMC= 'G4_GG-420 ' / bench filter c (in beam)
FILTNAMD= 'none ' / bench filter d (in beam)
FILTPOSA= 0.000 / bench filter number a
FILTPOSB= 0.000 / bench filter number b
FILTPOSC= 1.000 / bench filter number c
FILTPOSD= 0.000 / bench filter number d
FIBCABLE= 'RED ' / Hydra Fibre Cable Name
ICSSEQ = 45 / ICSINFO-sequence
ICSDATE = 'Feb 8 03:45:22 ' / ICSINFO-timestamp
=
WCSDIM = 1 /2
LTM1_1 = 1.
WAT0_001= 'system=equispec'
WAT1_001= 'wtype=linear label=Wavelength units=angstroms'
TRIM = 'Feb 14 14:14 Trim data section is [1:929,1:2048]'
OVERSCAN= 'Feb 14 14:14 Overscan section is [934:963,1:2048] with
mean=953.3094'
ZEROCOR = 'Feb 14 14:14 Zero level correction image is Zero_n1.fits'
CCDPROC = 'Feb 14 14:14 CCD processing done'
LJD = 2453774.
APNUM1 = '85 1 165.09 170.60'
CTYPE1 = 'LINEAR '
CRVAL1 = 3.64443859422733
CRPIX1 = -127.
CDELT1 = 1.05884953326994E-4
CD1_1 = 1.05884953326994E-4
DC-FLAG = 1 /1 /1
DCLOG1 = 'REFSPEC1 = a1030a1032.ms'
SKYSUB = 'yes '
MSTITLE = 'HD 107328'
BANDID1 = 'spectrum - background none, weights none, clean no'
LTV1 = -128.
HJD = 2453774.94213742
VHELIO = +35.7
VLSR = 22.8858023106138
VSUN = ' 20. 18. 30. 1900.'
POST Splot spectral header:test2.fits[1884][real]: HD107328
No bad pixels, min=0., max=0. (old)
Line storage mode, physdim [1884], length of user area 4172 s.u.
Created Mon 08:31:57 27-Nov-2006, Last modified Mon 08:31:57 27-Nov-2006
Pixel file "test2.fits" [ok]
EXTEND = F / File may contain extensions
ORIGIN = 'NOAO-IRAF FITS Image Kernel July 2003' / FITS file originator
DATE = '2006-11-27T13:31:57' / Date FITS file was generated
IRAF-TLM= '08:31:57 (27/11/2006)' / Time of last modification
OBJECT = 'HD107328' / Name of the
OBSERVAT= 'KPNO ' / observatory
OBSERVER= 'R.Munoz ' / observers
PROPID = '06A-0347 ' / observing proposal ID
EXPTIME = 56. / actual integration time, seconds
DARKTIME= 740. / total elapsed time, seconds
IMAGETYP= 'object ' / object, flat, bias, etc.
TELESCOP= 'kpno ' / telescope name
COMMENT
DATE-OBS= '2006-02-08 ' / Y2K compliant (yyyy-mm-ddThh:mm:ss)
UT = '10:30:47.802 ' / universal time (start of exposure)
JD = '2453774.93837729' / julian date of observation
INITTIME= 0.182 / software/hardware initialization (sec)
PREPTIME= 8.144 / CCD preparation time (sec)
OBSTIME = 739.384 / active observing time (sec)
READTIME= 121.392 / CCD readout time (sec)
TOTTIME = 869.106 / total time spent in task (sec)
COMMENT
DETECTOR= 't2ka ' / detector name
PREFLASH= 0 / preflash time, seconds
GAIN = 2.1 / gain, electrons per adu
DWELL = 5 / sample integration time, microseconds
RDNOISE = 4. / read noise, electrons
CAMTEMP = -87 / camera temperature, Celsius
DEWTEMP = -155 / dewar temperature, Celsius
CCDSEC = '[600:1528,1:2048]' / orientation to full frame
ORIGSEC = '[1:2048,1:2048] ' / original size full frame
CCDSUM = '1 1 ' / on chip summation
INSTRUME= 'test ' / instrument name
MJDHDR = 2.4537749481550E+06 / MJD
LSTHDR = '12:32:33.719 ' / local sidereal time
ELMAP = '61:10:58.848 ' / mapper elevation
AZMAP = '186:07:36.399 ' / mapper azimuth
TRACK = 'Tracking ' / tracking
COORDS = ' Helio. mean FK5 ' / coordinate system
EPOCH = 2006.100 / same as EQUINOX (for back compat.)
TARGRA = '12:20:40.798 ' / Telescope target RA
TARGDEC = '03:16:35.047 ' / Telescope target DEC
EQUINOX = 2006.100 / equinox of position
RA = '12:20:40.8' / Telescope RA
DEC = '+03:16:35' / Telescope DEC
RAOFFST = '-23:59:59.939 ' / Telescope RA offset
DECOFFST= '-00:00:00.23 ' / Telescope DEC offset
ZD = '28:48:27.472 ' / zenith distance
AIRMASS = 1.141 / airmass at end of exposure
ROTANGLE= 0.000 / MNIR rotator angle, degrees CW from N
ROTOFF = 0.000 / MNIR rotator angle offset
FOCUS = 8143.062 / Telescope focus (microns)
ROTPORT = 'mnir ' / Tertiary fold position
FOLDPOS = 'On ' / Tertiary port selection
GRATNAME= '600@10.1 ' / grating name
GRATANGL= 25.288 / grating angle (degrees)
DEWARANG= 0.200 / dewar azimuth (degrees)
CAMANGLE= 30.000 / camera-collimator angle (degrees)
FIBFOCUS= -0.212 / fiber focus (inches)
FIBRNAME= 'Red ' / fiber name
CAMFOCUS= -0.010 / camera focus (inches)
CAMNAME = 'Bench Spect (Red) Camera' / camera name
FILTNAMA= 'none ' / bench filter a (in beam)
FILTNAMB= 'none ' / bench filter b (in beam)
FILTNAMC= 'G4_GG-420 ' / bench filter c (in beam)
FILTNAMD= 'none ' / bench filter d (in beam)
FILTPOSA= 0.000 / bench filter number a
FILTPOSB= 0.000 / bench filter number b
FILTPOSC= 1.000 / bench filter number c
FILTPOSD= 0.000 / bench filter number d
FIBCABLE= 'RED ' / Hydra Fibre Cable Name
ICSSEQ = 45 / ICSINFO-sequence
ICSDATE = 'Feb 8 03:45:22 ' / ICSINFO-timestamp
=
WCSDIM = 1 /2
LTM1_1 = 1.
WAT0_001= 'system=equispec'
WAT1_001= 'wtype=linear label=Wavelength units=angstroms
units_display=angstrom'
TRIM = 'Feb 14 14:14 Trim data section is [1:929,1:2048]'
OVERSCAN= 'Feb 14 14:14 Overscan section is [934:963,1:2048] with
mean=953.3094'
ZEROCOR = 'Feb 14 14:14 Zero level correction image is Zero_n1.fits'
CCDPROC = 'Feb 14 14:14 CCD processing done'
LJD = 2453774.
APNUM1 = '85 1 165.09 170.60'
CTYPE1 = 'LINEAR '
CRVAL1 = 3.64044143551894
CRPIX1 = -127.
CDELT1 = 1.07872600621072E-4
CD1_1 = 1.07872600621072E-4
DC-FLAG = 1 /1 /1
DCLOG1 = 'REFSPEC1 = a1030a1032.ms'
SKYSUB = 'yes '
MSTITLE = 'HD 107328'
BANDID1 = 'spectrum - background none, weights none, clean no'
LTV1 = -128.
HJD = 2453774.94213742
VHELIO = +35.7
VLSR = 22.8858023106138
VSUN = ' 20. 18. 30. 1900.'
WAT1_002= 's '

This is a bug when the spectrum is log sampled and an image section has been extracted. The bug log is below. I don't see an easy workaround. The user should find an alternative means for modifying the spectrum than using SPLOT.Frank Valdes
[quote:683874f6b5]
NUMBER: 557
MODULE: splot
SYSTEM: -V2.13beta
DATE: Mon Dec 4 15:29:44 MST 2006
FROM: valdesBUG: When the dispersion function for a spectrum is log sampled,
e.g. DC-FLAG=1, and there is also a logical/physical pixel offset,
i.e. LTV1 is not zero, then writing out a new spectrum with the
'i' key the results in a wrong WCS. There is no workaround.STATUS: Fixed for a future release.
[/quote:683874f6b5]