分享一套适用大气专业的Python编程经验 (最近更新：2021/10/21)

Penny_Kuan

biubiubiu123 发表于 2020-3-5 10:52
报错信息是什么呢？
```
a = np.array([-999.0,10,11])

我这会再试了一下，a[a==udf] = np.nan 没有问题了，可以替换缺测值，之前的具体报错我找不到了。
但是出现了新的问题，AttributeError: 'module' object has no attribute 'nanmean'。

lleoiu

biubiubiu123 发表于 2020-3-5 13:57
这个函数我并没有使用过，多文件处理我通常使用CDO提前合并，或者使用os.listdir()获取全部文件名然后排 ...

谢谢回复。
简单写了一个计算入偏差的脚本，正常运行。但是结果不对。也不知道哪里出错。

1. 
   
2. import netCDF4 as nc
   
3. import numpy as np
   
4. import xarray as xr           
   
5.    
   
6. dyst = 1
   
7. dyed = 7
   
8. ndy  = 7
   
9.                                     
   
10. dir2 = "/mnt/d/data/jra55/daily/surface/"                                                        
    
11. dir4 = "/mnt/d/data/jra55/long_term_mean/daily/surface/"                     
    
12.                                                                                                                                 
    
13. tmfile = dir2 + "air.*.nc"                                                                                                                                                                           
    
14. t2mfile = dir4 + "t2m.day.ltm.nc"                                            
    
15.                                                                                    
    
16. ds  = xr.open_mfdataset(tmfile,concat_dim="time",combine="by_coords", engine="netcdf4")
    
17. lon = ds.variables['lon'][:]      
    
18. lat = ds.variables['lat'][:]   
    
19. 
    
20. data = ds['air']-273.15
    
21. ds5  = xr.open_dataset(t2mfile)  
    
22. tc   =  ds5.variables['t2m'][dyst:dyed,:,:] - 273.15
    
23. 
    
24. airf = data.loc[data.time.dt.month==7]
    
25. air  = airf.loc[airf.time.dt.dayofyear<=dyed]
    
26. 
    
27. time  = air.time                 
    
28. airdt = air.values
    
29.                  
    
30. nlon = np.shape(air)[2]
    
31. nlat = np.shape(air)[1]
    
32. ntim = np.shape(air)[0]
    
33. nyr  = int(float(ntim)/ndy)
    
34. 
    
35. dtm  = np.zeros(np.shape(airdt))   
    
36. cdtm = np.zeros(np.shape(airdt))   
    
37. 
    
38. for i in range (0,ndy):                 
    
39.      cdtm[i:ndy:ntim,:,:] = tc[i,:,:]
    
40. 
    
41. dtm = airdt - cdtm   

复制代码

biubiubiu123

Penny_Kuan 发表于 2020-3-6 13:27
我这会再试了一下，a[a==udf] = np.nan 没有问题了，可以替换缺测值，之前的具体报错我找不到了。
但是 ...

nanmean的话我感觉只能用a = np.nanmean(b) ,用a = b.nanmean就会报错，具体原因我还没找到

wdd1992

有心了，学习楼主，感谢分享

杰杰的任

感谢楼主分享

Penny_Kuan

biubiubiu123 发表于 2020-3-6 14:15
nanmean的话我感觉只能用a = np.nanmean(b) ,用a = b.nanmean就会报错，具体原因我还没找到

就是用的np模块里的，可能我系统上这个模块代码有点问题，没有nanmean函数【捂脸】

量子侠

谢谢

一大碗年糕

给楼主点个赞虽然还在用NCL，想必将来还是要用python的哈哈哈先收藏着

WULI杨飘呀

码住，等学完python基本语法了就来看

jimolengshan

正在学习气象统计方法，这个帖子给力，学习ing