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

Go to the documentation of this file.

subroutine sourceterm_HaC_drot_SFC(sou)
  use phys_constant, only : long, pi
  use grid_parameter, only : nrf, ntf, npf
  use def_metric_on_SFC_CF, only : psif, alphf
  use def_matter, only : emd, jomef_int
  use def_matter_parameter, only : ber, radi
  use interface_interpo_linear_type0
  implicit none
  real(long), pointer :: sou(:,:,:)
  real(long) :: emdfc, rhofc, prefc, hhfc, utfc, rhoHc
  real(long) :: psifc, alpfc, ene, jomef_intfc
  integer    :: irf, itf, ipf
!
! --- Source of the Hamiltonian constraint to compute 
!     the conformal factor psi.
!
  do ipf = 0, npf
    do itf = 0, ntf
      do irf = 0, nrf
        psifc =  psif(irf,itf,ipf)
        alpfc = alphf(irf,itf,ipf)
        emdfc =   emd(irf,itf,ipf)
        jomef_intfc = jomef_int(irf,itf,ipf)
        if (irf.eq.nrf) then 
!          emdfc = 1.0d-15
          emdfc = 0.0d0
        end if
        call peos_q2hprho(emdfc, hhfc, prefc, rhofc, ene)
        utfc  = hhfc/ber*exp(jomef_intfc)
        rhoHc = hhfc*rhofc*(alpfc*utfc)**2 - prefc
!
        sou(irf,itf,ipf) = - radi**2*2.0d0*pi*psifc**5*rhoHc
      end do
    end do
  end do
end subroutine sourceterm_HaC_drot_SFC

---