ncl站点插值 数组超长

团团

用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,71
    olon(i)=70+i
end do  
do l=0,41
    olat(l)=15+l
end do
olat!0="zlat"
olat@long_name="zlat"
olat@units="degrees_north"
olat&lat=olat
olon!0="zlon"
olon@long_name="zlon"
olon@units="degrees_east"
olon&lon=olon
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

矩阵长度 71 的话，那下标就是从 0到70. 循环0到71会报错。

团团

packard 发表于 2017-8-26 13:04
矩阵长度 71 的话，那下标就是从 0到70. 循环0到71会报错。

我改了只有又有错误了
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="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



data1=new((/41,71/), "float")

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

团团

packard 发表于 2017-8-26 13:04
矩阵长度 71 的话，那下标就是从 0到70. 循环0到71会报错。

错误是
cannot create child process:[errno=11]
error occurred at or near line 11955 in file load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"

packard

团团 发表于 2017-8-26 13:12
错误是
cannot create child process:[errno=11]
error occurred at or near line 11955 in file load  ...

load 那些系统的script不要在 begin 和 end 之间load。应该在之前。希望问题是这个。