24小时热门版块排行榜

返回列表

贺仪

铁杆木虫 (著名写手)

应助: 93 (初中生)
贵宾: 0.636
金币: 5937
帖子: 2499
在线: 588.7小时
虫号: 995125

[交流] fortran 文件读取与写入问题。

大家好！最近写程序遇到一个比较奇怪的问题：由read读取完之后，在do循环内write写出的时候格式和文字显示是正确的。但是若在end do
后面再write输出就有问题了。以下是代码以及输出结果。write放在do循环里面得输出是S19，是正确的，而循环外面的那个只输出了S1，没有9了。这是怎么回事呢？期待高手解答！谢谢！

引用回帖:

do i= 1,natom
      read(301,332) atom,j,atomn(i),rname,nr,&
      & (coord(j,i),j=1,3),tmp,tmp,elemt(i)
      write(*,332) atom,j,trim(atomn(i)),rname,nr,&
      & (coord(j,i),j=1,3),tmp,tmp,elemt(i)
end do
  close(301)

输出的是：
ATOM    4  S19 C91    1    0.600  -2.293  -0.171  0.00  0.00
ATOM    4 C15  C91    1    -0.515  -9.800  -6.394  0.00  0.00
ATOM    4  S19 C91    1    -3.613  -0.506  -1.445  0.00  0.00
ATOM    4 C15  C91    1    2.933 0.963 5.719  0.00  0.00
ATOM    4  S19 C91    1    -4.966 4.043  -1.252  0.00  0.00
ATOM    4 C15  C91    1    -6.223  -3.564  -7.324  0.00  0.00
ATOM    4  S19 C91    1    -9.197 5.822  -2.474  0.00  0.00
ATOM    4 C15  C91    1    -2.510 7.390 4.536  0.00  0.00
ATOM    4  S19 C91    1    -10.533  10.378  -2.333  0.00  0.00
ATOM    4 C15  C91    1    -11.929 2.673  -8.250  0.00  0.00
ATOM    4  S19 C91    1    -14.781  12.150  -3.504  0.00  0.00

而

引用回帖:

do i= 1,natom
read(301,332) atom,j,atomn(i),rname,nr,&
      & (coord(j,i),j=1,3),tmp,tmp,elemt(i)
end do
   close(301)
do i=1,natom
   write(*,332) atom,j,trim(atomn(i)),rname,nr,&
      & (coord(j,i),j=1,3),tmp,tmp,elemt(i)
  end do

输出的是
ATOM    4  S1  C91    1    0.600  -2.293  -0.171  0.00  0.00
ATOM    4 C15  C91    1    -0.515  -9.800  -6.394  0.00  0.00
ATOM    4  S1  C91    1    -3.613  -0.506  -1.445  0.00  0.00
ATOM    4 C15  C91    1    2.933 0.963 5.719  0.00  0.00
ATOM    4  S1  C91    1    -4.966 4.043  -1.252  0.00  0.00
ATOM    4 C15  C91    1    -6.223  -3.564  -7.324  0.00  0.00
ATOM    4  S1  C91    1    -9.197 5.822  -2.474  0.00  0.00
ATOM    4 C15  C91    1    -2.510 7.390 4.536  0.00  0.00
ATOM    4  S1  C91    1    -10.533  10.378  -2.333  0.00  0.00
ATOM    4 C15  C91    1    -11.929 2.673  -8.250  0.00  0.00
ATOM    4  S1  C91    1    -14.781  12.150  -3.504  0.00  0.00
ATOM    4 C15  C91    1    -7.955  13.814 3.349  0.00  0.00

回复此楼