- 积分
- 78
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-11-26
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
用ERA5数据画了个时间剖面图,想把时间轴改成北京时,求教该怎么转换
f = addfile("k:/MI/618.nc")
udata = f['u'][:,::-1,"36.65","117.15"]
vdata = f['v'][:,::-1,"36.65","117.15"]
rmdata=f['r'][:,::-1,"36.65","117.15"]
tmdata=f['t'][:,::-1,"36.65","117.15"]-273.15
Tmdata=f['time'][:]
wmdata=f['w'][:,::-1,"36.65","117.15"]
u=udata.T*2.5
v=vdata.T*2.5
r=rmdata.T
t=tmdata.T
w=wmdata.T
level = w.dimvalue(0)
height = meteolib.pressure_to_height_std(level)
tt = w.dimvalue(1)
levs = arange(-0.8,1.0,0.1)
layer1 = contourf(tt, height, w, levs)
colorbar(layer1, fontsize=16, bold="Ture")
ttm, heightm = meshgrid(tt, height)
#layer2 = barbs(ttm, heightm, u, v, color='b', size=12)
xaxis(axistype='time', timetickformat='ddHH')
xticks(fontsize=15, bold=True)
xlabel('Time(UTC)', fontsize=15, bold="Ture")
ylim(height[0], height[-1])
yticks(height, ['{:.0f}'.format(l) for l in level], bold="Ture")
ylabel('Pressure (hPa)', fontsize=15, bold="Ture")
title(u'垂直速度(Pa/s)', fontsize=20,fontname=u'黑体') |
|