请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 8186|回复: 17

ERA5数据meteoinfo剖面图

[复制链接]

新浪微博达人勋

发表于 2021-1-27 09:57:42 | 显示全部楼层 |阅读模式

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

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

x
f = addfile("F:/download.nc")
gh= f['z'][:,::-1,"40","115"]
udata = f['u'][:,::-1,"40","115"]
vdata = f['v'][:,::-1,"40","115"]
rmdata=f['r'][:,::-1,"40","115"]
tmdata=f['t'][:,::-1,"40","115"]-273.15
Tmdata=f['time'][:]
wmdata=f['w'][:,::-1,"40","115"]
u=udata.T*2.5
v=vdata.T*2.5
r=rmdata.T
t=tmdata.T
w=wmdata.T
levs1 = arange(60,100,10)
cols = [(255,255,255),(0,229,238),(0,238,0),(0,205,0),(0,160,0),(0,110,0)]
layer1 = contourf(r,levs1,colors=cols)
colorbar(layer1,fontsize=16,bold="Ture")
layer2 = barbs(u,v,color='k',size=12)
levs = arange(-60,36,2)
bss2 = []
for i in range(len(levs)):#温度样式
    bs2={}
    bs2['value'] = levs[i]
    bs2['color'] = 'r'
    bs2['size'] = 2
    if levs[i] < 0:
        bs2['linestyle'] = '--'
    bss2.append(bs2)
ls2 = plt.makesymbolspec('line',bss2)
layer3 = contour(t,symbolspec=ls2)
clabel(layer3)

#levw=arange(-10,10,0.1)
#layer4=contour(w,levw,color='k',size=2)
#clabel(layer4)

lev1 = tmdata.dimvalue(0)
lev1=lev1[::-1]
lev2 = meteo.pressure_to_height_std(lev1)
lev2 = lev2[:]/1000
tdata.setdimvalue(0, lev2)

xaxis(axistype='time',timetickformat='ddHH')
xticks(fontsize=15,bold=True)
xlabel('Time(UTC)',fontsize=15,bold="Ture")
ylim(480,1020)
levy=array([1000,975,950,925,900,875,850,825,800,775,750,700,650,600,550,500])
yticks(t.dimvalue(0),levy,bold="Ture")
ylabel('Pressure (hPa)',fontsize=15,bold="Ture")

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2021-7-30 17:25:20 | 显示全部楼层
请老师有时间帮忙看看,谢谢!
f = addfile("D:/nc/20210621sjcz.nc")
udata = f['u'][:,::-1,"26.68","100.25"]
vdata = f['v'][:,::-1,"26.68","100.25"]
rmdata=f['r'][:,::-1,"26.68","100.25"]
tmdata=f['t'][:,::-1,"26.68","100.25"]-273.15
Tmdata=f['time'][:]
wmdata=f['w'][:,::-1,"26.68","100.25"]
u=udata.T*2.5
v=vdata.T*2.5
r=rmdata.T
t=tmdata.T
w=wmdata.T
levs1 = arange(30,100,10)
cols = [(255,255,255),(0,229,238),(0,238,0),(0,205,0),(0,160,0),(0,110,0)]
layer1 = contourf(r,levs1,colors=cols)
colorbar(layer1,fontsize=16,bold="Ture")
layer2 = barbs(u,v,color='k',size=12)

levw=arange(-10,10,0.1)
layer4=contour(w,levw,color='k',size=2)
clabel(layer4)

lev1 = wmdata.dimvalue(0)
lev1=lev1[::-1]
lev2 = meteo.pressure_to_height_std(lev1)
lev2 = lev2[:]/1000
wmdata.setdimvalue(0, lev2)
rmdata.setdimvalue(0, lev2)
udata.setdimvalue(0, lev2)
vdata.setdimvalue(0, lev2)
xaxis(axistype='time',timetickformat='ddHH')
xticks(fontsize=15,bold=True)
xlabel('Time(UTC)',fontsize=15,bold="Ture")
#ylim(480,1020)
levy=array([750,700,650,600,550,500,450,400,350,300,250,225,200,175,150,125,100])
yticks(w.dimvalue(0),levy,bold="Ture")
ylabel('Pressure (hPa)',fontsize=15,bold="Ture")

20210621sjcz.nc

12.92 MB, 下载次数: 26, 下载积分: 金钱 -5

