wrf.interplevel returns NaN when interpolating at a level where the sample points are
还没有人认领这个 Issue。
评估
调研方向
首先运行提供的 NumPy 和 wrf.interplevel 示例，然后跟踪 interplevel 的实现及其对精确样本层的处理。将 0.5 处的结果与 yg[5,:,:] 以及 0.5+1e-8 的情况进行比较。当精确层上的插值不再产生 NaN，且现有行为对于附近层仍然正确时，即表示完成。
由索引模型根据 Issue 内容生成。
描述
Hi, I recently used wrf.interplevel to interpolate a 3d field onto a single pressure level 50000 Pa. I found a very wired result that every time when the 50000 Pa appears in the sample points, the result is NaN. If I added a little value to 50000Pa, it returns the correct value (approximately). Below is just a simple example to reproduce the issue. Please help!
# construct vert array
xx = np.linspace(0,1,11)
xg = np.zeros((11,4,3))
for j in np.arange(4):
for i in np.arange(3):
xg[:,j,i] = xx
# assign field3d array
yg = np.random.rand(11,4,3)
# print results for comparison
print(' True answer: ', yg[5,:,:])
print(' NaN result: ', wrf.interplevel(yg,xg,0.5).to_numpy())
print('Corrected result: ', wrf.interplevel(yg,xg,0.5+1e-8).to_numpy())
- 主要语言
- Python
- 星标
- 498
- 派生
- 178
- PR 合并指标
- 30 天内没有已合并 PR
环境准备
从这里开始
- 先读完整个 Issue，再读项目的贡献指南。
- 在 Issue 下留言说明你要接手 —— 这能避免两个人做同样的事。
- Fork 仓库，在一个分支上完成修改。
- 提交 Pull Request，并在描述里引用这个 Issue 编号。
NCAR/wrf-python 的其他 Issue
-
Add netCDF4 dependency and/or make it truly optional可能重新可做 @kafitzgerald 于 221 天前认领，目前没有进行中的 PR。 未关闭
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可能重新可做 @kafitzgerald 于 437 天前认领，目前没有进行中的 PR。 未关闭
NCAR/wrf-python#274 · 5 条评论 · 1 个 reaction · 已指派 1 人 ·
相似的 Issue
-
难度 2/5 1-3 小时 新手友好度 76/100
PedestrianDynamics/pyFDS-Evac#199 ·
维护者通常 1 天内回复
-
难度 2/5 1-3 小时 新手友好度 65/100
521xueweihan/HelloGitHub#3790 ·
-
难度 2/5 1-3 小时 新手友好度 78/100
sandialabs/atlas-ui-3#978 ·
维护者通常 1 天内回复
-
area: tests perceived difficulty: 2
难度 2/5 1-3 小时 新手友好度 72/100
Nitjsefnie-Harness-Commons/daedalus#1255 ·
维护者通常 1 天内回复
-
难度 2/5 1-3 小时 新手友好度 86/100
EleutherAI/lm-evaluation-harness#4256 ·
维护者通常 1 天内回复