- 积分
- 720
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2014-12-29
- 最后登录
- 1970-1-1
|
楼主 |
发表于 2022-1-17 23:39:12
|
显示全部楼层
#read data
fn = r'D:/MeteoInfo_2.1.4/MeteoInfo/map/dixing-china.nc'
f = addfile(fn)
elev = f['z']['26:30','100:103']/1000.0
lchina = shaperead('D:/MeteoInfo_2.1.4/MeteoInfo/map/bou2_4p.shp')
#Plot
figure(figsize=[1400,850],newfig=False)
ax = axes3d(tickfontsize=20,position=[0.1, 0.13, 0.75, 0.85],bbox=False)
levs = arange(0,6.1,0.1)
cols = makecolors(len(levs) + 1, cmap='temp_diff_18lev',reverse=False)
ls = ax.plot_surface(elev, levs, colors=cols, edge=False)
colorbar(ls,ticks=arange(0,6.1,1),orientation='horizontal',label='Altitude(km)',labelloc='right',fontsize=24,shrink=0.8,aspect=40,yshift=35)
x=[102.2181,101.6]
y=[27.875,27.5]
z=[1.529,3]
points = ax.scatter(x,y, z)
colorbar(points,shrink=0.8)
ylabel('Latitude (degrees)',fontsize=20)
xlabel('Longitude (degrees)',fontsize=20)
zlabel('Altitude (km)',fontsize=20)
zlim(0, 6)
xlim(101.3,103)
ylim(27.3,28.2) |
|