關於forall和where在程序裡使用的問題

max533

下面是我的一個程序，但我不清楚為什麼forall和where在這個案例上的應用。為什麼forall可以，但where不行。
分別在line57和line60，如果使用forall的話，在螢幕上就只能顯示4個字的字串。如果使用where的話，在螢幕上就只能顯示3個字的字串。
程式的相關附件在壓縮檔中。
先在這邊感謝各位高手。

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

查看代码

打印?

01	program add_gamit_station_info

02	implicit none

03

04	integer :: count,stat,err,err2,i,stat2,stat3,k,stat4

05	character :: example_line*300,filename*30,line*300,site*4

06	character (len=300),allocatable :: example_line_array(:)

07	character (len=4),allocatable :: example_site_array(:)

08

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

10	close (11,status='delete')

11	call system ('dir /b *.info > list.txt')

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

13	open (22,file='station_cwb.info.db')

14	open (33,file='station_expt.info.db')

15	!---------------------------------------------> to check array size

16	do k=1,5,1

17	read (22,*,iostat=stat4) example_line

18

end do

19

count=0

20

stat=0

21	do while (stat==0)

22	read (22,'(a300)',iostat=stat) example_line

23	if (stat/=0) exit

24	count=count+1

25

end do

26	!---------------------------------------------> to check array size

27	rewind (22)

28	!---------------------------------------------> to allocate the example_line_array & example_site_array

29	allocate (example_line_array(count),stat=err)

30

31	if (err==0) then

32	write (*,*) 'The example_line_array allocated sucessfully.'

33	else if (err/=0) then

34	write (*,*) 'The example_line_array failed to allocate.'

35

end if

36

37	allocate (example_site_array(count),stat=err2)

38

39	if (err2==0) then

40	write (*,*) 'The example_site_array allocated sucessfully.'

41	else if (err/=0) then

42	write (*,*) 'The example_site_array failed to allocate.'

43

end if

44	!---------------------------------------------> to allocate the example_line_array & example_site_array

45

46	!---------------------------------------------> to assign array value from station_cwb.info into example_line_array

47	do k=1,5,1

48	read (22,*,iostat=stat4) example_line

49

end do

50

51	do i=1,count,1

52	read (22,'(a300)') example_line_array(i)

53

end do

54	!---------------------------------------------> to assign array value from file into example_line_array

55

56	!---------------------------------------------> to assign array value from example_line_array(2:5) into example_site_array

57	!where (example_site_array/=' ') example_site_array=example_line_array(2:5)

58

59

60	forall (i=1:count:1) example_site_array(i)=example_line_array(i)(2:5)

61	!---------------------------------------------> to assign array value from example_line_array(2:5) into example_site_array

62	do i=1,count,1

63	write (*,*) example_site_array(i)

64

end do

65	rewind (22)

66	!---------------------------------------------> to check whether is new site.

67	!---------------------------------------------> if it is a new site , we write down the name of site.

68	!---------------------------------------------> if it isn't a new site , we skip site.

69

stat2=0

70	do while (stat2==0)

71	read (11,*,iostat=stat2) filename

72	if (stat2/=0) exit

73	open (44,file=trim(filename))

74

75	do k=1,5,1

76	read (44,*,iostat=stat4) line

77

end do

78

79

stat3=0

80	do while (stat3==0)

81	read (44,'(a300)',iostat=stat3) line

82	if (stat3/=0) exit

83	site=line(2:5)

84	if (any(example_site_array==site)) cycle

85	write (33,'(a4)') site

86

end do

87

end do

88	!---------------------------------------------> to check whether is new site.

89	!---------------------------------------------> if it is a new site , we write down the name of site.

90	!---------------------------------------------> if it isn't a new site , we skip site.

91	close (11)

92	close (22)

93	close (33)

94	close (44)

95	end program

fcode

改为
where (example_site_array/=' ') example_site_array(:)=example_line_array(:)(2:5)
你那样写的话，编译器会认为 (2:5) 是数组片段，而不是子字符串。

max533

fcode 发表于 2015-8-24 21:21
改为
where (example_site_array/=' ') example_site_array(:)=example_line_array(:)(2:5)
你那样写的话， ...

感謝fcode的說明，居然犯了這種錯誤，感謝fcode提醒。