00001 subroutine calc_AH_BBH_CF
00002 use phys_constant, only : long, pi
00003 use grid_parameter, only : ntg, npg
00004 use def_quantities_bh, only : AHmass, AHarea
00005 use interface_source_AHarea_CF
00006 use interface_surf_int_grav_rg
00007 use make_array_2d
00008 implicit none
00009 real(long) :: surf
00010 real(long), pointer :: sousf(:,:)
00011
00012 call alloc_array2d(sousf,0,ntg,0,npg)
00013
00014 call source_AHarea_CF(sousf)
00015 call surf_int_grav_rg(sousf,surf,0)
00016 AHarea = surf
00017 AHmass = dsqrt(AHarea/(16.0d0*pi))
00018
00019 deallocate(sousf)
00020 end subroutine calc_AH_BBH_CF
