【已解决】求季节平均，出来的结果和实际不符，为什么？

shirly

我的程序：
'set gxout fwrite'
'set fwrite d:\lw\dat\season\sum\sum.p.ave.dat'
'set x 1 127'
'set y 1 79'
n=1
while(n<=714)
'set t 'n''
'define pj=ave(p,t='n',t='n'+2)'
'd pj'
n=n+12
endwhile
'disable fwrite'
这是只求夏天6-8月平均的程序，但是出来有些点的值明显不对，比如第一个点第一年6-8月的数为1、2、3，按理说按这个程序出来第一个点第一年夏天的值为2，但是这个程序出来的是0.为什么呢？

做个霸气的木头

貌似原理是对的，但是有两个问题存在：1、pj这个定义的变量每次循环其实时间只有一个点，只是平均值不同，所以'set t 'n''应当是'set t 1'，2、ave函数的，时间截止'n'+2不会识别成你要的n+2的值，只能识别变量不识别算式；你可以新设一个变量，在while后面加一行j=n+2，然后把t='n'+2换成'j'。
总之就是：
'set gxout fwrite'
'set fwrite d:\lw\dat\season\sum\sum.p.ave.dat'
'set x 1 127'
'set y 1 79'
n=1
while(n<=714)
j=n+2
'set t 1'
'define pj=ave(p,t='n',t='j')'
'd pj'
n=n+12
endwhile
'disable fwrite'
;

shirly

做个霸气的木头 发表于 2013-3-23 21:35

                               
登录/注册后可看大图

貌似原理是对的，但是有两个问题存在：1、pj这个定义的变量每次循环其实时间只有一个点，只是平均值不同，所 ...

当结果出来对的那一刻，真是……一个下午的困惑终于解开了。
感谢感谢感谢！！！！！！！！！！！！！！
还有个问题，我求完平均之后，缺省值变成-9.99e+8（好像是grads默认缺省）了，之前的99999缺省没有了，这正常吗？

做个霸气的木头

shirly 发表于 2013-3-23 22:07

                               
登录/注册后可看大图

当结果出来对的那一刻，真是……一个下午的困惑终于解开了。
感谢感谢感谢！！！！！！！！！！！！！！ ...

正常，grads会自动忽略缺省值进行计算，默认的计算方法很合理，应该没事的