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

program overlap

real, dimension(88,88) :: ol
character(len=80) :: line,lfm
integer :: i, j, k, l, 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
      l=0
      do j=i,88
         if (l < 5) l=l+1
         read(idat,'(a)') line
         write(lfm, '(a,i0,a)') '(tr1,',l,'(tr2,f9.6))'
!        print *, lfm, min(i+4,j)
         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

Originally posted by snoopyzhao at 2011-05-02 15:42:17:
¸Õ²ÅÄǸöÀïÃæ¹ØÓÚ lfm µÄ¸³ÖµÓеã¶ùÎÊÌ⣬¾¡¹Ü²»Ó°ÏìʹÓ㬺Ǻǡ­¡­ÏÂÃæÕâ¸öÊǸÄÕýºóµÄ¡£

[code]
program overlap

real, dimension(88,88) :: ol
character(len=80) :: line,lfm
integer :: i, j, k, l, ...

Originally posted by kathy2008 at 2011-05-02 15:57:54:
¼ÌÐøÇë½Ì: idat = 12£¬ÊÇÔõôÑù
È·¶¨°Ñ12¸¶¸øidatÕâ¸ö±äÁ¿µÄÄØ£¿Ò²¾ÍÊÇ˵ΪʲôҪѡÔñ12£¬¶ø²»ÊÇÆäËûÊý×Ö£¬Ð»Ð»¡£