calc_mass_WL.f90

Go to the documentation of this file.
00001 subroutine calc_mass_WL
00002   use phys_constant, only  :   long, pi
00003   use grid_parameter, only :   nrg, ntg, npg, nrf, ntf, npf
00004   use def_matter_parameter, only : radi
00005   use make_array_3d
00006   use def_quantities, only : admmass, komarmass, komarmass_nc
00007   use interface_source_adm_mass_WL
00008   use interface_source_komar_mass_peos
00009   use interface_source_komar_mass_compact_WL
00010   use interface_vol_int_grav
00011   use interface_vol_int_fluid
00012   implicit none
00013   real(long)       ::     fac2pi, fac4pi
00014   real(long)       ::     volg, volf
00015   real(long),pointer :: soug(:,:,:), souf(:,:,:)
00016 !
00017 !
00018   call alloc_array3d(soug, 0, nrg, 0, ntg, 0, npg)
00019   call alloc_array3d(souf, 0, nrf, 0, ntf, 0, npf)
00020 !      
00021   call source_adm_mass_WL(soug,souf)
00022   call vol_int_grav(soug,volg)
00023   call vol_int_fluid(souf,volf)
00024 !
00025   fac2pi = 0.5d0/pi
00026   admmass = fac2pi*(radi*volg + radi**3*volf)
00027 !
00028   call source_komar_mass_peos(soug,souf)
00029   call vol_int_grav(soug,volg) 
00030   call vol_int_fluid(souf,volf)
00031 !
00032   fac4pi = 0.25d0/pi
00033   komarmass_nc = fac4pi*(radi*volg + radi**3*volf)
00034 !
00035   call source_komar_mass_compact_WL(souf)
00036   call vol_int_fluid(souf,volf)
00037 !
00038   komarmass = radi**3*volf
00039 !
00040 !  write (6,'(a20,1p,e14.6)') ' ADM   mass =       ', admmass
00041 !  write (6,'(a20,1p,e14.6)') ' Komar mass =       ', komarmass
00042 !
00043   deallocate(soug)
00044   deallocate(souf)
00045 end subroutine calc_mass_WL