MeteoInfoLab脚本汇总贴

redersun

MeteoInfo 发表于 2022-3-11 08:47
单步调试功能目前没有，可以选中部分代码，点击鼠标邮件，在弹出菜单中选择 Evaluate Selection 运行被选 ...

好的，谢谢老师

ta383045387

谢谢老师分享

gaosrn

谢谢王老师分享{:5_213:}

石榴之爱

麻烦请教王老师，在MeteoInfo网站上有contourfm用法的介绍及举例，下图一。在模仿这个例子的过程中，layer=contourfm(ps,20)这里出错，出错原因如下图二，不知道什么问题也不会解决。

MeteoInfo

石榴之爱 发表于 2022-10-5 22:21
麻烦请教王老师，在MeteoInfo网站上有contourfm用法的介绍及举例，下图一。在模仿这个例子的过程中，layer= ...

用MeteoInfo最新版本试试

石榴之爱

MeteoInfo 发表于 2022-10-6 22:13
用MeteoInfo最新版本试试

好的，谢谢王老师！

Beckyeahyeah

老师有FY4a的例子吗？

MeteoInfo

Beckyeahyeah 发表于 2022-11-7 13:31
老师有FY4a的例子吗？

参考这里：http://bbs.06climate.com/forum.p ... &extra=page%3D1

guo066600

王老师，有MI对高维数据聚类的例子？

风格而才

麻烦王老师给看一下时间垂直剖面图，画出来的图两边有空白是怎么回事，单独画相对湿度没问题，加了风场就有白边了，数据太大传不上去。

1. f = addfile(r'C:\Users\Administrator\Desktop\618_Geo.nc')
   
2. lev1 = f['level'][1:10]
   
3. tdata = f['r'][0:44,'1000:200','38.47','106.2']
   
4. tdata = tdata[0:44,::-1]
   
5. data = tdata.T
   
6. #高度转化      
   
7. lev1=lev1[::-1]
   
8. lev2= meteo.p2h(lev1)
   
9. levels = []
   
10. for j in range(0, len(lev1)):
    
11.     levels.append('%i' % lev1[j])  
    
12. #时间轴   
    
13. tt = f.gettimes()[0:44]
    
14. timstrs = []
    
15. for tim in tt:
    
16.     tim = tim + datetime.timedelta(hours=8)
    
17.     timstrs.append(tim.strftime('%m/%d/%H'))
    
18. #风场
    
19. u = f['u'][0:44,'1000:200','38.47','106.2']
    
20. u = u[0:44,::-1]
    
21. u = u.T
    
22. v = f['v'][0:44,'1000:200','38.47','106.2']
    
23. v = v[0:44,::-1]
    
24. v = v.T
    
25. speed = sqrt(u*u+v*v)
    
26. 
    
27. levs = arange(70, 100, 10)
    
28. cols = makecolors(len(levs)+1,cmap='suyang')
    
29. cols[0]='w'
    
30. layer = contourf(array(tt), lev2, data, levs, colors=cols)
    
31. colorbar(layer)
    
32. xx, yy = meshgrid(array(tt), lev2)
    
33. layer= barbs(xx, yy, u, v, speed,color='k', size = 10)
    
34. xticks(tt, timstrs, rotation=0)
    
35. yticks(lev2,levels)
    
36. antialias(True)

复制代码