ncl 站点插值出错

团团

中国160站站点插值的程序和错误如图，求解答
begin
fr="e:/cyg/cygdrive/tsummer_d.txt"
r=asciiread(fr, 160, "float")
;printVarSummary(r)
;flat="e:/cyg/cygdrive/lat.txt"
;lat=asciiread(flat, 41, "float")
;flon="e:/cyg/cygdrive/lon.txt"
;lon=asciiread(flon, 71, "float")
flatlon="e:/cyg/cygdrive/lat_lon.txt"
latlon=asciiread(flatlon, (/160,2/), "float")
zlat=latlon(:,0)
zlon=latlon(:,1)
;print(zlat)
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
olon=new(71, "float")
olat=new(41, "float")
do i=0,70
    olon(i)=70+i
end do  
do l=0,40
    olat(l)=15+l
end do
olon!0="zlon"
olon@long_name="lon"
olon@units="degrees_east"
olon&lon=olon
olat!0="zlat"
olat@long_name="lat"
olat@units="degrees_north"
olat&lat=olat

data1=new((/41,71/), "float")
;cha=new((/41,71/), float)
r@_FillValue=9999
rscan=(/10,5,3/)
data1=obj_anal_ic_deprecated(zlon, zlat, r, olon ,olat, rscan, False)

packard

olon!0="zlon"
上句，你定义了一个叫 "zlon"的 dimension.
下句，你在存取一个叫 "lon" 的dimension。这两个名字要一样
olon&lon=olon

yeah...

应该是olon&zlon =olon吧