ERA5数据meteoinfo剖面图

lixiliang2008

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")

920daibingbing

请老师有时间帮忙看看，谢谢！
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")

大漠汪洋

杠杠滴，优秀

华仔_beyond

小白请教个问题啊，请问大佬ERA5 single levels，pressure levels和ERA5-Land
这三种数据有什么区别啊？

loftydongshi

华仔_beyond 发表于 2021-1-27 12:51
小白请教个问题啊，请问大佬ERA5 single levels，pressure levels和ERA5-Land
这三种数据有什么区别啊？

我了解的，single leves都是地面或近地面单层的要素，pressure是高空多层的，具体的要素和分层你可以直接去网站看一下就清楚了，land没有接触过。

lixiliang2008

已附上剖面图！

920daibingbing

你好，我按照您的程序改了一下路径，为什么报错，能帮忙看看吗，谢谢！
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

MeteoInfo

920daibingbing 发表于 2021-7-29 21:47
你好，我按照您的程序改了一下路径，为什么报错，能帮忙看看吗，谢谢！
f = addfile("D:/nc/20210621.nc") ...

你看看是不是第二行代码缩进的问题，Python最重要的语法就是代码缩进。

920daibingbing

老师你好，确实是缩进问题，但坐标轴反了，要怎么调整，谢谢！

MeteoInfo

920daibingbing 发表于 2021-7-30 13:30
老师你好，确实是缩进问题，但坐标轴反了，要怎么调整，谢谢！

哪个轴反了？把图贴出来看看

920daibingbing

纵坐标反了