fRE Function

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)


Calls

proc~~fre~2~~CallsGraph proc~fre~2 fRE proc~deg2rad~2 deg2rad proc~fre~2->proc~deg2rad~2

Called by

proc~~fre~2~~CalledByGraph proc~fre~2 fRE proc~sample_distribution sample_distribution proc~sample_distribution->proc~fre~2 proc~get_equilibrium_distribution get_equilibrium_distribution proc~get_equilibrium_distribution->proc~sample_distribution

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