怎么做TRMM数据提取每年３-５月份的降水平均

┏_(゜ω゜)=

　　　首先是批处理　 'reinit' 'set gxout fwrite'
'set fwrite G:\jiangshui\jiangshuiall.dat'
i=1
while(i<=192)
'sdfopen G:\trmm\3a12.'i'.nc'  
'set  x 540 965'    * 3a12数据是有179.75度开始，每隔0.5度，所以这样设置，做的是太平洋的范围’
'set  y 1 160'
'set z 1'
'd sfcp'
'close 1'
i=i+1

endwhile

'disable fwrite'
'reinit'
可以得到数据，在编写一个描述他的ctl
dset G:\jiangshui\jiangshuiall.dat
title  surface precipitation
undef -9999.9
xdef   426 linear 89.75 0.5
ydef   160 linear  -39.75   0.5
zdef   1 levels  0      
tdef  192 linear 01jan98 1mo
vars 1
sfcp 1 t,y,x  surface precipitation
endvars
'reinit'
最后画图
'reinit'
'open G:\jiangshui\all.ctl'
'enable print G:\jiangshui\222.gmf'

'set lon 89.75 302.25'
'set lat -39.75 39.75'
'set z 1'

'set gxout shaded'
'define a=ave(sfcp,t=6,t=186,12)+ave(sfcp,t=7,t=187,12)+ave(sfcp,t=8,t=188,12)'  *三个月的再平均
'define b=a/3'
'd b'
'cbarn.gs'
'print'
'disable print'

'reinit'

可是怎么也画不出来

lqouc

没有print自然没有图

┏_(゜ω゜)=

根据上面的程序，画出来的图是这样的 ，感觉十分的奇怪，求指教啊。。。。

┏_(゜ω゜)=

lqouc 发表于 2014-4-24 15:56
没有print自然没有图

你好，我想问下我为什么画的的很奇怪呢？

┏_(゜ω゜)=

已解决