目的:把一个叫poro的文件的数据读按顺序读取出来成一维,以后就可以直接用poro(i)来表示第i个数据
问题:写完了运行后没看见数据,是空的。求指导
[Fortran] 纯文本查看 复制代码 program shujv1
implicit none
integer::NxD,NyD,NzD
integer,allocatable::poro(:)
! Local Variables
integer::i,j,k,l ! Counter
integer::stat
character::line*160
character::key*12
NxD=20
NyD=30
NzD=1
!allocate(poro(NxD*NyD*NzD))
! Execution Part
k=0
open(unit=4,file='.\shujv2\poro.dat',iostat=stat)
if(stat==0)then
do
read(4,'(a)',iostat=stat)line
! if (stat/=0)exit
if ((line(1:2)/='--').and.(line/=' ').and.(line(1:6)/='poro'))then
k=k+1
end if
end do
!rewind(4)
allocate(poro(k))
! if(stat==0)
do i=1,k
read(4,*)poro(i)
end do
! end if
do i=1,k
write(*,*)poro(i)
end do
endif
end program
|