python求月异常？

是冉冉升起的冉

ncl中求月异常有这个函数calcMonAnomLLLT，最近学习python好像没有这种函数(或者是我没有发现?)，不知道python下面这样做有没有什么问题，分解时间维度；分解之后变量变成year*mon*lev*lat*lon,变量命名有什么办法呢？个人感觉好像有一点复杂化了。

1. # =============================================================================
   
2. 
   
3. time,lat,lon,lev=var_clim.time,var_clim.lat,var_clim.lon,var_clim.level
   
4. 
   
5. '''
   
6.     算月平均气候态 time:year*mon
   
7.     将数据从 timexlatxlon 转化为 year x mon x lat x lon
   
8.     对转换后的year维度求平均得到月平均的结果
   
9. '''
   
10. nmon=3
    
11. ntime,levs,lats,lons=var_clim.shape
    
12. 
    
13. #year mon lev lat lon
    
14. 
    
15. varym  = var_clim.data.reshape((int(ntime/nmon),nmon,levs,lats,lons))
    
16. varclm = varym.mean(axis=0) # along the year axis
    
17. 
    
18. var_pos_ym=var_pos.data.reshape((pos_year.shape[0],nmon,levs,lats,lons))                                          
    
19. var_pos_ym_ano=np.zeros((int(ntime/nmon),nmon,levs,lats,lons))
    
20.                            
    
21. 
    
22. '''
    
23. '''
    
24. for i in range(0, pos_year.shape[0]):
    
25.     var_pos_ym_ano[i,:,:,:,:]=var_pos_ym[i,:,:,:,:]-varclm[:,:,:,:]
    
26. 
    

复制代码

李知恩

月异常是什么意思呢？我的理解是不是这个月的值和前三十年这个月的距平呢？其实还是很好实现的，简单举个例子用python实现位势高度的距平：

李知恩

z500_month = z500_ds.sel(time=slice('1990','2022'),lon=lons,lat=lats,level=500.0).groupby('time.month').mean(dim='time')
z500_jun_clim = z500_month.sel(month=9)  # June climatology
z500_jul_22 = z500_ds.sel(time='2023-09', lon=lons, lat=lats, level=500.0).mean(dim=['time'])

是冉冉升起的冉

李知恩 发表于 2024-8-10 21:14
z500_month = z500_ds.sel(time=slice('1990','2022'),lon=lons,lat=lats,level=500.0).groupby('time.mont ...

哇呜 谢谢！！！