Originally posted by snoopyzhao at 2011-05-02 11:13:46:
下载不下来，呵呵……

如果你知道矩阵的行数或列数，自己计算一下其中的数据个数，然后把它们全部读到一个一维数组中，或者直接到二维数组也可以……

类似于这样……

read(*,*) (a(i), i = 1, 100)

...

program overlap

real, dimension(88,88) :: ol
character(len=80) :: line,lfm
integer :: i, j, k, ios
integer :: idat = 12

   open(unit=idat, file='overlap.txt')

   do i=1,88
      do j=1,88
         ol(i,j) = 0.0
      end do
   end do

   do i=1,88,5
      read(idat,'(a)', iostat = ios ) line
      if (ios < 0) exit
      read(idat,'(a)') line
      do j=i,88
         read(idat,'(a)') line
         write(lfm, '(a,i0,a)') '(tr1,',min(i+4,j),'(tr2,f9.6))'
         read(line,lfm) (ol(j,k),k=i,min(i+4,j))
!        write(*,lfm) (ol(j,k),k=i,min(i+4,j))  
      end do
      read(idat,'(a)', iostat = ios) line
      if (ios < 0) exit
   end do

end program overlap