00001 subroutine source_trfreeG_WL_BH(souten)
00002 use grid_parameter, only : nrg, ntg, npg
00003 use interface_sourceterm_trfreeG_WL
00004 use make_array_4d
00005 implicit none
00006 real(8), pointer :: souten(:,:,:,:)
00007 real(8), pointer :: sou2(:,:,:,:)
00008 integer :: ia
00009
00010 call alloc_array4d(sou2,0,nrg,0,ntg,0,npg,1,6)
00011
00012 call sourceterm_trfreeG_WL(sou2)
00013 souten(0:nrg,0:ntg,0:npg,1:6) = sou2(0:nrg,0:ntg,0:npg,1:6)
00014
00015 deallocate(sou2)
00016
00017 end subroutine source_trfreeG_WL_BH
