Update neutral population
Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
type(NeutralPopulation), | intent(inout) | :: | pop | Neutral Population |
subroutine parallel_merge_populations(pop)
!+Update neutral population
type(NeutralPopulation), intent(inout) :: pop
!+ Neutral Population
integer :: i, j, k, ic, ind(3)
call parallel_sum(pop%dens)
do ic=1,beam_grid%ngrid
call ind2sub(beam_grid%dims,ic,ind)
i = ind(1) ; j = ind(2) ; k = ind(3)
call parallel_merge_reservoirs(pop%res(i,j,k))
enddo
end subroutine parallel_merge_populations