pro contrast
file=pickfile(filter="*.fts")
a=readfits(file)
ss=size(a)
erase
tvscl,a
print,"Input Y down"
cursor,x,y,3,/device
y1=y
print,"Input Y up"
cursor,x,y,3,/device
y2=y
print,y1,y2
erase
tvscl,a(0:400, y1:y2)
profiles, a(0:400, y1:y2)
;tvscl,a(1135:185, y1:y2)
;profiles, a(1135:185, y1:y2)
end