求助！用grads画平均降水量时遇到了问题

林国英

'reinit'
'enable print d:\LJ\1507jpp.gmf'
'open d:\LJ\grid.ctl'站点描述文件
'open d:\LJ\r1607.ctl'降水量描述文件
'set grads off'
'set grid off'
'set lon 73 135'
'set lat 17  55'
'set lev 1000'
'set t 65'
'set mpdset cnworld'
'define rr=ave(r.2,t=31,t=60,1)'
'a1=oacres(g,rr,10,7,4,2,1)'
'aa1=maskout(a1,g-0.5)'
'aaa1=smth9(aa1)'
'a=oacres(g,r.2,10,7,4,2,1)'
'aa=maskout(a,g-0.5)'
'aaa=smth9(aa)'
'define jp=aaa-aaa1'
'define jpp=jp/aaa'
'set xlopts 1 10 0.18'
'set gxout shaded'
'd jpp'
'cbarn'
'draw title jp percent for July 2015'
'print'
'disable print'
;

HSY

我发现我好low啊，这种问题我都是先fortran算好，再画图

lcqxzsl

好好学习 天天向上

river

ave函数不能直接做站点的平均，应该先插值成格点的然后再求平均

surelee

刚好遇到了同样的问题学习了

林国英

river 发表于 2016-4-28 08:01
ave函数不能直接做站点的平均，应该先插值成格点的然后再求平均

已解决，谢谢~

林国英

surelee 发表于 2016-4-28 09:19
刚好遇到了同样的问题学习了

'reinit'
'enable print d:\LJ\1507jp.gmf'
'open d:\LJ\grid.ctl'
'open d:\LJ\r1607.ctl'
'set grads off'
'set grid off'
'set lon 73 135'
'set lat 17  55'
'set lev 1000'
'set mpdset cnworld'
'define b=0'
i=31
while(i<=60)
'set t 'i
'a=oacres(g,r.2,10,7,4,2,1)'
'aa=maskout(a,g-0.5)'
'aaa=smth9(aa)'
'define b=b+aaa'
i=i+1
endwhile
'define pj=b/30'
'set t 65'
'a1=oacres(g,r.2,10,7,4,2,1)'
'aa1=maskout(a,g-0.5)'
'aaa1=smth9(aa)'
'define jp=aaa1-pj'

'set xlopts 1 10 0.18'
'set gxout shaded'
'd jp'
'cbarn'
'draw title jp for July 2015'
'print'
'disable print'
;

林国英

HSY 发表于 2016-4-28 12:22
我发现我好low啊，这种问题我都是先fortran算好，再画图

不，是我low，我用不惯Fortran

洪hong

以前都是一直用FORTRAN计算中国160站的物理量才用GRADS画图，今天长见识了。