MeteoInfoLab脚本示例：HYSPLIT模式浓度数据

MeteoInfo

HYSPLIT模式在核应急扩散模拟中有广泛的应用，这里给出一个脚本绘制多年前用HYSPLIT模拟的福岛核电站泄露扩散情况的图形。

1. fn = r'D:\Working\MIScript\JapanNuclear\cdump15'
   
2. f = addfile_hyconc(fn)
   
3. tidx = 5
   
4. data = f['C137'][tidx,0,:,:]
   
5. 
   
6. #Plot
   
7. axesm(outerposition=[0,0,1,1])
   
8. geoshow('cn_province', edgecolor='b')
   
9. geoshow('country', edgecolor='b')
   
10. levs = [1e-14,1e-13,1e-12,1e-11]
    
11. cols = ['w',(0,255,255),'g','y','r']
    
12. layer = contourfm(data, levs, colors=cols, smooth=False)
    
13. scatterm(141.0, 37.5, size=10, marker='S', color='r')
    
14. text(142.5, 39, 'Fukushima', bold=True)
    
15. xlim(110, 200)
    
16. ylim(15, 70)
    
17. yticks(arange(20, 61, 20))
    
18. colorbar(layer, label='C137', labelloc='top', shrink=0.8)
    
19. t = f.gettime(tidx)
    
20. title('Simulation of cesium-137 emitted from Fukushima (%s)' % t.strftime('%Y-%m-%d'))

fn = r'D:\Working\MIScript\JapanNuclear\cdump15'<br /> f = addfile_hyconc(fn)<br /> tidx = 5<br /> data = f['C137'][tidx,0,:,:]<br /> <br /> #Plot<br /> axesm(outerposition=[0,0,1,1])<br /> geoshow('cn_province', edgecolor='b')<br /> geoshow('country', edgecolor='b')<br /> levs = [1e-14,1e-13,1e-12,1e-11]<br /> cols = ['w',(0,255,255),'g','y','r']<br /> layer = contourfm(data, levs, colors=cols, smooth=False)<br /> scatterm(141.0, 37.5, size=10, marker='S', color='r')<br /> text(142.5, 39, 'Fukushima', bold=True)<br /> xlim(110, 200)<br /> ylim(15, 70)<br /> yticks(arange(20, 61, 20))<br /> colorbar(layer, label='C137', labelloc='top', shrink=0.8)<br /> t = f.gettime(tidx)<br /> title('Simulation of cesium-137 emitted from Fukushima (%s)' % t.strftime('%Y-%m-%d'))

独孤酒见

王老师，威武。我辈又可以好好学习了。

lulu

王老师威武！最近在学习hysplit浓度扩散模块 研究中！

lulu

王老师您好，我想问一下您的第十四行是什么意思啊？两个值是如何定的？还有我出现以下错误是什么意思呢
>>> run script...
Index out of range!
Traceback (most recent call last):
  File "<iostream>", line 12, in <module>
  File "D:\MeteoInfo\pylib\mipylib\plotlib\miplot.py", line 4804, in contourfm
    ls = plotutil.getlegendscheme(args, gdata.min(), gdata.max(), **kwargs)
AttributeError: 'NoneType' object has no attribute 'min'

MeteoInfo

lulu 发表于 2018-3-9 10:11
王老师您好，我想问一下您的第十四行是什么意思啊？两个值是如何定的？还有我出现以下错误是什么意思呢
>> ...

14行是福岛的经纬度，

仅凭一点错误信息很难看出什么。

可以加入MeteoInfo QQ群更方便交流

lulu

MeteoInfo 发表于 2018-3-9 12:38
14行是福岛的经纬度，

仅凭一点错误信息很难看出什么。

我以为13行是经纬度王老师，那13行是什么呢

MeteoInfo

lulu 发表于 2018-3-9 14:47
我以为13行是经纬度王老师，那13行是什么呢

也是经纬度，标注城市名的时候要有移动，避免压盖。

lulu

MeteoInfo 发表于 2018-3-9 15:29
也是经纬度，标注城市名的时候要有移动，避免压盖。

原来是这样啊，谢谢王老师

KIMO23

王老師，想請問您第一行要讀取的資料(Cdump15)，是使用HYSPLIT模式中Concentraction選項所執行出來的OUTPUT檔案嗎????

MeteoInfo

KIMO23 发表于 2018-3-14 16:20
王老師，想請問您第一行要讀取的資料(Cdump15)，是使用HYSPLIT模式中Concentraction選項所執行出來的OUTPUT ...

是的