有关用指数时间序列画直方图的问题

nuistzhou

是用来求一个指数40年的时间序列图的。指数值为两个地区风场的距平值相减。

'reinit'
'sdfopen e:\uwnd.mon.mean.nc '
i=282
'set lev 850'
while(i<=750)
'set t 'i''
'summer=ave(uwnd,t='i',t='i'+2)'    提取夏季资料
i=i+12
endwhile
'a1=aave(summer,x=41,x=61,y=41,y=45)'
'a2=aave(summer,x=41,x=61,y=47,y=51)'    做区域平均均
'aa1=ave(a1,t=1,t=40)'           求40年平均值
'aa2=ave(a2,t=1,t=40)'
'set t 1 40'
'juping1=a1-aa1'     分别求2个区域的距平值
'set t 1 40'
'juping2=a2-aa2'
'set gxout bar'     画直方图
'set t 1 40'
'd  juping2-juping1'    距平差，也就是指数
'printim e:\123.png'
;

运行后不出图，一片空白。
求释疑

做个霸气的木头

     你文件求夏季平均的时候就出问题了，你这样不会得到多个年份夏季平均的summer，每次的summer会被替换掉的。我出问题时的解决办法是用循环把每年夏季平均输出出来，写了个ctl，然后这个文件里是每年夏季平均的值，接下来就好做了~
     不知道其他高手有没有高明的办法分享个~

nuistzhou

做个霸气的木头 发表于 2012-5-18 12:37

                               
登录/注册后可看大图

你文件求夏季平均的时候就出问题了，你这样不会得到多个年份夏季平均的summer，每次的summer会被替换掉 ...

确实有这个问题。。。早上没辙，只好提出来再用f90处理，再excel 的。。。好低端

做个霸气的木头

nuistzhou 发表于 2012-5-18 13:01

                               
登录/注册后可看大图

确实有这个问题。。。早上没辙，只好提出来再用f90处理，再excel 的。。。好低端

我们老师讲过，凡是能解决问题的方法都可以。数据提出来后，直接写个对应ctl，就可以接着用了，不用fortan，后面好做很多了！

易小凯

本身资料是每月一个值，而你现在要把时间间隔变成年，不太方便实现。

mofangbao

你一步步的分别出图看看，看看是到哪里变成空白的

nuistzhou

mofangbao 发表于 2012-5-18 16:23

                               
登录/注册后可看大图

你一步步的分别出图看看，看看是到哪里变成空白的

如木头所说，就是只得到了一年的夏季平均。。。
从那之后就有问题了。
我有个疑问，后边的那些变量计算时，我并没有set t 1 40 ，那系统会不会默认这个变量计算后也包含了40个时次的数据？

polestar7

做个霸气的木头 发表于 2012-5-18 12:37

                               
登录/注册后可看大图

你文件求夏季平均的时候就出问题了，你这样不会得到多个年份夏季平均的summer，每次的summer会被替换掉 ...

你好，我想问一下，你提取的夏季平均可以多个年份的放在同一个nc里么？因为我用循环的话，提取出的值只有一个时刻的值，是不是需要用多个nc来放多年的值？

mofangbao

做出来是这样的么，两个变量的距平？

做个霸气的木头

polestar7 发表于 2012-5-18 17:09

                               
登录/注册后可看大图

你好，我想问一下，你提取的夏季平均可以多个年份的放在同一个nc里么？因为我用循环的话，提取出的值只有 ...

你没有明白我的意思，我是说年份是循环，每一次循环求 玩夏季平均直接输出出来，这样循环几次，就会输出几年的夏季值了，这个数据文件里放的就是每年的夏季平均值：给你个我以前的参考下，求夏季平均输出的gs：
'reinit'
'sdfopen c:\data\sst.mnmean.nc'
'set gxout fwrite'
'set fwrite c:\data\asrsst.grd'
i=1409
while(i<=1757)
j=i+3
'asrt=ave(sst,t+'i',t='j',1)'
'set x 1 180'
'set y 1 89'
'd asrt'
i=i+12
endwhile
'disable fwrite'
;