procedure teluric (inlist, outpref, flat)

string inlist {prompt="list of standard images"}
string outpref {prompt="output image name"}
string flat { prompt="flat field name"}
string geomap{ prompt="geometrical transformation image"}

struct *instr

begin

	string infile, outfile
	
	instr=inlist
#	imdel ('t*.imh')

	imreplace('@'//inlist,0.0, lower=32000)

        imcomb ('@'//inlist,"temp.med",
plfile="", sigma="", logfile="STDOUT", combine="median", reject="pclip",
project=no, outtype="real", offsets="", masktype="none", maskvalue=0.,
blank=0., scale="none", zero="none", weight="none", statsec="", expname="",
lthreshold=INDEF, hthreshold=INDEF, nlow=1, nhigh=1, mclip=yes, lsigma=3.,
hsigma=3., rdnoise="0.", gain="1.", sigscale=0.1, pclip=-0.5, grow=0)

	
	while(fscan(instr, infile) !=EOF) {
		outfile='t'//infile
		imarith(infile, "-" , "temp.med.imh", outfile)
		imarith(outfile, "/" , flat, outfile)
#		display( outfile, 1)
		hedit(outfile, "dispaxis", "2", ver-, add+)

transform (outfile,
outfile, geomap , database="database", interptype="poly5",
x1=INDEF, x2=INDEF, dx=INDEF, nx=INDEF, xlog=no, y1=INDEF, y2=INDEF, dy=INDEF,
ny=INDEF, ylog=no, flux=yes, logfiles="STDOUT,logfile")


apall (outfile ,
1, output="", references="", profiles="", interactive=no ,format="oned",
find=yes,
 recenter=yes, resize=no, edit=yes, trace=yes, fittrace=no, extract=yes,
extras=no, review=no, line=INDEF, nsum=1, lower=-5.0,upper=5.0,
apidtable="", b_function="chebyshev", b_order=1, b_sample="-28:-8,8:28",
b_naverage=1, b_niterate=5, b_low_reject=2., b_high_rejec=2., b_grow=0.,
width=5., radius=10., threshold=0., minsep=5., maxsep=1000.,
order="increasing", apertures="", npeaks=INDEF, shift=yes, llimit=INDEF,
ulimit=INDEF, ylevel=0.1, peak=yes, bkg=yes, r_grow=0., avglimits=no,
t_nsum=10, t_step=10, t_function="cheb", t_order=2, t_sample="*",
t_naverage=1, t_niterate=5, t_low_reject=2., t_high_rejec=2., t_grow=0.,
background="fit", skybox=1, weights="variance", clean=yes, saturation=INDEF,
readnoise="300", gain="30", lsigma=4., usigma=4., nsubaps=1)

	}
	imdel('temp.med.imh')

        imcomb ("t*.0001.imh",outpref,
plfile="", sigma="", logfile="STDOUT", combine="median", reject="pclip",
project=no, outtype="real", offsets="", masktype="none", maskvalue=0.,
blank=0., scale="median", zero="none", weight="none", statsec="", expname="",
lthreshold=INDEF, hthreshold=INDEF, nlow=1, nhigh=1, mclip=yes, lsigma=3.,
hsigma=3., rdnoise="0.", gain="1.", sigscale=0.1, pclip=-0.5, grow=0)

#	imdel('t*.imh')

	end