数组输出问题
目下，我有一组数据，存入（143，86）数组中，读取文件写入文件13数据正常，在把相同数组写入文件14，就不正常了，不知什么情况。文件如下：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 13,14行
应该为
read(11,*)freq(l,i),f(l,i)
write(13,*)l,f(l,i)
页:
[1]