PROGRAM	near2

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

1001	FORMAT (3f8.3)
	OPEN (unit=1, file='PCDall.pdb', status='OLD', readonly)
	OPEN (unit=2, file='wat.pdb', status='OLD', readonly)
	DO WHILE (.TRUE.)
	  READ (2, 1002, end=99) line
	  READ (line(31:54), 1001) x
	  WRITE (6, 1002)
	  WRITE (6, 1002) line
	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. 16.0 .AND. dd .GT. 1.) WRITE (6, 1002) line, SQRT(dd)
	endDO
9	REWIND (unit=1)
	endDO
99	CLOSE (unit=1)
	CLOSE (unit=2)
	end