; @notes-stars


input = fltarr(9,348-3)
s=strarr(1)

openr, 1, 'dfsc4761.alf'
readf, 1, s
print, s
readf, 1, s
print, s
readf, 1, s
print, s

readf, 1, input
close, 1

ynum = input(0,*)
ymag = input(3,*)


input = fltarr(9,346-3)

openr, 1, 'dfsc4763.alf'
readf, 1, s
print, s
readf, 1, s
print, s
readf, 1, s
print, s

readf, 1, input
close, 1

vnum = input(0,*)
vmag = input(3,*)
vxpos = input(1,*)
vypos = input(2,*)

; convert to 1D
ynum=ynum[*]
vnum=vnum[*]
ymag=ymag[*]
vmag=vmag[*]

; Select stars in both lists
;
; match will contain the index in ynum of stars that appear in both lists

match = vnum*0
for i=0, n_elements(vnum)-1 do   match[i] = where(ynum EQ vnum[i])
match = match(where(match NE -1))
doubles = ynum[match]


y = fltarr(n_elements(doubles))
v = fltarr(n_elements(doubles))
vx = fltarr(n_elements(doubles))
vy = fltarr(n_elements(doubles))

for i=0, n_elements(doubles)-1 do y[i] = ymag(where(ynum EQ doubles[i]))
for i=0, n_elements(doubles)-1 do v[i] = vmag(where(vnum EQ doubles[i]))

for i=0, n_elements(doubles)-1 do vx[i] = vxpos(where(vnum EQ doubles[i]))
for i=0, n_elements(doubles)-1 do vy[i] = vypos(where(vnum EQ doubles[i]))

;plot, vx, vy, psym=3
;;plot, v-y, y, psym=3, yr=[20,8], xr=[-1,3]
;;select = where((y lt 12.5) and (y gt 11.5))
;select = where((v-y lt 0) and (y lt 11.5))
;vx=vx(select)
;vy=vy(select)
;v=v(select)
;y=y(select)
;oplot, vx, vy, psym=4
;;oplot, v-y, y, psym=4

plot, v-y, y, psym=3, yr=[20,8], xr=[-1,3]
select = where((vx gt 101) and (vx lt 1400) and (vy gt 401) and (vy lt 1600) )
vx=vx(select)
vy=vy(select)
v=v(select)
y=y(select)
oplot, v-y, y, psym=4