attenuate Subroutine

  1. fidasim.f90
  2. libfida
  3. attenuate

public subroutine attenuate(ri, rf, vi, states, dstep_in)

Attenuate states along a trajectory

Arguments

Type IntentOptional AttributesName
real(kind=Float64), intent(in), dimension(3):: ri

Initial position in beam grid coordinates
real(kind=Float64), intent(in), dimension(3):: rf

Final position in beam grid coordinates
real(kind=Float64), intent(in), dimension(3):: vi

Initial velocity of neutral
real(kind=Float64), intent(inout), dimension(nlevs):: states

Density of states
real(kind=Float64), intent(in), optional :: dstep_in

Step length [cm]

Calls

proc~~attenuate~~CallsGraph proc~attenuate attenuate proc~colrad colrad proc~attenuate->proc~colrad proc~get_plasma get_plasma proc~attenuate->proc~get_plasma proc~eigen eigen proc~colrad->proc~eigen proc~get_rate_matrix get_rate_matrix proc~colrad->proc~get_rate_matrix proc~matinv matinv proc~colrad->proc~matinv proc~in_plasma in_plasma proc~get_plasma->proc~in_plasma proc~hqr2 hqr2 proc~eigen->proc~hqr2 proc~balback balback proc~eigen->proc~balback proc~elmtrans elmtrans proc~eigen->proc~elmtrans proc~elmhes elmhes proc~eigen->proc~elmhes proc~balance balance proc~eigen->proc~balance interface~interpol_coeff interpol_coeff proc~get_rate_matrix->interface~interpol_coeff proc~ludcmp ludcmp proc~matinv->proc~ludcmp proc~in_plasma->interface~interpol_coeff proc~xyz_to_uvw xyz_to_uvw proc~in_plasma->proc~xyz_to_uvw proc~interpol1d_coeff interpol1D_coeff interface~interpol_coeff->proc~interpol1d_coeff proc~interpol2d_coeff_arr interpol2D_coeff_arr interface~interpol_coeff->proc~interpol2d_coeff_arr proc~interpol2d_coeff interpol2D_coeff interface~interpol_coeff->proc~interpol2d_coeff proc~interpol1d_coeff_arr interpol1D_coeff_arr interface~interpol_coeff->proc~interpol1d_coeff_arr proc~hqrvec hqrvec proc~hqr2->proc~hqrvec proc~outerprod outerprod proc~ludcmp->proc~outerprod proc~interpol2d_coeff_arr->proc~interpol2d_coeff proc~interpol1d_coeff_arr->proc~interpol1d_coeff

Called by

proc~~attenuate~~CalledByGraph proc~attenuate attenuate proc~npa_weights npa_weights proc~npa_weights->proc~attenuate proc~npa_f npa_f proc~npa_f->proc~attenuate proc~npa_mc npa_mc proc~npa_mc->proc~attenuate program~fidasim fidasim program~fidasim->proc~npa_weights program~fidasim->proc~npa_f program~fidasim->proc~npa_mc

Contents

Source Code

attenuate

Source Code

subroutine attenuate(ri, rf, vi, states, dstep_in)
    !+ Attenuate `states` along a trajectory
    real(Float64), dimension(3), intent(in)        :: ri
        !+ Initial position in beam grid coordinates
    real(Float64), dimension(3), intent(in)        :: rf
        !+ Final position in beam grid coordinates
    real(Float64), dimension(3), intent(in)        :: vi
        !+ Initial velocity of neutral
    real(Float64), dimension(nlevs), intent(inout) :: states
        !+ Density of states
    real(Float64), intent(in), optional            :: dstep_in
        !+ Step length [cm]

    type(LocalProfiles) :: plasma
    real(Float64) :: photons, vabs, dt, dstep, dis,max_dis
    real(Float64), dimension(3) :: r0
    real(Float64), dimension(nlevs) :: dens

    if(present(dstep_in)) then
        dstep=dstep_in
    else
        dstep = sqrt(inter_grid%da) !cm
    endif

    max_dis = norm2(rf-ri)

    vabs = norm2(vi)
    dt = dstep/vabs

    call get_plasma(plasma,pos=ri)
    r0 = ri
    dis = 0.d0
    do while (dis.le.max_dis)
        call colrad(plasma,beam_ion,vi,dt,states,dens,photons)
        r0 = r0 + vi*dt
        dis = dis + dstep
        call get_plasma(plasma,pos=r0)
    enddo

end subroutine attenuate