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

Go to the documentation of this file.

subroutine source_virial_WL_MHD(sou_Tkin,sou_Pint,sou_Memf,sou_Wgra)
  use phys_constant,  only : long
  use grid_parameter, only : nrg, ntg, npg
  use def_metric,     only : psi
  use def_SEM_tensor_EMF, only : trsm_EMF
  use interface_source_virial_WL
  use interface_interpo_linear_type0
  implicit none
  real(long), pointer :: sou_Tkin(:,:,:), sou_Pint(:,:,:), 
                        sou_Memf(:,:,:), sou_Wgra(:,:,:)
  real(long) :: psigc, psigc6
  integer :: irg, itg, ipg
!
! --- Compute source terms for virial relations.
! --- sources for kinetic and internal energies are defined on SCF grid points.
! --- sources for EM, and gravitational energies are defined on mid points.
!
  call source_virial_WL(sou_Tkin,sou_Pint,sou_Wgra)
!
  do ipg = 1, npg
    do itg = 1, ntg
      do irg = 1, nrg
        call interpo_linear_type0(psigc,psi,irg,itg,ipg)
        psigc6 = psigc**6
        sou_Memf(irg,itg,ipg)= trsm_EMF(irg,itg,ipg)*psigc6
      end do
    end do
  end do
!
end subroutine source_virial_WL_MHD

---