<meta http-equiv="Content-Type" content="text/html; charset=GB18030"><div>Dear NCL Community,</div><div>I want to interpolate the zonal wind in the wrfout-data to a height of 100m from the ground. </div><div>I use the function 'wrf_user_interp_level'. </div><div>Before interpolation, the data of all grid points have a value. After the interpolation, the data of one grid point is the default value. After checking, I found that the data at the height of 100m from the grid point to the ground can be obtained without interpolation, so I would like to ask you if there is a problem with this function. </div><div>My NCL version is 6.6.2. The script runs without errors.</div><div><br></div><div>My script is:</div><div>begin</div><div>a=addfile("wrfout_d01_2005-12-20_00:00:00.nc","r")</div><div>u = wrf_user_getvar(a,"ua",-1)</div><div>z = wrf_user_getvar(a,"z",-1)</div><div>hgt = a->HGT</div><div>dimll=dimsizes(z)</div><div>hfg=new((/dimll(0),dimll(1),dimll(2),dimll(3)/),float)</div><div>copy_VarMeta(z, hfg)</div><div><br></div><div>do time=0,dimll(0)-1</div><div>    do lev=0,dimll(1)-1</div><div>        do jj=0,dimll(3)-1</div><div>            do ii=0,dimll(2)-1</div><div>                hfg(time,lev,ii,jj)=z(time,lev,ii,jj)-hgt(time,ii,jj)</div><div>            end do</div><div>        end do</div><div>    end do</div><div>end do</div><div><br></div><div>u_plane = wrf_user_interp_level(u,hfg,100.,False)</div><div>loc = wrf_user_ll_to_ij(a,95.5,33,True)</div><div>locX = loc(0) - 1</div><div>locY = loc(1) - 1</div><div>print(u(3,:,locY,locX-1))</div><div>print(u_plane(3,locY,locX-1))</div><div>print(hfg(3,:,locY,locX-1))</div><div>end</div><div><br></div><div>And the result is:</div><div><font editorwarp__="1" style="display: inline; font-size: 14px; font-family: Verdana; color: rgb(0, 0, 0); background-color: rgba(0, 0, 0, 0); font-weight: 400; font-style: normal;"><img src="cid:BF07F122@88BD9C44.E16A4461" style="" id="img_insert_163187244682407148741469316555"><br><img src="cid:CE7AF370@E817BD29.E16A4461" style="" id="img_insert_163187246967907048646088301835"><br><img src="cid:37072071@55F57078.E16A4461" style="" id="img_insert_163187248717104969652558518789"><br>You see, hfg(3)=100, and i</font>t is the _FillValue<span style="background-color: rgba(0, 0, 0, 0);"> when the zonal wind is interpolated to 100m£¬but u(3)=9.758717. And u(3) is the</span> value obtained by interpolating the zonal wind to a height of 100m above the ground. So I would like to ask if this is the problem with the interpolation function?</div><div>I want to upload my wrfout data(1.8G) and scripts. Can I upload them?</div><div><br></div><font editorwarp__="1" style="display: inline; font-size: 14px; font-family: Verdana; color: rgb(0, 0, 0); background-color: rgba(0, 0, 0, 0); font-weight: 400; font-style: normal;"><div style="display: block; font-size: 14px; font-family: Verdana; color: rgb(0, 0, 0); background-color: rgba(0, 0, 0, 0); font-weight: 400; font-style: normal;">Any hint is much appreciated!</div></font><div></div><font color="#333333" face="verdana, sans-serif"><span style="font-size: 13.3333px;"></span></font>