GrADS和NCL求出来区域平均不一样，求解！！！

qq469015280

aa22556688 发表于 2015-4-13 09:09
那在NCL里面哪个求平均的函数是加权的呢？能否告诉一下。
是用dim_avg_wgt_n_Wrap 这个函数么？

请问我用latRegWgt算wgty，再用wgt_areaave_Wrap（q,wgty,wgtx,opt）算区域平均，这样对吗？，
还有wgtx 该怎么算啊？？？

藕茭白空心菜

vivianal 发表于 2015-11-15 17:49
个人认为，如果计算的是中高纬度则需要加上权重，如果是较低纬则不需要，保险起见的话还是都用吧

请问为什么呢

vivianal

藕茭白空心菜 发表于 2017-3-17 13:10
请问为什么呢

因为中高纬度的地图投影变形比低纬度的严重。

藕茭白空心菜

vivianal 发表于 2017-3-24 09:11
因为中高纬度的地图投影变形比低纬度的严重。

谢谢啦~有没有数据每个格点都是等距的呢，就是不存在中高纬度的地图投影变窄的问题

东风破浪

求平均一般都是要进行面积加权的！

藕茭白空心菜

东风破浪 发表于 2017-3-24 09:24
求平均一般都是要进行面积加权的！

请问加权的话ncl用哪个函数

vivianal

藕茭白空心菜 发表于 2017-3-24 09:23
谢谢啦~有没有数据每个格点都是等距的呢，就是不存在中高纬度的地图投影变窄的问题

按我的理解，一般等距的格点数据应该都有变形，高斯格点不清楚。

藕茭白空心菜

vivianal 发表于 2017-3-26 22:22
按我的理解，一般等距的格点数据应该都有变形，高斯格点不清楚。

谢谢啦~~~~

vivianal

藕茭白空心菜 发表于 2017-3-24 09:55
请问加权的话ncl用哪个函数

  面积加权
   rad    = 4.*atan(1.)/180.
   clat   = input_sst->lat({-10:10})
   clat   = sqrt( cos(rad*clat) )


   wVAR   = ssta                                  ; copy meta data
   wVAR   = ssta*conform(ssta, clat, 1)
   wVAR@long_name = "Wgt: "+wVAR@long_name
    ssta_avg=dim_avg_n_Wrap(wVAR(:,:,:),(/1,2/))  ;get nino3 index,area mean

藕茭白空心菜

vivianal 发表于 2017-4-10 10:00
面积加权
   rad    = 4.*atan(1.)/180.
   clat   = input_sst->lat({-10:10})

非常感谢~~