有人使用的python代码(这里(作者盖尔·阿恩Waagbo)][1]绘制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 = 60.2 lon_focuspoint = 11.08不幸的是,我不知道焦点指的是,当我使用作者设定的值,显然我得到一个错误,并且错误- * IndexError:索引越界。*我发现这[链接][2]有用提供定义的属性问这个问题,但是我不完全确定如果我也应该得到细胞指数的纬度和经度WRF模式网格。你的想法是非常感谢![1]:https://code.google.com/p/pywrfplot/source/browse/ svn % 2 ftrunk [2]: //www.hoelymoley.com/questions/2732/how-to-get-cell-indices-from-latitude-and-longitude-in-wrf-model-grids