爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8301|回复: 11

[脚本编辑] GRADS绘制台风路径出错,请帮忙指点!

[复制链接]

新浪微博达人勋

发表于 2018-7-13 17:58:21 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

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'
;


密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-13 18:00:56 | 显示全部楼层
添加的图片居然没有显示出来,我把错误提示在这写一下:
ENABLE error: Missing file name
Error in read: NULL file name
ERROR occured on line 33
IN file D:\fnl\20180711-12\path.gs
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-14 09:54:13 | 显示全部楼层
第一句不是说enable 错误吗,你先看看那一句,改完了再看其他的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-14 11:11:02 | 显示全部楼层
检查一下outpng这个变量,path()函数的参数传递进去的,一步一步的倒着去查。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-14 16:49:53 | 显示全部楼层
wangke166 发表于 2018-7-14 11:11
检查一下outpng这个变量,path()函数的参数传递进去的,一步一步的倒着去查。

大神,请问一下 function path(args)函数怎么使用的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-14 16:50:00 | 显示全部楼层
wangke166 发表于 2018-7-14 11:11
检查一下outpng这个变量,path()函数的参数传递进去的,一步一步的倒着去查。

大神,请问一下 function path(args)函数怎么使用的
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-7-14 16:51:16 | 显示全部楼层
river 发表于 2018-7-14 09:54
第一句不是说enable 错误吗,你先看看那一句,改完了再看其他的

enable 出错,我想大概是因为我没有具体的文件名,是用PATH函数引用的,这个函数该怎么用啊,我不太会,这个GS是从论坛下载的别人的,还请多多指教啊
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2018-7-16 15:10:37 | 显示全部楼层
用法你的说明里已经写了:
ga>路径\path.gs 台风文件路径 输出台风图片路径
“台风文件路径”,是一个输入文件,看程序,是文本格式,格式如下:
时间 强度 纬度 经度 未知(第五列未用到)风速
所以你首先需要一个台风路径文件。
“输出台风图片路径”是你电脑上画图的路径和文件名。路径尽量用“/”代替“\”,可以减少很多出错几率。
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-8-17 10:54:40 | 显示全部楼层
wangke166 发表于 2018-7-16 15:10
用法你的说明里已经写了:
ga>路径\path.gs 台风文件路径 输出台风图片路径
“台风文件路径”,是一个输 ...

感谢指导,现在看懂了,我再试试
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2018-8-17 10:54:47 | 显示全部楼层
wangke166 发表于 2018-7-16 15:10
用法你的说明里已经写了:
ga>路径\path.gs 台风文件路径 输出台风图片路径
“台风文件路径”,是一个输 ...

感谢指导,现在看懂了,我再试试
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表