dim_avg_n函数求平均

hnqx

有一组grib2数据，利用systemfunc，以及ListSetType等批量读入，需要求某一段时间的海平面气压平均，利用dim_avg_n(p,0)可以得到所读取的所有grib2数据的平均，但是有个问题，比如我读入了9个时间段的grib2数据,时间是从2014年6月1日00时开始,6月3日06时结束，但是我只想求其中某一段时间，譬如6月1日00时-2日00时，这时候问题就出来了： dim_avg_n(p,0)是对p的第一维度求平均，但是如何选择性的对所读取的数据（譬如第一个时间段到第四个时间段）求平均呢，第一次发帖，问题表达的不是很清楚，请见谅！




-------------------------------------------------------
begin
;---Open some netCDF files.
diri = "/home/tl/grib2/"   ; input directory
fils = systemfunc ("ls " + diri + "*.grib2" ) ; file paths
f    = addfiles (fils, "r")   
ListSetType (f, "join")       ; join
; print(f)
;---Read variables off the file and do some conversions.
p    = f[:]->PRES_P0_L101_GLL0
p_ave=dim_avg_n(p,0)
;printVarSummary(p_ave)  
xwks=gsn_open_wks("ps","prs")
plot=gsn_contour(xwks,p_ave,False)

end

-----------------------------------------------------------
p变量如下，他是一个三维数组，第一个维度可以理解为时间，

Variable: p
Type: float
Total Size: 2345760 bytes
            586440 values
Number of Dimensions: 3
Dimensions and sizes: [ncl_join | 9] x [lat_0 | 181] x [lon_0 | 360]
Coordinates:
            lat_0: [90..-90]
            lon_0: [ 0..359]

hnqx

有同学帮忙解答一下吗？？

clarmy

dim_avg_n(p(a:b,:,:),0)
a、b是时间的起终点。

hnqx

clarmy 发表于 2015-8-31 20:20
dim_avg_n(p(a:b,:,:),0)
a、b是时间的起终点。

可以了，感谢，思维被套到笼子里去了，哈哈

clarmy

hnqx 发表于 2015-8-31 20:46
可以了，感谢，思维被套到笼子里去了，哈哈

不客气

当时年少春衫薄

clarmy 发表于 2015-8-31 20:20
dim_avg_n(p(a:b,:,:),0)
a、b是时间的起终点。

你好如果求那个经纬度区域的平均，但是分辨率是0.5*0.5的应该怎么办呢

clarmy

当时年少春衫薄 发表于 2016-5-29 15:50
你好如果求那个经纬度区域的平均，但是分辨率是0.5*0.5的应该怎么办呢

按角标计算起点和终点。
比如，0.5*0.5的分辨率，经纬范围是10*10。那么实际上的格点数就是21*21，那么经纬坐标中（1.5,2.5）的角标实际上就是（1.5/0.5-1=2,2.5/0.5-1=4）。

当时年少春衫薄

clarmy 发表于 2016-5-31 23:14
按角标计算起点和终点。
比如，0.5*0.5的分辨率，经纬范围是10*10。那么实际上的格点数就是21*21，那么 ...

哦哦是这样，谢谢啦

好久不见!

clarmy 发表于 2015-8-31 20:20
dim_avg_n(p(a:b,:,:),0)
a、b是时间的起终点。

您好，如果想要计算的是指定高度范围内的平均应该怎么做呢？请问您知道吗?我只会插值到某一高度层，但是不会求指定高度范围内的平均呢

冰糖雪梨Anais

hnqx 发表于 2015-8-31 20:46
可以了，感谢，思维被套到笼子里去了，哈哈

dim_avg_n(p(a:b,:,:),0)
a、b是时间的起终点。
请问时间的起终点是需要用索引表示吗？如果是，怎么得到需要求平均的时间起终点的索引？