关于ave函数求时间平均的问题

爱知天气 · 发表于 2014-4-26 18:04:58

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

ave中的时间设置可不可以是表达式？比如下边的语句。运行的时候说这一行出错，应该是t=表达式这里出错了吧
'd ave(omega,t='(yr-1)*5+1',t='(yr-1)*5+5')'

御风者 · 发表于 2014-4-26 19:16:06

这么写好像确实不行的。

穿云搅浪 · 发表于 2014-4-26 19:18:03

御风者发表于 2014-4-26 19:16
这么写好像确实不行的。

要应该分开写才行？

刷牙 · 发表于 2014-4-26 21:36:04

能不能先令t1=(yr-1)*5+1，t2=(yr-1)*5+5
然后'd ave(omega,t=t1,t=t2)'这样呢

sun_shine_Xia · 发表于 2014-4-26 21:36:31

要写成'（yr-1）'*5+1吧

爱知天气 · 发表于 2014-4-26 23:29:44

刷牙发表于 2014-4-26 21:36
能不能先令t1=(yr-1)*5+1，t2=(yr-1)*5+5
然后'd ave(omega,t=t1,t=t2)'这样呢

这样做可以运行，但是如果前边有循环的话，出来的数据不正确，某一年的资料会覆盖掉所有。就象这样
yr=1
while(yr<=20)
tt=(yr-1)*5+1
ttt=(yr-1)*5+5
'd ave(omega,t='tt',t='ttt')'
出来的数据只有一个值，画出的20个时次的图是一样的

爱知天气 · 发表于 2014-4-26 23:35:09

御风者发表于 2014-4-26 19:16
这么写好像确实不行的。

有没有其他解决办法如果想让它们可以变动的话

刷牙 · 发表于 2014-4-27 10:46:33

爱知天气发表于 2014-4-26 23:29
这样做可以运行，但是如果前边有循环的话，出来的数据不正确，某一年的资料会覆盖掉所有。就象这样
yr=1 ...

能不能把整个gs贴出来

爱知天气 · 发表于 2014-4-27 11:23:54

本帖最后由爱知天气于 2014-4-27 14:05 编辑

刷牙发表于 2014-4-27 10:46
能不能把整个gs贴出来

应该是这个
'reinit'
'open E:\BY\figure2\data6to10\omega1958-1977JtoO.ctl'
'set fwrite e:\BY\figure2\20yrmean\omega1958-1977.grd'
'set gxout fwrite'
yr=1
while(yr<=20)
tt=(yr-1)*5+1
ttt=(yr-1)*5+5
'd ave(omega,t='tt',t='ttt')'
yr=yr+1
endwhile
'disable fwrite'
;
其实不同tt、ttt表示，直接放在里面也是一样的

刷牙 · 发表于 2014-4-27 13:11:51

爱知天气发表于 2014-4-27 11:23
应该是这个
'reinit'
'open E:\BY\figure2\data6to10\omega1958-1977JtoO.ctl'

怎么只给出了ttt的表达式，没给出tt的，我以前做过和这个差不多的，这么写应该没有问题吧，你把tt的表达式加上，再看看

		自动登录	找回密码
密码			立即注册

[图形美化] 关于ave函数求时间平均的问题

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