问题是这样的:
数据样本为:
1 5.402000 -4.100000 6.486000
2 5.106000 10.426000 14.167000
3 8.123000 20.117000 -1.005000
4 -11.142000 17.711000 3.485000
5 -17.692000 8.982000 -10.812000
.....
160 ...
每个文件包含160个ID, 一共有几万个那么多。但是不需要全部使用,我单挑选那些文件号为10的倍数的来用,这样需要处理的文件个数为几千个。我希望结果是以ID号作为RESULT_(ID号)文件名,然后里面包含提取出几千个的数据。
我参考论坛里另外一个帖子(数据分类问题求助),仿照着fcode同学给的例子写了一个:
[Fortran] 纯文本查看 复制代码 Program main
Implicit None
Character(160) :: Filename
Character(160) :: fn1,fn2,fn3,tmp
Integer :: I, J,num
real :: A(4), B(4), C(4)
integer ios
fn1='c6c8_ID_'
fn2='c6_'
fn3='c8_'
num=160 ! total number of result files
print*, 'Running'
do J = 1, num
write (tmp,*)J
C print*, J
open (j,File="D:\c6c8\"//trim(adjustl(fn1))//trim(adjustl(tmp))//'.dat',status='unknown') ! save result here
write (j,'(3f8.3)') 0.0,0.0,1.0
enddo
open (20, File='IDlist.dat',status='old') ! 里面存放的是10的倍数,共计9200个。
do I = 1, 9200
read (20, *,iostat=ios) Filename
if (ios /=0) then
exit
endif
c print*, filename
open (200,File="D:\C6\"//trim(adjustl(fn2)) //trim(adjustl(filename))//'.dat')
open (300,File="D:\C8\"//trim(adjustl(fn3)) //trim(adjustl(filename))//'.dat')
Do J = 1, num
Read (200, *) A(:) ! first row=ID number, then x,y,z
Read (300, *) B(:) ! first row=ID number
C(:)=B(:)-A(:)
c print*, c(:)
Write (c(1), '(3f8.3)') C(2:4) ! first row=ID number
Close (300)
Close (200)
enddo
enddo
Close (20)
Do J = 1, num
Close (J)
End Do
End Program
请您指导如何修改,谢谢! |