用ncl手动求区域平均，得到的都是缺省值

李小毛123

用ncl手动求区域平均，得到的都是缺省值，代码如下
ty=new(51,double)
ty=0
do tt=1,51
do i=0,nlat-1
  do j=0,nlon-1
   ty(tt)=ty(tt)+dtdy1(tt,i,j)
  end do
end do
ty=ty(tt)/2529  ;2529是总的格点数
print(ty)
end do
最后得到的结果是(0)     9.969209968386869e+36
(1)     9.969209968386869e+36
(2)     9.969209968386869e+36
(3)     9.969209968386869e+36
(4)     9.969209968386869e+36
(5)     9.969209968386869e+36
(6)     9.969209968386869e+36
(7)     9.969209968386869e+36
(8)     9.969209968386869e+36
(9)     9.969209968386869e+36
(10)    9.969209968386869e+36
(11)    9.969209968386869e+36
(12)    9.969209968386869e+36
(13)    9.969209968386869e+36
(14)    9.969209968386869e+36
(15)    9.969209968386869e+36
(16)    9.969209968386869e+36
(17)    9.969209968386869e+36
(18)    9.969209968386869e+36
(19)    9.969209968386869e+36
(20)    9.969209968386869e+36
(21)    9.969209968386869e+36
(22)    9.969209968386869e+36
(23)    9.969209968386869e+36
(24)    9.969209968386869e+36
(25)    9.969209968386869e+36
(26)    9.969209968386869e+36
(27)    9.969209968386869e+36
(28)    9.969209968386869e+36
(29)    9.969209968386869e+36
(30)    9.969209968386869e+36
(31)    9.969209968386869e+36
(32)    9.969209968386869e+36
(33)    9.969209968386869e+36
(34)    9.969209968386869e+36
(35)    9.969209968386869e+36
(36)    9.969209968386869e+36
(37)    9.969209968386869e+36
(38)    9.969209968386869e+36
(39)    9.969209968386869e+36
(40)    9.969209968386869e+36
(41)    9.969209968386869e+36
(42)    9.969209968386869e+36
(43)    9.969209968386869e+36
(44)    9.969209968386869e+36
(45)    9.969209968386869e+36
(46)    9.969209968386869e+36
(47)    9.969209968386869e+36
(48)    9.969209968386869e+36
(49)    9.969209968386869e+36
(50)    9.969209968386869e+36
不太明白是怎么回事，望各位大神不吝赐教。

songwei

ty是新定义的数组，都是缺省值，在循环里即使加上别的东西，等号左端的东西永远是缺省。

李小毛123

songwei 发表于 2017-10-10 18:19
ty是新定义的数组，都是缺省值，在循环里即使加上别的东西，等号左端的东西永远是缺省。

ty=new(51,double)
ty=0
m=2529
do tt=1,51
do i=0,nlat-1
  do j=0,nlon-1
    if(any(ismissing(dtdy1(tt,i,j)))) then
   dtdy1(tt,i,j)=0
   m=m-1
else
  ty(tt)=ty(tt)+dtdy1(tt,i,j)
    end if
  end do
end do
ty=ty(tt)*(m^(-1))  ;2529是总的格点数
print(ty)
end do我把ty=0赋值为0不行吗？我把缺省值单拎出来之后，出来的结果是(0)     -9.530448159410801
(1)     -9.530448159410801
(2)     -9.530448159410801
(3)     -9.530448159410801
(4)     -9.530448159410801
(5)     -9.530448159410801
(6)     -9.530448159410801
(7)     -9.530448159410801
(8)     -9.530448159410801
(9)     -9.530448159410801
(10)    -9.530448159410801
(11)    -9.530448159410801
(12)    -9.530448159410801
(13)    -9.530448159410801
(14)    -9.530448159410801
(15)    -9.530448159410801
(16)    -9.530448159410801

packard

对于一个时间，如果水平区域内任一地点有缺省值，那这个点就会把整个求和弄成缺省值。
建议直接用函数吧，简单一些。
ty = dim_avg_n(dtdy1,(/1,2/))

橙子鲜

建议同上楼，ncl还是用函数吧！！

李小毛123

packard 发表于 2017-10-11 06:31
对于一个时间，如果水平区域内任一地点有缺省值，那这个点就会把整个求和弄成缺省值。
建议直接用函数吧， ...

谢谢谢谢，我开始就是用的那个函数，计算结果有问题，上午我重新计算了，结果出来了，还是谢谢啦。

李小毛123

橙子鲜 发表于 2017-10-11 10:57
建议同上楼，ncl还是用函数吧！！

恩恩，都感谢啦！

九曜.

我也遇到这样的问题 请问如何解决的？

雨里1998

packard 发表于 2017-10-11 06:31
对于一个时间，如果水平区域内任一地点有缺省值，那这个点就会把整个求和弄成缺省值。
建议直接用函数吧， ...

大佬，如果要求区域的加权平均值怎么做呢，dim_avg_wgt_n_Wrap中的权重只能是一维的

packard

雨里1998 发表于 2020-12-10 16:42
大佬，如果要求区域的加权平均值怎么做呢，dim_avg_wgt_n_Wrap中的权重只能是一维的

那可以先用reshape函数, 把空间那两维合并成一个维度, 然后再调 dim_avg_wgt_n