fzero Function

private function fzero(r, a, ko, P) result(f)

@param f Value of function. @param r Guess value of radial position of the particles. @param a Minor radius of the toroidal distribution @param ko Decay rate of radial distribution, see of \ref korc_spatial_distribution.exponential_torus. @param P Deviate of a random uniform distribution in the interval .

Arguments

Type IntentOptional AttributesName
real(kind=rp), intent(in) :: r
real(kind=rp), intent(in) :: a
real(kind=rp), intent(in) :: ko
real(kind=rp), intent(in) :: P

Return Value real(kind=rp)


Called by

proc~~fzero~~CalledByGraph proc~fzero fzero proc~exponential_elliptic_torus exponential_elliptic_torus proc~exponential_elliptic_torus->proc~fzero proc~exponential_torus exponential_torus proc~exponential_torus->proc~fzero

Contents

Source Code


Source Code

FUNCTION fzero(r,a,ko,P) RESULT(f)
	REAL(rp) 				:: f
	REAL(rp), INTENT(IN) 	:: r
	REAL(rp), INTENT(IN) 	:: a
	REAL(rp), INTENT(IN) 	:: ko
	REAL(rp), INTENT(IN) 	:: P

	f = EXP(-ko*r)*(1.0_rp + r*ko) + ( 1.0_rp - EXP(-ko*a)*(1.0_rp + a*ko) )*P - 1.0_rp
END FUNCTION fzero