[已解决]程序设计求助

muggle

我最近在编写一个用于筛选逐日最高温超过阈值的程序，遇到了一个我怎么都找不到错误的问题。

我的目的是打算，按照年份，把每天的192站点挑选出来，超过阈值的站点显示最高温，没超过的用0代替，按道理这样每天都是输出一个192行的txt。
但奇怪的是，有些天是对的，有些就莫名其妙的跑错位置了。比如5月26日第一个站跑到5月25日最后一行，6月12日一半跑到6月11日文件中。
我检查过，读取原始资料没有问题。。一下子还真想不到哪里出了错，有没有高人指点一下

除了读取资料部分，值得研究的就是那个判断语句，原始资料中是按照每年每站每天排列，我用a1,a2,a3分别代表年月日来判断。

程序是很简单，可能有哪些小地方我没注意到，导致了错误。急求解答！！



解决方法：71和72行掉个顺序就ok了

1. program zhandian
   
2. implicit none
   
3. 
   
4. 
   
5. integer,parameter::sta_num=192,ns=35,day=153
   
6. real sta(day,sta_num,ns),yr(day,sta_num,ns),mon(day,sta_num,ns),da(day,sta_num,ns),tmax(day,sta_num,ns)
   
7. real yz(sta_num,ns),a1(ns),a2(5),a3(31),a4(30)
   
8. integer i,j,k,ii,jj,kk
   
9. character*4 year(ns)
   
10. character*3 dy(153)
    
11. 
    
12. open(21,file='F:\192\dat\yr.txt',status='old')
    
13. read(21,*) year
    
14. close(21)
    
15. 
    
16. open(22,file='F:\192\dat\dy.txt',status='old')
    
17. read(22,*) dy
    
18. close(22)
    
19. 
    
20. 
    
21. 
    
22. do i=1,ns
    
23.   a1(i)=1978+i
    
24. enddo
    
25. 
    
26. do i=1,5
    
27.   a2(i)=4+i
    
28. enddo
    
29. 
    
30. do i=1,31
    
31.   a3(i)=i
    
32. enddo
    
33. 
    
34. 
    
35. 
    
36. !读取观测数据
    
37. 
    
38.   open(11,file='F:\192\sta_fenbu\35192153obs.txt',status='old')
    
39.     do k=1,ns
    
40.           do j=1,sta_num
    
41.         do i=1,day
    
42.               read(11,*) sta(i,j,k),yr(i,j,k),mon(i,j,k),da(i,j,k),tmax(i,j,k)
    
43.         enddo
    
44.           enddo
    
45.     enddo
    
46. 
    
47. 
    
48. !读取阈值文件
    
49. 
    
50. 
    
51. open(22,file='F:\192\yuzhi\r95t35.txt',status='old')
    
52. 
    
53. read(22,"(192F8.2)") ((yz(i,j),i=1,sta_num),j=1,ns)
    
54. 
    
55. !检验读取正确与否
    
56. !write(*,*) yz(1,1),yz(1,35)
    
57. !write(*,*) sta(25,1,1),yr(25,1,1),mon(25,1,1),da(25,1,1),tmax(25,1,1)
    
58. !write(*,*) a1(1),a2(1),a3(25)
    
59. !write(*,*) sta(2,1,35),yr(2,1,35),mon(2,1,35),da(2,1,35),tmax(2,1,35)
    
60. 
    
61. 
    
62. 
    
63. !挑选出大于阈值的站点，并按照年份输出
    
64. do k=1,ns
    
65.     do i=1,day
    
66.        do j=1,sta_num
    
67.            do ii=1,ns
    
68.                  do jj=1,5   
    
69.                         do kk=1,31
    
70.                            if (yr(i,j,k)==a1(ii) .and. (mon(i,j,k)==a2(jj)) .and. (da(i,j,k)==a3(kk))) then
    
71.                        if (tmax(i,j,k) <= yz(i,k)) then
    
72.                           open(33,file='f:\192\sta_fenbu\'//year(k)//'\'//dy(i)//'.txt')
    
73.                           write(33,*) sta(i,j,k),yr(i,j,k),mon(i,j,k),da(i,j,k),0.0
    
74.                        else
    
75.                                   write(33,*) sta(i,j,k),yr(i,j,k),mon(i,j,k),da(i,j,k),tmax(i,j,k)
    
76.                        endif
    
77.                    endif
    
78.                         enddo
    
79.                  enddo
    
80.            enddo
    
81.             enddo
    
82.         enddo
    
83. enddo
    
84. 
    
85. 
    
86. close(11)
    
87. close(22)
    
88. close(33)
    
89. end program
    
90.   
    
91. 
    

复制代码

muggle

因为错误是每年都存在的，所以应该是程序的问题？大神帮忙解答的时候只要运行1979年就可以了，省去了很多循环，谢谢

muggle

虽然不知道为什么，我把71和72行调换了个顺序就OK了~！！！{:lxm_24:}