关于do while下iostat的问题

Zhuming

首先，这个问题我已经通过修改代码解决了。但还是没有搞懂为什么是这样，故来问一下。
问题描述：从do while 开始，我循环读取了两个文本，使用status1、status2来查询文件读取状态。接着是在do循环下，使用了两个并列的if语句。程序执行完，结果只有第二个if语句判断后的结果。当我把两个文件中任何一个文件的读取给注释掉，执行结果都正确。
最后我的解决办法是，写两个do while。第一个do while下，读取并判断第一个文件；第二个do while 读取判断第二个文件。结果经过我对比检查是正确的。
so，我还是不知道第一种方法为什么错，请求赐教！谢谢！！

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

查看代码

打印?

01	program flash_density_test

02	implicit none

03

04

!声明变量

05	character(len=512)::cRead,dRead

06	character(len=10)::datee

07	character(len=15)::time

08	real::lat,lon

09	character(len=4)::timingerror

10	integer::station

11	integer::status1=0,status2=0

12	integer::i,j

13	integer,parameter::nx=150,ny=140

14	real,dimension(nx,ny)::md

15	character(len=10)::year,month,day,hour,minn,sec,ti,station1,ex1,ex2,ex3

16	real::lat1,lon1

17

18

do i=1,nx

19

do j=1,ny

20	md(i,j)=0.0

21

end do

22

end do

23

24

!读取文件

25	open(11,file='2006.txt',status='old')

26	OPEN(12,file='output1.txt',status='replace')

27	OPEN(13,file='xinan.grd',form='unformatted',status='replace')

28	OPEN(14,file='20170210_27.txt',status='old')

29

30	do while(.true.)

31	read(11,"(A512)",iostat=status1)cRead

32	if(status1/=0) exit

33	CALL parserRead(cRead)

34	!write(*,*)lat,lon

35	lon=lon-97

36	lat=lat-21

37

38	READ(14,"(A512)",iostat=status2)dRead

39	if(status2/=0) exit

40	CALL parserRead1(dRead)

41	!write(*,*)lat1,lon1

42	lon1=lon1-97

43	lat1=lat1-21

44

45

do i=1,nx

46

do j=1,ny

47	if(lon>=(i-1)*0.1.and.lon<i*0.1.AND.lat>=(j-1)*0.1.AND.lat<j*0.1) then

48	md(i,j)=md(i,j)+1.0

49

end if

50	if(lon1>=(i-1)*0.1.and.lon1<i*0.1.AND.lat1>=(j-1)*0.1.AND.lat1<j*0.1) then

51	md(i,j)=md(i,j)+1.0

52

end if

53

end do

54

end do

55

end do

56

57

close(11)

58	write(12,'(1x,150f7.1)')((md(i,j),i=1,nx),j=1,ny)

59

CLOSE(12)

60	write(13)((md(i,j),i=1,nx),j=1,ny)

61

CLOSE(13)

62

CLOSE(14)

63

64

contains

65

!调用子程序1

66	subroutine parserRead(c)

67	character(len=*)::c

68	integer::k

69

70	do k=1,LEN_TRIM(c)

71	if (c(k:k)=="/") c(k:k)="|"

72

end do

73	READ(c,*)datee,time,lat,lon,timingerror,station

74	do k=1,LEN_TRIM(datee)

75	if (datee(k:k)=="|") datee(k:k)="/"

76

end do

77	end subroutine parserRead

78

!调用子程序2

79	subroutine parserRead1(d)

80	character(len=*)::d

81	integer::h

82

83	do h=1,LEN_TRIM(d)

84	if(d(h:h)=="/") d(h:h)="|"

85

end do

86	READ(d,*)year,month,day,hour,minn,sec,lat1,lon1,ti,station1,ex1,ex2,ex3

87	do h=1,LEN_TRIM(year)

88	if(year(h:h)=="|") year(h:h)="/"

89

end do

90	end subroutine parserRead1

91

92	end program

li913

很大可能是两个文件的行数不同，常见于末尾多了空行。

Zhuming

li913 发表于 2020-6-18 09:36
很大可能是两个文件的行数不同，常见于末尾多了空行。

点醒我！！第二个文件行数比第一行少太多了，差几个数量级。写的时候没考虑到这个问题。感谢！！