grads 用ave 做5天滑动平均出错

baobaoniou2006

我的ctl 文件

dset osez1981-2010.grd
title qv
undef -9999.0
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 8  levels 1000 925 850 700 600 500 400 300
tdef 365  linear 00z01jan1983 1440mn
vars 1
ose 8 -999 Monthly Mean q*u
endvars


我的gs 文件

'reinit'
'open osez.ctl'
'set grads off'
*'set grid off'
* 'set lat -90 90'
* 'set lon 0 357.5'
           
          nt=1
          while(nt<=10)
          nz=1
          while(nz<=2)
         'set t 'nt''
         'set z 'nz''

          say nt
*          'd ose'
         'd ave(ose,t,t+4)'

         'printim  'time%nt%'-level'%nz%'.jpg white'
         'c'

          nz=nz+1
          endwhile
          nt=nt+1
          endwhile            

        'reinit'   
  错误提示为

Syntax Error:  Invalid dimension expression
  Expecting +/-/= operator, found Error from AVE:  1st dimension expression inva
lid
Operation Error:  Error from ave function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = ave(ose,t,t+4)


我觉得自己没有错啊，不知道ave 函数那个地方弄错了，请大家指教，谢谢。

baobaoniou2006

终于明白了，非常感谢lqouc 的回答。

lqouc

怎么就觉得没有错呢，报错已经很清楚了，自己看看就知道啊，人家都告诉你了ave在设置维数的时候要用等号，也就是t=.....

mofangbao

lqouc 发表于 2013-4-19 17:56

                               
登录/注册后可看大图

怎么就觉得没有错呢，报错已经很清楚了，自己看看就知道啊，人家都告诉你了ave在设置维数的时候要用等号，也 ...

要不是你回答了我就准备删帖子了

lqouc

mofangbao 发表于 2013-4-19 19:36

                               
登录/注册后可看大图

要不是你回答了我就准备删帖子了

哈哈，多谢手下留情啦，其实我一直觉得吧，这些错误你我看来可能不应该，也很简单，但是对他们来说真的有点困难，总要允许有一个进步的过程，这样才留的住人。不过确实很多人发帖问完就再也不来了，很可惜啊。

river

你是不是问完得到结果就不来了?那我不回答了,错误提示说清楚了,自己去看ave的用法去吧

蔚蓝的天

应该就是那个问题了

baobaoniou2006

抱歉，可能问的问题太简单了，但是实在是不知道错在哪里了，我按照帖子上面的说法，把gs 文件改为

   'd ave(ose,t=t,t=t+4)'


但是还是报错啊，结果如下

Syntax Error:  Invalid dimension expression
  Dimension value missing or invalid
Error from AVE:  1st dimension expression invalid
Operation Error:  Error from ave function
  Error ocurred at column 1
DISPLAY error:  Invalid expression
  Expression = ave(ose,t=t,t=t+4)



我找了一本grads 教程，把所有ave的用法都试了一遍，
比如把gs 变成'd ave(ose,t=1,t=5)'  
发现如果我不做循环，语句就能通过，但是一旦做 t 的循环过不了。

我想这个可能是牵扯到grads 纬度设置的问题，实在不知道怎么改，麻烦大家给看看

liuruixin

谢谢了，遇到同样的问题

liuruixin

谢谢了，遇到同样的问题，'define v=ave(ave(vwnd,t='t0',t='t0+11'),x=1,x=144)'
是不是前面的t是系统默认的时间维数，而t0是自己定义的，就要加引号