关于grads中插值报错的求助

zane

我最近几天刚上手grads，我的目的是想用smth9命令，对某流场散度进行9点区域平均的。在上网也搜索到了别人做的关于降雨场的9点区域平均代码。我自己重现的时候，遇到了报错。关于报错，它说我oacres的使用中，输入的第二个表达式必须是站点数据，我找不出毛病。不能能解决，所以在此求助。
得到的grid格点数据都是通过matlab生成，参照方法http://bbs.06climate.com/forum.php?mod=viewthread&tid=20308。
下面是我的ctl，gs文件以及报错内容。希望大家给我指导帮助，谢谢！
下面是ctl文件：
dset D:\interdiv02win.grd
title Phases of divergence in case winter 2002
undef -9.96921e+36
xdef 81 linear 111 0.05
ydef 101 linear 13 0.05
zdef 1  levels 1000
tdef 6  linear 09feb2002 9dy
vars 1
grid 0 99 prepare for interpolation
endvars

dset D:\div02win2.grd
title Phases of divergence in case winter 2002
undef -9.96921e+36
xdef 21 linear 112 0.1
ydef 31 linear 14 0.1
zdef 1  levels 1000
tdef 6  linear 09feb2002 9dy
vars 1
div 0 99 divergence in each phase
endvars

gs文件如下：
'reinit'
'open d:\interdiv02win.ctl'
'open d:\div02win.ctl'
'set grid off'
'set grads off'
'set map 3 1 5'
'set mpdset cnworld'
'set xlopts 3 5 0.18'
'set ylopts 3 5 0.18'
'set parea 0.7 10.4 0.6 7.8'
'set time 09feb2002'
'set lon 111 115'
'set lat 13 18'
'set lev 1000'
'define a1=oacres(grid,div.2)'
'define b1=maskout(a1,g-0.5)'
'define bb1=smth9(b1)'
'd bb1'
;

zane

请管理员给我删帖，我知道错误在哪里了。我编写站点数据错误。