R/sph.taildown.R

In SSN: Spatial Modeling on Stream Networks

sph.taildown <-
function(dist.hydro, a.mat, b.mat, parsil = parsil, range1 = range1,
	useTailDownWeight, weight = NULL)
{
	flow.connect <- b.mat == 0
	V <- parsil*(1 - 1.5*(dist.hydro/range1) +
		0.5*(dist.hydro/range1)^3)*(a.mat < range1)*flow.connect +
		parsil*(1 - 1.5*b.mat/range1 +
		0.5*a.mat/range1)*(1 - a.mat/range1)^2*
		(a.mat < range1)*(1 - flow.connect)
	V
}