为什么调用aave总是出错！！

酒国

在做垂直风切变相关的程序，要用到区域平均这个函数，但在调用时总是出错，是怎么回事呢？
下面是程序和grads出错截图。
'reinit'
'open d:\test\20110923\fnl_20110923_00_00.ctl'
'set lon 100 130'
'set lat 5 30'
'set lev 850'
'set t 23'
'set grads off'
'set grid off'
'set gxout contour'
'set csmooth on'
'set cthick 10'
'set ccolor 1'
'set xlopts 1 7 0.17'
'set ylopts 1 7 0.17'
'set map 1 1 6'
'draw map'
'set clopts 0 6 0.17'
'set xlint 2.0'
'set ylint 2.0'
'set lon 100 130'
'set lat 5 30'
'define u1=ugrdprs'
'define v1=vgrdprs'


'open d:\test\20110923\fnl_20110923_00_00.ctl'
'set lon 100 130'
'set lat 5 30'
'set t 23'
'set lev 200'
'set gxout contour'
'set csmooth on'
'set cthick 10'
'set ccolor 1'
'set xlopts 1 8 0.16'
'set ylopts 1 8 0.16'
'set map 1 1 6'
'set xlint 2.0'
'set ylint 2.0'
'set clopts 0 6 0.17'
'define u2=ugrdprs'
'define v2=vgrdprs'

'define a=u2-u1'
'define b=v2-v1'
'define c=mag(a,b)'
'define ac=aave(c,lon=109,lon=120,lat=13,lat=23)'
'd ac'

f117hqh

貌似应该是
'define u2=ugrdprs.2'
'define v2=vgrdprs.2'

酒国

                               
登录/注册后可看大图

谢谢你的回答。。。。现在的错误提示变成了这个，我的维度也不多啊，就2维度变化。

酒国

f117hqh 发表于 2016-2-29 16:23
貌似应该是
'define u2=ugrdprs.2'
'define v2=vgrdprs.2'

谢谢你的回答。。。。现在的错误提示变成了这个，我的维度也不多啊，就2维度变化。

星雨

你看你提取数据的区域
'set lon 100 130'
'set lat 5 30'
和计算的区域
'define ac=aave(c,lon=109,lon=120,lat=13,lat=23)'
配得上吗？
单从这一点，程序就不对，细心点吧。

f117hqh

酒国 发表于 2016-2-29 16:33
谢谢你的回答。。。。现在的错误提示变成了这个，我的维度也不多啊，就2维度变化。

实际上你在计算aave的时候，lat和lon依然是一个范围（100-130E，5-30N）
你需要把它定义到一个点上
在'define ac=aave(c,lon=109,lon=120,lat=13,lat=23)'之前，加上
'set lat 1'
'set lon 1'

以及ls说的范围不对应的问题，我不知道你是不是有意为之，就不多做评论了

酒国

星雨 发表于 2016-2-29 18:47
你看你提取数据的区域
'set lon 100 130'
'set lat 5 30'

求区域平均，求的区域和画图的区域不用一样大吧。

酒国

f117hqh 发表于 2016-2-29 19:36
实际上你在计算aave的时候，lat和lon依然是一个范围（100-130E，5-30N）
你需要把它定义到一个点上
在' ...

谢谢。在之前加了，还是不行。
'set lon 105'
'set lat 18'
'define ac=aave(c,lon=109,lon=119,lat=13,lat=23)'
'd ac'

星雨

'set lon 105'
'set lat 18'
'define ac=aave(c,lon=109,lon=119,lat=13,lat=23)'
'd ac'
你已经把lat定义到18°N了，后面c就只有18°N上的值，怎么求13-23°N内的平均。

星雨

还有，你的fctl文件里对应的FNL资料中确实有23个时次吗？
我猜你这gs文件是拷贝别人的，你把绘图语句去掉，先把计算语句运行通过在添加绘图语句。
你把ctl文件也贴出来，我看着给你编gs文件。