task
jsmoker wrote on May 08, 2017
Hi,
I have a problem with scombine in sptables. I want to combine ESO advanced data product spectra. The normal scombine does not work (c.f. http://archive.eso.org/cms/eso-data/help/1dspectra.html#IRAFSPTABLE )
0) I can run scombine as normal within onedspec for example
1) The parameter file seems to be in the correct place within sptable - xonedspec (?)
cd /iraf/iraf
[jsmoker@guest172-232 iraf]$ find . -name 'splot*par' -print
./noao/onedspec/splot.par
./extern/rvsao/Templates/splot.par
./extern/rvsao/Util/splot.par
./extern/guiapps/xapphot/splotpars.par
./extern/sptable/src/xonedspec/splot.par
[jsmoker@guest172-232 iraf]$ find . -name 'scombine*par' -print
./noao/onedspec/scombine/scombine.par
./extern/mscred/src/ccdred/src/combine/scombine.par
./extern/sptable/src/xonedspec/scombine.par
./extern/nfextern/newfirm/scombine.par
./extern/nfextern/src/combine/scombine.par
2) I have done a mkiraf
[jsmoker@guest172-232 CRIRES_OpticalData_NIRDIBS]$ mkiraf
-- creating a new uparm directory
Terminal types: xgterm,xterm,gterm,vt640,vt100,etc.
Enter terminal type: xgterm
A new LOGIN.CL file has been created in the current directory.
You may wish to review and edit this file to change the defaults.
3) I can run splot from sptables ok
NOAO/IRAF PC-IRAF Revision 2.16.1 EXPORT Mon Oct 14 21:40:13 MST 2013
This is the EXPORT version of IRAF V2.16 supporting PC systems.
Welcome to IRAF. To list the available commands, type ? or ??. To get
detailed information about a command, type `help '. To run a
command or load a package, type its name. Type `bye' to exit a
package, or `logout' to get out of the CL. Type `news' to find out
what is new in the version of the system you are using.
Visit http://iraf.net if you have questions or to report problems.
*** Checking update status... Your IRAF system is up to date
*** Initializing SAMP .... No Hub Available
The following commands or packages are currently defined:
adccdrom. dataio. esowfi. guiapps. lists. mtools. obsolete. proto. song. stecf. system. upsqiid. xdimsum.
cfh12k. dbms. euv. images. mem0. nfextern. optic. rvsao. sptable. stella. tables. utilities.
ctio. deitab. fitsutil. language. mscred. noao. plot. softools. sqiid. stsdas. ucsclris. vo.
vocl> sptable
xonedspec. xrv.
sptable> xonedspec
aidpars@ continuum dopcor mkspec refspectra sarith scopy sfit slist standard
autoidentify deredden fitprofs names reidentify sbands sensfunc sflip specplot telluric
bplot dispcor identify ndprep rspectext scombine setairmass sinterp specshift wspectable
calibrate disptrans lcalib odcombine sapertures scoords setjd skytweak splot wspectext
xonedspec> splot ADP.2014-10-01T10:21:08.647.fits - works fine
4) I cannot seem to get scombine to function.
xonedspec> epar scombine
ERROR: task `scombine' has no param file
called as: `eparam (scombine)'
xonedspec> ?
aidpars@ continuum dopcor mkspec refspectra sarith scopy sfit slist standard
autoidentify deredden fitprofs names reidentify sbands sensfunc sflip specplot telluric
bplot dispcor identify ndprep rspectext scombine setairmass sinterp specshift wspectable
calibrate disptrans lcalib odcombine sapertures scoords setjd skytweak splot wspectext
Does anybody have any any ideas?
Thanks!
Jonathan
Mike Fitzpatrick wrote on May 09, 2017
Thanks for the report, the task declaration is incorrect. The solution is to do
cl> sptable # load packages
cl> xonedspec
cl> copy xonedspec$scombine.par xonedspec$scombine/scombine.par
You should then be able to do "lpar scombine" or run the task normally.
jsmoker wrote on May 09, 2017
Hello,
Thanks a lot for having a look at that so quickly. I applied the fix you came up with and:
- I can now see and edit the parameter file OK (thanks!). It is cut and pasted at the end.
- Unfortunately when I try and run scombine within xonedspec I get a:
ERROR: segmentation violation
(I have done mkiraf a few times, logged out and logged into IRAF etc).
- I have tried with two different spectra (making a copy of one of them in order to combine it with itself), both of which will splot using the xonedspec splot. One of the spectra was a non ADP spectra and I could combine using the normal onedspec (which of course does not work on the ADP which is the reason I tried xonedspec in the first place).
- The file is an ESO advanced data product and I put an example here:
wget http://www.sc.eso.org/~jsmoker/ADP.2015-06-09T07:54:50.160.fits.gz
I realise that this may not be easy to track down or spend time on. In any case thanks a lot for your initial suggestion. If in the end this has problems I can always just write out the spectra to ascii and then into fits again and combine them with the normal scombine.
Best regards,
Jonathan
NOAO/IRAF PC-IRAF Revision 2.16.1 EXPORT Mon Oct 14 21:40:13 MST 2013
This is the EXPORT version of IRAF V2.16 supporting PC systems.
PACKAGE = xonedspec
TASK = scombine
input = ADP.2015-06-09T07:54:50.160.fits,ADP.2015-06-09T07:54:50.160_copy.fits List of i
output = nicetest.fits List of output spectra
(noutput= ) List of output number combined spectra
(logfile= STDOUT) Log file
(apertur= ) Apertures to combine
(group = apertures) Grouping option
(combine= average) Type of combine operation
(reject = none) Type of rejection
(first = no) Use first spectrum for dispersion?
(w1 = INDEF) Starting wavelength of output spectra
(w2 = INDEF) Ending wavelength of output spectra
(dw = INDEF) 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
(mode = ql)
ERROR: segmentation violation
Last post on May 09, 2017