Hi all,
There are 365 days of soil moisture NC files and I am trying to find out how
many days the values are below and above this certain threshold are repeated by
R. However, I couldn't reach exactly what I wanted. For example, Daily soil
moisture is below 0.3 without interrupting how many days in 365 days. NC file
contains annual soil moisture values daily
nctoarray <- function(ncfname, varid = NA) { nc <- nc_open(ncfname)
a <- aperm(ncvar_get(nc), c(2,1,3)) nc_close(nc) a }
function(x, threshold = 0.28, below = TRUE) {
if (below) {
y <- ifelse(x < threshold,1,0)
} else y <- ifelse(x > threshold,1,0)
y2 <- rle(y)
sel <- which(y2$values == 1)
max(y2$lengths[sel])
}
m1 <- suppressWarnings(apply(a,c(1,2), consechours, 0.3, TRUE))
m2 <- suppressWarnings(apply(a,c(1,2), consechours, 0.4, FALSE))
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
