NCL双线性插值后居然没变化

zpxie

path_dir = ../TRMM/3B42_daily.2012.06.25.7.nc"
path_nc=./wrfout_d03_2012-06-24_18:00:00.nc"
a=addfile(path_dir, "r")
b=addfile(path_nc, "r")

lat= a->latitude(321:336)                    ;提取一个区域的数据
lon= a->longitude(361:378)
rain=a->r(0,321:336,361:378)                 ;latitude:[30.375..34.125]
                                                             ;longitude:[90.375..94.625]

lat2d=b->XLAT(0,:,:)            
lon2d=b->XLONG(0,:,:)

;using bilinear interpolation
x0=new((/121/),float)          ;构造需要插值的经纬度格点
y0=new((/121/),float)

y0(0)=30.375
x0(0)=90.375

do nn=1,120
  y0(nn)=nn*0.03125+30.375
  x0(nn)=nn*0.03541667+90.375
end do


linint2_rain=linint2_Wrap(lon, lat, rain, False, x0, y0, 0) ;双线性插值

;Change the map-projection
xgrid=rgrid2rcm(y0, x0, linint2_rain, lat2d, lon2d, 1)

经过双线性插值后的结果跟没有插值，直接xgrid=rgrid2rcm(lat, lon, rain, lat2d, lon2d, 1)的画图结果一致。
不知道什么原因，请大神赐教。

风往北吹

我觉得没必要用双线性插值啊，rgrid2rcm里面应该也用到了插值吧。

zpxie

风往北吹 发表于 2014-11-1 08:59
我觉得没必要用双线性插值啊，rgrid2rcm里面应该也用到了插值吧。

照这样的结果 我觉得应该也是用到了插值。但是我很好奇的是 插值前的数据格点多少并不影响插值的结果，是不是说明双线性插值在这种情况下并不具备优势呢？您一般用什么插值 对于降水、或者土壤温湿度来说。

天生我才

请问问题解决了吗？