爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 877|回复: 2

求助ERA5数据如何把世界时转换为北京时

[复制链接]

新浪微博达人勋

发表于 2023-6-26 15:30:57 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

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'黑体')
垂直速度-时间.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2023-6-26 15:51:01 | 显示全部楼层
tt = tt + 3.6 * 8
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2023-6-26 16:25:55 | 显示全部楼层
万分感谢王老师!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表