wrf.interplevel and wrf.vinterp results in very high values for vertical temperature gradients
还没有人认领这个 Issue。
评估
调研方向
首先复现 issue 中报告的 interplevel 和 vinterp 计算，然后阅读链接的 wrf.vinterp 文档，并比较高度和温度输入。完成的标准是确定高梯度是否符合预期，或确定一个具体的插值缺陷以及可复现的修正路径。
由索引模型根据 Issue 内容生成。
描述
Hello, I am trying to calculate the vertical temperature gradient between the surface and 1 km AGL.
I did the following:
temp2m= wrf.getvar(wrfnc,'T2') #wrfnc is the WRF output file for a specific hour
tk= wrf.getvar(wrfnc,'tk') # Temperature in Kelvin
height_agl= wrf.getvar(wrfnc,'height_agl') # Model Height for Mass Grid (AGL) in m
temp_1km = wrf.interplevel(tk, height_agl, 1000) # tempertaure interpolated to 1000 m above AGL
Then I caculate dthe temperature gradient between surace and 1 km AGL as follows:
del_T1= temp_1km-temp2m # kelvin
del_h1= 998 #1000 m- 2 m
Temp_gradient= -(del_T1/del_h1)*1000 # -dT/dz (K/km)
I am getting very high values for temperature gradient between surface and 1 km above ground level, and I am not sure if they are meteorologically correct. I have attached a screenshot of the values.
I also tried using wrf.vintep as follows:
temp_1km = wrf.vinterp(wrfnc, field=tk, vert_coord="ght_agl", interp_levels=[1], extrapolate=True,field_type="tk")
The document (https://wrf-python.readthedocs.io/en/main/user_api/generated/wrf.vinterp.html) says ’ght_agl’= grid point height agl [km], so I set the interp_levels as 1 km. But I similar results with slight difference.
- 主要语言
- Python
- 星标
- 498
- 派生
- 178
- PR 合并指标
- 30 天内没有已合并 PR
贡献指南
从这里开始
- 先读完整个 Issue，再读项目的贡献指南。
- 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
- Fork 仓库，在一个分支上完成修改。
- 提交 Pull Request，并在描述里引用这个 Issue 编号。
NCAR/wrf-python 的其他 Issue
-
NCAR/wrf-python#309 · 已指派 1 人 ·
-
难度 4/5 3-5 天 新手友好度 38/100
NCAR/wrf-python#307 · 5 条评论 ·
-
难度 4/5 3-5 天 新手友好度 35/100
NCAR/wrf-python#286 · 3 条评论 ·
-
难度 5/5 一周以上 新手友好度 25/100
NCAR/wrf-python#282 · 1 条评论 ·
-
PyPI releases 未关闭
NCAR/wrf-python#274 · 5 条评论 · 1 个 reaction · 已指派 1 人 ·
相似的 Issue
-
难度 1/5 1 小时以内 新手友好度 75/100
-
hcocena 未关闭policies-accepted pre-review precheck-passed
难度 1/5 1 小时以内 新手友好度 88/100
Bioconductor/BiocContributions#214 · 5 条评论 ·
-
难度 1/5 1 小时以内 新手友好度 92/100
TencentCloud/Octop#1169 · 1 条评论 ·
-
难度 2/5 1-3 小时 新手友好度 70/100
521xueweihan/HelloGitHub#3778 ·
-
The version checker's trailing attribute region has no control for a less-than inside a quoted value 未关闭area: dashboard area: tests bug perceived difficulty: 2 python
难度 2/5 1-3 小时 新手友好度 84/100
Nitjsefnie-Harness-Commons/daedalus#1105 · 1 条评论 ·