00001 subroutine source_trG_WL(sou)
00002 use grid_parameter, only : nrg, ntg, npg, nrf, ntf, npf
00003 use interface_sourceterm_trG_CF
00004
00005 use interface_sourceterm_trG_drot_SFC
00006 use interface_sourceterm_trG_WL
00007 use interface_interpo_fl2gr_midpoint
00008 use interface_correct_matter_source_midpoint
00009 use make_array_3d
00010 implicit none
00011 real(8), pointer :: sou(:,:,:), sou1(:,:,:), sou2(:,:,:),
00012 sou3(:,:,:), souf(:,:,:)
00013 call alloc_array3d(sou1,0,nrg,0,ntg,0,npg)
00014 call alloc_array3d(sou2,0,nrg,0,ntg,0,npg)
00015 call alloc_array3d(sou3,0,nrg,0,ntg,0,npg)
00016 call alloc_array3d(souf,0,nrf,0,ntf,0,npf)
00017 call sourceterm_trG_CF(sou1)
00018
00019 call sourceterm_trG_drot_SFC(souf)
00020 call interpo_fl2gr_midpoint(souf,sou2)
00021 call correct_matter_source_midpoint(sou2)
00022 call sourceterm_trG_WL(sou3)
00023 sou(0:nrg,0:ntg,0:npg) = sou1(0:nrg,0:ntg,0:npg) &
00024 & + sou2(0:nrg,0:ntg,0:npg) &
00025 & + sou3(0:nrg,0:ntg,0:npg)
00026 deallocate(sou1)
00027 deallocate(sou2)
00028 deallocate(sou3)
00029 deallocate(souf)
00030
00031 end subroutine source_trG_WL
