procedure lplot_feiii (patch, field, mask, num, z)

# This procedure plots a CNOC2 object spectrum in the REST frame
# for the given redshift, and overplots vertical lines at the positions
# of FeIII lines provided by Katya Verner.

string patch	{prompt="4-digit patch code"}
string field	{prompt="2-digit ALPHANUMERIC field code"}
string mask	{prompt="1-digit mask code"}
string num	{prompt="Object number in mask"}
real   z	{prompt="Redshift"}

begin

string ptch, fld, msk, nm, igicommands
real zed, xmin, xmax, ymin, ymax, junk

ptch = patch
fld = field
msk = mask
nm = num
zed = z

del ("oog", ver-, >& "dev$null")
del ("goo", ver-, >& "dev$null")
del ("goog", ver-, >& "dev$null")

listpix ("/d/larkya/lin/specarchive/spec"//ptch//"/"//ptch//fld//msk//".cleancal[*,"//nm//"]", wcs="world", >> "oog")

filecalc ("oog", "$1/"//(1+zed)//";$2", >> "goo")

head ("goo", nlines=1) | scan (xmin, junk)
tail ("goo", nlines=1) | scan (xmax, junk)
sort ("goo", col=2, num+, rev+) | head (junk, nlines=1) | scan (junk, ymax)

igicommands = "goog"
print ("limits ",xmin,xmax,0.,ymax, >igicommands)
print ("data goo; xcol 1; ycol 2", >>igicommands)
print ("box; connect", >>igicommands)
print ("xlabel "//ptch//fld//msk//nm//" z="//zed, >>igicommands)
print ("ltype 1", >>igicommands)
print ("move 4079.70 0; draw 4079.70 "//ymax, >>igicommands)
print ("move 4658.05 0; draw 4658.05 "//ymax, >>igicommands)
print ("move 4667.01 0; draw 4667.01 "//ymax, >>igicommands)
print ("move 4701.53 0; draw 4701.53 "//ymax, >>igicommands)
print ("move 4733.91 0; draw 4733.91 "//ymax, >>igicommands)
print ("move 4754.69 0; draw 4754.69 "//ymax, >>igicommands)
print ("move 4769.43 0; draw 4769.43 "//ymax, >>igicommands)
print ("move 4777.68 0; draw 4777.68 "//ymax, >>igicommands)
print ("move 4881.00 0; draw 4881.00 "//ymax, >>igicommands)
print ("move 4930.54 0; draw 4930.54 "//ymax, >>igicommands)
print ("move 5011.26 0; draw 5011.26 "//ymax, >>igicommands)
print ("move 5084.77 0; draw 5084.77 "//ymax, >>igicommands)
print ("move 5270.40 0; draw 5270.40 "//ymax, >>igicommands)
print ("move 5411.98 0; draw 5411.98 "//ymax, >>igicommands)
print ("ltype 0", >>igicommands)
print ("move "//5577/(1+zed)//" 0; draw "//5577/(1+zed)//" "//ymax, >>igicommands)
print ("move "//5892/(1+zed)//" 0; draw "//5892/(1+zed)//" "//ymax, >>igicommands)

igi (initcmd="input "//igicommands//" ; end", >"dev$null")

del ("oog", ver-)

end