用Grads直接求NC资料减去气候场的距平

duoduoduo

大家好：

      前些日子翻到10多年前飞鸟回答的关于如何直接用grads求多年距平的帖子，突然想到是不是也可以照此法来求减去气候平均场的距平。不知我的想法是否正确，请高手指点。
    飞鸟当时的原文：'set z 1 17'
                                 'set x 1 144'
                                 'set y 1 73'
                                 'set t 1 12'
                                 'hgtave=ave(hgt,t+0,t=672,12)'
                                 'modify hgtave seasonal'
                                 'set t 1 676'
                                 'hgta=hgt-hgtave'
                                 hgta就是所有的距平值
                                 'set t 1 12'
                                 'hgtave=ave(hgt,t+0,t=672,12)'
                                 set 一句等于是 do it=1,12的循环控制，t+0则是grads中t=it+0的意思。最后我加的12是运算的“步长”，就是：对it＝1时，从时间序                                  中选取1,13,25,.....求ave；对it=2时，从时间序列中取2,14,26,.....求ave。
                                 modify一句的含义，手册中有。是把hgtave这个你求的变量作为气候均值。

那么以此推断，假如一个数据中1981年1月为t=13,2010年12月为t=372,则
                            'set t 1 12'
                             'hgtave=ave(hgt,t+12,t=372,12)'就意味着求出了1981-2010年逐月平均值，再用 'hgta=hgt-hgtave'就是求出了减去气候场的距平值。


这样的理解是否正确？
      

river

是的哇，就是这个意思

午后阳光

10多年前？家园不是2011年建的吗？

栗子

午后阳光 发表于 2015-4-22 18:50
10多年前？家园不是2011年建的吗？

他说的是动力论坛

duoduoduo

栗子 发表于 2015-4-22 20:19
他说的是动力论坛

{:lol:}{:lol:}是的是的 哈哈

limoxi

学习了，之前还笨笨的分月份做平均呢

duoduoduo

river 发表于 2015-4-22 18:48
是的哇，就是这个意思

谢谢 有你肯定我就放心用啦 哈哈

duoduoduo

limoxi 发表于 2015-4-25 22:35
学习了，之前还笨笨的分月份做平均呢

我之前也是呢

limoxi

duoduoduo 发表于 2015-4-26 09:31
我之前也是呢

但是这种方法好像用在矢量场中不对啊。
'reinit''
'sdfopen f:\data\u20.nc'
'sdfopen f:\data\v20.nc'

'define u20ave=ave(uwnd.1,t+0,t=360,12)'
'define v20ave=ave(vwnd.2,t+0,t=360,12)'
'set t 1 420'
'define u20=(uwnd.1-u20ave)'
'define v20=(vwnd.2-v20ave)'
'set grads off'
'set grid off'
'set lon 60 200'
'set lat 0 70'
'set lev 200'
'set t 414'
'set gxout vector'
'd skip(u20,2,2);skip(v20,2,2)'
'gxprint d:/pic/wind206.png'
'c'
'set t 415'
'set gxout vector'
'd skip(u20,2,2);skip(v20,2,2)'
'gxprint d:/pic/wind207.png'
'c'
'set t 416'
'set gxout vector'
'd skip(u20,2,2);skip(v20,2,2)'
'gxprint d:/pic/wind208.png'
'c'

;

limoxi

结果画出的矢量风不对，显然这距平值太大了怎么会这样？求解惑！