00001 subroutine compute_alps2wmeN_mpt(impt)
00002 use phys_constant, only : long, nmpt
00003 use grid_parameter, only : nrg, ntg, npg
00004 use def_metric, only : psi, alph
00005 use def_metric_pBH, only : wme, log_wme, log_N, index_wme
00006 use interface_compute_fnc_inversion_index
00007 implicit none
00008 real(long) :: index, fac2
00009 integer :: irg, itg, ipg, impt
00010
00011 index = - dble(index_wme)
00012 fac2 = 1.0d0/dsqrt(2.0d0)
00013 call compute_fnc_inversion_index(psi,wme,index)
00014 log_wme(1:nrg,0:ntg,0:npg) = dlog(wme(1:nrg,0:ntg,0:npg))
00015 log_N( 1:nrg,0:ntg,0:npg) = fac2*dlog(alph(1:nrg,0:ntg,0:npg))
00016 if (impt.ne.nmpt) then
00017 wme(0,0:ntg,0:npg) = 0.0d0
00018 log_wme(0,0:ntg,0:npg) = -40.0d0
00019 log_N( 0,0:ntg,0:npg) = -40.0d0
00020 else
00021 log_wme(0,0:ntg,0:npg) = dlog(wme(0,0:ntg,0:npg))
00022 log_N( 0,0:ntg,0:npg) = fac2*dlog(alph(0,0:ntg,0:npg))
00023 end if
00024
00025 end subroutine compute_alps2wmeN_mpt
