求区域为70-145E之间的位势高度（55N减去35°N）的平均值

...该...

请求各位Grads高手，有个关于求区域为70-145E之间的位势高度（55N减去35°N）的平均值的gs的问题
gs如下，我这样写的不能出结果，不知道哪出问题了
望指正。
'reinit'
'sdfopen f:/data/hgt.mon.mean.nc'
'set lon 70 145'
'set lat 35'
'set lev 500'
'set t 1 779'
'define west1=ave(hgt,lon=70,lon=145)'
'set lon 70 145'
'set lat 55'
'set lev 500'
'set t 1 779'
'define west2=ave(hgt,lon=70,lon=145)'
'define west3'reinit'
'sdfopen f:/data/hgt.mon.mean.nc'
'set lon 70 145'
'set lat 35'
'set lev 500'
'set t 1 779'
'define west1=ave(hgt,lon=70,lon=145)'
'set lon 70 145'
'set lat 55'
'set lev 500'
'set t 1 779'
'define west2=ave(hgt,lon=70,lon=145)'
'define west3'
'd west3=west2-west1'
'd west3'

兰溪之水

'reinit'
'sdfopen f:/data/hgt.mon.mean.nc'
'set lon 70 145'
'set lat 35'
'set lev 500'
'set t 1 779'
'define west1=ave(hgt,lon=70,lon=145)'
'set lon 70 145'
'set lat 55'
'set lev 500'
'set t 1 779'
'define west2=ave(hgt,lon=70,lon=145)'
'define west3=west2(lat=55)-west1(lat=35)'
'd west3'

...该...

兰溪按照你这样改的，没有任何结果出来。是不是时间变量要给确定不变的啊！

qxtlyf

先做一个时次的，看看结果
'reinit'
'sdfopen f:/data/hgt.mon.mean.nc'
'set t 1'
'set lev 500'
'set lat 35'
'define west1=ave(hgt,lon=70,lon=145)'
'set lat 55'
'define west2=ave(hgt,lon=70,lon=145)'
'define west3=west2-west1'
'd west3'

river

貌似没什么问题吧，gs还可以进一步简化，我用hgt.2012.nc资料试了一下：
'reinit'
'sdfopen f:/ncep/hgt.2012.nc'
'set lon 70 145'
'set lat 35'
'set lev 500'
'set t 1 779'
'define west1=ave(hgt,lon=70,lon=145)'
'set lat 55'
'define west2=ave(hgt,lon=70,lon=145)'
'd west2-west1'
'printim f:\test\1\hgt.gif white'
;