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

aa22556688

求解为什么GrADS（图1）和NCL（图2） 求出来全球SST平均温度不一样，而且差异很大？NCL初学者跪求高手指点，已经求气候态和异常值20多天了，始终发现有问题，跪求指点{:soso_e101:}！！！


Grads脚本：
'sdfopen E:\Data\sst.mnmean.nc'
'set x 1'
'set y 1'
'set z 1'
'set t 1 1932'
'd tloop(aave(sst,lon=0,lon=360,lat=-90,lat=90))'


NCL脚本：
begin
fsst      = addfile(dir+fil,"r")
st        = fsst->sst
global  = dim_avg_n_Wrap(st(:,:,:),(/1,2/))
print(global)
exit
end

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

niuda

同样求问，发现好多的函数ncl的结果和grads或fortran处理出来的结果不一样的

尽头的尽头

看你的数据是不是short型的，如果是，在ncl里面要转换数据类型

longlivehj

aave()      Calculates an areal average over an X-Y region (latitude-weighted)

aa22556688

niuda 发表于 2015-4-13 06:46
同样求问，发现好多的函数ncl的结果和grads或fortran处理出来的结果不一样的

可能是对NCL函数没有深入理解，不是很了解，所以就会出现应用错误的情况。

aa22556688

尽头的尽头 发表于 2015-4-13 08:35
看你的数据是不是short型的，如果是，在ncl里面要转换数据类型

数据是float型的，应该不用转吧

aa22556688

longlivehj 发表于 2015-4-13 09:02
aave()      Calculates an areal average over an X-Y region (latitude-weighted)

那在NCL里面哪个求平均的函数是加权的呢？能否告诉一下。
是用dim_avg_wgt_n_Wrap 这个函数么？

longlivehj

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

The aave function assumes that the world coordinates are longitude in the X dimension and latitude in the Y dimension, and does weighting in the latitude dimension by the difference between the sines of the latitude at the northern and southern edges of the grid box. For areal averaging without latitude weighting, use the amean function.

根据上面红色部分提示，跟ncl中latRegWgt函数（http://www.ncl.ucar.edu/Document ... ted/latRegWgt.shtml）求取权重的方式一样。楼主试试先用latRegWgt获得权重，然后用dim_avg_wgt_Wrap平均。

King@Control

grads里amean为等权

niuda

longlivehj 发表于 2015-4-13 09:21
The aave function assumes that the world coordinates are longitude in the X dimension and latitude ...

仅仅是求区域平均的话，不用加权重的吧，直接用dim_avg_n_Wrap（*, (/1,2/)）不可以的吗？