5
\ begingroup美元

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

你的想法是非常感谢!

\ endgroup美元

    1回答1

    2
    \ begingroup美元

    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概要文件。

    \ endgroup美元
    5
    • \ begingroup美元 谢谢你的澄清@casey♦。我一直想获得经度和纬度的值的域自很久以前,但从未成功地做到这一点。你知道任何方法吗?然而,我所知道的位置数据,在明尼苏达州。如果我对你的答案是正确的理解,明尼苏达州的坐标是44.9 n和93.2 w。所以我试着使用它作为我的焦点,但我仍然得到IndexError。 \ endgroup美元
      - - - - - -cwmwl
      2015年6月23日14:47
    • \ begingroup美元 @cwmwl你可以问一个新的问题。情节域应该设置在你运行这个模型的所以你应该知道你的域的界限。内陆点,我会做一个阴谋与纬度/经度的轴(所以你知道在哪里你想要skew-t数据)。 \ endgroup美元
      - - - - - -凯西
      2015年6月23日,在52。丢下
    • \ begingroup美元 @cwmwl我假设您输入92.3 w-92.3吗? \ endgroup美元
      - - - - - -凯西
      2015年6月23日,53
    • \ begingroup美元 ♦好了,我会问一个新的问题。通知你,我没有运行模型。我只是使用的输出模型。不知道,我需要知道这些信息模型。 \ endgroup美元
      - - - - - -cwmwl
      2015年6月23日14:56
    • \ begingroup美元 @cwmwl好吧,如果你知道你想要的网格索引的情节,你可以修改程序直接使用它们,而不是计算出来的焦点。这是另一个选项如果你已经知道你想要的情节。 \ endgroup美元
      - - - - - -凯西
      2015年6月23日,在58。曾

    你的答案

    通过点击“发布你的答案”,你同意我们服务条款,隐私政策饼干的政策

    不是你要找的答案?浏览其他问题标记问你自己的问题