ncl站点插值后成缺测

我不喜欢吃生姜

目前查看了这两个网页的内容，发现他们都是lat ,lon赋值出错，但我检查自己的这方面没有发现错误，请各位大神帮我看一看，是哪个地方出错了。。
http://bbs.06climate.com/forum.p ... hread&tid=45458



脚本如下：

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
path = "/public/home/data/Z9370/"
filenames = systemfunc("cd "+path+";ls CAPPI*.txt")
nfiles = dimsizes(filenames)
print(nfiles)
olon = new(250,"float")
olat = new(250,"float")
ref = new((/250,250/),"float")
do i = 1,250
        olon(i-1) = 118.2 + (i-1)*0.02
end do
do j = 1,250
        olat(j-1) = 25.0 + (j-1)*0.02
end do

     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 = (/0.05,0.03,0.02/)

i = 1
filename = filenames(i-1)
print("reading   "+filename)
dummy = asciiread(path+filename,-1,"float")
ncol = 3
nsta = dimsizes(dummy)/ncol
stationdata = onedtond(dummy,(/nsta,ncol/))
lat = stationdata(:,0)
lon = stationdata(:,1)
stn = stationdata(:,2)

ref = obj_anal_ic(lon,lat,stn,olon,olat,rscan,False)
ref!1 = "lon"
ref!0 = "lat"
ref&lat = olat
ref&lon = olon
ref&lat@units = "degrees_north"
ref&lon@units = "degrees_east"
print(ref)
printVarSummary(ref)
end

我不喜欢吃生姜

另一个帖子
http://bbs.06climate.com/forum.php?mod=viewthread&tid=35503

tyq嗯

请问楼主是怎么解决的

bgs3

我怎么没有看到怎去txt文件中的数据函数啊，不清楚systemfunc是否能取代这个特点

bgs3

bgs3 发表于 2017-12-6 13:42
我怎么没有看到怎去txt文件中的数据函数啊，不清楚systemfunc是否能取代这个特点

抱歉，没有看到后面的，不过好像 olon&lon=olon 最好再看下