密码修改失败请联系微信:mofangbao
回复 支持 0 反对 1

使用道具 举报

新浪微博达人勋

发表于 2021-1-27 10:15:26 | 显示全部楼层
杠杠滴,优秀
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-1-27 12:51:08 | 显示全部楼层
本帖最后由 华仔_beyond 于 2021-1-27 12:52 编辑

小白请教个问题啊,请问大佬ERA5 single levels,pressure levels和ERA5-Land
这三种数据有什么区别啊?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-1-27 14:14:50 | 显示全部楼层
华仔_beyond 发表于 2021-1-27 12:51
小白请教个问题啊,请问大佬ERA5 single levels,pressure levels和ERA5-Land
这三种数据有什么区别啊?

我了解的,single leves都是地面或近地面单层的要素,pressure是高空多层的,具体的要素和分层你可以直接去网站看一下就清楚了,land没有接触过。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2021-1-30 16:23:06 | 显示全部楼层
已附上剖面图!
poumiantu.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-29 21:47:33 | 显示全部楼层
你好,我按照您的程序改了一下路径,为什么报错,能帮忙看看吗,谢谢!
f = addfile("D:/nc/20210621.nc")
gh= f['z'][:,::-1,"26","100"]
udata = f['u'][:,::-1,"26","100"]
vdata = f['v'][:,::-1,"26","100"]
rmdata=f['r'][:,::-1,"26","100"]
tmdata=f['t'][:,::-1,"26","100"]-273.15
Tmdata=f['time'][:]
wmdata=f['w'][:,::-1,"26","100"]
u=udata.T*2.5
v=vdata.T*2.5
r=rmdata.T
t=tmdata.T
w=wmdata.T
levs1 = arange(60,100,10)
cols = [(255,255,255),(0,229,238),(0,238,0),(0,205,0),(0,160,0),(0,110,0)]
layer1 = contourf(r,levs1,colors=cols)
colorbar(layer1,fontsize=16,bold="Ture")
layer2 = barbs(u,v,color='k',size=12)
levs = arange(-60,36,2)
bss2 = []
for i in range(len(levs)):#温度样式
    bs2={}
     bs2['value'] = levs[i]
     bs2['color'] = 'r'
     bs2['size'] = 2
     if levs[i] < 0:
         bs2['linestyle'] = '--'
     bss2.append(bs2)
ls2 = plt.makesymbolspec('line',bss2)
layer3 = contour(t,symbolspec=ls2)
clabel(layer3)

#levw=arange(-10,10,0.1)
#layer4=contour(w,levw,color='k',size=2)
#clabel(layer4)

lev1 = tmdata.dimvalue(0)
lev1=lev1[::-1]
lev2 = meteo.pressure_to_height_std(lev1)
lev2 = lev2[:]/1000
tdata.setdimvalue(0, lev2)

xaxis(axistype='time',timetickformat='ddHH')
xticks(fontsize=15,bold=True)
xlabel('Time(UTC)',fontsize=15,bold="Ture")
ylim(480,1020)
levy=array([750,700,650,600,550,500,450,400,350,300,250,225,200,175,150,125,100])
yticks(t.dimvalue(0),levy,bold="Ture")
ylabel('Pressure (hPa)',fontsize=15,bold="Ture")


报错:
>>> run script...
  File "<iostream>", line 2
    gh= f['z'][:,::-1,"26","100"]
    ^
SyntaxError: mismatched input '' expecting EOF
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-30 08:38:18 | 显示全部楼层
920daibingbing 发表于 2021-7-29 21:47
你好,我按照您的程序改了一下路径,为什么报错,能帮忙看看吗,谢谢!
f = addfile("D:/nc/20210621.nc") ...

你看看是不是第二行代码缩进的问题,Python最重要的语法就是代码缩进。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-30 13:30:10 | 显示全部楼层
老师你好,确实是缩进问题,但坐标轴反了,要怎么调整,谢谢!
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-30 14:36:00 | 显示全部楼层
920daibingbing 发表于 2021-7-30 13:30
老师你好,确实是缩进问题,但坐标轴反了,要怎么调整,谢谢!

哪个轴反了?把图贴出来看看
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2021-7-30 14:49:29 | 显示全部楼层
纵坐标反了
1.png
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

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

本版积分规则

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

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

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