绘制任意两点流场剖面提示错误

松子茶

根据家园里的帖子编写的任意两点间的流场剖面，总是提示错误，请问哪位大神知道应该怎么修改呢？
'reinit'
'open d:/82/2/2.ctl'
'open d:/82/2/orog1.ctl'
'set grid on'
'set grads off'
'set zlog on'
lon1=114.105
lon2=114.636
lat1=36.342
lat2=36.563
'define alfa=atan2('lat2-lat1','lon2-lon1')'
'set x 1'
'set y 1'
'set lev 1000 500'
'set t 5'
lon=lon1
'collect 1 free'
'collect 2 free'
'collect 3 free'
while(lon<=lon2)
lat=lat1+(lat2-lat1)*(lon-lon1)/(lon2-lon1)
'collect 1 gr2stn(vvelprs,'lon','lat')'
'collect 2 gr2stn(ugrdprs*cos(alfa)+vgrdprs*sin(alfa),'lon','lat')'
'collect 3 gr2stn(hgtprs,'lon','lat')'
lon=lon+0.01
endwhile
*'set lat 34 39'
'set lon 114.28 114.40'
'set ylevs 1000 925 850 700 600 500'

'set gxout shaded'
'set csmooth on'
'set clab off'
'set cmin 10'
*'set ccols 0 5 0'
'd coll2gr(3,-ugrdprs)'

'set gxout stream'
'set cmin 0'
*'set cint 4'
'set clab on'
'd maskout(coll2gr(2,-ugrdprs),-coll2gr(3,-ugrdprs));maskout(coll2gr(1,-ugrdprs),-coll2gr
(3,-ugrdprs));mag(coll2gr(2,-ugrdprs),coll2gr(1,-ugrdprs))'
'draw title along (36.40,114.28)to(36.46,114.40)'
'printim d:/82/thetase2.png white'
‘print’

river

请仔细阅读原贴，好好看看报错信息

小其其格

'd coll2gr(3,-ugrdprs)'这里有问题，请根据river大神的帖子好好去看看coll2gr函数用法

松子茶

river 发表于 2017-9-16 20:12
请仔细阅读原贴，好好看看报错信息

好的，谢谢您~

松子茶

river 发表于 2017-9-16 20:12
请仔细阅读原贴，好好看看报错信息

我会认真看的

松子茶

小其其格 发表于 2017-9-17 16:31
'd coll2gr(3,-ugrdprs)'这里有问题，请根据river大神的帖子好好去看看coll2gr函数用法

非常感谢