00001 subroutine calc_proper_mass_peos_spin
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 : propermass
00007 use interface_source_proper_mass_peos_spin
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_proper_mass_peos_spin(souf)
00018 call vol_int_fluid(souf,volf)
00019
00020 propermass = radi**3*volf
00021
00022 write (6,'(a20,1p,e14.6)') ' Proper mass = ', propermass
00023
00024 deallocate(souf)
00025 end subroutine calc_proper_mass_peos_spin
