求不同年份、季节的水汽通量和水汽通量散度叠加的gs，请问哪里有问题么？因为我不太肯

baggio

求不同年份、季节的水汽通量和水汽通量散度叠加的gs，请问哪里有问题么？因为我不太肯定是不是对的！而且不管时间怎么设置，出的图像都是一样的！求助啊！
'reinit'
'sdfopen e:\xl\uwnd.mon.mean.nc'
'sdfopen e:\xl\vwnd.mon.mean.nc'
'sdfopen e:\xl\rhum.mon.mean.nc'
'sdfopen e:\xl\air.mon.mean.nc'
'set parea 1.5 10 1.5 7'
'set lon 69.999 120.001'
'set lat 14.999 50.001'
'set xlint 5'
'set ylint 5'
'set lev 500'
'set grads off'
'set grid off'
'set mpdset cnworld'
'set map 1 3 8'
'define tt=air.4+273.15'
'define es=6.11*exp((17.26*(tt-273.15))/(tt-35.86))'
'define qs=0.622*es/(lev-0.378*es)'
'define qq=rhum.3/100.*qs'
'define xflux=qq*uwnd.1'
'define yflux=qq*vwnd.2'
'define sqdiv=hdivg(xflux,yflux)*1e8'
'define aa=ave(sqdiv,t=33,t=141,12)'
'define ab=ave(sqdiv,t=34,t=142.12)'
'define ac=ave(sqdiv,t=35,t=143,12)'
'define a=(aa+ab+ac)/3'
'set cthick 6'
'd a'
'qu=(uwnd.1*qq/9.8)'
'qv=(vwnd.2*qq/9.8)'
'uvq=mag(qu,qv)*1000'
'set cthick 5'
'define ba=ave(uvq,t=33,t=141,12)'
'define bb=ave(uvq,t=34,t=142,12)'
'define bc=ave(uvq,t=35,t=143,12)'
'define b=(ba+bb+bc)/3'
'd b'
'draw title 1950-1959.autumn.vapourflux.500hPa'
'enable print e:\xl\2.gmf '
'print'
'disable print'
;

baggio

@mofangbao

baggio

@mofangbao

天边的那边

我来学习的

river

你求的不是平均么，难道一组数据的平均值还能有好几个？

baggio

river 发表于 2014-1-2 15:44

                               
登录/注册后可看大图

你求的不是平均么，难道一组数据的平均值还能有好几个？

先是对单月求，然后对季节求一下！不明白river你的意思！

baggio

@做个霸气的木头

river

baggio 发表于 2014-1-2 15:52

                               
登录/注册后可看大图

先是对单月求，然后对季节求一下！不明白river你的意思！

意思就是1 2 3 4 5的平均值是3，就这一个，不会有第二个平均值

baggio

river 发表于 2014-1-2 17:45

                               
登录/注册后可看大图

意思就是1 2 3 4 5的平均值是3，就这一个，不会有第二个平均值

在这gs中，他t=33、34、35是1950年9月、10月、11月，然后t<143(1959年11月)，这样算是50 年代的。以此类推，可以算60年代、70年代...但后面时间t就算是改了，但图还是和50年代图是一样的！不知道我这样说river你还能明白么？

river

baggio 发表于 2014-1-2 21:52

                               
登录/注册后可看大图

在这gs中，他t=33、34、35是1950年9月、10月、11月，然后t

不明白，如果有人明白的话肯定都会来回答的。问题再仔细描述一下，你说的改t是个怎么改法？