; File: silver.pro
; Author: Erik Brisson

; Plot of data in the file 'silver.dat'
; The values of the third column are plotted against the values of the second column

; Get data
n = 58
allin = fltarr(4,n)
openr, 1, 'silver_mod.dat'
readf, 1, allin
close, 1
x = allin(1,*)
y = allin(2,*)
for i=0, n-1 do begin
    print, x(i), y(i)
endfor

; Set up some specific stuff for the graphics system
DEVICE, DECOMPOSED=0
LOADCT,  5
PRE_PCOLOR = !P.COLOR
!P.COLOR=0

; Do the plot (first everything in black, then overwrite the plot in red)
plot, psym=1, background=255, x, y
oplot, psym=1, color=100, x, y
xyouts, 400, 200, color=100, '+'
xyouts, 420, 200, 'silver.dat'

end