MeteoInfoLab脚本示例：micaps数据绘制垂直剖面图

MeteoInfo

micaps第4类数据是格点数据，每个文件只有一个变量一个时次一个高度。这里示例读取不同高度的micaps数据绘制垂直剖面图。

脚本程序：

1. levels = ['1000','925','850','700','500']
   
2. levs = []
   
3. for i in range(0, len(levels)):
   
4.     levs.append(i)
   
5. datadir = 'D:/Temp/test/micaps'
   
6. i = 0
   
7. data1 = []
   
8. for lev in levels:
   
9.     indir = os.path.join(datadir, lev)
   
10.     fn = os.path.join(indir, '15100808.000')
    
11.     f = addfile_micaps(fn)
    
12.     var = f['var']
    
13.     data = var[0,0,[40],:]
    
14.     if i == 0:
    
15.         lon = data.dimvalue(0)
    
16.     data1.append(data.aslist())
    
17.     i += 1
    
18. 
    
19. #Plot
    
20. levs = array(levs)
    
21. data1 = array(data1)
    
22. layer = contourf(lon, levs, data1,20)
    
23. yticks(levs, levels)
    
24. xlabel('Degree east')
    
25. ylabel('Pressure (hPa)')
    
26. colorbar(layer)
    
27. title('Vertical speed')

levels = ['1000','925','850','700','500']<br /> levs = []<br /> for i in range(0, len(levels)):<br />     levs.append(i)<br /> datadir = 'D:/Temp/test/micaps'<br /> i = 0<br /> data1 = []<br /> for lev in levels:<br />     indir = os.path.join(datadir, lev)<br />     fn = os.path.join(indir, '15100808.000')<br />     f = addfile_micaps(fn)<br />     var = f['var']<br />     data = var[0,0,[40],:]<br />     if i == 0:<br />         lon = data.dimvalue(0)<br />     data1.append(data.aslist())<br />     i += 1<br /> <br /> #Plot<br /> levs = array(levs)<br /> data1 = array(data1)<br /> layer = contourf(lon, levs, data1,20)<br /> yticks(levs, levels)<br /> xlabel('Degree east')<br /> ylabel('Pressure (hPa)')<br /> colorbar(layer)<br /> title('Vertical speed')

╯①路霺笑

Traceback (most recent call last):
  File "<iostream>", line 9, in <module>
NameError: name 'os' is not defined
原样复制，请问显示的这个是什么意思啊？
我是该软件的初学者，麻烦您不吝指教一下

╯①路霺笑

raceback (most recent call last):
  File "<iostream>", line 15, in <module>
AttributeError: 'NoneType' object has no attribute 'dimvalue'

MeteoInfo

╯①路霺笑 发表于 2016-9-2 19:34
raceback (most recent call last):
  File "", line 15, in
AttributeError: 'NoneType' object has no ...

要有相应的数据文件才可以i

liutong

如何做时间高度的剖面图呢？

半颗麦芽糖

王老师，最近看到您的这个脚本，想试着画一画，也找了相关的数据，但是按照您的这个脚本一直提示错误，貌似是维度的问题，想试着自己解决问题，但是刚接触这门语言有些语句还是不懂，比如data = var[[40],:]，我不明白后面的40以及中括号里面内容的意思，所以想请王老师指导一下，附上报错和文件，万分感谢！
indices must be 4 dimensions!
Traceback (most recent call last):
  File "<iostream>", line 15, in <module>
AttributeError: 'NoneType' object has no attribute 'dimvalue'
test.rar (8.39 KB, 下载次数: 14)

2016-11-4 21:53 上传

点击文件名下载附件

MeteoInfo

半颗麦芽糖 发表于 2016-11-4 21:43
王老师，最近看到您的这个脚本，想试着画一画，也找了相关的数据，但是按照您的这个脚本一直提示错误，貌似 ...

参考这里：http://www.meteothinker.com/docs ... /data_tutorial.html

半颗麦芽糖

MeteoInfo 发表于 2016-11-4 22:07
参考这里：http://www.meteothinker.com/docs/meteoinfolab/userguide/data_tutorial.html

王老师，根据您提供的参考网站我看懂了这个脚本并针对我的数据进行了修改，成功地出了图，非常感谢！我还有一个问题就是横坐标的纬度能否进行设定，比如设置成30N—40N，因为这个脚本中纬度是作为数据提取出来的，所以不知道该在哪里添加设置语言，前后试了lat = [30,40]都出了错，请王老师指导，谢谢！附上脚本和图片

1. levels = ['1000','925','850','700','500']
   
2. levs = []
   
3. for i in range(0, len(levels)):
   
4. levs.append(i)
   
5. datadir = 'F:/micaps/test'
   
6. i = 0
   
7. data1 = []
   
8. for lev in levels:
   
9. indir = os.path.join(datadir, lev)
   
10. fn = os.path.join(indir, '16102108.000')
    
11. f = addfile_micaps(fn)
    
12. var = f['var']
    
13. data = var[0,0,:,[90]]
    
14. if i == 0:
    
15. lat = data.dimvalue(0)
    
16. data1.append(data.aslist())
    
17. i += 1
    
18. 
    
19. #Plot
    
20. levs = array(levs)
    
21. data1 = array(data1)
    
22. layer = contour(lat, levs, data1,30)
    
23. yticks(levs, levels)
    
24. xlabel('Lat')
    
25. ylabel('Pressure (hPa)')
    
26. colorbar(layer)
    
27. title('Vertical vort')

复制代码

MeteoInfo

半颗麦芽糖 发表于 2016-11-5 10:43
王老师，根据您提供的参考网站我看懂了这个脚本并针对我的数据进行了修改，成功地出了图，非常感谢！我还 ...

用xlim函数。

半颗麦芽糖

MeteoInfo 发表于 2016-11-5 19:14
用xlim函数。

谢谢王老师的回复，终于把图画出来了，获益匪浅，以后学习meteoinfo的过程中恳请老师多多指教