这个需求可以很简单的实现,楼主您写复杂了。
如果后面不再需要,ns 和 qa 都可以不用数组。
[Fortran] 纯文本查看 复制代码 program mascentre
implicit none
Integer , parameter :: N = 10000
Real(8) :: qa(N,3)
Integer :: ns(N) , i
character::u
open(11,file='tmp100.txt')
open(12,file='mass_centre.txt')
do i=1,N
read(11,*) u ,u , ns(i),qa(i,:)
write(12,'(i6,3(f13.6,1x))') ns(i) , qa(i,:)
enddo
close(11)
close(12)
end program mascentre |