[交流] 【求助】如何写FORTRAN程序实现求平均最近邻距离已有2人参与

Originally posted by snoopyzhao at 2011-02-22 02:55:14:
你的问题是如何把这些数都读入，呵呵……

另外，你的每个原子就一行，还是有很多行？

你的数据有相应的格式规范吗？

Originally posted by snoopyzhao at 2011-02-22 02:55:14:
你的问题是如何把这些数都读入，呵呵……

另外，你的每个原子就一行，还是有很多行？

你的数据有相应的格式规范吗？

Originally posted by snoopyzhao at 2011-02-22 03:17:12:
你没有明白我的意思

我的问题有二：

1）你的每个原子是一行，还是有换行？你的这些数据是从哪里来的，别的软件输出的，还是你自己算的？

2）你的数据有固定的格式吗？比如 I2F6.2 之类的……

最好给 ...

Originally posted by snoopyzhao at 2011-02-22 06:19:48:
可能比较麻烦，回头写一个看看吧，呵呵

program ex

implicit none
character(len = 128) :: line, fm
integer :: ios, len_line, i, j
real, dimension(5000) :: table

do
   read (*,'(a)', iostat = ios) line
   if (ios < 0) exit
   if (index(line, 'nearest neighbor table') /= 0) exit
end do

j = 0

do
   read (*,'(a)', iostat = ios) line
   if (ios < 0) exit
   if (index(line, 'LATTYP') /= 0) exit
   if (len(trim(line)) == 0) exit
   len_line = len(trim(line(27:)))
   write(fm,'(a,i0,a)') '(', len_line/9,'(tr5,f4.2))'
   read (line(27:), fm) (table(i), i=j+1,j+len_line/9)
   j = j + len_line/9
end do

write (*,'(f4.2)') (table(i), i= 1, j)

write(*,*) sum(table(1:j))/j

end program ex