grads画台风路径出现问题

帽帽

循环之后只有第一张有路径，从第二张开始只有标题，没有路径显示
[url=]台风资料用的是中国台风网(www.typhoon.gov.cn)  [/url]“CMA-STI 热带气旋最佳路径数据集”
我也是之前找到别人分享的gs，修改之后写的
'reinit'
'enable print e:\1957.gmf'
'open e:\model.le.ctl'
'set lat 5 35'
'set lon 110 140'
'set lev 500'
'set xlopts 1 4 0.15'
'set ylopts 1 4 0.15'
'set mpdset cnworld'
'set grads off'
'set grid off'
'set cmax 0'
'd z'
line=1
while(line<=952)
parap=read('e:/CH1957BST.txt')
parap=sublin(parap,2)
lno=subwrd(parap,3)
name=subwrd(parap,8)
tno=subwrd(parap,4)
say lno
say name
lp=1
while(lp<=lno)
parap=read('e:/CH1957BST.txt')
parap=sublin(parap,2)
*逐行读取
*get day time latitude lontitude*
*
ymd1=subwrd(parap,1)
pres=subwrd(parap,5)
*获取时间信息
year1=substr(ymd1,1,4)
day1=substr(ymd1,7,2)
*获取日期
time1=substr(ymd1,9,2)
*获取时次
say day1
say time1
lat1=subwrd(parap,3)
*分解得到经纬度值
lon1=subwrd(parap,4)
latv=lat1/10
lonv=lon1/10
say lonv
say latv
'q w2xy 'lonv' 'latv''
*转成直角坐标
lonx=subwrd(result,3)
laty=subwrd(result,6)
say laty
say lonx
'draw wxsym 40 'lonx' 'laty' 0.1 2 5'
*在当前点位置画台风符号
if(time1=00)
*如果是00时标出日期
'set strsiz 0.08 0.08'
'set string 1 bc'
*'draw string 'lonx' 'laty-0.1' 'day1'.08h'
*'draw string 'lonx' 'laty+0.1' 'pres''
endif
if(lp>=2)
*如果不是第一个点，当前点与前一点连线
'draw line 'lonxb' 'latyb' 'lonx' 'laty
'draw title 'year' No. 'tno' typhoon 'name''
Endif
lp=lp+1
latyb=laty
lonxb=lonx
endwhile
'print'
'c'
line=line+1
endwhile
'disable print'
;