提取想要月份做平均

杰西卡·芬

ncl有month_to_annual、month_to_season等内置函数，但只能计算12个月或者某三个月的逐年平均值。若想要挑选其中几个月进行分析则没有内置函数，相关帖子回复得少(可能是简单，别人都不愿意讲吧)。想要提取每年4-9月(共有150年)降水数据，例子如下。有更简便的方法可以多多交流。
f=addfile("D:\pr_185001-200012.nc","r")
pr=f->pr
time=f->time
lat=f->lat
lon=f->lon
ntim=dimsizes(time)
nlat=dimsizes(lat)
nlon=dimsizes(lon)
nyears=ntim/12
pr4d=reshape(pr,(/nyears,12,nlat,nlon/))                ;用reshape重新定义数据维度。提取前后维度的乘积要相等。
pr_4_9=reshape(pr4d(:,3:8,:,:),(/nyears*6,nlat,nlon/))
j=0
pr150=new((/150,nlat,nlon/),float,"No_FillValue")
do i =0,149
   pr150(i,:,:)=(pr_4_9(j,:,:)+pr_4_9(j+1,:,:)+pr_4_9(j+2,:,:)+pr_4_9(j+3,:,:)+pr_4_9(j+4,:,:)+pr_4_9(j+5,:,:))/6.0
   j=j+6
end do
;printVarSummary(pr150)
pr150_ave=dim_avg_n_Wrap(pr150,(/1,2/))
printVarSummary(pr150_ave)

lt5201314

楼主，我现在遇到了这个问题，有点类似，做的多年逐月资料算逐月气候态，不是12的整数倍没办法直接用函数求，请问楼主有什么解决办法吗

杰西卡·芬

lt5201314 发表于 2020-9-17 19:41
楼主，我现在遇到了这个问题，有点类似，做的多年逐月资料算逐月气候态，不是12的整数倍没办法直接用函数求 ...

可能你的数据中有一年的数据只有几个月。你把你的数据切开两份来算（没满12月份的和满12月的分开算）。请我我讲清楚了吗？

lt5201314

分开算具体是怎么操作呢，我之前打算在ncl里写循环，把不同年份的同一月份提取出来，但后来实现不了

杰西卡·芬

lt5201314 发表于 2020-9-18 14:09
分开算具体是怎么操作呢，我之前打算在ncl里写循环，把不同年份的同一月份提取出来，但后来实现不了

就是先算前面没有满12个月那年的平均。后面都有12个月的年份再平均。把不同年的同一月份提取出来算平均好像有函数的。你可以翻一下那些大神的帖子，我记得是有的。

lt5201314

嗯嗯，你说的那个可以把前面第一年的先算，再算后面的，但第一年的数据怎么加进去就是个问题了，我翻了两天了，都没看到

杰西卡·芬

lt5201314 发表于 2020-9-18 14:53
嗯嗯，你说的那个可以把前面第一年的先算，再算后面的，但第一年的数据怎么加进去就是个问题了，我翻了两天 ...

画图的时候加进去就可以了。就相当于一张图用了两组数据

lt5201314

嗯嗯，我平均做出来了，正在出图了，谢谢指点

niuyingli

lt5201314 发表于 2020-9-19 06:44
嗯嗯，我平均做出来了，正在出图了，谢谢指点

怎么做的呢，求指点

lt5201314

niuyingli 发表于 2020-11-15 11:24
怎么做的呢，求指点

我用的纯数学计算，先算12个月整的，再把多的月份加进去算平均，不过气候态一般取30年，后来我直接截取30年的数据进行计算的。