procedure mkref  (inflat )
 
#
# Procedure to trace a stellar spectra to make reference images 
# by Richard Elston 8/17/93
#

string inflat {prompt="input stellar specrtra to trace"}


begin

	string inimage, outname
	inimage=inflat
	outname='ref'

#
# make the subimages for each band
#

	hedit (inimage, "dispaxis", 1, add+, ver-)
	imdel(outname//".?.imh")
	imcopy(inimage//"[1:256,23:87]", outname//".k")
	imcopy(inimage//"[1:256,79:165]", outname//".h")
	imcopy(inimage//"[1:256,129:209]", outname//".j")
	imcopy(inimage//"[1:256,165:253]", outname//".i")
#
# magnify the input specta
#
	magnify(outname//".?.imh",outname//".?.imh",1.0,2.0,flux-,
		interpol="nearest")
	
#
# Aprnomalize the flats
#
	apnormal(outname//".k", outname//".ktemp",
references="", interactive=yes, find=yes, recenter=yes, resize=no,
edit=yes, trace=yes, fittrace=yes, normalize=no, fitspec=yes, line=INDEF,
nsum=2, cennorm=yes, threshold=100., background="none", weights="none",
pfit="fit1d", clean=no, skybox=1, saturation=INDEF, readnoise="0.", gain="1.",
lsigma=4., usigma=4., function="legendre", order=7, sample="*", naverage=1,
niterate=3, low_reject=2.5, high_reject=2.5, grow=0.)

	apnormal(outname//".h", outname//".htemp",
references="", interactive=yes, find=yes, recenter=yes, resize=no,
edit=yes, trace=yes, fittrace=yes, normalize=no, fitspec=yes, line=INDEF,
nsum=10, cennorm=no, threshold=1000., background="none", weights="none",
pfit="fit1d", clean=no, skybox=1, saturation=INDEF, readnoise="0.", gain="1.",
lsigma=4., usigma=4., function="legendre", order=7, sample="*", naverage=1,
niterate=3, low_reject=2.5, high_reject=2.5, grow=0.)
	
	apnormal(outname//".j", outname//".jtemp",
references="", interactive=yes, find=yes, recenter=yes, resize=no,
edit=yes, trace=yes, fittrace=yes, normalize=no, fitspec=yes, line=INDEF,
nsum=10, cennorm=no, threshold=1000., background="none", weights="none",
pfit="fit1d", clean=no, skybox=1, saturation=INDEF, readnoise="0.", gain="1.",
lsigma=4., usigma=4., function="legendre", order=7, sample="*", naverage=1,
niterate=3, low_reject=2.5, high_reject=2.5, grow=0.)
	
	apnormal(outname//".i", outname//".itemp",
references="", interactive=yes, find=yes, recenter=yes, resize=no,
edit=yes, trace=yes, fittrace=yes, normalize=no, fitspec=yes, line=INDEF,
nsum=10, cennorm=no, threshold=1000., background="none", weights="none",
pfit="fit1d", clean=no, skybox=1, saturation=INDEF, readnoise="0.", gain="1.",
lsigma=4., usigma=4., function="legendre", order=7, sample="*", naverage=1,
niterate=3, low_reject=2.5, high_reject=2.5, grow=0.)


	imdel(outname//".?temp.imh")


	end