月平均求季节平均权重

一大碗年糕

试了下geocat的month_to_season（原ncl month_to_season)函数，与自己计算的两种重考虑权重的月平均均有差异。说明该函数确实没考虑不同月天数的权重，与考虑权重时的误差在E-2，如果不在意这点误差可以考虑使用

1. ### gc ###
   
2. gc_ds = gc.climatologies.month_to_season(ds.sel(time=slice("1940","2022")),"JJA")
   
3. 
   
4. 
   
5. ### weighted ###
   
6. ds_season = ds.groupby('time.season')['JJA']
   
7. month_length = ds_season.time.dt.days_in_month
   
8. month_length
   
9. # 按“ time.season”分组来计算权重。
   
10. weights = month_length.groupby('time.season')['JJA'].groupby('time.year') / month_length.groupby('time.season')['JJA'].groupby('time.year').sum()
    
11. weights
    
12. # 计算加权平均值
    
13. ds_season_weighted = (ds_season * weights).groupby('time.year').sum(dim='time')

复制代码

内马尔

好的，谢谢楼主

一大碗年糕

内马尔 发表于 2023-7-25 12:53
好的，谢谢楼主