Nothing
sinkfill <- function(DEM,res,degree) {
## data preparation
stopifnot(is(DEM,"matrix"))
if(min(as.vector(DEM[!is.na(DEM)])) < -9000)
stop("DEM contains unrealistic values (< -9000)")
DEM[is.na(DEM)] <- -9999
nrow <- dim(DEM)[1]
ncol <- dim(DEM)[2]
## exclude an area of the DEM from the calculations
## still to implement in R (using NA's)
result <- .C("c_sinkfill",
PACKAGE = "topmodel",
as.double(DEM),
result = double(nrow*ncol + 2),
as.integer(nrow),
as.integer(ncol),
as.double(res),
as.double(degree))$result
## formatting of the results
## the first position in result gives the number of iterations that have been performed
result[result > 999998] <- NA
if(result[1] == -1) cat("incomplete sink removal\n")
else {
if(result[1] == 100) cat("Maximum number of iterations reached (100).\nSink removal is probably not complete; please run again\n")
else cat("number of iterations = ", result[1], "\n")
}
cat("number of sinks left = ", result[2], "\n")
return(matrix(result[3:(nrow*ncol+2)], nrow = nrow))
}
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.