gyro_radius Function

  1. fidasim.f90
  2. libfida
  3. gyro_radius

public function gyro_radius(fields, energy, pitch, Ai) result(gyro_rad)

Calculates mean gyro-radius

Arguments

TypeIntentOptionalAttributesName
type(LocalEMFields), intent(in) :: fields

Electromagnetic fields at guiding center
real(kind=Float64), intent(in) :: energy

Energy of particle
real(kind=Float64), intent(in) :: pitch

Particle pitch w.r.t the magnetic field
real(kind=Float64), intent(in) :: Ai

Atomic Mass of ion

Return Value real(kind=Float64)

Mean gyro-radius


Calls

proc~~gyro_radius~~CallsGraph proc~gyro_radius gyro_radius proc~pitch_to_vec pitch_to_vec proc~gyro_radius->proc~pitch_to_vec proc~gyro_step gyro_step proc~gyro_radius->proc~gyro_step proc~cross_product cross_product proc~gyro_step->proc~cross_product

Contents

Source Code

gyro_radius

Source Code

function gyro_radius(fields, energy, pitch, Ai) result (gyro_rad)
    !+ Calculates mean gyro-radius
    type(LocalEMFields), intent(in)          :: fields
        !+ Electromagnetic fields at guiding center
    real(Float64), intent(in)                :: energy
        !+ Energy of particle
    real(Float64), intent(in)                :: pitch
        !+ Particle pitch w.r.t the magnetic field
    real(Float64), intent(in)                :: Ai
        !+ Atomic Mass of ion
    real(Float64) :: gyro_rad
        !+ Mean gyro-radius

    real(Float64), dimension(3) :: vi_norm, r_step
    real(Float64) :: vabs, phi
    integer :: i,n

    vabs  = sqrt(energy/(v2_to_E_per_amu*Ai))

    gyro_rad = 0.d0
    n = 6
    do i=1,n
        phi = i*2*pi/n
        call pitch_to_vec(pitch, phi, fields, vi_norm)
        call gyro_step(vabs*vi_norm, fields, Ai, r_step)
        gyro_rad = gyro_rad + norm2(r_step)/n
    enddo

end function gyro_radius