在NCL中手动求区域平均得到的都是缺省值

Castalia

do tt =0,dimsizes(t_ec&initial_time0_hours)-1
  do jj = 0,dimsizes(t_ec&lat_0)-1
    do ii = 0,dimsizes(t_ec&lon_0)-1
    spread_area_ave(tt) = spread_area_ave(tt)+spread_area(tt,jj,ii)
   end do
  end do
  spread_area_ave(tt)=spread_area_ave(tt)*1.0/(dimsizes(t_ec&lat_0)*dimsizes(t_ec&lon_0))
end do
print(spread_area_ave)

想要用这个来跟wgt_areaave_Wrap函数对比一下

longlivehj

可能是spread_area_ave没定义好。
如果是用new定义的，可能没有初始化，默认全为缺省值，导致后面的累加不起作用。

longlivehj

fealbai 发表于 2014-9-21 09:56
没注意看问题，我都是做两次算区域平均，

你是说两次调用dim_avg_n么？
这个函数的第二个参数可以是数组，比如a(time, level, lat, lon)这样的一个变量，求区域平均的话，可以是dim_avg_n(a, (/2, 3/))。

fealbai

区域平均，dim_avg_n求区域平均好像不会，试试

Castalia

longlivehj 发表于 2014-9-20 17:35
可能是spread_area_ave没定义好。
如果是用new定义的，可能没有初始化，默认全为缺省值，导致后面的累加不 ...

对的~~谢啦！！

Castalia

fealbai 发表于 2014-9-20 17:56
区域平均，dim_avg_n求区域平均好像不会，试试

这个不是求区域平均的应该。。。

fealbai

Castalia 发表于 2014-9-20 22:47
这个不是求区域平均的应该。。。

没注意看问题，我都是做两次算区域平均，

fealbai

longlivehj 发表于 2014-9-21 15:31
你是说两次调用dim_avg_n么？
这个函数的第二个参数可以是数组，比如a(time, level, lat, lon)这样的一 ...

学到了！

LiuJiawei

longlivehj 发表于 2014-9-21 15:31
你是说两次调用dim_avg_n么？
这个函数的第二个参数可以是数组，比如a(time, level, lat, lon)这样的一 ...

如果数组是两维的比如a( lat, lon)， b=dim_avg_n(a, (/0, 1/))，运行报错得不到结果，该如何解决呢？

longlivehj

LiuJiawei 发表于 2015-2-13 15:42
如果数组是两维的比如a( lat, lon)， b=dim_avg_n(a, (/0, 1/))，运行报错得不到结果，该如何解决呢？

什么错误？
这种情况，直接avg就行了。