- 积分
- 8049
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-5-29
- 最后登录
- 1970-1-1
|

楼主 |
发表于 2013-12-5 09:44:27
|
显示全部楼层
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'
|
|