00001 subroutine surf_int_grav(sousf,surf,irg)
00002 use phys_constant, only : long
00003 use grid_parameter, only : ntg, npg
00004 use coordinate_grav_r, only : hrg
00005 use weight_midpoint_grav, only : hwtpgsf
00006 implicit none
00007 real(long), pointer :: sousf(:,:)
00008 real(long), intent(out) :: surf
00009 integer, intent(in) :: irg
00010 integer :: itg, ipg
00011 surf = 0.0d0
00012 do ipg = 1, npg
00013 do itg = 1, ntg
00014 surf = surf + sousf(itg,ipg)*hwtpgsf(itg,ipg)
00015 end do
00016 end do
00017 surf = surf*hrg(irg)**2
00018 end subroutine surf_int_grav
