procedure xcomb  (inmage)
 
#
# Procedure to combine ABBA spectra from OSIRIS x-disp spectra
# by Richard Elston 9/3/93
#

string inmage  {prompt="input raw image prefix"}
string offset  {"none", prompt="offset files", mode="h"}

begin

	string inimage, rfsky,ps
	inimage=inmage

#
# put mask frames into the headers 
#

	hedit(inimage//"*.ksub.imh", "BPM", "k.mask.pl", add+, ver-)
	hedit(inimage//"*.hsub.imh", "BPM", "h.mask.pl", add+, ver-)
	hedit(inimage//"*.jsub.imh", "BPM", "j.mask.pl", add+, ver-)
	hedit(inimage//"*.isub.imh", "BPM", "i.mask.pl", add+, ver-)
#
# imcomb the images with out any shifts
#


	imdel(inimage//".kavg.imh", ver-)
	
imcomb (inimage//"*.ksub.imh",inimage//".kavg.imh",
plfile="", sigma="", logfile="STDOUT", combine="average",
reject="pclip", project=no, outtype="real", offsets=offset,
masktype="badvalue", maskvalue=0., blank=0., scale="none", zero="none",
weight="none", statsec="", expname="", lthreshold=INDEF, hthreshold=INDEF,
nlow=1, nhigh=1, nkeep=1, mclip=yes, lsigma=3., hsigma=3., rdnoise="0.",
gain="1.", snoise="0.", sigscale=0.1, pclip=-0.5, grow=0)

	imdel(inimage//".havg.imh", ver-)
	
imcomb (inimage//"*.hsub.imh",inimage//".havg.imh",
plfile="", sigma="", logfile="STDOUT", combine="average",
reject="pclip", project=no, outtype="real", offsets=offset,
masktype="badvalue", maskvalue=0., blank=0., scale="none", zero="none",
weight="none", statsec="", expname="", lthreshold=INDEF, hthreshold=INDEF,
nlow=1, nhigh=1, nkeep=1, mclip=yes, lsigma=3., hsigma=3., rdnoise="0.",
gain="1.", snoise="0.", sigscale=0.1, pclip=-0.5, grow=0)

	imdel(inimage//".javg.imh", ver-)
	
imcomb (inimage//"*.jsub.imh",inimage//".javg.imh",
plfile="", sigma="", logfile="STDOUT", combine="average",
reject="pclip", project=no, outtype="real", offsets=offset,
masktype="badvalue", maskvalue=0., blank=0., scale="none", zero="none",
weight="none", statsec="", expname="", lthreshold=INDEF, hthreshold=INDEF,
nlow=1, nhigh=1, nkeep=1, mclip=yes, lsigma=3., hsigma=3., rdnoise="0.",
gain="1.", snoise="0.", sigscale=0.1, pclip=-0.5, grow=0)

	imdel(inimage//".iavg.imh", ver-)
	
imcomb (inimage//"*.isub.imh",inimage//".iavg.imh",
plfile="", sigma="", logfile="STDOUT", combine="average",
reject="pclip", project=no, outtype="real", offsets=offset,
masktype="badvalue", maskvalue=0., blank=0., scale="none", zero="none",
weight="none", statsec="", expname="", lthreshold=INDEF, hthreshold=INDEF,
nlow=1, nhigh=1, nkeep=1, mclip=yes, lsigma=3., hsigma=3., rdnoise="0.",
gain="1.", snoise="0.", sigscale=0.1, pclip=-0.5, grow=0)


#
# display the averaged spectra
#


display (inimage//".kavg.imh",
1, erase=yes, border_erase=no, select_frame=yes, repeat=no, fill=no,
zscale=yes, contrast=0.25, zrange=yes, nsample_line=5, xcenter=0.5,
ycenter=0.8, xsize=1., ysize=1., xmag=2., ymag=1., order=0, z1=INDEF,
z2=INDEF, ztrans="linear", lutfile="")

display (inimage//".havg.imh",
1, erase=no, border_erase=no, select_frame=yes, repeat=no, fill=no,
zscale=yes, contrast=0.25, zrange=yes, nsample_line=5, xcenter=0.5,
ycenter=0.6, xsize=1., ysize=1., xmag=2., ymag=1., order=0, z1=INDEF,
z2=INDEF, ztrans="linear", lutfile="")

display (inimage//".javg.imh",
1, erase=no, border_erase=no, select_frame=yes, repeat=no, fill=no,
zscale=yes, contrast=0.25, zrange=yes, nsample_line=5, xcenter=0.5,
ycenter=0.4, xsize=1., ysize=1., xmag=2., ymag=1., order=0, z1=INDEF,
z2=INDEF, ztrans="linear", lutfile="")

display (inimage//".iavg.imh",
1, erase=no, border_erase=no, select_frame=yes, repeat=no, fill=no,
zscale=yes, contrast=0.25, zrange=yes, nsample_line=5, xcenter=0.5,
ycenter=0.2, xsize=1., ysize=1., xmag=2., ymag=1., order=0, z1=INDEF,
z2=INDEF, ztrans="linear", lutfile="")

end