求多年月平均值【已解决】

刷牙

要求34年1-12月的月平均值
不知道哪里有问题
奇怪的的是1-5月的结果都有问题
而6-11月结果正常，12月的又有错
请大家帮忙看看问题出在哪里了
谢谢了

ctl描述如下：
dset  D:\study\era\data\hgt\month\hgtmonth.dat   
title 4x daily NMC reanalysis            
UNDEF  -99.0                                 

XDEF 43 LINEAR  72 1.5           
YDEF 28 LINEAR  15 1.5         
ZDEF 1 linear 0 1
TDEF 408 LINEAR 00Z01jan1979 1mo  
vars 1                                                     
hgt     0   99   tropopause height
ENDVARS

gs文件如下：
'reinit'
'open D:\study\era\data\hgt\month\month.ctl'
'set fwrite D:\study\era\data\hgt\month\hgtlongtermmonthmean.dat'
'set gxout fwrite'
'set lat 15 55.5'
'set lon 72 135'

*求1-12月的平均值
*1月34年平均值
'define month=ave(hgt,t=1,t=408,12)'
'd month'
*2月34年平均值
'define month=ave(hgt,t=2,t=408,12)'
'd month'
*3月34年平均值
'define month=ave(hgt,t=3,t=408,12)'
'd month'
'define month=ave(hgt,t=4,t=408,12)'
'd month'
'define month=ave(hgt,t=5,t=408,12)'
'd month'
'define month=ave(hgt,t=6,t=408,12)'
'd month'
'define month=ave(hgt,t=7,t=408,12)'
'd month'
'define month=ave(hgt,t=8,t=408,12)'
'd month'
'define month=ave(hgt,t=9,t=408,12)'
'd month'
'define month=ave(hgt,t=10,t=408,12)'
'd month'
'define month=ave(hgt,t=11,t=408,12)'
'd month'
'define month=ave(hgt,t=12,t=408,12)'
'd month'
'disable fwrite'
;

好吧，自己找到原因了
原来是计算以后我之前定义的缺测值被grads改成了默认的缺测值-9.99e8
ctl文件中改一下就好了{:soso_e129:}

刷牙

逛论坛那是相当的长知识

老男孩

拿回家消化一下。

海阔天又空

很有帮助谢谢楼主

乱目

很有帮助谢谢楼主分享

halida0101

谢谢了楼主