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

藕茭白空心菜

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

为什么只加权这么少部分的纬度呢{-10:10}

vivianal

藕茭白空心菜 发表于 2017-4-10 11:04
为什么只加权这么少部分的纬度呢{-10:10}

看你算的区域啦，我当时只算了-10到10度

藕茭白空心菜

vivianal 发表于 2017-4-10 18:18
看你算的区域啦，我当时只算了-10到10度

嗯，谢谢你的耐心解答

MARZOLINI

麦田_smile 发表于 2015-9-19 17:51
我试了下 的确wgt加上权重就和grads差不多了 那么到底一般的区域平均要不要加权重呢？例如北大西洋的海温 ...

请问一下，你用ncl求得区域权重平均怎么算的，我用ncl求得区域平均加了权重怎么还是和grads的结果差很多

麦田_smile

MARZOLINI 发表于 2018-4-14 14:26
请问一下，你用ncl求得区域权重平均怎么算的，我用ncl求得区域平均加了权重怎么还是和grads的结果差很多

    lat    = f3->lat({0:60})
    printVarSummary(lat)
   wgtx   = latRegWgt(lat, "float", 0)
    lat01  = 0
    lat02  = 60
    lon01  = 280
    lon02  = 360           
    index  = wgt_areaave_Wrap(x3(time|:,{lat|lat01:lat02},{lon|lon01:lon02}),wgtx,1.0,0)

MARZOLINI

麦田_smile 发表于 2018-4-15 23:21
lat    = f3->lat({0:60})
    printVarSummary(lat)
   wgtx   = latRegWgt(lat, "float", 0)

嗯嗯，谢谢，我再试试

莎朗嘿哟YY

麦田_smile 发表于 2018-4-15 23:21
lat    = f3->lat({0:60})
    printVarSummary(lat)
   wgtx   = latRegWgt(lat, "float", 0)

您好，我想问，如果是模式数据网格做区域平均还需要先做面积加权吗？

莎朗嘿哟YY

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

您好，我想问，如果是模式数据网格用ncl做区域平均还需要先做面积加权吗？

卢小傻子

个人觉得，GrADS中aave是加权的。ncl那样算应该没加权。见GrADS官网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

Xiaoooooooo-

麦田_smile 发表于 2015-9-19 17:51
我试了下 的确wgt加上权重就和grads差不多了 那么到底一般的区域平均要不要加权重呢？例如北大西洋的海温 ...

您好！请问对站点数据进行区域平均，数据是（station，year），我直接用dim_avg_n_wrap对station这维求平均可以吗？谢谢您！