求助，利用ncl怎么将风云2f卫星CTT数据插值到wrf网格

油麦面疙瘩

rt。想要把风云2f卫星CTT数据插值到wrf网格上再输出为nc格式，目前是使用了ncl自带的ESMF网格转换，但是画出来图很奇怪，就好像是丢失了一部分数据一样。如下图：
我的脚本如下：
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/wrf/WRFUserARW.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/esmf/ESMF_regridding.ncl"

begin
   file_path="fy-2f/CTT/FY2F_CTT_MLT_NOM_20180521_0600.hdf"
   fi=addfile(file_path,"r")
   CTT=fi->CTT_Hourly_Product
   LVR = CTT@LowerValidRange ; The valid range of data is [LVR,UVR].
   UVR = CTT@UpperValidRange

   file2_path="/fy-2f/CTT/NOM_ITG_2288_2288(0E0N)_LE.dat"
   latlon2d=fbindirread(file2_path, 0,(/2,2288,2288/), "float")
   src_lat = latlon2d(1,:,:)
   src_lon = latlon2d(0,:,:)
   src_lat@_FillValue=300.0
   src_lon@_FillValue=300.0
   src_lon = src_lon+112.0
   src_lat@units="degrees_north"
   src_lon@units="degrees_east"

   CTT@lat2d      = src_lat
   CTT@lon2d      = src_lon
   CTT@units      = "degree Kelvin"
   CTT@long_name  = "Temperature of Cloud Top"
   CTT@coordinates = "lat2d lon2d"
   CTT@_FillValue  = -999.

   CTT_1d               = ndtooned(CTT)
   ind_notValid         = ind(CTT_1d.le.LVR.or.CTT_1d.gt.UVR)
   CTT_1d(ind_notValid)    = CTT@_FillValue
   CTT               = onedtond(CTT_1d,dimsizes(CTT))

   CTT = CTT-273.15  ;(convert degree Kelvin to degree Celsius)
   CTT@units = "degree Celsius"
   printMinMax(CTT,0)
   printVarSummary(CTT)

   ;插值
   file3_path="WRF/wrfout_d01_2018-05-21_05:00:00.nc"
   fi3=addfile(file3_path, "r")
   dst_lat = fi3->XLAT(0,:,:)
   dst_lon = fi3->XLONG(0,:,:)

   Opt =True
   Opt@Interpmethod = "bilinear"

   Opt@SrcGridLat        = src_lat           ; source grid
   Opt@SrcGridLon        = src_lon
   Opt@SrcInputFileName = file2_path
   Opt@SrcRegional       = True

   Opt@DstGridLat        = dst_lat           ; destination grid
   Opt@DstGridLon        = dst_lon
   Opt@DstRegional       = True

   Opt@ForceOverwrite   = True
   Opt@Debug            = True
   Opt@PrintTimings     = True

   CTT_regrid  = ESMF_regrid(CTT,Opt)     ; Do the regridding
   printVarSummary(CTT_regrid)

   ;设置绘图属性

   wks=gsn_open_wks("png", "ctt-1-0")

   res=True

    res@cnFillOn=True

   res@mpMinLatF = min(dst_lat)
   res@mpMaxLatF = max(dst_lat)
   res@mpMinLonF = min(dst_lon)
   res@mpMaxLonF = max(dst_lon)

   plot_regrid = gsn_csm_contour_map(wks, CTT_regrid, res)

end

我用到的数据也贴上来了，求助一下大家看我这个是哪里的问题？假如是插值的问题，那么用哪个插值函数比较合适呢？先谢谢各位大佬了（鞠躬）
wrfout文件太大了，传不上来(꒦_꒦)
FY2F_CTT_MLT_NOM_20180521_0600.hdf (6.75 MB, 下载次数: 1)

2020-2-11 16:11 上传

点击文件名下载附件



NOM_ITG_2288_2288(0E0N)_LE.dat (39.94 MB, 下载次数: 37)

2020-2-11 16:14 上传

点击文件名下载附件

油麦面疙瘩

自己顶一下

denghao

楼主解决了吗？是不是CTT变量的经纬度属性设置问题?