00001 subroutine calc_ang_mom_peos
00002
00003 use phys_constant, only : long, pi
00004 use grid_parameter, only : nrf, ntf, npf
00005 use def_matter_parameter, only : radi
00006 use def_quantities, only : angmom
00007 use make_array_3d
00008 use interface_source_ang_mom_peos
00009 use interface_vol_int_fluid
00010 implicit none
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_ang_mom_peos(souf)
00018 call vol_int_fluid(souf,volf)
00019 angmom = radi**4*volf
00020
00021
00022
00023 deallocate(souf)
00024
00025 end subroutine calc_ang_mom_peos
