想用ncl把一个数组求和，出现了错误

q863510557

出错的那一截如下所示，刚开始有一个三维数组，时间（每小时），经度，纬度。时间跨度为一个月，我想求每小时的月平均，于是就先把noo2（：，：，：）赋值给一个二维的求和数组su（：，：），然后再用这个二维数组除以总天数求平均。循环24次。但是每次都有报错如下：
su=new((/50,50/),"float")
ts=1
do ts=1,24
ttt=ts
do while(ttt.le.696)
su(:,:)=su(:,:)+noo2(ttt,:,:)
ttt=ttt+24
end do
noo2ave(ts,:,:)=su/28
su(:,:)=0
end do
报错是这样的：
fatal:Plus: Dimension size, for dimension number 0, of operands does not match, can't continue
会的同志们麻烦教一下我，不胜感激。

longlivehj

noo2_4d = reshape(noo2, (/24, 696 / 24, 50, 50/))
noo2avg = dim_avg_n(noo2_4d, 1)

longlivehj

至于你程序里的错误，请告知：
1、出错的具体行（ncl错误提示里有）；
2、noo2的维数信息。696小时是29天，为何平均时除以28？

另外，ncl下标都是从0开始，ts=1这种操作很容易下标越界！

q863510557

longlivehj 发表于 2014-4-11 18:03
至于你程序里的错误，请告知：
1、出错的具体行（ncl错误提示里有）；
2、noo2的维数信息。696小时是29天 ...

谢谢大神，我知道错误了，我的数据只有28天的，写成696个小时，数据超出了。

longlivehj

q863510557 发表于 2014-4-11 18:31
谢谢大神，我知道错误了，我的数据只有28天的，写成696个小时，数据超出了。

哦，好的！
没必要用循环，两行就够了。
noo2_4d = reshape(noo2, (/24, 28, 50, 50/))
noo2avg = dim_avg_n(noo2_4d, 1)

q863510557

longlivehj 发表于 2014-4-11 18:33
哦，好的！
没必要用循环，两行就够了。
noo2_4d = reshape(noo2, (/24, 28, 50, 50/))

恩恩  这两行我要好好研究研究