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

LiuJiawei · 发表于 2015-2-13 16:29:46

longlivehj 发表于 2015-2-13 15:48
什么错误？
这种情况，直接avg就行了。

报错是incomplete information to create return value form builtin function
avg果然有用，谢谢！
另外，请教一下，a(lat,Lon),b(lat,Lon)求它们的相关系数用哪个函数？

1luochen1 · 发表于 2015-4-9 21:37:28

楼主最后对比两个函数的结果如何？

happyaliyun · 发表于 2015-12-26 14:12:21

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

什么叫没有初始化呐~~~~~~

我还是叫小确幸 · 发表于 2016-4-26 18:12:25

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

在计算区域平均~很有帮助

当时年少春衫薄 · 发表于 2016-5-29 15:47:14

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

你好如果是某一经纬度范围内的情况呢

MARZOLINI · 发表于 2018-4-14 14:32:47

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

你好，我想问一下就是求区域平均的时候，需要加权重吗，我用ncl加了权重算出来的区域平均为什么会和grads求出来的区域平均相差很多，这是我用ncl求区域平均的脚本，能帮我看一下吗，谢谢。
rad = 4.*atan(1.)/180.
clat = fspan(30,18,25)
clat = sqrt(cos(rad*clat))
wVAR =tp1
wVAR =tp1*conform(tp1, clat, 1)
wVAR@long_name = "Wgt: "+wVAR@long_name
tp1_avg=dim_avg_n_Wrap(wVAR(:,:,:),(/1,2/))

x^2 · 发表于 2018-10-23 08:40:17

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

这个是不是等权重的平均，对于大的区域的格点资料是不是要做球面平均？用这个函数是个lat*cos(纬度)吗？

CTCTCTCT · 发表于 2020-5-31 22:33:34

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

你好，我也遇到这个问题，我是要用new定义的是四维数组，请问这个该怎么初始化呢？

		自动登录	找回密码
密码			立即注册