00001 subroutine source_HaC_CF_qeos(sou)
00002 use grid_parameter, only : nrg, ntg, npg, nrf, ntf, npf
00003 use interface_sourceterm_HaC_CF
00004
00005 use interface_sourceterm_HaC_drot_SFC_qeos
00006 use interface_interpo_fl2gr_midpoint
00007 use interface_correct_matter_source_midpoint
00008 use make_array_3d
00009 implicit none
00010 real(8), pointer :: sou(:,:,:), sou1(:,:,:), sou2(:,:,:), souf(:,:,:)
00011
00012 call alloc_array3d(sou1,0,nrg,0,ntg,0,npg)
00013 call alloc_array3d(sou2,0,nrg,0,ntg,0,npg)
00014 call alloc_array3d(souf,0,nrf,0,ntf,0,npf)
00015 call sourceterm_HaC_CF(sou1)
00016
00017 call sourceterm_HaC_drot_SFC_qeos(souf)
00018 call interpo_fl2gr_midpoint(souf,sou2)
00019 call correct_matter_source_midpoint(sou2)
00020 sou(0:nrg,0:ntg,0:npg) = sou1(0:nrg,0:ntg,0:npg) &
00021 & + sou2(0:nrg,0:ntg,0:npg)
00022 deallocate(sou1)
00023 deallocate(sou2)
00024 deallocate(souf)
00025
00026 end subroutine source_HaC_CF_qeos
