提取数据后按年存放出现问题，不知道问题在哪

如若

1. program tiqu
   
2. integer,parameter::m=52,day=22647,yr=62,dd=92*62,yrd=92  
   
3. character*5 id(m)
   
4. character*4 yrc
   
5. integer year(m,day),mon(m,day),dy(m,day),temp,tmax,tmin,rain(m,day)
   
6. integer year1(m,dd),mon1(m,dd),dy1(m,dd),rain1(m,dd)
   
7. real rainyr(1951:2012)
   
8. integer a(1951:2012)
   
9. integer i,yrn
   
10. 
    
11. !!!!!!读入站点名!!!!!!
    
12. open(10,file='f:\fortran\shuju\52sta.txt')
    
13. read(10,*) (id(i),i=1,m)
    
14. do i=1,m
    
15.    j=1
    
16.   !print *,id(i)
    
17. !!!!!!读入51-13年数据!!!!!!
    
18.     open(33,file='f:\fortran\shuju\51-13\'//id(i)//'.txt')
    
19. 100 read(33,*)        id(i),year(i,j),mon(i,j),dy(i,j),temp,tmax,tmin,rain(i,j)
    
20.    !  print *, id(i),year(i,j),mon(i,j),dy(i,j),rain(i,j)
    
21.             j=j+1
    
22.         if(.not.eof(33))        goto 100
    
23.         write(33,*)
    
24.         close(33)
    
25. enddo
    
26. !!!!!!提取3-5月降水!!!!!!
    
27. do i=1,m
    
28.    open(44,file='f:\fortran\shuju\51-13\'//id(i)//'.txt')
    
29.    !open(55,file='f:\fortran\shuju\chunji\'//id(i)//'.txt')
    
30.    do j=1,day
    
31.           if(mon(i,j).ge.6.and.mon(i,j).le.8) then
    
32.                  if(rain(i,j)>=30000) then
    
33.                         rain(i,j)=0.0
    
34.                  end if
    
35.                  if(rain(i,j)==-32744.or.rain(i,j)==32744.or.rain(i,j)==32766.or.rain(i,j)==-32766) then
    
36.                            rain(i,j)=0.0
    
37.                  end if   
    
38.                  ! print *, id(i),year(i,j),mon(i,j),dy(i,j),rain(i,j)
    
39.                 write(44,*) id(i),year(i,j),mon(i,j),dy(i,j),rain(i,j)  !按台站存放!
    
40.           endif
    
41.         enddo
    
42.         close(44)
    
43.         !close(55)
    
44. enddo
    
45. !print *,'ok'
    
46. !!!!!!按年存放!!!!!!
    
47. do i=1,m
    
48.         open(100,file='f:\fortran\shuju\51-13\'//id(i)//'.txt')
    
49.         do j=1,day
    
50.            read(100,*) id(i),year1(i,j),mon1(i,j),dy1(i,j),rain1(i,j)
    
51.            ! print *,id(i),year1(i,j),mon1(i,j),dy1(i,j),rain1(i,j)
    
52.         enddo
    
53. print *,id(i),year1(i,j),mon1(i,j),dy1(i,j),rain1(i,j)
    
54. enddo
    
55. print *,'ok'
    
56. j=1
    
57. do yrn=1951,2012
    
58.         rainyr(yrn)=0.0
    
59.         write(yrc,'(i4)') yrn
    
60.         open(200,file='f:\fortran\shuju\yearly\'//yrc//'.txt')
    
61.         jj=j
    
62.         do i=1,m
    
63.            j=jj
    
64. 400           if(year1(i,j)==yrn)        then
    
65.           !print *,id(i),year1(i,j),mon1(i,j),dy1(i,j),rain1(i,j)
    
66.                   write(200,*) id(i),year1(i,j),mon1(i,j),dy1(i,j),rain1(i,j)
    
67.                   rainyr(yrn)=rainyr(yrn)+rain1(i,j)
    
68.                   j=j+1
    
69.                   if(j<=day)  goto 400
    
70.            endif
    
71.         enddo
    
72. enddo
    
73. open(200,file='f:\fortran\shuju\yr_ave.txt')
    
74. write(200,*) (rainyr(i),i=1951,2012)
    
75. 
    
76. end
    
77. 
    

复制代码

运行出来如下图
求大神告知问题出在那里了！！！！

weinihou

都感觉有点儿同情lz了。。。

貌似几天下来还是再搞一个东西

weinihou

lz跪了，是在17~26行之间么？
lz一点儿一点儿调诶，用debug模式，不要一下就跑了
既然编译连接什么的没出错，那不就可以一句一句慢慢走了嘛，看看变量怎么变的
为了方便，可以把没用的、暂时没用的先全都注释掉

如若

weinihou 发表于 2016-4-22 20:42
都感觉有点儿同情lz了。。。

貌似几天下来还是再搞一个东西

{:sad:}{:sad:}{:sad:}

如若

weinihou 发表于 2016-4-22 20:42
都感觉有点儿同情lz了。。。

貌似几天下来还是再搞一个东西

{:sad:}{:sad:}{:sad:}