有人使用的python代码吗(盖尔·阿恩Waagbo写的)去把skewT WRF输出?

我有一个问题关于代码中使用的一些功能pywrfplotUtils中找到。py——这是getXY(经度、纬度)函数中被调用getDimensions (nc)的功能。

下面的代码:

def getDimensions (nc): Nx = nc.getncattr (WEST-EAST_GRID_DIMENSION) 1纽约= nc.getncattr (SOUTH-NORTH_GRID_DIMENTION) 1新西兰= nc.getncattr (BOTTOM-TOP_GRID_DIMENSION) 1 dx = nc.getncattr (dx) dy = nc.getncattr (dy)朗=数控。变量(' XLAT '][0]背阔肌=数控。变量(' XLONG '][0] #寻找焦点坐标x, y = getXY(经度(Ny / 2:),背阔肌(:,Nx / 2))返回Nx, Ny,新西兰朗,背阔肌,dx, dy, x, y def getXY(经度、纬度):x_nr = 0 y_nr = 0时(经度(x_nr) < lon_focuspoint): x_nr + = 1时(lat [y_nr] < lat_focuspoint): y_nr + = 1打印“x_nr:“str (x_nr),“朗:“str(朗[x_nr])打印“y_nr:“str (y_nr),“纬度:“str(经度[x_nr])返回x_nr y_nr

pywrfplotParams.py的焦点是一组变量:

#这定义的网格框skewT-plots等基础lat_focuspoint lon_focuspoint = 11.08 = 60.2

不幸的是,我不知道焦点指的是,当我使用作者设定的值,显然,我得到一个错误,错误是-IndexError:索引越界。

我发现这链接提供有用的定义属性问这个问题,但是我不完全确定如果我也应该得到细胞指数的纬度和经度WRF模式网格。

你的想法是非常感谢!