WRF_out文件绘制降水空间填色图时出现白线

taiyang

使用ERA5作为初始场，运行WRF模式，对输出的文件进行处理绘制，使用python降水填色空间场。出现以下问题：
- 当指定levels时，在赤道以北3°左右会出现一条白线，已经检测了该范围内的数据结果，并没有异常值
- 只有使用python提取数据绘制填色图时出现此问题；
- ncl、ncview、xarray.plot()都没有该问题



绘图code如下：

1. from wrf import getvar,get_cartopy
   
2. import xarray as xr
   
3. import numpy as np
   
4. import matplotlib.pyplot as plt
   
5. import cartopy.crs as ccrs
   
6. import cartopy.feature as cfeature
   
7. from cartopy.mpl.ticker import LongitudeFormatter,LatitudeFormatter
   
8. from netCDF4 import Dataset as nc
   
9. import cmaps
   
10. 
    
11. ncfile=nc(r'D:/wrf/wrfout_d01_2004-06-25_04_00_00.nc')

复制代码

希望大家提提建议

taiyang

1. from wrf import getvar,get_cartopy
   
2. import xarray as xr
   
3. import numpy as np
   
4. import matplotlib.pyplot as plt
   
5. import cartopy.crs as ccrs
   
6. import cartopy.feature as cfeature
   
7. from cartopy.mpl.ticker import LongitudeFormatter,LatitudeFormatter
   
8. from netCDF4 import Dataset as nc
   
9. import cmaps
   
10. 
    
11. ncfile=nc(r'D:/wrf/wrfout_d01_2004-06-25_04_00_00.nc')
    
12. 
    
13. 
    
14. slp=getvar(ncfile, "slp")
    
15. 
    
16. rain_exp=np.array(getvar(ncfile, "RAINNC"))
    
17. rain_con=np.array(getvar(ncfile, "RAINC"))
    
18. rain_tot=np.array(rain_exp+rain_con)
    
19. rain1=getvar(ncfile, "RAINNC")
    
20. rain2=getvar(ncfile, "RAINC")
    
21. lon=getvar(ncfile, "XLONG")
    
22. lat=getvar(ncfile, "XLAT")
    
23. rain=(rain1+rain2)
    
24. 
    
25. fig = plt.figure(figsize=(10, 14),dpi=100)
    
26. 
    
27. 
    
28. box=[0,180,-30,30.5]
    
29. scale='50m'
    
30. xstep, ystep = 40, 15
    
31. 
    
32. #######plot####################
    
33. def make_map(ax, title):
    
34.     # set_extent  set crs
    
35.     ax.set_extent(box, crs=ccrs.PlateCarree())
    
36.     land = cfeature.NaturalEarthFeature('physical',
    
37.                                         'land',
    
38.                                         scale,
    
39.                                         edgecolor='grey',
    
40.                                         facecolor='none',zorder=1
    
41.                                         )
    
42.     ax.add_feature(land)  # set land color
    
43.     ax.coastlines(scale)  # set coastline resolution
    
44.     # set coordinate axis
    
45.     ax.set_xticks(np.arange(box[0],box[1], xstep),crs=ccrs.PlateCarree())
    
46.     ax.set_yticks(np.arange(box[2], box[3]+ystep, ystep),crs=ccrs.PlateCarree(central_longitude=0))
    
47.     ax.xaxis.set_major_formatter(LongitudeFormatter(zero_direction_label =False))#经度0不加标识
    
48.     ax.yaxis.set_major_formatter(LatitudeFormatter())
    
49. 
    
50.     ax.set_title(title, fontsize=15, loc='left',pad=10)
    
51.     ax.tick_params(which='major',
    
52.                     direction='out',
    
53.                     length=8,
    
54.                     width=0.99,
    
55.                     pad=3,
    
56.                     labelsize=15,
    
57.                     bottom=True, left=True, right=False, top=False)
    
58.     return ax
    
59. ax0=fig.add_subplot(311,projection=ccrs.PlateCarree())
    
60. ax0=make_map(ax0,'2004-06-25_04--RAINC(set levels(0,5.1,0.1))')
    
61. 
    
62. cb0=ax0.contourf(lon,lat,rain2,transform=ccrs.PlateCarree(),levels=np.arange(0,5.1,0.1))
    
63. c0=fig.colorbar(cb0,ax=ax0,shrink=0.8,pad=0.05,aspect=13)
    
64. 
    
65. ax1=fig.add_subplot(312,projection=ccrs.PlateCarree())
    
66. ax1=make_map(ax1,'2004-06-25_04--RAINNC(set levels(0,5.1,0.1))')
    
67. 
    
68. cb1=ax1.contourf(lon,lat,rain1,transform=ccrs.PlateCarree(),levels=np.arange(0,5.1,0.1))
    
69. c1=fig.colorbar(cb1,ax=ax1,shrink=0.8,pad=0.05,aspect=13)
    
70. 
    
71. ax2=fig.add_subplot(313,projection=ccrs.PlateCarree())
    
72. ax2=make_map(ax2,'2004-06-25_04--RAINNC+RAINC(set levels(0,5.1,0.5))')
    
73. 
    
74. cb2=ax2.contourf(lon,lat,rain,transform=ccrs.PlateCarree(),levels=np.arange(0,5.1,0.5))
    
75. c2=fig.colorbar(cb2,ax=ax2,shrink=0.8,pad=0.05,aspect=13)
    
76. 
    
77. 
    
78. plt.show()

复制代码

taiyang

问题解决了，原因是，模式中的lon的排列是（0-180,-180-0）导致的，将负的那一部分加上360.使得lon的排列范围是0-360，就解决该问题啦~

努力学python

taiyang 发表于 2022-11-11 10:12
问题解决了，原因是，模式中的lon的排列是（0-180,-180-0）导致的，将负的那一部分加上360.使得lon的排列范 ...

谢谢楼主！帮了我大忙