parameter (maxstr=50000) integer chip(maxstr), cc real v(maxstr), i(maxstr), x(maxstr), y(maxstr), ii open(1,file='phot-1.16/lmc_u2c501.phot',status='old') n=1 1 read(1,*,end=99) chip(n), x(n), y(n), a, a, a, a, a, a, a, a, v(n), a, a, a, a, a, a, a, a, i(n) n=n+1 goto 1 99 n=n-1 open(2,file='phot-1.15/lmc_u2c501.phot',status='old') 2 read(2,*,end=199) cc, xx, yy, a, a, a, a, a, a, a, a, vv, a, a, a, a, a, a, a, a, ii do j=1,n if (cc .eq. chip(j)) then if (abs(xx-x(j)) .lt. 0.2) then if (abs(yy-y(j)) .lt. 0.2) then print 101, cc, xx, yy, vv, ii, vv-v(j), ii-i(j) 101 format(i4,2f9.2,4f8.3) end if end if end if end do goto 2 199 continue stop end