R/infiltration.R

In topmodel: Implementation of the Hydrological Model TOPMODEL in R

infiltration <- function(rain, parameters) {

  result <- .C("c_infiltration",
               PACKAGE = "topmodel",
               as.double(rain),
               as.double(parameters),
               as.integer(length(rain)),
               result = double(length(rain)))$result

  result[result < -9000] <- NA
  
  return(result)

}