我从ECAD（欧洲气候评估和数据集）下载了一个数据库，其中包含降水数据。
摘要（对照） rr（变量1）具有形状464 201 23376
所以它有23376次时间步，因为这对于一次运行来说太多了，我每次都单独运行。但是当试图运行12000以上的时间步时，我会得到一个致命错误（11000不是问题）。我除了内存问题，但我的笔记本电脑没有显示内存问题。此步骤导致问题：cntlprcp1 = get.var.ncdf（nc = control，varid =“rr”，start，count）。
control = open.ncdf(con="rr_0.25deg_reg_v9.0.nc", write=FALSE, readunlim=FALSE)
paste(control$filename,"has",control$nvars,"variables"), fill=TRUE)
lonmat = matrix(get.var.ncdf(nc=control,varid="longitude"),nrow=464,ncol=201)
latmat = matrix(get.var.ncdf(nc=control,varid="latitude"),nrow=201,ncol=464)
timearr = get.var.ncdf(nc=control,varid="time")
t1 = julian(x=1, d=1, y=1950, origin=c(month = 1, day = 1, year = 1601))
targettime = julian(x=7,d=4,y=1980,origin=c(month = 1, day = 1, year = 1601))
inds = (1:dim(timearr))
tind = inds[targettime-t1 == timearr]
ndims = control$var[['rr']]$ndims
varsize = control$var[['rr']]$varsize
start = c( 1, 1, **16000**)
count = c(varsize[1], varsize[2], 1)
cntlprcp1 = get.var.ncdf(nc=control,varid="rr",start,count)
cntlprcp1[cntlprcp1<0] = NA
x = 1:nrow(cntlprcp1) # R plots rows along X axis!
y = 1:ncol(cntlprcp1)
image.plot(x,y,cntlprcp1,col=tim.colors())