- 积分
- 541
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-10-19
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
python读取:
uwnd_925 = nc_u.variables['uwnd'][:, 1, lat_index_B:lat_index_E + 1, lon_index_B:lon_index_E + 1]
uwnd_850 = nc_u.variables['uwnd'][:, 2, lat_index_B:lat_index_E + 1, lon_index_B:lon_index_E + 1]
uwnd_700 = nc_u.variables['uwnd'][:, 3, lat_index_B:lat_index_E + 1, lon_index_B:lon_index_E + 1]
vwnd_925 = nc_v.variables['vwnd'][:, 1, lat_index_B:lat_index_E + 1, lon_index_B:lon_index_E + 1]
vwnd_850 = nc_v.variables['vwnd'][:, 2, lat_index_B:lat_index_E + 1, lon_index_B:lon_index_E + 1]
vwnd_700 = nc_v.variables['vwnd'][:, 3, lat_index_B:lat_index_E + 1, lon_index_B:lon_index_E + 1]
# 求平均风速
wnd_925 = np.sqrt(uwnd_925 ** 2 + vwnd_925 ** 2)
wnd_850 = np.sqrt(uwnd_850 ** 2 + vwnd_850 ** 2)
wnd_700 = np.sqrt(uwnd_700 ** 2 + vwnd_700 ** 2)
# plotplt.plot(year_list, wnd_925_y, 'r', lw=2.2, ls='-', c='red', label='925hPa')
plt.plot(year_list, wnd_850_y, 'r', lw=2.2, ls='-', c='blue', label='850hPa')
plt.plot(year_list, wnd_700_y, 'r', lw=2.2, ls='-', c='green', label='700hPa')
感觉也没有错
|
-
|