NCL将wrf格点数据插值到站点，除了一两个站点都是缺测值，求求大家帮我看看！

海棠未眠

file_wrf=addfile(dir_wrf+"wrfout_d01_2011-03-15_06_00_00","r")
f=addfile("wrfchemi_d01_2011-03-15_06_00_00","r")
f_obs=asciiread(dir_obs+file_obs, (/32,8/), "float");(:,3)
lat2d=file_wrf->XLAT(0,:,:)
lon2d=file_wrf->XLONG(0,:,:)
lat=lat2d(:,0)
lon=lon2d(0,:)
e_ces1=f->E_CES(:,0,:,:)
cesnd=e_ces1
;;e_ces=3600*(3000*3000)*e_ces1 ;;单位换算
;copy_VarMeta(e_ces1,e_ces)
cesnd!0="lat"
cesnd!1="lon"
cesnd&lat=lat
cesnd&lon=lon
;printVarSummary(f_obs)
;站点的经纬度-------------------------
LAT=f_obs(:,0)
LON=f_obs(:,1)
;-----------------------------------------
  ;fo = rcm2points_Wrap (lat2d,lon2d,cesnd,LON,LAT, 0)   ; => fo(ntim,nlvl,5)
   fo=linint2_points_Wrap(lon,lat,cesnd, False, LON,LAT, 0)