grads计算均方根误差

3120101008

现在在做wrfout与站点资料的均方根误差的计算，参照论坛之前一篇帖子的脚本，计算出的结果是7.225e+07，感觉不对啊，不知道是脚本的问题，还是我数据有问题？我的站点资料是中国气象局七百多个站的值，为了计算方便wrfout经过处理之后，与站点资料插值后的格点数是一致的，麻烦大家帮我看看！

'reinit'
'open E:\grads\grid.ctl'
'open E:\grads\rain\r07.ctl'
'open E:\grads\wrfout71\r07w.ctl'
'define a=oacres(rain,rain.2)'
'define aa=aave(a,x=1,x=121,y=1,y=79)'
'set t 1'
'set z 1'
i=1
j=1
'define cha1=0'
while(i<121)
'set x ' i
while(j<79)
'set y 'j
'define quan=rain.3'
'define cha1=(aa-quan)*(aa-quan)+cha1'
j=j+1
endwhile
i=i+1
endwhile
'define cha2=cha1/(121*79)'
'define cha=sqrt(cha2)'
;

边缘人糊涂乐

顶~我不懂，我来顶下！加油！加油！

00吴广振

grads的功能很强大的

lichangcan

楼主问题解决了吗遇到了同样的问题