computeLubForce_tangential

public function computeLubForce_tangential(h, hc, mu, R, ut, utr) result(Flub)

Type	Intent		Name
real(kind=rk),	intent(in)	::	h	Gap between particles (physical units)
real(kind=rk),	intent(in)	::	hc	Critical gap (physical units) If h < hc we use hc to compute the lubrication forces to prevent them from diverging as h -> 0.
real(kind=rk),	intent(in)	::	mu	Dynamic viscosity (physical units)
real(kind=rk),	intent(in)	::	R	Radius of the particles
real(kind=rk),	intent(in)	::	ut(3)	Relative velocity in tangential direction between the two particles
real(kind=rk),	intent(in)	::	utr(3)	Relative velocity in tangential direction caused by rotational motion of particles