R/sph.tailup.R

In jayverhoef/SSN: Spatial Modeling on Stream Networks

sph.tailup <- function(dist.hydro, weight, parsil = parsil, range1 = range1)
{
	no <- length(dist.hydro[,1])
	np <- length(dist.hydro[1,])
	parsil*(matrix(rep(1, times = no*np), nrow = no) - 1.5*(dist.hydro/range1) +
		0.5*(dist.hydro/range1)^3)*(dist.hydro < range1)*weight
}