implicit double precision (a-h,o-z)
	character*1 sign
	character*120 inline1,inline2
*  SECONDS TO RADIANS
      PARAMETER (S2R=0.72722052166430398D-4)
*  ARC SECONDS TO RADIANS
      PARAMETER (AS2R=0.48481368110953599D-5)
C
*  SECONDS TO RADIANS
      STOR(RH,RM,RS)=S2R*(60D0*(60D0*RH+RM)+RS)
*  ARC SECONDS TO RADIANS
      ASTOR(DD,DM,DS)=AS2R*(60D0*(60D0*DD+DM)+DS)

C matchrd.f matches the n346vubr.dat with the n346sch.cat...
	nv=0
	nb=0
	nr=0
	nu=0
	open(unit=1,file='n346vubr.dat',status='old')
	open(unit=2,file='n346sch.cat',status='old')
1	read(1,2,end=99) inline1
2	format(a120)
	if(inline1(1:1).eq.'#') go to 1
	read(inline1(39:61),3) ihr,imin,sec,sign,ideg,iamin,asec
3	format(i2,1x,i2,1x,f5.2,1x,a1,i2,1x,i2,1x,f4.1)
	if(sign.eq.'-') then
	ideg=-ideg
	iamin=-iamin
	asec=-asec
	end if
	ra=stor(ihr,imin,sec)
	dec=astor(ideg,iamin,asec)
	write(6,4) ihr,imin,sec,ideg,iamin,asec
	rewind(2)
C Now go hunting for it....
22	read(2,2,end=1) inline2
	read(inline2(1:23),3) ihr,imin,sec,sign,ideg,iamin,asec
	if(sign.eq.'-') then
        ideg=-ideg
        iamin=-iamin
        asec=-asec
        end if
        rasch=stor(ihr,imin,sec)
        decsch=astor(ideg,iamin,asec)
	call apa(ra,dec,rasch,decsch,dist,pa)
	d=dist/as2r
	if(d.gt.2.) go to 22
4	format(i4,1x,i4,1x,f6.3,1x,i4,1x,i4,1x,f6.2)
	write(6,4) ihr,imin,sec,ideg,iamin,asec
	write(6,*) d
	write(6,2) inline1
	write(6,2) inline2
C We have a winner...maybe...now what?
	read(inline2(24:120),*) v,bmv,umb,vmr,it,ve,v1,bmve,bmv1,
     &      umbe,umb1,vmre
	read(inline1(18:32),*) vi,vie
	read(inline1(63:120),*) ui,uie,n1,bi,bie,n2,ri,rie,n3
C if the errors are too great, don't bother...
C First, V:
	if(vie.gt.0.05) go to 12
	if(ve.gt.0.05) go to 12
	write(7,133) bmv,vi-v,d,vi,vie,v,ve
133	format(7f9.3)
C Now, U:
12	if(n1.eq.0) go to 13
	if(uie.gt.0.1) go to 13
	if(dsqrt(ve**2+bmve**2+umbe**2).gt.0.1) go to 13
	u=v+bmv+umb
	write(8,133) umb,ui-u,d,ui,u
C Now, B
13	if(n2.eq.0) go to 14
	if(bie.gt.0.05) go to 14
	if(dsqrt(ve**2+bmve**2).gt.0.1) go to 14
	b=v+bmv
	write(9,133) bmv,bi-b,d,bi,b
C finally, R:
14	if(n3.eq.0) go to 15
	if(rie.gt.0.05) go to 15
	if(dsqrt(ve**2+vmre**2).gt.0.1) go to 15
	r=v-vmr
	write(10,133) vmr,ri-r,d,ri,r
15	continue
	go to 1
99	continue
	stop
	end
        subroutine apa(rac,dcc,rao,dco,dist,pa)
c computes the distance (delt) and position (pa) along the
c sky FROM (rac,dec) to (rao,dco)
c written aug 14, 1987
C rewritten August 21, 1989, to do using "standard coords", rather
C than the old and peculiar (but possibly equivalent!) way
c plm
      implicit double precision (a-h,o-z)
        parameter (PI=3.141592653589793238462643)
C CONVERT TO PROPER COORDS
        CALL DS2TP(RAO,DCO,RAC,DCC,XI,ETA)
        pa=0.
        if(xi*eta.ne.0.)PA=DATAN2(XI,ETA)
        if(pa.le.0.) pa=pa+2.*pi
        DIST=DSQRT(XI**2+ETA**2)
        return
        end
        include '/tofu/home/massey/survey/astromsub.f'