Welcome to iraf.net Friday, March 29 2024 @ 05:23 AM GMT


 Forum Index > Help Desk > General IRAF New Topic Post Reply
  task
   
jsmoker
 05/08/2017 06:57PM (Read 1910 times)  
+----
Newbie

Status: offline


Registered: 09/15/2008
Posts: 8
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 https://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\$this->_split2($m[0]) sptable
xonedspec. xrv.

sptable\$this->_split2($m[0]) 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\$this->_split2($m[0]) splot ADP.2014-10-01T10:21:08.647.fits - works fine

4) I cannot seem to get scombine to function.

xonedspec\$this->_split2($m[0]) epar scombine
ERROR: task `scombine' has no param file
called as: `eparam (scombine)'
xonedspec\$this->_split2($m[0]) ?
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



 
Profile Email
 Quote
fitz
 05/09/2017 04:32PM  
AAAAA
Admin

Status: offline


Registered: 09/30/2005
Posts: 4040
Thanks for the report, the task declaration is incorrect. The solution is to do

cl\$this->_split2($m[0]) sptable # load packages
cl\$this->_split2($m[0]) xonedspec
cl\$this->_split2($m[0]) copy xonedspec$scombine.par xonedspec$scombine/scombine.par

You should then be able to do "lpar scombine" or run the task normally.

 
Profile Email
 Quote
jsmoker
 05/09/2017 06:24PM  
+----
Newbie

Status: offline


Registered: 09/15/2008
Posts: 8
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


 
Profile Email
 Quote
   
Content generated in: 0.09 seconds
New Topic Post Reply

Normal Topic Normal Topic
Sticky Topic Sticky Topic
Locked Topic Locked Topic
New Post New Post
Sticky Topic W/ New Post Sticky Topic W/ New Post
Locked Topic W/ New Post Locked Topic W/ New Post
View Anonymous Posts 
Anonymous users can post 
Filtered HTML Allowed 
Censored Content 
dog allergies remedies cialis 20 mg chilblain remedies


Privacy Policy
Terms of Use

User Functions

Login