求海陆热力差、、求帮助

appy!!!

气候态=纬向均差+海陆热力差   （高度场的）
'reinit'
'sdfopen f:/heigh/hgt.2003.nc'
'enable print f:/heigh/hgt.gif'
'set lat -90 90'
'set lon 0 360'
t1=1
'set z 6'
while(t1<=73)
'set grads off'
t2=(t1-1)*5+1
t3=t2+4
'define a1=ave(hgt,t='t2',t='t3')'
'd a1-ave(a1,lon=0,lon=360)'
'draw title 't1'pented hgt '
'printim ps'%t1%'.gif  white'
'c'
'disable print'
t1=t1+1
endwhile
'reinit'
这个画不出图、、还有、如果我想用气候态减去纬向均差的时候怎么弄啊？