什么是“焦点”在网格WRF-plots箱子吗?(python代码)

Question

有人使用的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模式网格。

你的想法是非常感谢!

Answer 1

Skew-T / logp情节是垂直廓线从一个点,它出现在你们的重点策划Skew-T的代码是60.2纬度,经度11.08。的代码片段:

def getXY(经度、纬度):x_nr = 0 y_nr = 0时(朗[x_nr] < lon_focuspoint): x_nr + = 1时(lat [y_nr] < lat_focuspoint): y_nr + = 1

发现最近的网格点西部和南部的焦点由x和y现设置为0和递增索引只要网格点的纬度/经度值小于焦点。这是焦点是什么,你想要的纬度/经度位置skew-t / logp情节。

一个例外IndexError:索引越界。意味着上述代码试图访问一个值经度(x_nr)或lat [y_nr]这是禁止入内的。这意味着你的注意力的纬度/经度坐标点并不是在你的域。

成功运行代码,编辑焦点是点在你的域模型,特别是经度和纬度的值,你想skew-t概要文件。