如何通过wrfout文件计算不同边界层高度内的平均风速

AllenLu

向大家请教，如何通过wrfout文件计算不同边界层高度内的平均风速，或者说，怎样计算不同高度的风速，有大神知道嘛，望不吝赐教[url=]可怜[/url]

AllenLu

或者说用ncl插值怎么弄
WRF-Python也行

AllenLu

在官网找到这个，但是运行出错
Interpolate Relative Humidity to Boundary Layer Heights

1. from netCDF4 import Dataset
   
2. from wrf import getvar, interplevel
   
3. 
   
4. wrfin = Dataset("wrfout_d02_2010-06-13_21:00:00")
   
5. 
   
6. rh = getvar(wrfin, "rh")
   
7. height = getvar(wrfin, "height_agl")
   
8. pblh = getvar(wrfin, "PBLH")
   
9. 
   
10. rh_pblh = interplevel(rh, height, pblh)

from netCDF4 import Dataset<br /> from wrf import getvar, interplevel<br /> <br /> wrfin = Dataset("wrfout_d02_2010-06-13_21:00:00")<br /> <br /> rh = getvar(wrfin, "rh")<br /> height = getvar(wrfin, "height_agl")<br /> pblh = getvar(wrfin, "PBLH")<br /> <br /> rh_pblh = interplevel(rh, height, pblh)

错误：
Traceback (most recent call last):
  File "222.py", line 7, in <module>
    height = getvar(wrfin, "height_agl")
  File "/python3.7/site-packages/wrf/routines.py", line 343, in getvar
    raise ValueError("'%s' is not a valid variable name" % (varname))
ValueError: 'height_agl' is not a valid variable name
大家知道为什么嘛

AllenLu

AllenLu 发表于 2021-7-21 21:24
在官网找到这个，但是运行出错
Interpolate Relative Humidity to Boundary Layer Heights

继续自问自答：我用z替换了height_agl，相对湿度rh可以出来结果，但是风速要怎么弄，把rh替换为U、V好像不太行哎会报错：[code:python]  raise ValueError("arguments 0 and 1 must have the same shape")[/code]
而且要怎么把插值结果输出来呀，求解答

AllenLu

AllenLu 发表于 2021-7-22 09:54
继续自问自答：我用z替换了height_agl，相对湿度rh可以出来结果，但是风速要怎么弄，把rh替换为U、V好像 ...

不要替换，是因为我版本太低，才没有height_agl变量，更新到1.3.2就有了

johnLeee

AllenLu 发表于 2021-7-24 11:04
不要替换，是因为我版本太低，才没有height_agl变量，更新到1.3.2就有了

您好。请问您知道如何算了嘛，谢谢

AllenLu

johnLeee 发表于 2023-11-29 16:28
您好。请问您知道如何算了嘛，谢谢

这个好久了，现在毕业了，之前写的东西找不到了