用ncl做站点资料的插值

zyj0625

请教各位，我想用ncl做站点资料的插值，是画全国的pm2.5的分布，但是插值后没有报错，但是图不对，想问是哪里出现错误？这是我的代码
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"
begin
pm=asciiread("pm_new_daily.txt", -1, "float")      ;--读取雾霾的数据
latlon=asciiread("flatlon_new.txt", (/1579,2/), "float")    ;--读取经纬度数据
pmm=reshape(pm,(/1579,151/))         ;--共1579站，151天数据
pm1=pmm(:,0)      ;--用第一天的数据作图
;pm1=dim_avg_n_Wrap(pmm,1)     ;--做气候平均
pm1@_FillValue=99999.000
lonst=latlon(:,1)
latst=latlon(:,0)

delta=0.5    ;插值网格大小
lons=70
lone=140
lats=15
late=55
nlat=ceil((late-lats)/delta)+1
nlon=ceil((lone-lons)/delta)+1
glat=fspan(lats,late,toint(nlat))
glon=fspan(lons,lone,toint(nlon))
glon!0="lon"
glon@long_name="longitude"
glon@units="degree_east"
glat!0="lat"
glat@long_name="latitude"
glat@units="degree_north"
rscan=(/5.,2.5,1/)
pm1:=obj_anal_ic_Wrap(lonst, latst, pm1, glon, glat, rscan, False)    ;插值函数


wks=gsn_open_wks("png","pm1")
res=True
res@gsnDraw=False
res@gsnFrame=False
res@gsnAddCyclic=False
res@gsnMaximize=True
res@mpMinLatF=15
res@mpMaxLatF=55
res@mpMinLonF=70
res@mpMaxLonF=140
res@mpDataBaseVersion="MediumRes"
res@mpDataSetName="Earth..4"
res@mpOutlineOn=True
res@mpOutlineSpecifiers=(/"China"/)
; ; res@mpOutlineSpecifiers=(/"China:states"/)
res@mpNationalLineColor="black"
res@mpNationalLineThicknessF=4

res@mpAreaMaskingOn=True
res@mpMaskAreaSpecifiers=(/"China"/)
res@mpOceanFillColor=0
res@mpInlandWaterFillColor=0
res@mpLandFillColor=0
res@mpGeophysicalLineColor="black"
res@mpGeophysicalLineThicknessF=4

plot=gsn_csm_contour_map(wks, pm1, res)
;标记站点
; mres=True
; mres@gsMarkerSizeF=0.005
; mres@gsMarkerThicknessF=0.5
; dum=gsn_add_polymarker(wks, plot, lonst, latst, mres)
draw(plot)
frame(wks)

end

画出的图是这样的 就很奇怪

银狐

pm1:=obj_anal_ic_Wrap(lonst, latst, pm1, glon, glat, rscan, False)    ;插值函数
这句，pm1后面多了个冒号吧。。。

九十

请问楼主问题解决了吗？有脚本可以分享吗