bb_cx_rates Subroutine

  1. fidasim.f90
  2. libfida
  3. bb_cx_rates

public subroutine bb_cx_rates(denn, vn, vi, rates)

Get beam-beam neutralization/cx rates

Arguments

TypeIntentOptionalAttributesName
real(kind=Float64), intent(in), dimension(nlevs):: denn

Neutral density [cm^-3]
real(kind=Float64), intent(in), dimension(3):: vn

Neutral velocity [cm/s]
real(kind=Float64), intent(in), dimension(3):: vi

Ion velocity [cm/s]
real(kind=Float64), intent(out), dimension(nlevs):: rates

Reaction rates [1/s]

Calls

proc~~bb_cx_rates~~CallsGraph proc~bb_cx_rates bb_cx_rates interface~interpol_coeff interpol_coeff proc~bb_cx_rates->interface~interpol_coeff

Called by

proc~~bb_cx_rates~~CalledByGraph proc~bb_cx_rates bb_cx_rates proc~neutral_cx_rate neutral_cx_rate proc~neutral_cx_rate->proc~bb_cx_rates proc~fida_weights_los fida_weights_los proc~fida_weights_los->proc~bb_cx_rates proc~get_total_cx_rate get_total_cx_rate proc~get_total_cx_rate->proc~neutral_cx_rate program~fidasim fidasim program~fidasim->proc~fida_weights_los proc~halo halo program~fidasim->proc~halo proc~dcx dcx program~fidasim->proc~dcx proc~fida_f fida_f program~fidasim->proc~fida_f proc~npa_mc npa_mc program~fidasim->proc~npa_mc proc~fida_weights_mc fida_weights_mc program~fidasim->proc~fida_weights_mc proc~npa_f npa_f program~fidasim->proc~npa_f proc~npa_weights npa_weights program~fidasim->proc~npa_weights proc~fida_mc fida_mc program~fidasim->proc~fida_mc proc~halo->proc~neutral_cx_rate proc~dcx->proc~get_total_cx_rate proc~fida_f->proc~get_total_cx_rate proc~npa_mc->proc~get_total_cx_rate proc~fida_weights_mc->proc~get_total_cx_rate proc~npa_f->proc~get_total_cx_rate proc~npa_weights->proc~get_total_cx_rate proc~fida_mc->proc~get_total_cx_rate

Contents

Source Code

bb_cx_rates

Source Code

subroutine bb_cx_rates(denn, vn, vi, rates)
    !+ Get beam-beam neutralization/cx rates
    real(Float64), dimension(nlevs), intent(in)  :: denn
        !+ Neutral density [cm^-3]
    real(Float64), dimension(3),     intent(in)  :: vn
        !+ Neutral velocity [cm/s]
    real(Float64), dimension(3),     intent(in)  :: vi
        !+ Ion velocity [cm/s]
    real(Float64), dimension(nlevs), intent(out) :: rates
        !+ Reaction rates [1/s]

    real(Float64), dimension(nlevs,nlevs) :: neut  !!rate coeff
    real(Float64) :: eb !! relative Energy
    type(InterpolCoeffs1D) :: c
    real(Float64) :: dlogE, logEmin, logeb
    real(Float64) :: vrel !! relative velocity
    integer :: ebi, neb, err

    !Eeff
    vrel=norm2(vi-vn)
    eb=v2_to_E_per_amu*vrel**2  ! [kev/amu]
    logeb = log10(eb)
    logEmin = tables%H_H_cx_cross%logemin
    dlogE = tables%H_H_cx_cross%dlogE
    neb = tables%H_H_cx_cross%nenergy
    call interpol_coeff(logEmin,dlogE,neb,logeb,c,err)
    ebi = c%i
    if(err.eq.1) then
        if(inputs%verbose.ge.0) then
            write(*,'(a)') "BB_CX_RATES: Eb out of range of H_H_cx table. Using nearest energy value."
            write(*,'("eb = ",ES10.3," [keV]")') eb
        endif
        if(ebi.lt.1) then
            ebi=1
            c%b1=1.0 ; c%b2=0.0
        else
            ebi=neb-1
            c%b1=0.0 ; c%b2=1.0
        endif
    endif

    neut(:,:) = (c%b1*tables%H_H_cx_cross%log_cross(:,:,ebi) + &
                 c%b2*tables%H_H_cx_cross%log_cross(:,:,ebi+1))

    where (neut.lt.tables%H_H_cx_cross%minlog_cross)
        neut = 0.d0
    elsewhere
        neut = exp(neut*log_10)
    end where

    rates=matmul(neut,denn)*vrel

end subroutine bb_cx_rates