缺测数据

云中月影

在网上下载的160站降水资料，求夏季降水资料，但得出的数据有缺测值，怎么处理缺测值？？？program ex_106
integer,parameter::m=160,n=62 !160站，54年
integer i,j
real x(m,n),y(m,n),z(m,n),ave(m,n)
open(1,file='d:\160js\r16006.txt',status='old',form='formatted') 678月降水.zip (46.12 KB, 下载次数: 0)

2013-4-5 12:42 上传

点击文件名下载附件
open(2,file='d:\160js\r16007.txt',status='old',form='formatted')
open(3,file='d:\160js\r16008.txt',status='old',form='formatted')
do j=1,n
read(1,*)(x(i,j),i=1,m)
enddo
do j=1,n
read(2,*)(y(i,j),i=1,m)
enddo
do j=1,n
read(3,*)(z(i,j),i=1,m)
enddo
do j=11,60
do i=1,m
ave(i,j)=jun(i,j)+jul(i,j)+aug(i,j)
enddo
enddo
print*,ave
end

lqouc

不是很明白你到底想干什么。尤其是“求夏季降水资料”是个什么意思？简单说一下好了。希望有帮助。
你的程序就是求了个平均么？但是估计你应该是会出错的，假如像你注释的只有54年的数据那你的循环肯定会read end of file。最后你的数据做的那个求和最好不要用ave表示，那是平均的意思，求和是sum。
缺省值的问题，不知道你是想做统计还是要画图。做统计的话就用if判断下就好了。如果是作图的话那在ctl里面说明缺省值就好了。

尽头的尽头

如果仅仅是处理缺测值，if语句就可以了，然后如果用grads画图的话，在ctl里面写好缺测就行了

river

楼上两位都说了，我再补充一点，如果只是用grads来画图的话，首先要确定缺测值是不是都是一个值，如果有不同还要在程序中赋值成一个值

逸儿~~

原来grads画图只要ctl里标明了就可以不考虑缺测了，受教~~

中国队加油

还是看回复比较受教