lin.taildown <- function(dist.hydro, a.mat, b.mat, parsil = parsil,
range1 = range1, useTailDownWeight, weight = NULL)
{
flow.connect <- b.mat == 0
no <- length(a.mat[,1])
np <- length(a.mat[1,])
V <- parsil*(matrix(rep(1, times = no*np), nrow = no) - dist.hydro/range1)*
(a.mat < range1)*flow.connect +
parsil*(matrix(rep(1, times = no*np), nrow = no) - a.mat/range1)*
(a.mat < range1)*(1 - flow.connect)
if(useTailDownWeight == TRUE) V <- V*weight
V
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.