Ошибка в ncvar_get при переборе данных терраклимата - PullRequest
Новая
       83

Ошибка в ncvar_get при переборе данных терраклимата

0 голосов
/ 19 сентября 2019

Поэтому я пытаюсь перебрать весь проект 1001 генома с tarraclimate, но как только я получаю 12 записей, я получаю следующую ошибку:

Ошибка в ncvar_get_inner (ncid2use, varid2use, nc $ var [[li]] $ missval, addOffset,: Ошибка: переменная имеет 3 dim, но в start есть 2 записи. Они должны совпадать! Кроме того: Предупреждение: В varsize - start: длинная длина объекта не кратна короткой длине объекта

любые советы о том, что может быть причиной, которая будет принята с благодарностью. Координаты законны, должно появиться , и я не вижу ничего особенного между этими координатами и предыдущимив списке. Вот первые 15 или около того, это ошибки из цикла 12-го.

введите описание изображения здесь 

for (row in 1:nrow(coordlist)) {
y <- c(coordlist[row,"longitude"],coordlist[row,"latitude"])
baseurlagg <- paste0(paste0("http://thredds.northwestknowledge.net:8080/thredds/dodsC/agg_terraclimate_",var),"_1958_CurrentYear_GLOBE.nc")

nc <- nc_open(baseurlagg)
lon <- ncvar_get(nc, "lon")
lat <- ncvar_get(nc, "lat")
flat = match(abs(lat - y[2]) < 1/48, 1)
latindex = which(flat %in% 1)
flon = match(abs(lon - y[1]) < 1/48, 1)
lonindex = which(flon %in% 1)
start <- c(lonindex, latindex, 1)
count <- c(1, 1, -1)


# read in the full period of record using aggregated files

data <- as.data.frame(ncvar_get(nc, varid = var,start = start, count))
data$CS_Number<- coordlist[row,"CS_number"]
data2 <- cbind(data,dates)
outsideDF <- rbind(outsideDF,data2)
}
...