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

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 4702|回复: 9

MeteoInfoLab脚本示例:Sigma坐标插值为等压面坐标

[复制链接]

新浪微博达人勋

发表于 2019-5-16 15:26:33 | 显示全部楼层 |阅读模式

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

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

x
WRF等模式的输出数据垂直方向是Sigma坐标,可以用meteolib包里的log_interpolate_1d函数将数据转为等压面垂直坐标。

  1. f = addfile('D:/Temp/nc/wrf_example.nc')
  2. lat = f['lat'][:,0]
  3. lon = f['lon'][0,:]
  4. temperature = f['temperature'][:,:,:,:]
  5. pres = f['pressure'][:,:,:,:]
  6. hgt = f['height'][:,:,:,:]

  7. #Projection
  8. proj = projinfo(proj='lcc', lon_0=-100, lat_0=45, lat_1=33, lat_2=45)

  9. # Set the forecast hour
  10. FH = 1

  11. #Interpolate to a z level
  12. p = [70000., 80000.]
  13. height, temp = meteo.log_interpolate_1d(p, pres, hgt, temperature, axis=1)

  14. #Plot
  15. axesm(projinfo=proj, griddx=5, griddy=5)
  16. geoshow('D:/Temp/map/states.shp', edgecolor='b')
  17. geoshow('country', edgecolor=(0,0,255))
  18. levs = arange(-20, 20, 1)
  19. cols = makecolors(len(levs)+1, 'MPL_RdBu', reverse=True)
  20. layer = contourfm(lon, lat, temp[FH,0,:,:], levs, colors=cols, proj=f.proj)
  21. colorbar(layer, label='Celsius')
  22. cs = contourm(lon, lat, height[FH,0,:,:], arange(2880,4000,60), colors='k', proj=f.proj)
  23. clabel(cs, fontsize=10, drawshadow=False)
  24. t = f.gettime(FH)
  25. title('WRF-ARW Forecast VALID: %s UTC' % t.strftime('%Y-%m-%d %H:00') + \
  26.     '\n700 hPa Heights (m) and Temperature (C)')


sigma2pres_1.png
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2019-5-16 15:46:53 | 显示全部楼层
{:eb502:}{:eb502:}{:eb502:}
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-12-4 15:55:30 | 显示全部楼层
运行出错,好像是所用版本没有这个功能?

>>> run script...
Traceback (most recent call last):
  File "<iostream>", line 16, in <module>
AttributeError: 'module' object has no attribute 'log_interpolate_1d'
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-12-5 17:39:28 | 显示全部楼层
bianpuxi 发表于 2019-12-4 15:55
运行出错,好像是所用版本没有这个功能?

>>> run script...

下载了最新版本的软件,成功运行了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-12-6 08:45:48 | 显示全部楼层
王老师,MeteoInfoLab支持在脚本里参数调用运行吗?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-12-6 16:58:12 | 显示全部楼层
随缘 发表于 2019-12-6 08:45
王老师,MeteoInfoLab支持在脚本里参数调用运行吗?

没看明白你的问题。
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-12-6 17:07:56 | 显示全部楼层
MeteoInfo 发表于 2019-12-6 16:58
没看明白你的问题。

不好意思,我的意思是:MeteoInfoLab是可以脚本运行的吧?就是我要自动批处理数据绘图,在cmd或是什么的可以运行脚本自动运行
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

 楼主| 发表于 2019-12-6 18:13:23 | 显示全部楼层
随缘 发表于 2019-12-6 17:07
不好意思,我的意思是:MeteoInfoLab是可以脚本运行的吧?就是我要自动批处理数据绘图,在cmd或是什么的 ...

参考这里:http://www.meteothink.org/docs/i ... nfomap-meteoinfolab
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-12-6 21:03:09 | 显示全部楼层
MeteoInfo 发表于 2019-12-6 18:13
参考这里:http://www.meteothink.org/docs/introduction.html#running-meteoinfomap-meteoinfolab

多谢王老师
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

新浪微博达人勋

发表于 2019-12-10 16:42:09 | 显示全部楼层
王老师,程序虽然运行成功了,但是转换后的数据都是NaN是为什么呢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

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

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

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