MeteoInfoLab脚本示例：双Y轴，纬度平均，经度剖面图（叠加地形）

︶ㄣ安定■丶

这里给出一个：利用欧洲中心下的w变量，指定纬度平均的经度剖面图，然后叠加该区域的地形剖面，最后使用类似NCL的双Y轴（气压和高度）绘制，
结果如下图所示：


具体脚本如下：

1. f = addfile(r'H:\test\data\test.nc')
   
2. tdata = f['w'][0,:,'35:50',:]
   
3. tdata = mean(tdata,axis=1)
   
4. tdata = tdata[::-1,:]
   
5. 
   
6. lev1 = tdata.dimvalue(0)
   
7. lev1=lev1[::-1]
   
8. lev2 = meteo.pressure_to_height_std(lev1)
   
9. lev2 = lev2[:]/1000
   
10. tdata.setdimvalue(0, lev2)
    
11. 
    
12. fn1 = addfile(r'H:\alldata\dixing/elev.0.5-deg.nc')
    
13. height=fn1['data']['35:50',:]
    
14. height = mean(height,axis=0)
    
15. height = height/1000
    
16. x= fn1['lon'][:]
    
17. 
    
18. ax1=axes()
    
19. levs = arange(-0.12,0.14,0.02)
    
20. layer = contourf(tdata,levs)
    
21. fill_between(x,height,color='gray')
    
22. plot(x,height,color='k')
    
23. yaxis(tickin=False,tickfontsize=17)
    
24. xaxis(tickin=False,tickfontsize=17)
    
25. yaxis(location='right',tickin=True,tickfontsize=18)
    
26. yticks(tdata.dimvalue(0), lev1)
    
27. xlim(70,140)
    
28. ylim(lev2.min(),12.001)
    
29. ylabel('Hpa',fontsize=18)
    
30. xlabel('')
    
31. colorbar(layer,orientation='horizontal',aspect=35)
    
32. 
    
33. ax2 = twinx(ax1)
    
34. ylim(lev2.min(),12.001)
    
35. yaxis(location='right',tickin=False,tickfontsize=18)
    
36. ylabel('Hight(km)',fontsize=18,bold=False)
    
37. antialias(True)
    
38. 
    

f = addfile(r'H:\test\data\test.nc')<br /> tdata = f['w'][0,:,'35:50',:]<br /> tdata = mean(tdata,axis=1)<br /> tdata = tdata[::-1,:]<br /> <br /> lev1 = tdata.dimvalue(0)<br /> lev1=lev1[::-1]<br /> lev2 = meteo.pressure_to_height_std(lev1)<br /> lev2 = lev2[:]/1000<br /> tdata.setdimvalue(0, lev2) <br /> <br /> fn1 = addfile(r'H:\alldata\dixing/elev.0.5-deg.nc')<br /> height=fn1['data']['35:50',:]<br /> height = mean(height,axis=0)<br /> height = height/1000<br /> x= fn1['lon'][:] <br /> <br /> ax1=axes() <br /> levs = arange(-0.12,0.14,0.02)<br /> layer = contourf(tdata,levs)<br /> fill_between(x,height,color='gray')<br /> plot(x,height,color='k')<br /> yaxis(tickin=False,tickfontsize=17)<br /> xaxis(tickin=False,tickfontsize=17)<br /> yaxis(location='right',tickin=True,tickfontsize=18)<br /> yticks(tdata.dimvalue(0), lev1)<br /> xlim(70,140)<br /> ylim(lev2.min(),12.001)<br /> ylabel('Hpa',fontsize=18)<br /> xlabel('')<br /> colorbar(layer,orientation='horizontal',aspect=35)<br /> <br /> ax2 = twinx(ax1)<br /> ylim(lev2.min(),12.001)<br /> yaxis(location='right',tickin=False,tickfontsize=18)<br /> ylabel('Hight(km)',fontsize=18,bold=False)<br /> antialias(True)<br /> <br />

MeteoInfo

已放入脚本汇总帖。

静言_GRMC

非常好  点赞

呼小喵喵

{:eb502:}{:eb502:}{:eb502:}{:eb502:}

独孤酒见

安定老师和王老师都是本版的高人，更是热心肠的人，本版的类似于我的有福气了。

小其其格

1. #Get file names
   
2. # coding=utf-8
   
3. fn = addfile('F:\wfl\MeteoInfoJS/wrf_allfields_d02_2013_07_29.ctl')
   
4. u = fn['U'][10,:,'28.0:28.0',:]
   
5. w = fn['W'][10,:,'28.0:28.0',:]
   
6. speed=sqrt(u*u+w*w)
   
7. #w=w*20
   
8. lev1 = u.dimvalue(0)
   
9. lev2= meteo.p2h(lev1)
   
10. #lev2 = meteo.pressure_to_height_std(lev1)
    
11. lev2=lev2/1000
    
12. u.setdimvalue(0, lev2)
    
13. w.setdimvalue(0, lev2)
    
14. layer=quiver(u,w,speed,15,size=10)
    
15. xlim(119,122.5)
    
16. ylim(0,12)
    
17. colorbar(layer)

复制代码

我这样绘图，流线没有显示呀，只有色标，怎么回事？

Jackluo

谢谢楼主分享~

river

谢谢分享！图很漂亮

goddess豆豆

为安老师打call

bianpuxi

请问MeteoInfoLab有能把sigma坐标直接转换为P坐标或者高度的函数吗