wrf数据插值出错

酒国

参考http://bbs.06climate.com/forum.php?mod=viewthread&tid=44335中的程序，将wrf模拟的数据再插值到特定的网格点中，并输出到nc文件中。
运行下面的程序，插值那一行提示有错误。
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"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"

f = addfile("wrfout_d01_2009-10-01_00_00_00.nc","r")
TZ=f->T2(0,:,:)
print(TZ)

lat2d = f->XLAT(0,:,0)
lat2d@units = "degrees_north"
lon2d = f->XLONG(0,0,:)
lon2d@units = "degrees_east"

olon =fspan(70.0041666666667,139.99583333327,8400)
olat =fspan(10.0041666684426,54.9958333350973,5400)
olon!0        ="lon"
olon@long_name="lon"
olon@units    ="degrees-east"
olon&lon      =olon
olat!0        ="lat"
olat@long_name="lat"
olat@units    ="degrees-north"
olat&lat      =olat

rscan = (/10,5,3/)
fo = obj_anal_ic_deprecated(lon2d,lat2d,TZ,olon,olat,rscan,False)
错误提示：
fatal:Number of dimensions in parameter (2) of (obj_anal_ic_deprecated) is (2), (1) dimensions were expected
fatal:["Execute.c":8637]:Execute: Error occurred at or near line 37 in file 07.ncl

printVarSummary(fo)
;将数据写入nc文件中

FO=fo
FO!0="lat"
FO!1="lon"
FO@units="K"
FO@long_name="Temperature at 2m"
FO@missing_value=-1e+10
FO@average_op_ncl="dim_avg_n over dimension(s): lev"

fout=addfile("HiRes.nc","c")
setfileoption(fout,"DefineMode",True)
;生成全域属性
fileAtt=True
fileAtt@creation_data="FO"
;事先定义坐标变量
;dimNames=(/"olat","olon"/)
;dimSizes=(/-1,-1/)
;dimUnlim=(/True,True/)
;filedimdef(fout,dimNames,dimSizes,dimUnlim)

filevardef(fout,"lat",typeof(olat),getvardims(olat))
filevardef(fout,"lon",typeof(olon),getvardims(olon))
filevardef(fout,"TZ",typeof(FO),getvardims(FO))

filevarattdef(fout,"TZ",FO)

fout->lat=olat
fout->lon=olon
fout->TZ=FO
setfileoption(fout,"DefineMode",True)
ff=addfile("HiRes.nc","r")
print(ff)
p=ff->TZ
print(p(0,1,:))
end