def calcLSTC(Tb, rho_s, D_s, T_sc):
	"""
	Calculates sediment transport capacity for points along perimeter.

	Parameters
	----------
	Tb : boundary shear stress at each point [Pa]
	rho_s : sediment density [kg/m^3]
	D_s : sediment diameter [m]
	T_sc : critical shields stress of sediment [1]

	Returns
	-------
	q_t : sediment transport capacity at each point along perimeter

	Notes
	-----
	Returns -1 if there is not enough shear stress to transport sediment
	"""

	R_b = rho_s/999.97 - 1
	T_s = Tb/((rho_s - 999.97)*9.81*D_s)

	# If shear stress is under critical shields stress force alluviation
	if any((T_s-T_sc)<0):
		return -1
	else:
		q_t = 5.7*((T_s - T_sc)**(3./2))*2650*((R_b*9.81*D_s**3)**(1./2))
		return q_t