Fortran读写文件报错

ackery · 发表于 2019-3-17 12:47:41

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

01program test
02 
03implicit none
04 
05integer,parameter::nx=144,ny=73,nz=17,nt=12
06 
07real,dimension(nx,ny,nz,nt)::uwnd,vwnd
08 
09integer::ii,jj,kk,tt,irec
10 
11open(21,file='D:\workspace\work2\uwnd.mon.ltm.dat',form='unformatted',access='direct',recl=nx*ny)
12 
13open(22,file='D:\workspace\work2\vwnd.mon.ltm.dat',form='unformatted',access='direct',recl=nx*ny)
14 
15open(23,file='D:\workspace\work2\uvwnd.dat',status='replace',form='unformatted',access='direct',recl=nx*ny)
16 
17irec=1
18 
19do tt=1,nt
20 
21    do kk=1,nz
22 
23        read(21,rec=irec) ((uwnd(ii,jj,kk,tt),ii=1,nx),jj=1,ny)
24 
25        irec=irec+1
26 
27    end do
28 
29end do
30 
31irec=1
32 
33do tt=1,nt
34 
35    do kk=1,nz
36 
37        read(22,rec=irec) ((vwnd(ii,jj,kk,tt),ii=1,nx),jj=1,ny)
38 
39        irec=irec+1
40 
41        print*,irec
42 
43    end do
44 
45end do
46 
47irec=1
48 
49do tt=1,nt
50 
51    do kk=1,nz
52 
53        write(23,rec=irec) ((uwnd(ii,jj,kk,tt),ii=1,nx),jj=1,ny)
54 
55        irec=irec+1
56 
57    end do
58 
59    do kk=1,nz
60 
61        write(23,rec=irec) ((vwnd(ii,jj,kk,tt),ii=1,nx),jj=1,ny)
62 
63        irec=irec+1
64 
65    end do
66 
67end do
68 
69 
70 
71 
72 
73close(21)
74 
75close(22)
76 
77close(23)
78 
79end program test

运行后报错：

青衣巷 · 发表于 2019-3-18 10:09:13

这应该是文件里的数据量不足，使得你读文件的时候，遇到了记录的尾部。
也可能是recl要乘以4，就像这样
recl=nx*ny*4

[文件读写] Fortran读写文件报错

规矩勋章

新人勋章

元老勋章