Store birth particles/density in birth
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
integer(kind=Int32), | intent(in), | dimension(3) | :: | ind | beam_grid indices |
|
integer(kind=Int32), | intent(in) | :: | neut_type | Neutral type |
||
real(kind=Float64), | intent(in) | :: | dflux | Deposited flux |
subroutine store_births(ind, neut_type, dflux)
!+ Store birth particles/density in [[libfida:birth]]
integer(Int32), dimension(3), intent(in) :: ind
!+ [[libfida:beam_grid]] indices
integer(Int32), intent(in) :: neut_type
!+ Neutral type
real(Float64), intent(in) :: dflux
!+ Deposited flux
!$OMP CRITICAL(store_births_1)
birth%dens( neut_type,ind(1),ind(2),ind(3))= &
birth%dens(neut_type,ind(1),ind(2),ind(3)) + dflux
!$OMP END CRITICAL(store_births_1)
end subroutine store_births