program hcompare real hdat1(12,138), hdat2(12,138) open(unit=1,file='hdrift.dat.47011',type='UNKNOWN') open(unit=2,file='hdrift.dat.47021',type='UNKNOWN') do i=1,12 c write(*,*)'Reading ... i=',i do j=1,14 c write(*,*)'Reading ... j=',j if(j.eq.1) then kmax=7 m=0 lmin=(j-1)*10-m+1 lmax=lmin+kmax read(1,*)(hdat1(i,l),l=lmin,lmax) c write(*,100)(hdat1(i,l),l=lmin,lmax) read(2,*)(hdat2(i,l),l=lmin,lmax) c write(*,100)(hdat2(i,l),l=lmin,lmax) else kmax=9 m=2 lmin=(j-1)*10-m+1 lmax=lmin+kmax read(1,*)(hdat1(i,l),l=lmin,lmax) c write(*,101)(hdat1(i,l),l=lmin,lmax) read(2,*)(hdat2(i,l),l=lmin,lmax) c write(*,101)(hdat2(i,l),l=lmin,lmax) endif enddo enddo do i=1,138 write(30,102)i,(hdat1(j,i),j=1,12) write(31,102)i,(hdat2(j,i),j=1,12) enddo 100 format(8(f6.4,1x)) 101 format(10(f6.4,1x)) 102 format(i4,1x,12(f6.4,1x)) stop end