00001 subroutine source_trG_WL_BH(sou)
00002 use grid_parameter, only : nrg, ntg, npg
00003 use interface_sourceterm_trG_CF
00004 use interface_sourceterm_trG_WL
00005 use make_array_3d
00006 implicit none
00007 real(8), pointer :: sou(:,:,:), sou1(:,:,:), sou3(:,:,:)
00008 call alloc_array3d(sou1,0,nrg,0,ntg,0,npg)
00009 call alloc_array3d(sou3,0,nrg,0,ntg,0,npg)
00010 call sourceterm_trG_CF(sou1)
00011 call sourceterm_trG_WL(sou3)
00012 sou(0:nrg,0:ntg,0:npg) = sou1(0:nrg,0:ntg,0:npg) &
00013 & + sou3(0:nrg,0:ntg,0:npg)
00014 deallocate(sou1)
00015 deallocate(sou3)
00016
00017 end subroutine source_trG_WL_BH
