求助！！我想要调用wrf_user_intrp2d和wrf_user_intrp3d函数

打破砂锅纹到底

我想要调用wrf_user_intrp2d和wrf_user_intrp3d这两个函数，但使用数据不是wrfout的数据。
xlat = wrf_user_getvar(a, "XLAT",0)
plane = (/100,150,149,120/)
angle = 90
opts   = True
lat_plane = wrf_user_intrp2d(xlat,plane,angle,opts)

使用wrfout里的XLAT变量可以调用成功；print（lat_plane)的结果也是没错的

但是，我自己定义的二维经纬度格点，调用wrf_user_intrp2d函数会出现问题。好像只能识别起始点的经纬度。这部分脚本如下：
olon = new(501,"float")
olat = new(501,"float")
lat2d = new((/501,501/),"float")
lon2d = new((/501,501/),"float")

do i=0,501-1
     olon(i) = {start_lon} + 0.01*(i)
end do
do l=0,501-1
     olat(l) = {start_lat} + 0.01*(l)
end do

do k = 0,501-1
    lon2d(:,k)=olon(k)
end do
do m = 0,501-1
    lat2d(m,:)=olat(m)
end do


lat_plane = wrf_user_intrp2d(lat2d,plane,45.,True)
lon_plane = wrf_user_intrp2d(lon2d,plane,45.,True)
print(lon_plane)
dimsX = dimsizes(lat_plane)
do ii = 0,dimsX（0）-2
     gsn_polyline(wks,plot,(/lon_plane(ii),lon_plane(ii+1)/),(/lat_plane(ii),lat_plane(ii+1)/),lnres)
end do

插值后的lon_plane输出倒是没报错，但是输出的值都相同，好像都是起始点的经纬度，没有识别终止点的经纬度，本应画出的直线也都叠加于一点。

已解决
file:///C:\Users\dell\AppData\Roaming\Tencent\Users\229288458\QQ\WinTemp\RichOle\I9]{1WD8@25Q]9RNQS}0OQG.png

file:///C:\Users\dell\AppData\Roaming\Tencent\Users\229288458\QQ\WinTemp\RichOle\I9]{1WD8@25Q]9RNQS}0OQG.png

打破砂锅纹到底

自己顶一顶

freekiller

参考http://bbs.06climate.com/forum.php?mod=viewthread&tid=34227
TRMM数据插值和wrf剖面比较

打破砂锅纹到底

freekiller 发表于 2016-5-29 11:58
参考http://bbs.06climate.com/forum.php?mod=viewthread&tid=34227
TRMM数据插值和wrf剖面比较

谢谢楼上帮助