用besttrack画台风路径

ˇBelieve_

前几天刚学会画台风路径，分享给跟我一样的新手，希望大家共同进步。
台风路径资料可以在附件中的地址下载
gs文件如下：
'reinit'
'open d:\data\model.le.ctl'
'set lat 20 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 map 1 1 2'
'set grid off'
'set cmax 0'
'd z'
i=1
aa='d:\data\Fung-wong.txt'
while(i<=32)
a1=read(aa)
b1=sublin(a1,2)
say b1
la=subwrd(b1,3)/10
lo=subwrd(b1,4)/10
'q w2xy 'lo' 'la''
x1=subwrd(result,3)
y1=subwrd(result,6)
say x1
say y1
if(i>1)
'draw line 'x2' 'y2' 'x1' 'y1''
endif
i=i+1
x2=x1
y2=y1
endwhile
ff=close(aa)
j=1
bb='d:\data\Fung-wong-dian.txt'
while(j<=10)
c1=read(bb)
d1=sublin(c1,2)
say d1
la=subwrd(d1,3)/10
lo=subwrd(d1,4)/10
'q w2xy 'lo' 'la''
x2=subwrd(result,3)
y2=subwrd(result,6)
say x2
say y2
j=j+1
k=subwrd(d1,2)
if(k=1)
'draw mark 2 'x2' 'y2' 0.1'
endif
if(k=2)
'draw mark 3 'x2' 'y2' 0.1'
endif
if(k=3)
'draw mark 3 'x2' 'y2' 0.15'
endif
if(k=4)
'draw mark 10 'x2' 'y2' 0.1'
endif
if(k=5)
'draw mark 10 'x2' 'y2' 0.15'
endif
if(k=0)
'draw mark 6 'x2' 'y2' 0.1'
endif
endwhile
ff=close(bb)
*'draw title Fung-wong'
'printim d:\data\Fung-wong.png white'
;

yysimida

都是高手啊

lusang

好东西，学习了。

river

感谢分享，贴个效果图更直观

Guarurrrr

{:eb302:}好厉害么么哒！

Guarurrrr

{:eb301:}

limengxia

{:eb502:}

IAMQWW1

好东西，学习一下

mm_cat

你好，可以分享一下besttrack的ctl资料么，不知道这个应该怎么写....

meimeideqiqi

楼主你好，我用了你的程序，显示第20行出错，la=subwrd(b1,3)/10，错误原因是non-numeric args to numeric operation，这是为什么呢，求解答