计算出月水汽通量散度之后如何求多年平均？？（50积分求助）

龙行天下

我用NCEP资料计算出了月水汽通量散度，还想计算多年平均值，用ave函数试了好多次都不成功，请问该怎么做？？
时间i=169代表1962年1月，769代表2012年12月。
'sdfopen E:\test\uwnd.mon.mean.nc'
'sdfopen E:\test\vwnd.mon.mean.nc'
'sdfopen E:\test\rhum.mon.mean.nc'
'sdfopen E:\test\air.mon.mean.nc'
'set lon 70 140'
'set lat 15 55'
'set lev 600'
'set mpdset cnworld'
i=169

while(i<=769)
'set t 'i
'set grads off'
'set grid off'
'set gxout shaded'
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(air.4-273.15))/(air.4-29.65)))'
'define qs=(0.62197*es/(prs-0.378*es))'
'define q=qs*rhum.3/100'
'd hdivg(q*uwnd.1/g,q*vwnd.2/g)'
'run cbarn'

'printim E:\test\rhdiv'i'-600hPa.png white'
'c'

i=i+1
endwhile
;

river

用系统自带的ave函数

mofangbao

不成功当然是你没做对啊  不能说你没做对就否定这个方法的  你应该继续尝试ave  不要失败就放弃

龙行天下

river 发表于 2013-12-5 07:47

                               
登录/注册后可看大图

用系统自带的ave函数

river你好，年平均我已经算出来了。但是我想对两个时间段的平均值作差（下面程序的红色字体部分），出来的图所有值都为0，不知道是什么原因，你能帮我看看吗？
'sdfopen E:\test\uwnd.mon.mean.nc'
'sdfopen E:\test\vwnd.mon.mean.nc'
'sdfopen E:\test\rhum.mon.mean.nc'
'sdfopen E:\test\air.mon.mean.nc'
'set lon 65 110'
'set lat 20 45'
'set lev 600'
'set mpdset cnworld'
i=169

while(i<=769)
'set t 'i
'set grads off'
'set grid off'
'set gxout shaded'
'define prs=lev'
'define g=9.8'
'define es=(6.112*exp((17.67*(air.4-273.15))/(air.4-29.65)))'
'define qs=(0.62197*es/(prs-0.378*es))'
'define q=qs*rhum.3/100'
'define sandu=hdivg(q*uwnd.1/g,q*vwnd.2/g)'
'run cbarn'

'c'

i=i+1
endwhile
;

'd ave(sandu,t=169,t=528,12)-ave(sandu,t=529,t=768,12)'


'printim E:\test\rhdiv-600hPa-1992年之前.png white'

river

龙行天下 发表于 2013-12-5 09:44

                               
登录/注册后可看大图

river你好，年平均我已经算出来了。但是我想对两个时间段的平均值作差（下面程序的红色字体部分），出来的 ...

你这俩句在分号后面，根本就没执行。你从哪看出来都等于0的？难道你不明白gs最后加的那个分号的意义？

龙行天下

river 发表于 2013-12-5 12:09

                               
登录/注册后可看大图

你这俩句在分号后面，根本就没执行。你从哪看出来都等于0的？难道你不明白gs最后加的那个分号的意义？

river，我已经把分号放到最后了，但是计算出的结果还是不对的，你能不能再帮我看看啊，不胜感激

river

龙行天下 发表于 2013-12-5 17:15

                               
登录/注册后可看大图

river，我已经把分号放到最后了，但是计算出的结果还是不对的，你能不能再帮我看看啊，不胜感激

首先我有点疑问，你的ave函数里为什么要加一个时间增量呢，你是单独只计算某一个月的差值吗？还有你贴的图根本就看不出来它是0啊。你需要的是把错误提示的截图贴上来，而不是一张什么都看不出来的图。
在d前面加一句set t 1 试试。

飞飞

路过学习一下

SeoYuria

路过也支持一下哦——

xrpamhf

769是2012年1月吧