请选择 进入手机版 | 继续访问电脑版
爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

新浪微博登陆

只需一步, 快速开始

QQ登录

只需一步,快速开始

搜索
查看: 437|回复: 9

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

[复制链接] |关注本帖

新浪微博达人勋

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

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


hy_conc_nuclear.png

密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-3-8 13:17:55 | 显示全部楼层 |取消关注该作者的回复
王老师,威武。我辈又可以好好学习了。
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-3-9 09:38:37 | 显示全部楼层 |取消关注该作者的回复
王老师威武!最近在学习hysplit浓度扩散模块 研究中!
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-3-9 10:11:35 | 显示全部楼层 |取消关注该作者的回复
王老师您好,我想问一下您的第十四行是什么意思啊?两个值是如何定的?还有我出现以下错误是什么意思呢
>>> 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'
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-3-9 12:38:29 | 显示全部楼层 |取消关注该作者的回复
lulu 发表于 2018-3-9 10:11
王老师您好,我想问一下您的第十四行是什么意思啊?两个值是如何定的?还有我出现以下错误是什么意思呢
>> ...

14行是福岛的经纬度,

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

可以加入MeteoInfo QQ群更方便交流
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-3-9 14:47:18 | 显示全部楼层 |取消关注该作者的回复
MeteoInfo 发表于 2018-3-9 12:38
14行是福岛的经纬度,

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

我以为13行是经纬度王老师,那13行是什么呢
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-3-9 15:29:25 | 显示全部楼层 |取消关注该作者的回复
lulu 发表于 2018-3-9 14:47
我以为13行是经纬度王老师,那13行是什么呢

也是经纬度,标注城市名的时候要有移动,避免压盖。
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-3-12 14:03:52 | 显示全部楼层 |取消关注该作者的回复
MeteoInfo 发表于 2018-3-9 15:29
也是经纬度,标注城市名的时候要有移动,避免压盖。

原来是这样啊,谢谢王老师
密码修改失败请联系qq:937062711

新浪微博达人勋

发表于 2018-3-14 16:20:28 | 显示全部楼层 |取消关注该作者的回复
王老師,想請問您第一行要讀取的資料(Cdump15),是使用HYSPLIT模式中Concentraction選項所執行出來的OUTPUT檔案嗎????
密码修改失败请联系qq:937062711

新浪微博达人勋

 楼主| 发表于 2018-3-14 16:48:51 | 显示全部楼层 |取消关注该作者的回复
KIMO23 发表于 2018-3-14 16:20
王老師,想請問您第一行要讀取的資料(Cdump15),是使用HYSPLIT模式中Concentraction選項所執行出來的OUTPUT ...

是的
密码修改失败请联系qq:937062711
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表