iraf.net

We are trying to set up IRAF and use ircombbin, but we are getting the following error:Command line:bastet [33] [aaron] test: ls
binir040706.0225.fits images inj ir040706.flatk.fits login.cl~
coordsK inh ink login.cl tempbastet [34] [aaron] test: cl
setting terminal type to xterm... NOAO PC-IRAF Revision 2.12.2-EXPORT Sun Jan 25 16:09:03 MST 2004
This is the EXPORT version of PC-IRAF V2.12 supporting most PC systems. Welcome to IRAF. To list the available commands, type ? or ??. To get
detailed information about a command, type `help command'. 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. The following
commands or packages are currently defined: Kpsfmag ircombbin proto. tables.
cleanup ircopy rvsao. tirsky
color. language. softools. tirskybin
dataio. lists. sortfilter tyaloirproc
dbms. mvfiles sortfilter2 tyaloirprocbadq
ftools. noao. stecf. utilities.
images. obsolete. stsdas. xray.
ircombbadq plot. system.
cl> !ds9&
[1] 20616
cl> ircombbin
ERROR: task `ircombbin' has no param file
cl> This is our login.cl file:bastet [36] [aaron] test: cat login.cl
# LOGIN.CL -- User login file for the IRAF command language.# Identify login.cl version (checked in images.cl).
if (defpar ("logver"))
logver = "IRAF V2.12.2 January 2004"set home = "/local/data/bastet4/aaron/Results/Analysis/Scripts/"
set imdir = "/d0/imdirs/aaron/"
set uparm = "home$uparm/"
set userid = "aaron"#-------------------------------------------------------------------------------------------
task ircopy = /local/data/bastet4/aaron/Results/Analysis/Scripts/ircopy.cl
task ircombbin = /local/data/bastet4/aaron/Results/Analysis/Scripts/ircombbin.v2.cl
task ircombbadq = /local/data/bastet4/aaron/Results/Analysis/Scripts/ircombbadq.cl
task tirsky = /local/data/bastet4/aaron/Results/Analysis/Scripts/tirsky.cl
task tirskybin = /local/data/bastet4/aaron/Results/Analysis/Scripts/tirskybin.cl
task tyaloirproc = /local/data/bastet4/aaron/Results/Analysis/Scripts/tyaloirproc.cl
task tyaloirprocbadq = /local/data/bastet4/aaron/Results/Analysis/Scripts/tyaloirprocbadq.cl
task sortfilter = /local/data/bastet4/aaron/Results/Analysis/Scripts/sortfilter.cl
task sortfilter2 = /local/data/bastet4/aaron/Results/Analysis/Scripts/sortfilter2.cl
task cleanup = /local/data/bastet4/aaron/Results/Analysis/Scripts/cleanup.cl
task Kpsfmag = /local/data/bastet4/aaron/Results/Analysis/Scripts/Kpsfmag.cl
task mvfiles = /local/data/bastet4/aaron/Results/Analysis/Scripts/mvfiles.cl
keep
#-------------------------------------------------------------------------------------------# Set the terminal type.
if (envget("TERM") == "xterm") {
if (!access (".hushiraf"))
print "setting terminal type to xgterm..."
stty xgterm
} else {
if (!access (".hushiraf"))
print "setting terminal type to xterm..."
stty xterm
}# Uncomment and edit to change the defaults.
#set editor = vi
#set printer = lp
#set pspage = "letter"
#set stdimage = imt800
#set stdimcur = stdimage
#set stdplot = lw
#set clobber = no
#set filewait = yes
#set cmbuflen = 512000
#set min_lenuserarea = 64000
#set imtype = "imh"
#set imextn = "oif:imh fxf:fits,fit plf:pl qpf:qp stf:hhh,??h"
# XIMTOOL/DISPLAY stuff. Set node to the name of your workstation to
# enable remote image display. The trailing "!" is required.
#set node = "my_workstation!"# CL parameters you mighth want to change.
#ehinit = "nostandout eol noverify"
#epinit = "standout showall"
showtype = yes# Load the default CL package. Doing so here allows us to override package
# paths and load personalized packages from our loginuser.cl.
clpackage
# Default USER package; extend or modify as you wish. Note that this can
# be used to call FORTRAN programs from IRAF.package usertask $adb $bc $cal $cat $comm $cp $csh $date $dbx $df $diff = "$foreign"
task $du $find $finger $ftp $grep $lpq $lprm $ls $mail $make = "$foreign"
task $man $mon $mv $nm $od $ps $rcp $rlogin $rsh $ruptime = "$foreign"
task $rwho $sh $spell $sps $strings $su $telnet $tip $top = "$foreign"
task $vi $emacs $w $wc $less $rusers $sync $pwd $gdb = "$foreign"task $xc $mkpkg $generic $rtar $wtar $buglog = "$foreign"
#task $fc = "$xc -h $* -limfort -lsys -lvops -los"
task $fc = ("$" // envget("iraf") // "unix/hlib/fc.csh" //
" -h $* -limfort -lsys -lvops -los")
task $nbugs = ("$(setenv EDITOR 'buglog -e';" //
"less -Cqm +G " // envget ("iraf") // "local/bugs.*)")
task $cls = "$clear;ls"
task $clw = "$clear;w"
task $pg = ("$(less -Cqm $*)")if (access ("home$loginuser.cl"))
cl < "home$loginuser.cl"
;keep
prcache directory
cache directory page type help# Print the message of the day.
if (access (".hushiraf"))
menus = no
else {
clear; type hlib$motd
}# Delete any old MTIO lock (magtape position) files.
if (deftask ("mtclean"))
mtclean
else
delete uparm$mt?.lok,uparm$*.wcs verify-# List any packages you want loaded at login time, ONE PER LINE.
images # general image operators
plot # graphics tasks
dataio # data conversions, import export
lists # list processing# The if(deftask...) is needed for V2.9 compatibility.
if (deftask ("proto"))
proto # prototype or ad hoc taskstv # image display
utilities # miscellaneous utilities
noao # optical astronomy packageskeepThis is our uparm directory:bastet [49] [aaron] Scripts: cd uparm
/local/data/bastet4/aaron/Results/Analysis/Scripts/uparm
bastet [50] [aaron] uparm: ls
clemvfils.par tvdisply.par
bastet [51] [aaron] uparm: Do you have any ideas why we would be getting the error message:ERROR: task `ircombbin' has no param filewhen we call the ircombbin routine?Any help is GREATLY appreciated.Thanks,
Aaron

The 'task has no param file' message usually means that either the path to the task/script is incorrect, or that the task really doesn't have any parameters but was declared as if it did. For scripts, parameters are anything declared before the 'begin' statement in a procedure script (command scripts don't have parameters by definition). If your variables are all declared after the 'begin' then the script has no parameters and should be declared with a '$' before the task name, i.e.[code:1:b150d3b125]task $ircombbin = /path/ircombbin.cl[/code:1:b150d3b125]-Mike

This is my ircombbin.v2.cl file:bastet [47] [aaron] Scripts: cat ircombbin.v2.cl
procedure ircombbin()string refimg, shiftfile, coordfile, outfile,continue
string ink, inj, inh, refk, refj, refh, date1, one
struct *imglistbeginstring shiftname, imgname, date, cont, which
struct infile
int i#------------------SORTING IMAGES INTO RIGHT FILTERS-------------
if(!access("ink")){
if(!access("inj")){
if(!access("inh")){sortfilter2}}}#-----------------GIVING THE OUTPUT FILES DISTINCT NAMES-------
print("Please input a number (e.g. 1) for the output file:")
which = one
#-----------------------------------------------------------------#-----------------------------------------------------------------
#--------------REDUCE AND COMBINE J FILTER----------------------------------if(!access("inj2")){
if(access("inj")){
if(!access("coordsJ")){
inj="inj"
infile = inj
imglist = infile
} #end of !access coordsJ#-------------------------------TIRSKYBIN.CL----------------------!/Users/buxton/iraf/scripts/makefile.out #make files in prep for tirsky.cl
tirskybin(irflt="J")
tyaloirproc(irflt = "J") #normalizes flat and prepares sky image,
#does sky subtraction,
#flat fielding, bad pixel correction
#---------SHIFT IMAGES--------------------- inj="ssjimages"
infile = inj
imglist = infile
while(fscan(imglist,imgname) != EOF){
imexamine(imgname, logfile="coordsJ",keeplog+,use_dis+)
}
!/Users/buxton/iraf/scripts/getpos2.out #find x and y shifts for images
inj = "inj"
infile = inj
imglist = infile
i = 1
while(fscan(imglist, imgname) !=EOF) {
if(i == 1){
refj = imgname
imgets(imgname,"FILENAME")
date1=substr(imgets.value,3,9)
} #end of if statement
i=i+1
print("r"//imgname, >> "inj2")
print("shiftr"//imgname, >> "outshiftj")
} #end of while loopimalign(input="@inj2",referenc=refj,coords="coordsJ", output="@outshiftj",shifts="shiftsJ",boxsize=7,bigbox=11,interp="linear",verbose-)date=date1#---------COMBINE IMAGES---------------------#If the number of images = 2, then you must use no rejection when combining images, else use minmax rejectionif((i-1) == 2){
imcomb(input="@outshiftj",output="shiftrir"//date//"ADDEDJ."//which,combine="average",reject="none",mclip+,lsigma=3., hsigma=3.,nlow=1,nhigh=1,nkeep=1)}
else{imcomb(input="@outshiftj",output="shiftrir"//date//"ADDEDJ."//which,combine="average",reject="minmax",mclip+,lsigma=3., hsigma=3.,nlow=1,nhigh=1,nkeep=1)}copy(input="inj",output="inj."//which)
copy(input="coordsJ",output="coordsJ."//which)
copy(input="shiftsJ",output="shiftsJ."//which)delete("inj")
delete("shiftsJ")
delete("coordsJ")
delete("ssjimages")
cleanup} #end access inj
} #end !access inj2#-----------------------------------------------------------------------#-----------------------------------------------------------------
#--------------REDUCE AND COMBINE H FILTER----------------------------------if(!access("inh2")){
if(access("inh")){
if(!access("coordsH")){
inh="inh"
infile = inh
imglist = infile
} #end of !access coordsH#-------------------------------TIRSKYBIN.CL----------------------!/Users/buxton/iraf/scripts/makefile.out #make files in prep for tirsky.cl
tirskybin(irflt="H")
tyaloirproc(irflt = "H") #prepares flat and sky image, does sky subtraction,
#flat fielding, bad pixel correction
#---------SHIFT IMAGES--------------------- inh="sshimages"
infile = inh
imglist = infile
while(fscan(imglist,imgname) != EOF){
imexamine(imgname, logfile="coordsH",keeplog+,use_dis+)
}
!/Users/buxton/iraf/scripts/getpos2.out #find x and y shifts for images
inh = "inh"
infile = inh
imglist = infile
i = 1
while(fscan(imglist, imgname) !=EOF) {
if(i == 1){
refh = imgname
imgets(imgname,"FILENAME")
date1=substr(imgets.value,3,9)
} #end of if statement
i=i+1
print("r"//imgname, >> "inh2")
print("shiftr"//imgname, >> "outshifth")
} #end of while loopimalign(input="@inh2",referenc=refh,coords="coordsH", output="@outshifth",shifts="shiftsH",boxsize=7,bigbox=11,interp="linear",verbose-)date=date1#---------COMBINE IMAGES---------------------#If the number of images = 2, then you must use no rejection when combining images, else use minmax rejectionif((i-1) == 2){
imcomb(input="@outshifth",output="shiftrir"//date//"ADDEDH."//which,combine="average",reject="none",mclip+,lsigma=3., hsigma=3.,nlow=1,nhigh=1,nkeep=1)}
else{imcomb(input="@outshifth",output="shiftrir"//date//"ADDEDH."//which,combine="average",reject="minmax",mclip+,lsigma=3., hsigma=3.,nlow=1,nhigh=1,nkeep=1)}copy(input="inh",output="inh."//which)
copy(input="coordsH",output="coordsH."//which)
copy(input="shiftsH",output="shiftsH."//which)delete("inh")
delete("shiftsH")
delete("coordsH")
delete("sshimages")
cleanup} #end access inh
} #end !access inh2#-----------------------------------------------------------------------#-----------------------------------------------------------------
#--------------REDUCE AND COMBINE K FILTER----------------------------------if(!access("ink2")){
if(access("ink")){
if(!access("coordsK")){
ink="ink"
infile = ink
imglist = infile
} #end of !access coordsK#-------------------------------TIRSKYBIN.CL----------------------!/Users/buxton/iraf/scripts/makefile.out #make files in prep for tirsky.cl
tirskybin(irflt="K")
tyaloirproc(irflt = "K") #prepares flat and sky image, does sky subtraction,
#flat fielding, bad pixel correction
#---------SHIFT IMAGES--------------------- ink="sskimages"
infile = ink
imglist = infile
while(fscan(imglist,imgname) != EOF){
imexamine(imgname, logfile="coordsK",keeplog+,use_dis+)
}
!/Users/buxton/iraf/scripts/getpos2.out #find x and y shifts for images
ink = "ink"
infile = ink
imglist = infile
i = 1
while(fscan(imglist, imgname) !=EOF) {
if(i == 1){
refk = imgname
imgets(imgname,"FILENAME")
date1=substr(imgets.value,3,9)
} #end of if statement
i=i+1
print("r"//imgname, >> "ink2")
print("shiftr"//imgname, >> "outshiftk")
} #end of while loopimalign(input="@ink2",referenc=refk,coords="coordsK", output="@outshiftk",shifts="shiftsK",boxsize=7,bigbox=11,interp="linear",verbose-)date=date1#---------COMBINE IMAGES---------------------#If the number of images = 2, then you must use no rejection when combining images, else use minmax rejectionif((i-1) == 2){
imcomb(input="@outshiftk",output="shiftrir"//date//"ADDEDK."//which,combine="average",reject="none",mclip+,lsigma=3., hsigma=3.,nlow=1,nhigh=1,nkeep=1)}
else{imcomb(input="@outshiftk",output="shiftrir"//date//"ADDEDK."//which,combine="average",reject="minmax",mclip+,lsigma=3., hsigma=3.,nlow=1,nhigh=1,nkeep=1)}copy(input="ink",output="ink."//which)
copy(input="coordsK",output="coordsK."//which)
copy(input="shiftsK",output="shiftsK."//which)delete("ink")
delete("shiftsK")
delete("coordsK")
delete("sskimages")
cleanup} #end access ink
} #end !access ink2#-----------------------------------------------------------------------endI don't think it should have the $ in the login.cl file. I checked the uparm path in the login.cl file and nothing seems to be wrong. Is there any other reason why we would be getting this error?Thanks again,
Aaron

You'll need to rename the script file to 'ircombbin.cl'. When looking for the parameter file a filename is being constructed internally based on the task name (with a '.par' or a '.cl' extension). An alternative is to create an 'ircombbin.par' file to explicitly hold the parameters, if you did this it would look like:[code:1:44b313c87d]
refimg,s,h,,,,
shiftfile,s,h,,,,
coordfile,s,h,,,,
outfile,s,h,,,,
continue,s,h,,,,
ink,s,h,,,,
inj,s,h,,,,
inh,s,h,,,,
refk,s,h,,,,
refj,s,h,,,,
refh,s,h,,,,
date1,s,h,,,,
one,s,h,,,,
imglist,*struct,h,"",,,
mode,s,h,"ql",,,
[/code:1:44b313c87d]-Mike

Thanks so much =], looks like it solved the problem.