MeteoInfoLab脚本示例：wrfout数据

MeteoInfo

WRF输出的数据是netCDF格式，但是有它自己的约定（conventions）。MeteoInfo支持此conventions，因此读取wrfout数据并绘图比较简单。

示例脚本：

1. fn = 'D:/Temp/nc/wrfout_d02_2010-09-01_12_00_00'
   
2. f = addfile(fn)
   
3. data = f['T'][0,0,:,:]
   
4. #Plot
   
5. axesm(projinfo=f.proj, griddx=1, griddy=1, gridline=True)
   
6. lchina = shaperead('D:/Temp/map/bou2_4p.shp')
   
7. lhainan = shaperead('D:/Temp/map/hainan.shp')
   
8. geoshow(lhainan, edgecolor='gray')
   
9. geoshow(lchina)
   
10. layer = contourfm(data, 20, proj=f.proj)
    
11. colorbar(layer, label='K')
    
12. title('Perturbation potential temperature (theta-t0)')

fn = 'D:/Temp/nc/wrfout_d02_2010-09-01_12_00_00'<br /> f = addfile(fn)<br /> data = f['T'][0,0,:,:]<br /> #Plot<br /> axesm(projinfo=f.proj, griddx=1, griddy=1, gridline=True)<br /> lchina = shaperead('D:/Temp/map/bou2_4p.shp')<br /> lhainan = shaperead('D:/Temp/map/hainan.shp')<br /> geoshow(lhainan, edgecolor='gray')<br /> geoshow(lchina)<br /> layer = contourfm(data, 20, proj=f.proj)<br /> colorbar(layer, label='K')<br /> title('Perturbation potential temperature (theta-t0)')

nuist2006

不错啊这个

放逐流年

王老师，我看了一下午这个MeteoInfolab的代码以及数据处理，只发现这一个和wrfout有关代码，或者说画图的，正像您在帖子中说的“WRF输出的数据是netCDF格式，但是有它自己的约定（conventions）。MeteoInfo支持此conventions，因此读取wrfout数据并绘图比较简单。”这样，meteoinfolab支持wrfout的约定，那我们处理wrfout数据时应该怎么办呢？

MeteoInfo

放逐流年 发表于 2017-4-10 17:44
王老师，我看了一下午这个MeteoInfolab的代码以及数据处理，只发现这一个和wrfout有关代码，或者说画图的， ...

不明白你想问什么

vuchiyuan

非常好~

abc258886@

这个和ncl 有什么区别！！！

蹦蹦不是笨笨

。shp文件怎么来的呀

药草

wrfout的nc文件，数据格式是二维的经度、二维的维度、二维的数据，出图的时候，经纬度和数据无法关联，用python绘图使用contourf成功了，用meteoinfo中的scatterm也成功了，但是想在meteoinfo中画等值面用contourfm方法失败，没有设置经纬度的地方，是有其他的方法吗?

MeteoInfo

药草 发表于 2021-3-12 14:23
wrfout的nc文件，数据格式是二维的经度、二维的维度、二维的数据，出图的时候，经纬度和数据无法关联，用py ...

MeteoInfo会自动读取文件中的投影和坐标信息。如果绘制某些数据文件有问题可以把文件传上来，方便调试。

药草

MeteoInfo 发表于 2021-3-12 15:56
MeteoInfo会自动读取文件中的投影和坐标信息。如果绘制某些数据文件有问题可以把文件传上来，方便调试。

这个文件，您给看一下，我用python是可以的，现在该用python写了，但是还是想知道为什么meteoinfo不行，可能是我对meteoinfo了解的不够全面