新手求教，用ncl求区域平均值

张丹瑜婷

我想用求ncl求华北平原（北纬36.5°-40.5°，东经114.5°-119.5°）区域相对湿度（格点资料）的平均值，有没有大神指个路，我在哪可以找到类似的程序，让我参考一下，我不太会，谢谢大家！！！

三岁最可爱

f=addfile(filename,"r")
; 如果你的资料是四维资料（t，lev，lat，lon），可以把该区域取出来再做个lat和lon维的平均就行了
var_1=f->var(:, :, {36.5:40.5}, {114.5:119.5})
var_2=dim_avg_n_Wrap(var_1, (/2,3/))

HDdd

http://www.ncl.ucar.edu/Document ... n/wgt_areaave.shtml

_﹎__╱/●勄

三岁最可爱 发表于 2020-3-25 12:37
f=addfile(filename,"r")
; 如果你的资料是四维资料（t，lev，lat，lon），可以把该区域取出来再做个lat和 ...

你好，看到了你的方法，我尝试了一下，规定经纬度范围就出现了如下错误：fatal:["FileSupport.c":148]:Dimension (west_east) of file (wrfout_d03_2020-06-20_12) does not have an associated coordinate variable。不规定经纬度范围qsnow=a->QSNOW(:,:,:,:)   qsnow_ave=dim_avg_n_Wrap(qsnow, (/2,3/))，就是正确的，请问是什么原因。

九十

_﹎__╱/●勄 发表于 2020-6-27 21:26
你好，看到了你的方法，我尝试了一下，规定经纬度范围就出现了如下错误：fatal:["FileSupport.c":148]:Di ...

请问你的问题解决了吗？

天边的那边

请问你的问题解决了吗？

sourpomelo

_﹎__╱/●勄 发表于 2020-6-27 21:26
你好，看到了你的方法，我尝试了一下，规定经纬度范围就出现了如下错误：fatal:["FileSupport.c":148]:Di ...

应该是你的变量没有定义坐标变量，你可以用printVarSummary(qsnow)输出看看，如果是这样自己定义一下经纬度变量就可以了