求助，程序编译没问题，但计算结果出错，程序如下

ledio

站点资料求降水量随海拔高度的分布，每300m一层，海拔属于同一层次的降水量求和，最后求平均，程序编译没问题，但计算结果出错，程序筛选分层失败，降水量之和出现在同一层，而选出的降水总次数是正确的。请各位大神看看问题出在哪里，感激不尽~~

代码如下：

1. 
   
2. 
   
3. 
   
4. !1000为站点数，31为文件个数
   
5.     real x(1000,31),y(1000,31),z(1000,31),averain(8)
   
6.         character*13 name
   
7.     integer error
   
8.     integer ip(8)   
   
9. !按namelist读k个文件
   
10.         open(10,file='e:\r24-8-p\namelist.txt',form='formatted')
    
11.         k=1
    
12.         n=0
    
13. 101        read(10,*,end=1001)name
    
14.       write(*,*)name
    
15.       open(9,file='e:\r24-8-p\result_8_98\'//name//'.txt',form='formatted')
    
16. 
    
17. !读数据
    
18.       do while (.true.)
    
19.         read(9,*,iostat=error)a,b,c,d,e
    
20.         if(error /= 0) exit
    
21.         n=n+1
    
22.     x(n,k)=b
    
23.         y(n,k)=c
    
24.         z(n,k)=e
    
25.         end do
    
26.         close(9)
    
27.       k=k+1
    
28.       goto 101
    
29. 1001  close(10)
    
30. 
    
31. 
    
32. !求满足经、纬度范围的站点的平均降水量
    
33.     averain=0
    
34.         ip=0
    
35.         do k=1,31
    
36.       do n=1,1000
    
37.         if(x(n,k).ge.103.and.x(n,k).le.110.and.y(n,k).ge.24.and.y(n,k).le.29)then
    
38.         if(d>=300.and.d<600)then
    
39.     ip(1)=ip(1)+1
    
40.         averain(1)=averain(1)+z(n,k)
    
41.     else if(d>=600.and.d<900)then
    
42.     ip(2)=ip(2)+1
    
43.         averain(2)=averain(2)+z(n,k)
    
44.     else if(d>=900.and.d<1200)then
    
45.     ip(3)=ip(3)+1
    
46.         averain(3)=averain(3)+z(n,k)
    
47.     else if(d>=1200.and.d<1500)then
    
48.     ip(4)=ip(4)+1
    
49.         averain(4)=averain(4)+z(n,k)
    
50.     else if(d>=1500.and.d<1800)then
    
51.     ip(5)=ip(5)+1
    
52.         averain(5)=averain(5)+z(n,k)   
    
53.     else if(d>=1800.and.d<2100)then
    
54.     ip(6)=ip(6)+1
    
55.         averain(6)=averain(6)+z(n,k)
    
56.     else if(d>=2100.and.d<2400)then
    
57.     ip(7)=ip(7)+1
    
58.         averain(7)=averain(7)+z(n,k)
    
59.     else if(d>=2400.and.d<2700)then
    
60.     ip(8)=ip(8)+1
    
61.         averain(8)=averain(8)+z(n,k)
    
62.         end if
    
63.     end if
    
64.         end do
    
65.         end do
    
66.         do i=1,8
    
67. 
    
68.     averain(i)=averain(i)/ip(i)
    
69.     end do
    
70.     write(*,*)ip,averain
    
71.     open(20,file='e:\r24-8-p\ave_8_98.txt',form='formatted')
    
72.         do i=1,8
    
73.     write(20,*) ip(i),averain(i)
    
74.     end do
    
75. 
    
76.         end
    

复制代码

namelist.txt如下
98080108.0001
98080208.0001
98080308.0001
98080408.0001
98080508.0001
98080608.0001
98080708.0001
98080808.0001
98080908.0001
98081008.0001
98081108.0001
98081208.0001
98081308.0001
98081408.0001
98081508.0001
98081608.0001
98081708.0001
98081808.0001
98081908.0001
98082008.0001
98082108.0001
98082208.0001
98082308.0001
98082408.0001
98082508.0001
98082608.0001
98082708.0001
98082808.0001
98083008.0001
98083108.0001

言深深

粗略看了一下，有个疑问，你的d是高程数据吧？难道不应该是经纬度的函数吗？为什么就是一个值···？

ledio

言深深 发表于 2015-5-23 09:23
粗略看了一下，有个疑问，你的d是高程数据吧？难道不应该是经纬度的函数吗？为什么就是一个值···？

对对对，问题就出在这里，d应该存放在数组中参与判断

biubiubiu

{:5_214:}