文件数据读写

xiaochen · 发表于 2022-9-19 16:38:13

目的：把一个叫poro的文件的数据读按顺序读取出来成一维，以后就可以直接用poro（i）来表示第i个数据
问题：写完了运行后没看见数据，是空的。求指导

[Fortran] 纯文本查看 复制代码

01program shujv1
02    implicit none
03   integer::NxD,NyD,NzD  
04   integer,allocatable::poro(:)                   
05 
06   ! Local Variables
07   integer::i,j,k,l         ! Counter
08   integer::stat
09   character::line*160
10   character::key*12
11    NxD=20
12   NyD=30
13   NzD=1
14!allocate(poro(NxD*NyD*NzD))
15   ! Execution Part
16k=0
17open(unit=4,file='.\shujv2\poro.dat',iostat=stat)
18if(stat==0)then
19do
20  read(4,'(a)',iostat=stat)line
21        ! if (stat/=0)exit
22         if ((line(1:2)/='--').and.(line/=' ').and.(line(1:6)/='poro'))then
23 
24k=k+1
25end if
26end do
27 
28!rewind(4)
29allocate(poro(k))
30 
31!  if(stat==0)
32    do i=1,k
33 
34     read(4,*)poro(i)
35 
36   end do
37 
38! end if
39 
40  do i=1,k
41 
42write(*,*)poro(i)
43 
44end do
45endif
46end program

fcode · 发表于 2022-9-19 17:32:46

取消 if (stat/=0)exit 和 rewind(4) 的注释

xiaochen · 发表于 2022-9-20 08:57:31

fcode 发表于 2022-9-19 17:32
取消 if (stat/=0)exit 和 rewind(4) 的注释

还是不太行，而且我发现我的k的值一直没变，是我的do有问题吗

fcode · 发表于 2022-9-20 09:00:57

给出你要读取的文件（或部分内容）

xiaochen · 发表于 2022-9-20 09:05:49

fcode 发表于 2022-9-20 09:00
给出你要读取的文件（或部分内容）

文件是dat格式的

poro
0.203
0.191
0.194
0.256
0.231
0.227
0.285
0.249
0.195
0.159
0.148
0.193
0.200
0.244
0.240
0.222
0.223
0.240

xiaochen · 发表于 2022-9-20 10:28:02

我试着改简单一点，现在有数据了，但是得出的是一个-842150451*600和文件中的完全不一样

[Fortran] 纯文本查看 复制代码

01open(unit=4,file='.\poro.dat')
02 
03allocate(poro(NxD*NyD))
04  read(4,*)line
05 
06 if ((line(1:2)/='/').and.(line/=' ').and.(line(1:6)/='poro'))then
07     k=k+1
08     
09     read(4,*)poro
10      
11 end if
12     close(4)
13     write(*,*)line
14     write(*,*)poro
15     end

伊芸贾 · 发表于 2022-9-20 11:59:54

[Fortran] 纯文本查看 复制代码

01program shujv1
02  implicit none
03  real,allocatable::poro(:)
04  integer::k,stat
05  character(len=160) :: line
06  open(unit=4,file='.\poro.dat',iostat=stat)
07  if(stat/=0) stop
08  k=0
09  Do
10    read(4,'(a)',iostat=stat) line
11    if (stat/=0)exit
12    if ((line(1:2)/='--').and.(line/=' ').and.(line(1:6)/='poro')) k=k+1
13  End Do
14  allocate(poro(k))
15  rewind(4)
16  k = 0
17  Do
18    read(4,'(a)',iostat=stat) line
19    if (stat/=0)exit
20    if ((line(1:2)/='--').and.(line/=' ').and.(line(1:6)/='poro')) then
21      k = k + 1
22      read(line,*,iostat=stat) poro(k)
23    end if
24  End Do
25  write(*,'(f10.4)') poro
26end program shujv1

xiaochen · 发表于 2022-9-20 15:08:38

伊芸贾发表于 2022-9-20 11:59
[mw_shl_code=fortran,true]program shujv1
implicit none
real,allocatable::poro(:)

感谢大佬，可以了

我改了一天，谢谢谢谢

[文件读写] 文件数据读写

美女勋章

热心勋章

星光勋章

新人勋章

贡献勋章

管理勋章

帅哥勋章

爱心勋章

规矩勋章

元老勋章

水王勋章