求助：怎样用grads画1月、2月和12月的平均

guoli15

       本人下载的是ncep的多年月平均资料(nc文件）。
       其中关于时间的内置ctl文件如下：tdef 12 linear 00z01jan0001  1mo
       我想对资料按照季节进行平均，其中3 4 5、6 7 8、9 10 11都能搞定，可是12 1 2这三个月的平均怎么用grads直接画出来。不想用fortran...
       对各位好心的朋友提前说声感谢。

river

shirly 发表于 2013-3-23 20:49

                               
登录/注册后可看大图

你好~我也遇到一样的问题，向你讨教。
我的程序：
'set gxout fwrite'

'define pj=ave(p,t='n',t='n'+2)'这句有问题，应该是'define pj=ave(p,t='n',t='n+2')'

lqouc

如果就一年的话就用笨办法好了
set t 12
define a=var
set t 1
define b=a+var
set t 2
d (b+var)/3
当然你如果是多年的话就需要用别的办法了，可以循环

guoli15

lqouc 发表于 2013-3-15 10:28

                               
登录/注册后可看大图

如果就一年的话就用笨办法好了
set t 12
define a=var

感谢。我现在去试下。

guoli15

lqouc 发表于 2013-3-15 10:28

                               
登录/注册后可看大图

如果就一年的话就用笨办法好了
set t 12
define a=var

已搞定，感谢您的帮助

shirly

guoli15 发表于 2013-3-15 11:06

                               
登录/注册后可看大图

已搞定，感谢您的帮助

你好~我也遇到一样的问题，向你讨教。
我的程序：
'set gxout fwrite'
'set fwrite d:\lw\dat\season\sum\sum.p.ave.dat'
'set x 1 127'
'set y 1 79'
n=1
while(n<=714)
'set t 'n''
'define pj=ave(p,t='n',t='n'+2)'
'd pj'
n=n+12
endwhile
'disable fwrite'
这是只求夏天6-8月平均的程序，但是出来有些点的值明显不对，比如第一个点第一年6-8月的数为1、2、3，按理说按这个程序出来第一个点第一年夏天的值为2，但是这个程序出来的是0.为什么呢？

七夕路

之后那个循环程序我正好要用 大家都是人才

小小人

嗯嗯，用grads作图还是需要注意些的

Green_Hunter

讨教，我计算1948-2013年每年6-8月三个月的平均值，我的数据室nc数据，但是我出来的数据文件是错的，我的脚本什么地方错了啊！！！请帮我改一下谢谢。
'reinit'
'sdfopen E:\1\hgt.mon.mean.nc'
'set gxout fwrite'
'set fwrite E:\1\xj.nc'
it=1
while(it<=65)
   t1=(it-1)*12+6
   t2=t1+2
   'd ave(hgt,t='t1',t='t2')'
   it=it+1
endwhile
'disable fwrite'
'reinit'

流浪的心

领教