import xarray as xr
import matplotlib.pyplot as plt
import numpy as np
file=r'D:\data\data\fnl_20210617_12_00.grib2'
RH=xr.open_dataset(file,engine='cfgrib',backend_kwargs={'filter_by_keys': {'typeOfLevel': 'isobaricInhPa','cfVarName': 'r'}})
lon=RH['longitude'].loc[70:140]
level=RH['isobaricInhPa'][:]
r=RH['r'].loc[:,30,70:140]
fig=plt.figure(figsize=(5,3),dpi=500)
ax=fig.add_axes([0,0,1,1])
ac=ax.contourf(lon,level,r,cmap='Greens',levels=np.arange(0,110,10))
cb=fig.colorbar(ac,ax=ax,pad=0.01)
cb.ax.set_ylabel('relative humidity(%)',fontsize=9)
import cartopy.mpl.ticker as cmt
ax.xaxis.set_major_formatter(cmt.LongitudeFormatter())#给横坐标加经度符号
ax.invert_yaxis()#颠倒y轴
ax.set_yticks([1000,925,850,700,500,300,200,100])
ax.set_xlabel('Longitude',fontsize=12)
ax.set_ylabel('Level(hPa)',fontsize=12)
ax.set_title('Lon-Lev Cross Section',loc='left',fontsize=15)
ax.set_title('Lat=30°N',loc='right',fontsize=7)
气象剖面图
最新推荐文章于 2024-07-21 23:25:11 发布