求助纠正错误！u，v风场叠加绘距平图！

月是故乡明

打开两个NC文件，将u，v风场叠加绘距平图，这是我的gs文件，不知道对不对，画出的图看起差异不是很大（我绘制的冬季降水特征小年份和特征大值年份，分析风场的影响），希望大家能抽出几分钟帮我看一下，谢谢啦。下面是我选的九个特征小年的距平，特征大年类似，就不贴出来了。不知道这种叠加方式和计算方式是否正确？
'reinit'
'sdfopen d:\pressure\uwnd.mon.mean.nc'
'sdfopen d:\pressure\vwnd.mon.mean.nc'
'set lev 850'
'set grid off'
'set grads off'
'set gxout vector'
'set arrscl 0.13'
'define a=ave(uwnd,t=72,t=74)'
'define b=ave(uwnd,t=132,t=134)'
'define c=ave(uwnd,t=252,t=254)'
'define d=ave(uwnd,t=420,t=422)'
'define e=ave(uwnd,t=444,t=446)'
'define f=ave(uwnd,t=504,t=506)'
'define g=ave(uwnd,t=564,t=566)'
'define h=ave(uwnd,t=600,t=602)'
'define i=ave(uwnd,t=660,t=662)'
'define j=ave(uwnd,t+47,t+779,12)'
'define k=ave(uwnd,t+48,t+780,12)'
'define l=ave(uwnd,t+49,t+781,12)'
'define m=ave(vwnd.2,t=72,t=74)'
'define n=ave(vwnd.2,t=132,t=134)'
'define o=ave(vwnd.2,t=252,t=254)'
'define p=ave(vwnd.2,t=420,t=422)'
'define q=ave(vwnd.2,t=444,t=446)'
'define r=ave(vwnd.2,t=504,t=506)'
'define s=ave(vwnd.2,t=564,t=566)'
'define t=ave(vwnd.2,t=600,t=602)'
'define w=ave(vwnd.2,t=660,t=662)'
'define u=((a+b+c+d+e+f+g+h+i)/9-(j+k+l)/3)'
'define v=((m+n+o+p+q+r+s+t+w)/9-(j+k+l)/3)'
'd skip(u,2);skip(v,2)'
'enable print d:\pic\uvwnd1(850).gmf'
'print'
'disable print'
;
uvwnd特征小年距平图：

uvwnd特征小年距平图：

/猪头哼哼

http://www.esrl.noaa.gov/psd/cgi-bin/data/composites/printpage.pl
楼主可以去这个网站在线合成试试~~~看看结果是不是和你画的一致

kongfeng0824

对。可以吧网上的例子的数据和你画出来的对比一下既可以知道你的程序是否有问题。

river

我觉得这三句有点儿怪'define j=ave(uwnd,t+47,t+779,12)'
'define k=ave(uwnd,t+48,t+780,12)'
'define l=ave(uwnd,t+49,t+781,12)'
这应该每一句都能求出来好几个值，最后面求平均的时候感觉就不对了

月是故乡明

/猪头哼哼 发表于 2013-5-20 00:24

                               
登录/注册后可看大图

http://www.esrl.noaa.gov/psd/cgi-bin/data/composites/printpage.pl
楼主可以去这个网站在线合成试试~~~ ...

嗯，好的，我去看看

月是故乡明

kongfeng0824 发表于 2013-5-20 05:50

                               
登录/注册后可看大图

对。可以吧网上的例子的数据和你画出来的对比一下既可以知道你的程序是否有问题。

嗯好……我去试一下

月是故乡明

river 发表于 2013-5-20 07:26

                               
登录/注册后可看大图

我觉得这三句有点儿怪'define j=ave(uwnd,t+47,t+779,12)'
'define k=ave(uwnd,t+48,t+780,12)'
'define  ...

我这三句主要是求从1951年到2013年冬季的平均降水，后面用(j+k+l)/3)又求了平均，所以应该只有一个值了……我现在再重新调试一下，谢谢你

river

月是故乡明 发表于 2013-5-20 10:25

                               
登录/注册后可看大图

我这三句主要是求从1951年到2013年冬季的平均降水，后面用(j+k+l)/3)又求了平均，所以应该只有一个值了… ...

是的呢，我把加号理所当然的当成等于号了。

流水2012

知道了两个文件同时打开后，第二个文件的识别方法。

limoxi

月是故乡明 发表于 2013-5-20 10:25
我这三句主要是求从1951年到2013年冬季的平均降水，后面用(j+k+l)/3)又求了平均，所以应该只有一个值了… ...

用这种方法是可以分别算出特定3个月的多年平均吗？那t的初始值要怎么设置呢？