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

MeteoInfo

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. 
   
8. #Projection
   
9. proj = projinfo(proj='lcc', lon_0=-100, lat_0=45, lat_1=33, lat_2=45)
   
10. 
    
11. # Set the forecast hour
    
12. FH = 1
    
13. 
    
14. #Interpolate to a z level
    
15. p = [70000., 80000.]
    
16. height, temp = meteo.log_interpolate_1d(p, pres, hgt, temperature, axis=1)
    
17. 
    
18. #Plot
    
19. axesm(projinfo=proj, griddx=5, griddy=5)
    
20. geoshow('D:/Temp/map/states.shp', edgecolor='b')
    
21. geoshow('country', edgecolor=(0,0,255))
    
22. levs = arange(-20, 20, 1)
    
23. cols = makecolors(len(levs)+1, 'MPL_RdBu', reverse=True)
    
24. layer = contourfm(lon, lat, temp[FH,0,:,:], levs, colors=cols, proj=f.proj)
    
25. colorbar(layer, label='Celsius')
    
26. cs = contourm(lon, lat, height[FH,0,:,:], arange(2880,4000,60), colors='k', proj=f.proj)
    
27. clabel(cs, fontsize=10, drawshadow=False)
    
28. t = f.gettime(FH)
    
29. title('WRF-ARW Forecast VALID: %s UTC' % t.strftime('%Y-%m-%d %H:00') + \
    
30.     '\n700 hPa Heights (m) and Temperature (C)')

f = addfile('D:/Temp/nc/wrf_example.nc')<br /> lat = f['lat'][:,0]<br /> lon = f['lon'][0,:]<br /> temperature = f['temperature'][:,:,:,:]<br /> pres = f['pressure'][:,:,:,:]<br /> hgt = f['height'][:,:,:,:]<br /> <br /> #Projection<br /> proj = projinfo(proj='lcc', lon_0=-100, lat_0=45, lat_1=33, lat_2=45)<br /> <br /> # Set the forecast hour<br /> FH = 1<br /> <br /> #Interpolate to a z level<br /> p = [70000., 80000.]<br /> height, temp = meteo.log_interpolate_1d(p, pres, hgt, temperature, axis=1)<br /> <br /> #Plot<br /> axesm(projinfo=proj, griddx=5, griddy=5)<br /> geoshow('D:/Temp/map/states.shp', edgecolor='b')<br /> geoshow('country', edgecolor=(0,0,255))<br /> levs = arange(-20, 20, 1)<br /> cols = makecolors(len(levs)+1, 'MPL_RdBu', reverse=True)<br /> layer = contourfm(lon, lat, temp[FH,0,:,:], levs, colors=cols, proj=f.proj)<br /> colorbar(layer, label='Celsius')<br /> cs = contourm(lon, lat, height[FH,0,:,:], arange(2880,4000,60), colors='k', proj=f.proj)<br /> clabel(cs, fontsize=10, drawshadow=False)<br /> t = f.gettime(FH)<br /> title('WRF-ARW Forecast VALID: %s UTC' % t.strftime('%Y-%m-%d %H:00') + \<br />     '\n700 hPa Heights (m) and Temperature (C)')

随缘

{:eb502:}{:eb502:}{:eb502:}

bianpuxi

运行出错，好像是所用版本没有这个功能？

>>> run script...
Traceback (most recent call last):
  File "<iostream>", line 16, in <module>
AttributeError: 'module' object has no attribute 'log_interpolate_1d'

bianpuxi

bianpuxi 发表于 2019-12-4 15:55
运行出错，好像是所用版本没有这个功能？

>>> run script...

下载了最新版本的软件，成功运行了

随缘

王老师，MeteoInfoLab支持在脚本里参数调用运行吗？

MeteoInfo

随缘 发表于 2019-12-6 08:45
王老师，MeteoInfoLab支持在脚本里参数调用运行吗？

没看明白你的问题。

随缘

MeteoInfo 发表于 2019-12-6 16:58
没看明白你的问题。

不好意思，我的意思是：MeteoInfoLab是可以脚本运行的吧？就是我要自动批处理数据绘图，在cmd或是什么的可以运行脚本自动运行

MeteoInfo

随缘 发表于 2019-12-6 17:07
不好意思，我的意思是：MeteoInfoLab是可以脚本运行的吧？就是我要自动批处理数据绘图，在cmd或是什么的 ...

参考这里：http://www.meteothink.org/docs/i ... nfomap-meteoinfolab

随缘

MeteoInfo 发表于 2019-12-6 18:13
参考这里：http://www.meteothink.org/docs/introduction.html#running-meteoinfomap-meteoinfolab

多谢王老师

bianpuxi

王老师，程序虽然运行成功了，但是转换后的数据都是NaN是为什么呢