關於循序讀檔的問題

max533

如題：在第2個迴圈內，我因為同一個檔案想一次讀2行再讓他返回前一行，並進行計算，檔案還沒讀完但卻在讀完2行之後就直接跳出，我該怎麼辦呢?!
stat2和stat3都在讀完一次之後就不再是0，是我哪裡有設計錯誤嗎?!

感謝各位高手幫忙解答！！！
以下為程序代碼：

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

查看代码

打印?

01	program radiosonde_PWV_B

02	implicit none

03

04	integer :: stat,stat2,stat3,height_m1,height_m2,yr1,doy1,hr1,yr2,doy2,hr2

05	character :: filename*13,filename_new*13

06	real*8 :: timeall1,part_vapor_weight1,timeall2,part_vapor_weight2,pwv

07

08	open (11,file='list.txt')

09	close (11,status='delete')

10	call system ('dir /b *.radio2 > list.txt')

11	open (11,file='list.txt')

12

13

stat=0

14	do while (stat==0)

15	read (11,*,iostat=stat) filename

16	if (stat/=0) exit

17	open (22,file=filename)

18	filename_new=filename(1:6)//'.radio3'

19	open (33,file=filename_new)

20

21

stat2=0

22

stat3=0

23	do while (stat2==0 .and. stat3==0)

24	read (22,*,iostat=stat2) timeall1,height_m1,part_vapor_weight1,yr1,doy1,hr1

25	read (22,*,iostat=stat3) timeall2,height_m2,part_vapor_weight2,yr2,doy2,hr2

26	if (stat2/=0 .or. stat3/=0) exit

27	backspace (22)

28	if ( hr1==hr2 ) then

29	pwv=(part_vapor_weight1+part_vapor_weight2)*(height_m2-height_m1)*1D3/2

30	write (33,50) timeall1,height_m1,height_m2,pwv,yr1,doy1,hr1

31	50    format (f12.7,1x,i5.5,1x,i5.5,1x,f12.7,1x,i4,1x,i3.3,1x,i2.2)

32	else if ( hr1/=hr2 ) then

33

cycle

34

end if

35

end do

36

end do

37	end program

以下為測試數據：

fcode

yr1,doy1,hr1,yr2,doy2,hr2 这些变量应为 real 类型

max533

fcode 发表于 2015-2-10 10:35
yr1,doy1,hr1,yr2,doy2,hr2 这些变量应为 real 类型

為什麼呢?!這些變數都在原始資料中都是整數耶！！！

fcode

1994.0009406    00011             1022.1 14.0    76   0.0000104187 1994 001 08
timeall1    height_m1,part_vapor_weight1, yr1, doy1,           hr1

1994.0009406 timeall1
00011         height_m1
1022.1        part_vapor_weight1
14.0         yr1   这是real
76           doy1
0.0000104187 hr1  这也是

max533

fcode 发表于 2015-2-10 13:15
1994.0009406    00011             1022.1 14.0    76   0.0000104187 1994 001 08
timeall1    height_m1 ...

感謝大大提醒，的確因為沒有讀好檔案，造成檔案判斷讓iostat回傳的值不等於0。