ERA5数据meteoinfo剖面图

MeteoInfo

920daibingbing 发表于 2021-7-30 14:49
纵坐标反了

能不能把脚本程序和数据文件贴出来，我有空了仔细看看。

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

920daibingbing

老师您能帮忙看看吗，我的纵坐标反了，谢谢！

edwardli

920daibingbing 发表于 2021-8-7 16:15
老师您能帮忙看看吗，我的纵坐标反了，谢谢！

如果是matplotlib，ax.invert_yaxis()可以反向，plt.yscale("log")可以对数坐标

GZH2021

920daibingbing 发表于 2021-8-7 16:15
老师您能帮忙看看吗，我的纵坐标反了，谢谢！

请问您解决了吗？我也是出现这个问题

MeteoInfo

参考这个：

1. f = addfile("D:/Temp/nc/20210621sjcz.nc")
   
2. udata = f['u'][:,::-1,"26.68","100.25"]
   
3. vdata = f['v'][:,::-1,"26.68","100.25"]
   
4. rmdata=f['r'][:,::-1,"26.68","100.25"]
   
5. tmdata=f['t'][:,::-1,"26.68","100.25"]-273.15
   
6. Tmdata=f['time'][:]
   
7. wmdata=f['w'][:,::-1,"26.68","100.25"]
   
8. u=udata.T*2.5
   
9. v=vdata.T*2.5
   
10. r=rmdata.T
    
11. t=tmdata.T
    
12. w=wmdata.T
    
13. 
    
14. level = r.dimvalue(0)
    
15. height = meteolib.pressure_to_height_std(level)
    
16. tt = r.dimvalue(1)
    
17. levs = arange(30,100,10)
    
18. cols = [(255,255,255),(0,229,238),(0,238,0),(0,205,0),(0,160,0),(0,110,0)]
    
19. layer1 = contourf(tt, height, r, levs, colors=cols)
    
20. colorbar(layer1, fontsize=16, bold="Ture")
    
21. ttm, heightm = meshgrid(tt, height)
    
22. layer2 = barbs(ttm, heightm, u, v, color='b', size=12)
    
23. 
    
24. levw=arange(-10,10,0.1)
    
25. layer4=contour(tt, height, w, levw, color='k', size=1)
    
26. clabel(layer4)
    
27. 
    
28. xaxis(axistype='time', timetickformat='ddHH')
    
29. xticks(fontsize=15, bold=True)
    
30. xlabel('Time(UTC)', fontsize=15, bold="Ture")
    
31. ylim(height[0], height[-1])
    
32. yticks(height, ['{:.0f}'.format(l) for l in level], bold="Ture")
    
33. ylabel('Pressure (hPa)', fontsize=15, bold="Ture")

f = addfile("D:/Temp/nc/20210621sjcz.nc")<br /> udata = f['u'][:,::-1,"26.68","100.25"]<br /> vdata = f['v'][:,::-1,"26.68","100.25"]<br /> rmdata=f['r'][:,::-1,"26.68","100.25"]<br /> tmdata=f['t'][:,::-1,"26.68","100.25"]-273.15<br /> Tmdata=f['time'][:]<br /> wmdata=f['w'][:,::-1,"26.68","100.25"]<br /> u=udata.T*2.5<br /> v=vdata.T*2.5<br /> r=rmdata.T<br /> t=tmdata.T<br /> w=wmdata.T<br /> <br /> level = r.dimvalue(0)<br /> height = meteolib.pressure_to_height_std(level)<br /> tt = r.dimvalue(1)<br /> levs = arange(30,100,10)<br /> cols = [(255,255,255),(0,229,238),(0,238,0),(0,205,0),(0,160,0),(0,110,0)]<br /> layer1 = contourf(tt, height, r, levs, colors=cols)<br /> colorbar(layer1, fontsize=16, bold="Ture")<br /> ttm, heightm = meshgrid(tt, height)<br /> layer2 = barbs(ttm, heightm, u, v, color='b', size=12)<br /> <br /> levw=arange(-10,10,0.1)<br /> layer4=contour(tt, height, w, levw, color='k', size=1)<br /> clabel(layer4)<br /> <br /> xaxis(axistype='time', timetickformat='ddHH')<br /> xticks(fontsize=15, bold=True)<br /> xlabel('Time(UTC)', fontsize=15, bold="Ture")<br /> ylim(height[0], height[-1])<br /> yticks(height, ['{:.0f}'.format(l) for l in level], bold="Ture")<br /> ylabel('Pressure (hPa)', fontsize=15, bold="Ture")

escaper

风应该不要乘2.5吧

escaper

MeteoInfo 发表于 2021-10-17 17:57
参考这个：

风应该不要乘2.5吧，单位是m/s