admmass_frg.f90
Go to the documentation of this file.00001 subroutine admmass_frg
00002
00003 use phys_constant, only : pi
00004 use def_metric_1D, only : alphf, psif
00005 use CB_fR_param_flphy, only : ahoadm
00006 use def_matter_1D, only : emd, ber, pinx, radi
00007 use grid_parameter_1D, only : nrf
00008 use weight_grav_1D, only : wgdr
00009 implicit none
00010 real(8) :: alpfc, emdfc, hhfc, prefc, psifc, rhofc, rhoHc, sou, utfc, weiflu
00011 integer :: ir
00012
00013 ahoadm = 0.0d0
00014 emd(nrf) = 1.0d-14
00015 do ir = 0, nrf
00016
00017 if (emd(ir) <= 0.0d0) emd(ir) = 1.0d-14
00018 emdfc = emd(ir)
00019 rhofc = emd(ir)**pinx
00020 prefc = rhofc*emdfc
00021 alpfc = alphf(ir)
00022 psifc = psif(ir)
00023 hhfc = 1.0d0 + (pinx+1.0d0)*emdfc
00024 utfc = hhfc/ber
00025
00026 rhoHc = hhfc*rhofc*(alpfc*utfc)**2 - prefc
00027 sou = rhoHc*psifc**5
00028
00029 weiflu = wgdr(ir)
00030 ahoadm = ahoadm + sou*weiflu
00031
00032 end do
00033
00034 ahoadm = radi**3*ahoadm*4.0d0*pi
00035
00036 end subroutine admmass_frg
