爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 15497|回复: 9

wrf-python提取WRF行星边界层高度PBLH时出错

[复制链接]
发表于 2021-7-23 17:53:11 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

x
  1. <span class="kn" style="box-sizing: border-box; color: rgb(0, 112, 32); font-weight: bold;">from</span> <span class="nn" style="box-sizing: border-box; color: rgb(14, 132, 181); font-weight: bold;">netCDF4</span> <span class="kn" style="box-sizing: border-box; color: rgb(0, 112, 32); font-weight: bold;">import</span> <span class="n" style="box-sizing: border-box;">Dataset</span>
  2. <span class="kn" style="box-sizing: border-box; color: rgb(0, 112, 32); font-weight: bold;">from</span> <span class="nn" style="box-sizing: border-box; color: rgb(14, 132, 181); font-weight: bold;">wrf</span> <span class="kn" style="box-sizing: border-box; color: rgb(0, 112, 32); font-weight: bold;">import</span> <span class="n" style="box-sizing: border-box;">getvar</span><span class="p" style="box-sizing: border-box;">,</span> <span class="n" style="box-sizing: border-box;">interplevel</span>

  3. <span class="n" style="box-sizing: border-box;">wrfin</span> <span class="o" style="box-sizing: border-box; color: rgb(102, 102, 102);">=</span> <span class="n" style="box-sizing: border-box;">Dataset</span><span class="p" style="box-sizing: border-box;">(</span><span class="s2" style="box-sizing: border-box; color: rgb(64, 112, 160);">"wrfout_d02_2010-06-13_21:00:00"</span><span class="p" style="box-sizing: border-box;">)</span>

  4. <span class="n" style="box-sizing: border-box;">pblh</span> <span class="o" style="box-sizing: border-box; color: rgb(102, 102, 102);">=</span> <span class="n" style="box-sizing: border-box;">getvar</span><span class="p" style="box-sizing: border-box;">(</span><span class="n" style="box-sizing: border-box;">wrfin</span><span class="p" style="box-sizing: border-box;">,</span> <span class="s2" style="box-sizing: border-box; color: rgb(64, 112, 160);">"PBLH"</span><span class="p" style="box-sizing: border-box;">)</span>

复制代码
输出结果全是0,大家知道为啥嘛,感谢
  1. <xarray.DataArray 'PBLH' (south_north: 21, west_east: 18)>
  2. array([[0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  3.         0., 0.],
  4.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  5.         0., 0.],
  6.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  7.         0., 0.],
  8.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  9.         0., 0.],
  10.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  11.         0., 0.],
  12.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  13.         0., 0.],
  14.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  15.         0., 0.],
  16.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  17.         0., 0.],
  18.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  19.         0., 0.],
  20.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  21.         0., 0.],
  22. ...
  23.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  24.         0., 0.],
  25.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  26.         0., 0.],
  27.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  28.         0., 0.],
  29.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  30.         0., 0.],
  31.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  32.         0., 0.],
  33.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  34.         0., 0.],
  35.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  36.         0., 0.],
  37.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  38.         0., 0.],
  39.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  40.         0., 0.],
  41.        [0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., 0.,
  42.         0., 0.]], dtype=float32)
  43. Coordinates:
  44.     XLONG    (south_north, west_east) float32 111.0 111.0 111.0 ... 111.6 111.6
  45.     XLAT     (south_north, west_east) float32 40.09 40.09 40.09 ... 40.63 40.63
  46.     XTIME    float32 0.0
  47.     Time     datetime64[ns] 2019-10-01
  48. Dimensions without coordinates: south_north, west_east
  49. Attributes:
  50.     FieldType:    104
  51.     MemoryOrder:  XY
  52.     description:  PBL HEIGHT
  53.     units:        m
  54.     stagger:      
  55.     coordinates:  XLONG XLAT XTIME
  56.     projection:   LambertConformal(stand_lon=108.09700012207031, moad_cen_lat...
复制代码



密码修改失败请联系微信:mofangbao
 楼主| 发表于 2021-7-24 11:02:15 | 显示全部楼层
结贴,是只提取了第一个小时的,改timeidx=就行,可是如果我想全读怎么办呢?
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2021-8-5 15:54:43 | 显示全部楼层
AllenLu 发表于 2021-7-24 11:02
结贴,是只提取了第一个小时的,改timeidx=就行,可是如果我想全读怎么办呢?

在哪里改呀
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2021-8-6 11:08:58 | 显示全部楼层

类似height = getvar(wrfin, "height_agl",timeidx=1)
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2021-8-6 12:48:16 | 显示全部楼层
AllenLu 发表于 2021-8-6 11:08
类似height = getvar(wrfin, "height_agl",timeidx=1)

能问问是哪个具体的文件里面吗,新手
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2021-8-9 22:40:13 | 显示全部楼层
你好我想问一下wrf-python能不能安装到Windows上呀,怎么安装,感谢
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2021-8-12 20:04:41 | 显示全部楼层
坨坨and便便 发表于 2021-8-9 22:40
你好我想问一下wrf-python能不能安装到Windows上呀,怎么安装,感谢

好像不能,这是linux里的包,你可以试试子系统
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2021-8-12 20:07:17 | 显示全部楼层
小吴吴 发表于 2021-8-6 12:48
能问问是哪个具体的文件里面吗,新手

啊这,也太新了,比我还新
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

发表于 2021-8-13 14:26:45 | 显示全部楼层
AllenLu 发表于 2021-8-12 20:07
啊这,也太新了,比我还新

没有办法,网络上教程太少了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

 楼主| 发表于 2021-8-13 18:25:26 | 显示全部楼层
小吴吴 发表于 2021-8-13 14:26
没有办法,网络上教程太少了

你可以去wrf-python官方看看,教程很详细,琢磨琢磨就会了
密码修改失败请联系微信:mofangbao
回复 支持 反对

使用道具 举报

您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表