WRF-plots网格框中的“焦点”是什么?(python代码)

Question

有人使用过来自的python代码吗此处(作者:盖尔·阿恩·瓦格博)从WRF输出中绘制skewT图?

我有一个关于pywrfplotUtils.py中代码中使用的一些函数的问题-这是getDimensions(nc)函数中调用的getXY(lon,lat)函数。

代码如下:

def getDimensions(nc): Nx = nc.getncattr('WEST-EAST_GRID_DIMENSION')-1 Ny = nc.getncattr(' south - north_grid_dimension ')-1 Nz = nc.getncattr('BOTTOM-TOP_GRID_DIMENSION')-1 dx = nc.getncattr(' dx ') dy = nc.getncattr(' dy ') lons = nc.getncattr。变量['XLAT'][0] lats = nc。变量['XLONG'][0] #查找焦点x的坐标，y = getXY(lons[Ny/2，:]，lats[:，Nx/2]) return Nx,Ny,Nz,lons,lats,dx,dy,x,y def getXY(lon,lat): x_nr = 0 y_nr = 0 while (lon[x_nr] < lon_focuspoint): x_nr += 1 while (lat[y_nr] < lat_focuspoint): y_nr += 1 print "x_nr:" str(x_nr)，" lat:"，str(lon[x_nr]) print "y_nr:" str(y_nr)，" lat:"，str(lon[x_nr]) return x_nr,y_nr

焦点是pywrfplotParams.py中设置的一个变量:

#这定义了网格框，其中skewT-plots等基于lat_focuspoint = 60.2 lon_focuspoint = 11.08

不幸的是，我不知道焦点指的是什么，所以当我使用作者设置的值时，显然我会得到一个错误，这个错误是-IndexError:索引越界。

我发现了这个链接在为该问题中问到的属性提供定义时很有用，但我不完全确定是否也应该从WRF模型网格中的纬度和经度获得单元格索引。

非常感谢你的想法!

Answer 1

Skew-T/Log-p图是从一个点开始的垂直剖面，在该代码中，您要绘制Skew-T的点是纬度60.2，经度11.08。代码片段:

当(lon[x_nr] < lon_focuspoint): x_nr += 1 while (lat[y_nr] < lat_focuspoint): y_nr += 1

通过将x和y网格点设置为0，并在该网格点上的纬度/经度值小于焦点的情况下递增索引，找到在焦点以西和以南最近的网格点。这就是焦点的作用——您想要skew-t/log-p绘图的仰角/仰角位置。

一个例外IndexError:索引越界。意思是上面的代码试图访问一个值经度(x_nr)或lat [y_nr]这是越界。这意味着您的焦点的纬度/经度坐标不在您的域中。

要成功地运行代码，请将焦点编辑为模型域中的一个点，特别是希望使用skew-t概要文件的纬度和经度的值。