本帖最后由 wxy 于 2017-3-4 17:41 编辑
请问我想提取文件中的指定行数据,应该怎么改
[Fortran] 纯文本查看 复制代码 program find_close
implicit none
integer m,i
integer :: sol , country
integer , allocatable :: site(:)
real , allocatable :: lat(:) , long(:)
open(8,file="soil.sol")
open(10,file="find_close.txt")
!read(8,*)m !//m无法再文件中获取,需手动输入或获取文件行数减去3,得到。
m=414503 !文件行数
allocate(site(m),lat(m),long(m))
do i = 1,m
if(mod(i-3,9)==0)then
read(8,*) site(i) , country , lat(i) , long(i) , sol
end if
enddo
write(10,*) site(i) , country , lat(i), long(i) , sol
end program find_close
|