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

Go to the documentation of this file.

subroutine correct_matter_source_midpoint(sou)
  use phys_constant, only : long, pi
  use grid_parameter, only : nrg, ntg, npg
  use coordinate_grav_r, only  :   drg, rg
  use def_matter, only : rs
  use interface_interpo_linear_type0_2Dsurf
  implicit none
  real(long), pointer :: sou(:,:,:)
  real(long) :: hsurf
  integer    :: irg, itg, ipg
!
! -- Correct the weight for a matter source in a volume integral
  do ipg = 1, npg
    do itg = 1, ntg
      call interpo_linear_type0_2Dsurf(hsurf,rs,itg,ipg)
      do irg = 0, nrg
        if (hsurf < rg(irg)) then 
          sou(irg,itg,ipg) = sou(irg,itg,ipg)*(hsurf-rg(irg-1))/drg(irg)
          exit
        end if
      end do
    end do
  end do
end subroutine correct_matter_source_midpoint

---