计算青高热源和降水的时间序列相关，gs文件有些问题请教各位

youhuiqi7455

'reinit' 'set grid off' 'set grads off' 'set mpdset cnworld' 'open E:\lunwen\data\JRA-55-Q1Q2\q1nwns.ctl' 'open E:\lunwen\jiangshui\1.ctl' 'set t 1 12' 'define q11=ave(q1,t+0,t=444,12)' 'define q12=ave(q1,t+1,t=444,12)' 'define q13=ave(q1,t+11,t=444,12)' 'define qq1=(q11+q12+q13)/3' 'set t 1 12' 'define precip1=ave(precip.2,t+0,t=444,12)' 'define precip2=ave(precip.2,t+1,t=444,12)' 'define precip3=ave(precip.2,t+11,t=444,12)' 'define ppre=(precip1+precip2+precip3)/3' 'set x 1' 'set y 1' 'set z 1' 'set t 1 444' 'define a=aave(qq1,lon=73.3,lon=104.8,lat=26,lat=39.8)' 'set dfile 2' 'set lon 72 137.5' 'set lat 15 57' 'set z 1' 'set t 1 444' 'd tcorr(a,ppre,t=1,t=444)' 代码如上，我求的是1980-2017平均冬季青高热源q1和冬季降水之间的时间序列相关，出图是可以出图，但是我不管怎么改变t+0,t+1,t+11这个求平均的地方图片都一样，比如改成夏季t+5,t+6,t+7图还是一样的，但是如果热源和降水不一样，比如一个是t+0,t+1,t+11，另一个t+5,t+6,t+7图就会改变，这是为什么呢？ctl是没有问题的

苏拉苏拉

有个疑问哦，冬季的话应该是连续的本年12月和次年1月，2月。你t+0,t+1,t+11不成了本年度1月，2月，12月。

youhuiqi7455

苏拉苏拉 发表于 2020-4-30 15:03
有个疑问哦，冬季的话应该是连续的本年12月和次年1月，2月。你t+0,t+1,t+11不成了本年度1月，2月，12月。

没关系的吧，因为需要的是37年整体的平均，只需要算出对应月份的均值，然后再求3个月均值，这样每一年的冬季月份都算进去就可以了吧？ 这里我的数据没有1979的12月，但是应该影响不大

youhuiqi7455

问题貌似解决了，图片是不一样的，只是我眼瞎了

youhuiqi7455

冬季和夏季的图很相似，只是一个是正值一个是负值，一时没发现区别，以为自己写错了

youhuiqi7455

不对，还是一样的，问题没解决啊，有没有大佬帮忙看看

youhuiqi7455

貌似找到问题所在了，因为我降水的ctl是自己写的，批量描述了1980-2017的nc逐月降水资料，而每个nc文件的起始时间不一样，是不是就不能参与grads的运算呢？