procedure lplot (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 common spectral features: dotted lines for emission, dashed lines
# for absorption, dot-dash lines for features in the spectra of
# early-type galaxies which empirically mimic emission lines,
# and solid lines for 5577 and 5892 A night sky emission.

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 911.75 0; draw 911.75 "//ymax, >>igicommands)
print ("move 972.54 0; draw 972.54 "//ymax, >>igicommands)
print ("move 1026 0; draw 1026 "//ymax, >>igicommands)
print ("move 1034 0; draw 1034 "//ymax, >>igicommands)
print ("move 1216 0; draw 1216 "//ymax, >>igicommands)
print ("move 1240 0; draw 1240 "//ymax, >>igicommands)
print ("move 1302 0; draw 1302 "//ymax, >>igicommands)
print ("move 1335 0; draw 1335 "//ymax, >>igicommands)
print ("move 1397 0; draw 1397 "//ymax, >>igicommands)
print ("move 1402 0; draw 1402 "//ymax, >>igicommands)
print ("move 1549 0; draw 1549 "//ymax, >>igicommands)
print ("move 1640 0; draw 1640 "//ymax, >>igicommands)
print ("move 1663 0; draw 1663 "//ymax, >>igicommands)
print ("move 1858 0; draw 1858 "//ymax, >>igicommands)
print ("move 1909 0; draw 1909 "//ymax, >>igicommands)
print ("move 2142 0; draw 2142 "//ymax, >>igicommands)
print ("move 2326 0; draw 2326 "//ymax, >>igicommands)
print ("move 2471 0; draw 2471 "//ymax, >>igicommands)
print ("move 2798 0; draw 2798 "//ymax, >>igicommands)
print ("move 3133 0; draw 3133 "//ymax, >>igicommands)
print ("move 3346 0; draw 3346 "//ymax, >>igicommands)
print ("move 3426 0; draw 3426 "//ymax, >>igicommands)
print ("move 3727 0; draw 3727 "//ymax, >>igicommands)
print ("move 3869 0; draw 3869 "//ymax, >>igicommands)
print ("move 3969 0; draw 3969 "//ymax, >>igicommands)
print ("move 4102 0; draw 4102 "//ymax, >>igicommands)
print ("move 4340 0; draw 4340 "//ymax, >>igicommands)
print ("move 4471 0; draw 4471 "//ymax, >>igicommands)
print ("move 4686 0; draw 4686 "//ymax, >>igicommands)
print ("move 4861 0; draw 4861 "//ymax, >>igicommands)
print ("move 4959 0; draw 4959 "//ymax, >>igicommands)
print ("move 5007 0; draw 5007 "//ymax, >>igicommands)
print ("move 5199 0; draw 5199 "//ymax, >>igicommands)
print ("move 5876 0; draw 5876 "//ymax, >>igicommands)
print ("move 5977.0 0; draw 5977.0 "//ymax, >>igicommands)
print ("move 6300.23 0; draw 6300.23 "//ymax, >>igicommands)
print ("move 6364 0; draw 6364 "//ymax, >>igicommands)
print ("move 6548.06 0; draw 6548.06 "//ymax, >>igicommands)
print ("move 6562.82 0; draw 6562.82 "//ymax, >>igicommands)
print ("move 6583.57 0; draw 6583.57 "//ymax, >>igicommands)
print ("move 6717.0 0; draw 6717.0 "//ymax, >>igicommands)
print ("move 6731.3 0; draw 6731.3 "//ymax, >>igicommands)
print ("move 7135.8 0; draw 7135.8 "//ymax, >>igicommands)
print ("move 7237.3 0; draw 7237.3 "//ymax, >>igicommands)
print ("move 7325 0; draw 7325 "//ymax, >>igicommands)
print ("move 9069 0; draw 9069 "//ymax, >>igicommands)
print ("ltype 2", >>igicommands)
print ("move 1199.9674 0; draw 1199.9674 "//ymax, >>igicommands)
print ("move 1206.5 0; draw 1206.5 "//ymax, >>igicommands)
print ("move 1259.519 0; draw 1259.519 "//ymax, >>igicommands)
print ("move 1260.4221 0; draw 1260.4221 "//ymax, >>igicommands)
print ("move 1302.1685 0; draw 1302.1685 "//ymax, >>igicommands)
print ("move 1304.3702 0; draw 1304.3702 "//ymax, >>igicommands)
print ("move 1334.5323 0; draw 1334.5323 "//ymax, >>igicommands)
print ("move 1393.755 0; draw 1393.755 "//ymax, >>igicommands)
print ("move 1402.77 0; draw 1402.77 "//ymax, >>igicommands)
print ("move 1526.7066 0; draw 1526.7066 "//ymax, >>igicommands)
print ("move 1608.4511 0; draw 1608.4511 "//ymax, >>igicommands)
print ("move 1670.7874 0; draw 1670.7874 "//ymax, >>igicommands)
print ("move 1854.7164 0; draw 1854.7164 "//ymax, >>igicommands)
print ("move 1862.7895 0; draw 1862.7895 "//ymax, >>igicommands)
print ("move 2249.88 0; draw 2249.88 "//ymax, >>igicommands)
print ("move 2260.78 0; draw 2260.78 "//ymax, >>igicommands)
print ("move 2344.21 0; draw 2344.21 "//ymax, >>igicommands)
print ("move 2374.46 0; draw 2374.46 "//ymax, >>igicommands)
print ("move 2382.77 0; draw 2382.77 "//ymax, >>igicommands)
print ("move 2586.65 0; draw 2586.65 "//ymax, >>igicommands)
print ("move 2600.17 0; draw 2600.17 "//ymax, >>igicommands)
print ("move 2796.35 0; draw 2796.35 "//ymax, >>igicommands)
print ("move 2803.53 0; draw 2803.53 "//ymax, >>igicommands)
print ("move 2852 0; draw 2852 "//ymax, >>igicommands)
print ("move 3742 0; draw 3742 "//ymax, >>igicommands)
print ("move 3799 0; draw 3799 "//ymax, >>igicommands)
print ("move 3835 0; draw 3835 "//ymax, >>igicommands)
print ("move 3889 0; draw 3889 "//ymax, >>igicommands)
print ("move 3934 0; draw 3934 "//ymax, >>igicommands)
print ("move 3968 0; draw 3968 "//ymax, >>igicommands)
print ("move 3970 0; draw 3970 "//ymax, >>igicommands)
print ("move 4227 0; draw 4227 "//ymax, >>igicommands)
print ("move 4304 0; draw 4304 "//ymax, >>igicommands)
print ("move 4384 0; draw 4384 "//ymax, >>igicommands)
print ("move 5175 0; draw 5175 "//ymax, >>igicommands)
print ("move 5269 0; draw 5269 "//ymax, >>igicommands)
print ("move 5332 0; draw 5332 "//ymax, >>igicommands)
print ("move 5401 0; draw 5401 "//ymax, >>igicommands)
print ("move 5783 0; draw 5783 "//ymax, >>igicommands)
print ("move 5893 0; draw 5893 "//ymax, >>igicommands)
print ("move 6235 0; draw 6235 "//ymax, >>igicommands)
print ("ltype 4", >>igicommands)
print ("move 4025 0; draw 4025 "//ymax, >>igicommands)
print ("move 4505 0; draw 4505 "//ymax, >>igicommands)
print ("move 4555 0; draw 4555 "//ymax, >>igicommands)
print ("move 4625 0; draw 4625 "//ymax, >>igicommands)
print ("move 5051 0; draw 5051 "//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