fRE Function

  1. korc_simple_equilibrium_pdf.f90
  2. korc_simple_equilibrium_pdf
  3. fRE

private function fRE(eta, p)

Arguments

Type IntentOptional AttributesName
real(kind=rp), intent(in) :: eta
real(kind=rp), intent(in) :: p

Return Value real(kind=rp)


Contents

Source Code

fRE

Source Code

FUNCTION fRE(eta,p)
	REAL(rp), INTENT(IN) :: eta ! pitch angle in degrees
	REAL(rp), INTENT(IN) :: p ! momentum in units of mc
	REAL(rp) :: fRE
	REAL(rp) :: A

	A = (2.0_rp*pdf_params%E/(pdf_params%Zeff + 1.0_rp))*(p**2/SQRT(p**2.0_rp + 1.0_rp))
	fRE = 0.5_rp*A*EXP(A*COS(deg2rad(eta)))/SINH(A)
!	fRE = fRE*PR(eta,p,pdf_params%Bo,pdf_params%lambda)
END FUNCTION fRE