PROGRAM	near

	CHARACTER*66	line
	REAL		x(3), z(3)

	READ (5, 1001) x
1001	FORMAT (3f8.3)
	OPEN (unit=1, file='PCD.pdb', status='OLD', readonly)
	DO WHILE (.TRUE.)
	  READ (1, 1002, end=9) line
1002	  FORMAT (a, f8.2)
	  READ (line(31:54), 1001) z
	  dd = (x(1)-z(1))**2 + (x(2)-z(2))**2 + (x(3)-z(3))**2
	  IF (dd .LE. 20.0) WRITE (6, 1002) line, SQRT(dd)
	endDO
9	CLOSE (unit=1)
	end