【已解决】grads循环中，每隔24个数相加之后求平均出错。

shirly

yuan/// 发表于 2013-9-2 08:47

                               
登录/注册后可看大图

亲，为什么不直接在FORTRAN里面直接算出来需要的东西啊，这样多方便啊。

grads也可以的啊，理解得深刻了比fortran简单多了~！

冬日情愫+

shirly 发表于 2013-9-2 08:52

                               
登录/注册后可看大图

我们当时上课用的“Grads使用手册”（LASG编的），我查的关于ave的函数只有下面的类型：ave（var,t=1,t=2 ...

你这里面：ave(mag(u,v),t-1,t+1)用的就是相对时间，我之前用地也是相对时间。我好想看的是论坛上的介绍，忘了是不是大气所的教材。

HUYAN

谢谢啦 也在纠结中 终于搞定

shirly

HUYAN 发表于 2014-2-27 14:31

                               
登录/注册后可看大图

谢谢啦 也在纠结中 终于搞定

共同学习啦

jiefeng16

{:soso_e179:}{:soso_e179:}

liuyichao629

thanks a lot ~~~~~~~

香草拿铁

借鉴了一下
感谢分享
用了一个t+0,t+3的while，endwhile循环

wavelet2008

{:5_196:}

AnneChen

学习了，感谢分享讨论。

OomaulenOo

香草拿铁 发表于 2015-10-20 09:50
借鉴了一下
感谢分享
用了一个t+0,t+3的while，endwhile循环

嗨，你好，我的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'

;