00001 subroutine calc_angmom_asympto(cobj)
00002 use phys_constant, only : long, pi
00003 use grid_parameter, only : nrg, ntg, npg, rgout
00004 use coordinate_grav_r, only : rg
00005 use def_matter_parameter, only : radi
00006 use make_array_2d
00007 use def_quantities, only : angmom_asymp
00008 use interface_source_angmom_asympto
00009 use interface_surf_int_grav_rg
00010 implicit none
00011 real(long) :: fac8pi
00012 real(long) :: surf, rg_asympt
00013 real(long), pointer :: sousf(:,:)
00014 integer :: irg, ir
00015 character(len=2), intent(in) :: cobj
00016
00017 call alloc_array2d(sousf, 0, ntg, 0, npg)
00018
00019 rg_asympt = 10.0**(dlog10(rgout)*2.0/3.0)
00020 do ir = 0, nrg
00021 irg = ir
00022 if (rg(ir).ge.rg_asympt) exit
00023 end do
00024
00025 call source_angmom_asympto(sousf,irg)
00026 call surf_int_grav_rg(sousf,surf,irg)
00027 fac8pi = 1.0d0/(8.0d0*pi)
00028 angmom_asymp = fac8pi*radi**2*surf
00029 if (cobj.eq.'bh') angmom_asymp = fac8pi*surf
00030
00031
00032
00033 deallocate(sousf)
00034 end subroutine calc_angmom_asympto
