- 积分
- 144
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-12-24
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
最近在学习画台风路径,使用的资料是FNL每6小时一次的,使用的GS内容如下。出现这样的错误提示,请各位大神帮忙指点一下,该怎么解决,不胜感激啊
**使用方法:在ga>路径\path.gs 台风文件路径 输出台风图片路径
**把初始点的经纬度写改写在lat0,lon0处
function path(args)
'reinit'
infile=subwrd(args,1)
outpng=subwrd(args,2)
say infile
say outpng
'enable print 'outpng
'open d:\fnl\20180711-12\fnl_20180711_00_00.grib2.ctl'
'set parea 1 11 1 8'
'set grads off'
'set lat 10 35'
'set lon 110 140'
'set cmin 10000'
'set xlopts 1 0.1 0.16'
'set ylopts 1 0.1 0.16'
'd hgtprs'
i=1
strong1=1
c1=1
c2=1
c3=1
c4=1
c5=1
c6=1
lat0=25.7
lon0=147.4
'q w2xy 'lon0' 'lat0
x0=subwrd(result,3)
y0=subwrd(result,6)
while(1)
info=read(infile)
code=sublin(info,1)
if(code!=0)
break
endif
info=sublin(info,2)
latty=subwrd(info,3)*0.1
lonty=subwrd(info,4)*0.1
tim=subwrd(info,1)
wind=subwrd(info,6)
strong0=subwrd(info,2)
'q w2xy 'lonty' 'latty
x=subwrd(result,3)
y=subwrd(result,6)
'set string 2 l 2 0'
'set strsiz 0.15 0.15'
'set line 1 1 5'
'draw line 'x0' 'y0' 'x' 'y
x0=x
y0=y
st=substr(tim,5,10)
if(wind>10&wind<=17.1&c1=1)
'draw wxsym 41 'x' 'y' 0.2 4 1'
c1=c1+1
endif
if(wind>17.1&wind<=24.4&c2=1)
'draw wxsym 41 'x' 'y' 0.2 9 1'
c2=c2+1
endif
if(wind>24.4&wind<=32.6&c3=1)
'draw wxsym 41 'x' 'y' 0.2 7 1'
c3=c3+1
endif
if(wind>32.6&wind<=41.4&c4=1)
'draw wxsym 41 'x' 'y' 0.2 8 1'
c4=c4+1
endif
if(wind>41.4&wind<=50.9&c5=1)
'draw wxsym 41 'x' 'y' 0.2 6 1'
c5=c5+1
endif
if(wind>50.9&c6=1)
'draw wxsym 41 'x' 'y' 0.2 2 1'
c6=c6+1
endif
*if(math_fmod(i,4)=0)
if(strong0!=strong1|i=1)
strong1=strong0
if(strong0=1)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--TD'
'draw wxsym 41 'x' 'y' 0.2 4 1'
endif
if(strong0=2)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--TS'
'draw wxsym 41 'x' 'y' 0.2 9 1'
endif
if(strong0=3)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--STS'
'draw wxsym 41 'x' 'y' 0.2 7 1'
endif
if(strong0=4)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--TY'
'draw wxsym 41 'x' 'y' 0.2 8 1'
endif
if(strong0=5)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--STY'
'draw wxsym 41 'x' 'y' 0.2 6 1'
endif
if(strong0=6)
'set string 1 l 1 1'
'draw string 'x' 'y' 'st'--SuperTY'
'draw wxsym 41 'x' 'y' 0.2 2 1'
endif
endif
i=i+1
endwhile
'print'
'disable print'
;
|
|