数组输出问题

alohomora100

目下，我有一组数据，存入（143，86）数组中，读取文件写入文件13数据正常，在把相同数组写入文件14，就不正常了，不知什么情况。文件如下：

[Fortran] 纯文本查看 复制代码

      program vacf 
      implicit none
      real(4),parameter::natom=143,nn=86
      integer::i,j,k,l
      real(4),dimension(natom,nn)::freq,f
      real(4)::fn
      open(11,file='power-spectrum-for-atom')
      open(12,file='fn.txt')
      open(13,file='1.txt')
      open(14,file='2.txt')
      do l=1,natom
        do i=1,nn
        read(11,*)freq(l,nn),f(l,nn)
        write(13,*)l,f(l,nn)
        enddo
      enddo
      do i=1,nn
        fn=0.d0
        do l=1,3
          fn=fn+f(l,i)
        enddo
        write(12,*)fn/3.d0
      enddo
      do i=1,natom
        do j=1,nn
      write(14,*)f(i,j)
        enddo
      enddo
      close(11)
      close(12)
      close(13)
      close(14)
      end

kyra

13,14行
应该为
  read(11,*)freq(l,i),f(l,i)
  write(13,*)l,f(l,i)