MeteoInfoLab脚本汇总贴

karstification

谢谢

风格而才

麻烦王老师帮忙看一下哪出问题了，相对湿度叠加风场的空间垂直剖面，只显示填色的相对湿度数据，不显示风场。

1. f = addfile(r'D:\script_and_date\date\yanglie\EC.nc')
   
2. t = f.gettime(0)
   
3. print(t)
   
4. u = f['u'][0,:,'31:41','70:140:2']
   
5. u = mean(u,axis=1)
   
6. v = f['v'][0,:,'31:41','70:140:2']
   
7. v = mean(v,axis=1)
   
8. speed = sqrt(u*u+v*v)
   
9. speed = speed[::-1,:]
   
10. 
    
11. tdata = f['r'][0,:,'31:41','70:140']
    
12. tdata = mean(tdata,axis=1)
    
13. tdata = tdata[::-1,:]
    
14. 
    
15. lev1 = tdata.dimvalue(0)
    
16. lev1=lev1[::-1]
    
17. lev2 = meteo.pressure_to_height_std(lev1)
    
18. lev2 = lev2[:]/1000
    
19. tdata.setdimvalue(0, lev2)
    
20. speed.setdimvalue(0, lev2)
    
21. 
    
22. fn1 = addfile(r'D:\script_and_date\date\yanglie\elev.0.25-deg.nc')
    
23. height=fn1['data'][0,'31:41',:]
    
24. height = mean(height,axis=0)
    
25. height = height/1000
    
26. x= fn1['lon'][:]
    
27. 
    
28. axes()
    
29. levs = arange(10,100,10)
    
30. cols = makecolors(len(levs)+1, cmap='RH', reverse=True)
    
31. layer = contourf(tdata,levs,colors=cols)
    
32. layer1 = barbs(u, v, speed,color='k', size = 10)
    
33. fill_between(x,height,color='gray')
    
34. plot(x,height,color='k')
    
35. yaxis(tickin=False,tickfontsize=17)
    
36. xaxis(tickin=False,tickfontsize=17)
    
37. yaxis(location='right',tickin=True,tickfontsize=18)
    
38. yticks(tdata.dimvalue(0), lev1)
    
39. ylim(lev2.min(),lev2.max())
    
40. xlim(90,120)
    
41. colorbar(layer,orientation='horizontal',aspect=35)
    
42. antialias(True)

复制代码

MeteoInfo

风格而才 发表于 2021-3-9 15:37
麻烦王老师帮忙看一下哪出问题了，相对湿度叠加风场的空间垂直剖面，只显示填色的相对湿度数据，不显示风场 ...

1. f = addfile(r'D:\Temp\nc\EC.nc')
   
2. t = f.gettime(0)
   
3. print(t)
   
4. u = f['u'][0,:,'31:41','70:140:2']
   
5. u = mean(u,axis=1)
   
6. v = f['v'][0,:,'31:41','70:140:2']
   
7. v = mean(v,axis=1)
   
8. speed = sqrt(u*u+v*v)
   
9. speed = speed[::-1,:]
   
10. 
    
11. tdata = f['r'][0,:,'31:41','70:140']
    
12. tdata = mean(tdata,axis=1)
    
13. tdata = tdata[::-1,:]
    
14. 
    
15. lev1 = tdata.dimvalue(0)
    
16. lev1=lev1[::-1]
    
17. lev2 = meteo.pressure_to_height_std(lev1)
    
18. lev2 = lev2[:]/1000
    
19. tdata.setdimvalue(0, lev2)
    
20. speed.setdimvalue(0, lev2)
    
21. 
    
22. fn1 = addfile(r'D:\Temp\nc\elev.0.25-deg.nc')
    
23. height=fn1['data'][0,'31:41',:]
    
24. height = mean(height,axis=0)
    
25. height = height/1000
    
26. x= fn1['lon'][:]
    
27. 
    
28. axes()
    
29. levs = arange(10,100,10)
    
30. cols = makecolors(len(levs)+1, reverse=True)
    
31. layer = contourf(tdata,levs,colors=cols)
    
32. xx, yy = meshgrid(u.dimvalue(1), lev2)
    
33. layer1 = barbs(xx, yy, u, v, speed,color='k', size = 10)
    
34. fill_between(x,height,color='gray')
    
35. plot(x,height,color='k')
    
36. yaxis(tickin=False,tickfontsize=17)
    
37. xaxis(tickin=False,tickfontsize=17)
    
38. yaxis(location='right',tickin=True,tickfontsize=18)
    
39. yticks(tdata.dimvalue(0), lev1)
    
40. ylim(lev2.min(),lev2.max())
    
41. xlim(90,120)
    
42. colorbar(layer,orientation='horizontal',aspect=35)
    
43. antialias(True)

f = addfile(r'D:\Temp\nc\EC.nc')<br /> t = f.gettime(0)<br /> print(t)<br /> u = f['u'][0,:,'31:41','70:140:2']<br /> u = mean(u,axis=1)<br /> v = f['v'][0,:,'31:41','70:140:2']<br /> v = mean(v,axis=1)<br /> speed = sqrt(u*u+v*v)<br /> speed = speed[::-1,:]<br /> <br /> tdata = f['r'][0,:,'31:41','70:140']<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 /> speed.setdimvalue(0, lev2)<br /> <br /> fn1 = addfile(r'D:\Temp\nc\elev.0.25-deg.nc')<br /> height=fn1['data'][0,'31:41',:]<br /> height = mean(height,axis=0)<br /> height = height/1000<br /> x= fn1['lon'][:]<br /> <br /> axes()<br /> levs = arange(10,100,10)<br /> cols = makecolors(len(levs)+1, reverse=True)<br /> layer = contourf(tdata,levs,colors=cols)<br /> xx, yy = meshgrid(u.dimvalue(1), lev2)<br /> layer1 = barbs(xx, yy, u, v, speed,color='k', size = 10)<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 /> ylim(lev2.min(),lev2.max())<br /> xlim(90,120)<br /> colorbar(layer,orientation='horizontal',aspect=35)<br /> antialias(True)

风格而才

MeteoInfo 发表于 2021-3-9 17:47

谢谢 王老师出来了！

YyyRT

谢谢老师的分享啦！！

Joan0125

老师太棒了！！！佩服佩服

gaowl

非常强大，学习。

非对称

老师，脚本示例数据哪里下载啊？

MeteoInfo

非对称 发表于 2021-7-23 14:54
老师，脚本示例数据哪里下载啊？

MeteoInfo软件里自带了一些，MeteoInfo网站上也有一些数据文件：http://www.meteothink.org/downloads/data.html

dianzimi

谢谢王老师