wrf.interplevel returns NaN when interpolating at a level where the sample points are

[qiuyang50]

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())

[erogluorhan]

Hi @qiuyang50 thanks very much for sharing this with us! We'll look into it when we have some time.