00001 subroutine calc_rest_mass
00002 use phys_constant, only : long, pi
00003 use grid_parameter, only : nrf, ntf, npf
00004 use def_matter_parameter, only : radi
00005 use make_array_3d
00006 use def_quantities, only : restmass
00007 use interface_source_rest_mass
00008 use interface_vol_int_fluid
00009 implicit none
00010 real(long) :: fac2pi, fac4pi
00011 real(long) :: volf
00012 real(long), pointer :: souf(:,:,:)
00013 integer :: ir, it, ip
00014
00015 call alloc_array3d(souf, 0, nrf, 0, ntf, 0, npf)
00016
00017 call source_rest_mass(souf)
00018 call vol_int_fluid(souf,volf)
00019
00020 restmass = radi**3*volf
00021
00022
00023
00024 deallocate(souf)
00025 end subroutine calc_rest_mass
