00001 subroutine calc_ang_mom_asymp
00002
00003 use phys_constant, only : long, pi
00004 use grid_parameter, only : nrg, ntg, npg
00005 use coordinate_grav_r, only : hrg
00006 use def_matter_parameter, only : radi
00007 use def_quantities, only : angmom_asymp
00008 use make_array_2d
00009 use interface_source_ang_mom_asymp
00010 use interface_surf_int_grav
00011 implicit none
00012 real(long) :: surf
00013 real(long), pointer :: sousf(:,:)
00014 integer :: irg
00015
00016 call alloc_array2d(sousf,1,ntg,1,npg)
00017
00018 irg = nrg - 4
00019 call source_ang_mom_asymp(sousf,irg)
00020 call surf_int_grav(sousf,surf,irg)
00021 angmom_asymp = radi**2*surf/(8.0d0*pi)
00022
00023
00024
00025 deallocate(sousf)
00026
00027 end subroutine calc_ang_mom_asymp
