00001 subroutine source_MoC_WL_BH(souvec)
00002 use grid_parameter, only : nrg, ntg, npg
00003 use interface_sourceterm_MoC_CF_with_divshift
00004 use interface_sourceterm_MoC_WL
00005 use make_array_4d
00006 implicit none
00007 real(8), pointer :: souvec(:,:,:,:)
00008 real(8), pointer :: souvec2(:,:,:,:), souvec3(:,:,:,:)
00009 integer :: ia
00010
00011 call alloc_array4d(souvec2,0,nrg,0,ntg,0,npg,1,3)
00012 call alloc_array4d(souvec3,0,nrg,0,ntg,0,npg,1,3)
00013
00014 call sourceterm_MoC_CF_with_divshift(souvec2)
00015 call sourceterm_MoC_WL(souvec3)
00016 souvec(0:nrg,0:ntg,0:npg,1:3) = souvec2(0:nrg,0:ntg,0:npg,1:3) &
00017 & + souvec3(0:nrg,0:ntg,0:npg,1:3)
00018
00019 deallocate(souvec2)
00020 deallocate(souvec3)
00021
00022 end subroutine source_MoC_WL_BH
