procedure plotcors (aperture)

# Remember to first do   task plotcors = plotcors.cl   in IRAF.

int aperture		{prompt="Aperture for which to plot xcor results"}
bool cleanup=yes	{prompt="Delete temporary xcor images when done?"}

begin

int	ap
bool	cu
string	opt

ap = aperture
cu = cleanup

del ("list_corn", ver-, >& "dev$null")
del ("plotcors.mc", ver-, >& "dev$null")
imdel ("cor"//ap//"*", ver-, >& "dev$null")

# get red of IRAF banner on plot
opt=splot.options
splot.options=opt//" nosysid"

# copy xcor results to individual images
imcopy ("cor[*,"//ap//",2]", "cor"//ap//"2")
imcopy ("cor[*,"//ap//",3]", "cor"//ap//"3")
imcopy ("cor[*,"//ap//",4]", "cor"//ap//"4")
imcopy ("cor[*,"//ap//",5]", "cor"//ap//"5")
files ("cor"//ap//"*", >> "list_corn")

print ("Plotting xcor results in format:  2 4  for aperture #"//ap)
print ("                                  3 5")
bplot ("@list_corn", band=1, >G "plotcors.mc")
gkimosaic ("plotcors.mc", dev="stdgraph", nx=2, ny=2, interact-, fill+)
print("Hardcopy: '=gcur' '=' 'q'  Replot: 'gkimosaic plotcors.mc dev=stdgraph'")

if (cu) {
	del ("list_corn", ver-, >& "dev$null")
	imdel ("cor"//ap//"*", ver-, >& "dev$null")
}

# return splot options to original state
splot.options=opt

end