if (FALSE) {
## get the flowpath from a head catchment to
## the sea
## Be patient, this might take a while.
system.time(fps <- 973842223 %up% rhineCCM2$PFAFSTETTER)
flowpath <- subset(rhineCCM2, subset = fps)
sp::plot(flowpath)
## does parallel threats increase the speed
## and provide the same results?
op.old <- options()
options("hydrocode.parallel" = TRUE)
system.time(fpp <- 973842223 %up% rhineCCM2$PFAFSTETTER)
all.equal(fps, fpp)
## reset previous settings
options(op.old)
}
Run the code above in your browser using DataLab