getan Subroutine

public subroutine getan(antot, antota, antotp, local_only)

Compute velocity space integrals over $h_s$: $$\mathrm{antot} = n_s q_s \int d^3 v f_0 J_0(Z_s) h_s$$ $$\mathrm{antota} = 2 \beta n_s q_s \sqrt{\tfrac{T_s}{m_s}} \int d^3 v f_0 J_0(Z_s) h_s$$ $$\mathrm{antotp} = n_s T_s \int d^3 v f_0 \frac{v_\perp^2 J_1(Z_s)}{Z_s} h_s$$

See Colin's field equation notes for more details

Always uses gnew as $h_s$

FIXME: Get these notes on the website

Arguments

  subroutine getan (antot, antota, antotp, local_only)
    use dist_fn_arrays, only: gnew
    use theta_grid, only: ntgrid
    implicit none
    complex, dimension (-ntgrid:,:,:), intent (in out) :: antot, antota, antotp
    logical, intent(in), optional :: local_only
    call getan_from_dfn(gnew, antot, antota, antotp, local_only)
  end subroutine getan