MeteoInfoLab脚本示例：micaps4类数据

MeteoInfo

将两个micaps 4类数据读为2个二维数组（addfile_micaps()函数），相减获得24小时降水量，然后绘图。

1. fn1 = 'D:/Temp/micaps/14050120.024'
   
2. fn2 = 'D:/Temp/micaps/14050120.048'
   
3. f1 = addfile_micaps(fn1)
   
4. f2 = addfile_micaps(fn2)
   
5. t1 = f1.gettime(0)
   
6. t2 = f2.gettime(0)
   
7. data1 = f1['var'][0,0,:,:]
   
8. data2 = f2['var'][0,0,:,:]
   
9. data = data2 - data1
   
10. #Plot
    
11. axesm()
    
12. geoshow('cn_province')
    
13. geoshow('country', edgecolor='k')
    
14. levs = [0.1,2,5,10,20,30,50,80,100,120,150]
    
15. cols = makecolors(len(levs)+1, cmap='MPL_rainbow')
    
16. cols[0] = 'w'
    
17. #layer = contourfm(data, levs, colors=cols)
    
18. layer = imshowm(data, levs, colors=cols)
    
19. colorbar(layer, label='mm', extendrect=False)
    
20. yticks(arange(10, 61, 10))
    
21. title(u'24小时降水（' + t1.strftime('%Y-%m-%d %H:00') + ' - ' +\
    
22.     t2.strftime('%Y-%m-%d %H:00') + u'）', fontname=u'黑体')

fn1 = 'D:/Temp/micaps/14050120.024'<br /> fn2 = 'D:/Temp/micaps/14050120.048'<br /> f1 = addfile_micaps(fn1)<br /> f2 = addfile_micaps(fn2)<br /> t1 = f1.gettime(0)<br /> t2 = f2.gettime(0)<br /> data1 = f1['var'][0,0,:,:]<br /> data2 = f2['var'][0,0,:,:]<br /> data = data2 - data1<br /> #Plot<br /> axesm()<br /> geoshow('cn_province')<br /> geoshow('country', edgecolor='k')<br /> levs = [0.1,2,5,10,20,30,50,80,100,120,150]<br /> cols = makecolors(len(levs)+1, cmap='MPL_rainbow')<br /> cols[0] = 'w'<br /> #layer = contourfm(data, levs, colors=cols)<br /> layer = imshowm(data, levs, colors=cols)<br /> colorbar(layer, label='mm', extendrect=False)<br /> yticks(arange(10, 61, 10))<br /> title(u'24小时降水（' + t1.strftime('%Y-%m-%d %H:00') + ' - ' +\<br />     t2.strftime('%Y-%m-%d %H:00') + u'）', fontname=u'黑体')

MeteoInfo

风格而才 发表于 2018-1-18 15:08
王老师不好意思，又来麻烦您了！之前提的问题不知道怎么回事没了。 我参考您的这个帖子想利用Micamps第四类 ...

专门修正了这个bug，你下载最新版本的MeteoInfo试试。

好多鱼

很好的帖子，留下来了。

liutong

请教老师，为什么会显示函数名错误
Traceback (most recent call last):
  File "<string>", line 3, in <module>
NameError: name 'addfile_micaps' is not defined

MeteoInfo

liutong 发表于 2016-8-26 16:42
请教老师，为什么会显示函数名错误
Traceback (most recent call last):
  File "", line 3, in

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

liutong

MeteoInfo 发表于 2016-8-26 17:02
参考此贴：http://bbs.06climate.com/forum.php?mod=viewthread&tid=36151&extra=page%3D1

老师，还是不太明白，adddile_micaps（）不就是函数吗？为什么会显示不存在？

MeteoInfo

liutong 发表于 2016-8-30 23:17
老师，还是不太明白，adddile_micaps（）不就是函数吗？为什么会显示不存在？

你是怎么运行脚本的？需要在MeteoInfoLab中运行。

liutong

MeteoInfo 发表于 2016-8-31 08:33
你是怎么运行脚本的？需要在MeteoInfoLab中运行。

[img]file:///C:/Users/Administrator/AppData/Roaming/Tencent/Users/862048002/QQ/WinTemp/RichOle/U6_S%603Z_(YD$Q%7B9G$U%7D8@B3.png[/img] 就是在meteinfolab中运行的，但是就是有错的

liutong

另外python零基础，非常想学习meteinfo，但知道要先熟悉python这门语言，所以老师能不能给推荐本学习python的书？

MeteoInfo

liutong 发表于 2016-8-31 23:10
就是在meteinfolab中运行的，但是就是有错的

这是C#版的MeteoInfo，不是MeteoInfoLab。建议你下载MeteoInfo Java版，MeteoInfoLab在里面。

MeteoInfo

liutong 发表于 2016-8-31 23:47
另外python零基础，非常想学习meteinfo，但知道要先熟悉python这门语言，所以老师能不能给推荐本学习python ...

这个博客写得不错，可以用于Python入门。最重要还是要多练习。