爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
123
返回列表 发新帖
楼主: shirly

[分享资料] 【已解决】grads循环中,每隔24个数相加之后求平均出错。

[复制链接]
发表于 2017-12-28 18:05:51 | 显示全部楼层
楼主您好,我的gs也总是报这样的错误能不能帮我看一下啊,没有搞清楚哪里错的。我这是个画低涡路径的gs。

function path(args)
'reinit'
infile = subwrd(args,1)
outpng = subwrd(args,2)
infile='E:\zydwnj\199607\lujing\967\199607.txt'

'D:\OpenGrADS\Classic\scripts\define_colors.gs'
'open E:\zydwnj\199607\lujing\967\topo10.ctl'
'set grads off'
'set grid off'
'set xlpos -10'
'set ylpos -10'
'set map 1 1 2'
'set mpdset cnworld'
'set lat 35 70'
'set lon 50 100'
'set gxout shaded'

'set rgb 16 250 250 250'
'set rgb 17 235 235 235'
'set rgb 18 220 220 220'
'set rgb 19 205 205 205'
'set rgb 20 190 190 190'
'set rgb 21 175 175 175'
'set rgb 22 160 160 160'
'set rgb 23 145 145 145'
'set rgb 24 130 130 130'
'set rgb 25 115 115 115'
'set rgb 26 100 100 100'
'set rgb 27  85  85  85'
'set rgb 28  70  70  70'
'set rgb 29  55  55  55'
'set rgb 30  40  40  40'
'set rgb 31  25  25  25'
'set rgb 32  10  10  10'
'set rgb 33   5   5   5'
'set rbcols 16 17 18 19 20 21 22 23 24 25'

'd tp'
'cbarn 1 1 9.94 4.20'

'run D:\OpenGrADS\Classic\scripts\axis.gs -type L -lable on -position o -start 40 -interval 10 -suffix `3.  -lsize 0.25  -lthick 3 -hoffset -0.15'
'run D:\OpenGrADS\Classic\scripts\axis.gs -type B -lable on -position o -start 50  -interval 10 -suffix `3.  -lsize 0.25 -lthick 3 -voffset -0.1'

i=1
strong1=1
c1=1
c2=1
c3=1
c4=1
c5=1
c6=1
lat0=52.5
lon0=72.5
'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.17 0.17'
'set line 1 1 5'
'arrow3 'x' 'y' 'x0' 'y0' 1 1 6'


x0=x
y0=y
st=substr(tim,5,10)

'set font 2'
'set string 1 l 4 0'
'draw string 'x+0.2' 'y' 'st''

'set font 2'
'set string 1 c 4 0'
'set strsiz 0.2 0.2'
'draw string  'x' 'y' D'

i=i+1

endwhile



'print E:\zydwnj\199607\lujing\967\960707.eps'

;
密码修改失败请联系微信:mofangbao
发表于 2017-12-29 09:45:41 | 显示全部楼层
OomaulenOo 发表于 2017-12-28 18:05
楼主您好,我的gs也总是报这样的错误能不能帮我看一下啊,没有搞清楚哪里错的。我这是个画低涡路径的gs。
...

说实话,看层主的gs文件有点头大。建议先对一个站点/数据进行处理,一个运行无误后,再批量处理
密码修改失败请联系微信:mofangbao
发表于 2017-12-29 09:45:38 | 显示全部楼层
OomaulenOo 发表于 2017-12-28 18:05
楼主您好,我的gs也总是报这样的错误能不能帮我看一下啊,没有搞清楚哪里错的。我这是个画低涡路径的gs。
...

说实话,看层主的gs文件有点头大。建议先对一个站点/数据进行处理,一个运行无误后,再批量处理
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

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

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

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