parameter(maxstars=10000)
	dimension x(maxstars),y(maxstars),a(maxstars),e(maxstars),
     &  n(maxstars),istat(maxstars)
	hit2=25.
	do 178 i=1,maxstars
	istat(i)=0
178	continue
	istar=0
	open(unit=1,file='n346vxy.dat',status='old')
1       read(1,*,end=99) x1,y1,a1,e1,n1
	istar=istar+1
	if(istar.gt.maxstars) stop 'too many stars'
	x(istar)=x1
	y(istar)=y1
	a(istar)=a1
	e(istar)=e1
	n(istar)=n1
	go to 1
99	close(unit=1)
	do 100 i=1,istar
	if (i.eq.144) then
	iflag=1
	else
	iflag=0
	end if
C We are surveying star "i"..
	do 100 j=1,istar
	if(i.eq.j) go to 100
	dm=a(i)-a(j)
	if(dm.lt.-1.0) go to 100
	d2=(x(j)-x(i))**2+(y(j)-y(i))**2
	if(d2.lt.hit2) then
	istat(i)=1
	if(iflag.eq.1) write(6,*) i,j,dm,d2
	end if
100	continue
	continue
	do 103 i=1,istar
	if(istat(i).eq.1) go to 103
	write(3,101) x(i),y(i),a(i),e(i),n(i)
101	format(2f8.2,f8.3,f8.3,i5)
103	continue
	stop
	end