/home/uryu/Cocal_svn_checkinout/Cocal/code/Subroutine/source_HaC_WL_BH.f90

Go to the documentation of this file.

subroutine source_HaC_WL_BH(sou)
  use grid_parameter, only : nrg, ntg, npg
  use interface_sourceterm_HaC_CF
  use interface_sourceterm_HaC_WL
  use make_array_3d
  use def_metric, only : psi
  use coordinate_grav_r, only : hrg
  use interface_dadbscalar_type3
  implicit none
  real(8), pointer :: sou(:,:,:), sou1(:,:,:), sou3(:,:,:)
  real(8) :: dabfnc(3,3), Lappsi
  integer :: irg, itg, ipg
!
  call alloc_array3d(sou1,0,nrg,0,ntg,0,npg)
  call alloc_array3d(sou3,0,nrg,0,ntg,0,npg)
  call sourceterm_HaC_CF(sou1)
  call sourceterm_HaC_WL(sou3)
  sou(0:nrg,0:ntg,0:npg) = sou1(0:nrg,0:ntg,0:npg) &
  &                      + sou3(0:nrg,0:ntg,0:npg)
!
!itg = ntg/2-2; ipg = 2
!open(16,file='test_vec_pot',status='unknown')
!  do irg = 1, nrg
!    call dadbscalar_type3(psi,dabfnc,irg,itg,ipg)
!    Lappsi = dabfnc(1,1) + dabfnc(2,2) + dabfnc(3,3)
!    write(16,'(1p,12e20.12)') hrg(irg), (sou(irg,itg,ipg) - Lappsi)
!!!  /(abs(sou(irg,itg,ipg)) + abs(Lappsi))
!  end do
!close(16)
!stop
!
  deallocate(sou1)
  deallocate(sou3)
!
end subroutine source_HaC_WL_BH

---