- 积分
- 3079
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2015-11-2
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
# add data
f = addfile('E:/20130627.nc')
data = f['Vo']['3',:,'31:37','110:117']*100000
lev = array([200, 300, 500, 700, 850, 925, 1000])
lev = meteo.pressure_to_height_std(lev)/1000.0
data.setdimvalue(0, lev)
figure(figsize=[1000,1000],newfig=False)
tdata_2km = data['1.5',:,:]
tdata_4km = data['3',:,:]
tdata_6km = data['5.5',:,:]
tdata_9km = data['9',:,:]
#plot
levs = arange(-10,10,1)
cols=makecolors(len(levs)+1,cmap='NCV_jet')
ax = axes3d(position=[0, 0.1, 0.85, 0.88],tickfontsize=18,bbox=True)
ax.set_elevation(-75)
#ax.set_rotation(180)
lprovince = shaperead('E:/map/city.shp')
layer=ax.contourf(tdata_2km,levs,colors=cols,offset=1.5,alpha=1,zdir='z',smooth=True)
ax.plot_layer(lprovince,edgecolor='gray',offset=1.5)
layer=ax.contourf(tdata_6km,levs,colors=cols,offset=5.5,alpha=1,zdir='z',smooth=True)
ax.plot_layer(lprovince,edgecolor='gray',offset=5.5)
layer=ax.contourf(tdata_9km,levs,colors=cols,offset=9,alpha=1,zdir='z',smooth=True)
ax.plot_layer(lprovince,edgecolor='gray',offset=9)
xlim(111,117)
ylim(31,36.5)
zlim(0,10)
ylabel('Latitude ',fontsize=18)
xlabel('Longitude ',fontsize=18)
zlabel('Altitude (km)',fontsize=18)
colorbar(layer,extendrect=False,ticks=arange(-10,10,2),shrink=0.8,label=u'涡度$(×10^{-5}·s^{-1})$',labelloc='right',fontsize=18,fontname=u'宋体',yshift=55)
savefig('E:/Vo_20130627.jpg', dpi = 300,width=800, height=800)
王老师,三维图colorbar的lable中公式无法正确显示,平面图是可以正常显示出来的,请王老师指导。另外3.4.11版本改为中文字体也无法正常显示u'中文'的用法,是不是三维图有所调整。 |
-
|