grads如何对月平均数据求逐年平均并做时间序列

LemonTree007

1.这里要对月平均数据求逐年平均，然后做出区域平均的时间序列图，求大神指教。
2.如果是三个时间连续的月平均nc文件，要求他们的逐年平均，再做出时间序列又该怎么写？
个人觉得要定义新变量，做每隔12个月的循环求逐年平均数据，但写不出脚本怎么办？急求指导。
下面是一个例子，grads显示数据太大（190101-201412的月平均数据），会很慢，我想知道我写的这个逻辑是否正确，有哪些不符合脚本语言的，跪求指正！

1. 'reinit'
   
2. 'sdfopen D:\gradsht\spei01.nc'
   
3. 'set lon 70 140'
   
4. 'set lat 10 60'
   
5. 'set lev 0'
   
6. 'set t 1 1368'
   
7. 'define speit=spei'
   
8. 'set t 1'
   
9. n=1
   
10. while(n<125)
    
11.     'speit.x.y.n=ave(spei,t+0,t=t+12).(lat=10,lat=60).(lon=70,lon=140)'
    
12.     n=n+1
    
13.     'set t t+12'
    
14. endwhile
    
15. 'set gxout line'   
    
16. 'd tloop(aave(speit,lon=70,lon=140,lat=10,lat=60))'
    
17. pull dummy
    
18. 'c'
    

复制代码

琼小爷

想知道楼主做出来了吗？也遇到了同样的问题

LemonTree007

琼小爷 发表于 2018-5-15 20:33
想知道楼主做出来了吗？也遇到了同样的问题

还没有哎，这几天没什么进展，做不出来就用这个图了，；论文挺急的

LemonTree007

琼小爷 发表于 2018-5-15 20:33
想知道楼主做出来了吗？也遇到了同样的问题

不好意思，看错帖子了，这个我做出来了。但比较复杂。附件是我的code。amon.gs是我把三个文件合并成一个grd文件。
amon-pdsi-1901-2014.ctl是我给grd写的ctl文件。
然后又写了一个clamon.gs将多年月平均数据处理成，逐月多年平均。

琼小爷

LemonTree007 发表于 2018-5-16 10:14
不好意思，看错帖子了，这个我做出来了。但比较复杂。附件是我的code。amon.gs是我把三个文件合并成一个g ...

好，谢谢楼主。我研究一下

琼小爷

LemonTree007 发表于 2018-5-16 10:14
不好意思，看错帖子了，这个我做出来了。但比较复杂。附件是我的code。amon.gs是我把三个文件合并成一个g ...

好，谢谢楼主。我研究一下

琼小爷

LemonTree007 发表于 2018-5-16 10:14
不好意思，看错帖子了，这个我做出来了。但比较复杂。附件是我的code。amon.gs是我把三个文件合并成一个g ...

楼主能打包成压缩包吗，不够金钱下

LemonTree007

琼小爷 发表于 2018-5-16 10:59
楼主能打包成压缩包吗，不够金钱下

金钱不够的话ctl你不用下载，自己写就行了，反正我们的数据不一样。