Fortran Coder

标题: 数组输出问题 [打印本页]
作者: alohomora100    时间: 2018-1-18 16:19
标题: 数组输出问题
目下，我有一组数据，存入（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

1.txt

336.27 KB, 下载次数: 1

2.txt

192.16 KB, 下载次数: 1
作者: kyra    时间: 2018-1-18 18:50
13,14行
应该为
  read(11,*)freq(l,i),f(l,i)
  write(13,*)l,f(l,i)




欢迎光临 Fortran Coder (http://bbs.fcode.cn/) Powered by Discuz! X3.2