- 积分
- 80
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-4-16
- 最后登录
- 1970-1-1

|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
import numpy as np
import matplotlib.pyplot as plt
import xarray as xr
import cartopy.crs as ccrs
import cartopy.feature as cf
import cartopy.io.shapereader as shpreader
from cartopy.mpl.ticker import LongitudeFormatter,LatitudeFormatter
from cartopy.mpl.gridliner import LONGITUDE_FORMATTER, LATITUDE_FORMATTER
import matplotlib.ticker as mticker
plt.rcParams['font.sans-serif']=['SimHei']
data=xr.open_dataset(r'C:\Users\Administrator\Desktop\UV\uwnd.mon.mean.nc')
data1=xr.open_dataset(r'C:\Users\Administrator\Desktop\UV\vwnd.mon.mean.nc')
uwnd1=data.groupby('time.month').mean() #计算多年的U风月平均
vwnd1=data1.groupby('time.month').mean() #计算多年的V风月平均
lon=uwnd1['lon'][:]#读取经度,一维的
lat=uwnd1['lat'][:]#读取纬度,一维的
UGRD=uwnd1['uwnd'][:][:][:][:]#读取纬向风,四维的
VGRD=vwnd1['vwnd'][:][:][:][:]#读取经向风,四维的
pres= uwnd1['level'][:]#读取气压层数,一维的
montn1=uwnd1['month'][:]#读取月份,一维的
fig=plt.figure(figsize=(5,3),dpi=700)#添加画布
ax=fig.add_axes([0,0,1,1])#添加子图
ax.invert_yaxis()#反转纵轴,使1000hPa作为起点
ax.set_xlabel('月份',fontsize=7)
ax.set_ylabel('层次(hPa)',fontsize=7)
ax.tick_params(axis='both',which='both',labelsize=7)
ax.barbs(montn1[:],pres[:],UGRD[:,:,26,102],VGRD[:,:,26,102],
barb_increments={'half':2,'full':4,'flag':20},length=5)
|
|