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

Go to the documentation of this file.

subroutine test_source
  use phys_constant, only     : long, pi
  use coordinate_grav_r, only : rg
  use grid_parameter, only    : nrf, ntf, npf
  use def_matter, only        : emd, rs
  implicit none
  integer     ::   irf, itf, ipf
  real(long)  ::   zfac, small = 1.0d-15
!
!
  do ipf = 0, npf
    do itf = 0, ntf
    rs(itf,ipf) = 1.0d0
        emd(0,itf,ipf) = 4.0d0*pi
!!        emd(0,itf,ipf) = 4.0d0*pi
      do irf = 1, nrf
        emd(irf,itf,ipf) = 4.0d0*pi
!!        emd(irf,itf,ipf) = 4.0d0*pi*sin(pi*rg(irf))/(pi*rg(irf))
      end do
    end do
  end do
!
!
end subroutine test_source

---