- 积分
- 1274
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-7-30
- 最后登录
- 1970-1-1

|
发表于 2016-11-10 15:02:21
|
显示全部楼层
cdo可以做,ncl里直接构成新数组也可以做,我做了一个30年11月到3月的平均,亲测有效,主要如下:
date = cd_calendar(uwnd&time, -1)
; print(date)
dStrt = 198111
iStrt = ind(date.eq.dStrt)
; print(iStrt)
dims = dimsizes(uwnd)
nlat = dims(1)
nlon = dims(2)
uwnd500 = new((/30,nlat,nlon/),float,-32767)
do i = 0,29
uwnd500(i,:,:) = uwnd(0,:,:)
end do
do i = 0,29
uwnd500(i,:,:) = (uwnd(iStrt,:,:)+uwnd(iStrt+1,:,:)+uwnd(iStrt+2,:,:)+uwnd(iStrt+3,:,:)+uwnd(iStrt+4,:,:))/5.0
iStrt = iStrt+12
end do
; print(iStrt)
; printVarSummary(uwnd500)
; delete(uwnd)
clm_u500 = dim_avg_n_Wrap(uwnd500, 0) |
|