micaps数据作图

lyh0205

最近想用meteoinfo画涡度散度等物理量图，本人用的是micaps资料，在吧里搜索发现很多帖子都是一些脚本的例子，很多都搞不懂，想请问各位meteoinfo能不能直接读取micaps数据计算涡度散度并绘图呢？

MeteoInfo

参考这个脚本：

1. fn = 'D:/Temp/micaps/m_11_12080920.000'
   
2. f = addfile_micaps(fn)
   
3. u = f['U'][0,0,:,:]
   
4. v = f['V'][0,0,:,:]
   
5. vort = hcurl(u, v)
   
6. divg = hdivg(u, v)
   
7. 
   
8. #Plot
   
9. subplot(2,1,1,axestype='map')
   
10. geoshow('country')
    
11. layer = contourfm(vort, 20, smooth=False)
    
12. title('Vorticity')
    
13. 
    
14. subplot(2,1,2,axestype='map')
    
15. geoshow('country')
    
16. layer = contourfm(divg, 20, smooth=False)
    
17. title('Divergence')

fn = 'D:/Temp/micaps/m_11_12080920.000'<br /> f = addfile_micaps(fn)<br /> u = f['U'][0,0,:,:]<br /> v = f['V'][0,0,:,:]<br /> vort = hcurl(u, v)<br /> divg = hdivg(u, v)<br /> <br /> #Plot<br /> subplot(2,1,1,axestype='map')<br /> geoshow('country')<br /> layer = contourfm(vort, 20, smooth=False)<br /> title('Vorticity')<br /> <br /> subplot(2,1,2,axestype='map')<br /> geoshow('country')<br /> layer = contourfm(divg, 20, smooth=False)<br /> title('Divergence')

lyh0205

MeteoInfo 发表于 2020-3-16 16:57
参考这个脚本：

            谢谢老师！！！！