MeteoInfoLab脚本示例：计算涡度、散度

MeteoInfo

用U/V分量数据计算涡度和散度，计算涡度的函数是hcurl，计算散度的函数是hdivg，参数都是U, V。

脚本程序：

1. f = addfile('D:/Temp/GrADS/model.ctl')
   
2. u = f['U'][0,0,:,:]
   
3. v = f['V'][0,0,:,:]
   
4. vort = hcurl(u, v)
   
5. divg = hdivg(u, v)
   
6. axesm()
   
7. mlayer = shaperead('D:/Temp/map/country1.shp')
   
8. geoshow(mlayer)
   
9. layer = contourfm(vort, 20)
   
10. #layer = contourfm(divg, 20)
    
11. title('Vorticity')
    
12. #title('Divergence')
    
13. colorbar(layer, orientation='horizontal', aspect=50)

f = addfile('D:/Temp/GrADS/model.ctl') <br /> u = f['U'][0,0,:,:] <br /> v = f['V'][0,0,:,:] <br /> vort = hcurl(u, v) <br /> divg = hdivg(u, v) <br /> axesm() <br /> mlayer = shaperead('D:/Temp/map/country1.shp') <br /> geoshow(mlayer) <br /> layer = contourfm(vort, 20) <br /> #layer = contourfm(divg, 20) <br /> title('Vorticity') <br /> #title('Divergence') <br /> colorbar(layer, orientation='horizontal', aspect=50)

Jason@Hello

版主,你好!我想问一下,你知道meteoinfo求散度的函数处理的是什么数据?我现在在知道他传值的时候是u,v都是二维数据,而里面的一维数组分别代表什么,我就不太清楚,如果版主知道的话.麻烦讲解一下,不胜感激!

MeteoInfo

Jason@Hello 发表于 2015-8-6 09:15
版主,你好!我想问一下,你知道meteoinfo求散度的函数处理的是什么数据?我现在在知道他传值的时候是u,v都是二 ...

U/V都是二维数组，没涉及到一维数组呀。

Jason@Hello

MeteoInfo 发表于 2015-8-6 09:49
U/V都是二维数组，没涉及到一维数组呀。

我的意思是二维数组里面的一维数组。我猜u/v是插值好之后的格点数据，然而在程序中它是利用二维数组的格式来表示这个网格的是吗？

MeteoInfo

Jason@Hello 发表于 2015-8-6 13:06
我的意思是二维数组里面的一维数组。我猜u/v是插值好之后的格点数据，然而在程序中它是利用二维数组的格 ...

没看明白到底想问什么？

棉花猪仔

MeteoInfo 发表于 2015-8-6 14:15
没看明白到底想问什么？

版主，请问   PyGridData   怎么转   DimArray？

MeteoInfo

棉花猪仔 发表于 2015-8-6 15:13
版主，请问   PyGridData   怎么转   DimArray？

看来你看了MeteoInfoLab封装的Jython脚本。没有做PyGridData转DimArray的功能，事实上通过DimVariable slice获取的数组就是DimArray，似乎没有需要从PyGridData转为DimArray。

棉花猪仔

MeteoInfo 发表于 2015-8-6 15:24
看来你看了MeteoInfoLab封装的Jython脚本。没有做PyGridData转DimArray的功能，事实上通过DimVariable sl ...

感谢版主热心解答，
可是还没解决问题，我另开一贴把情况说一下，麻烦版主帮帮忙

Jason@Hello

MeteoInfo 发表于 2015-8-6 14:15
没看明白到底想问什么？

可能是我提问方式不对！但是还是感谢版主热心回答！

无尽的兰博

老师您好，我用您的脚本用grib2的数据画了涡度图，但是我想选择200hPa那层的，我的设置如下：

f = addfile('F:/data/study/fnl/fnl_20160622_18_00.grib2')
u = f['u-component_of_wind_pressure_difference_layer'][0,200,:,:]
v = f['v-component_of_wind_pressure_difference_layer'][0,200,:,:]
vort = hcurl(u, v)
axesm()
mlayer = shaperead('F:/MeteoInfo/map/country1.shp')
geoshow(mlayer)
layer = contourfm(vort, 20)
title('Vorticity')
colorbar(layer, orientation='horizontal', aspect=50)

但是运行后报错如下：

Index out of range!
Index out of range!
Traceback (most recent call last):
  File "<iostream>", line 8, in <module>
  File "F:\MeteoInfo\pylib\mipylib\miplot.py", line 3984, in contourfm
    ls = __getlegendscheme(args, gdata.min(), gdata.max(), **kwargs)
AttributeError: 'NoneType' object has no attribute 'min'

还麻烦老师帮我看一下需要怎么修改。谢谢您了！