00001 subroutine compute_psi2wme
00002 use phys_constant, only : long
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
00010
00011 index = - dble(index_wme)
00012 call compute_fnc_inversion_index(psi,wme,index)
00013 wme(0,0:ntg,0:npg) = 0.0d0
00014 log_wme(1:nrg,0:ntg,0:npg) = dlog(wme(1:nrg,0:ntg,0:npg))
00015 log_wme(0, 0:ntg,0:npg) = -40.0d0
00016 fac2 = 1.0d0/dsqrt(2.0d0)
00017 log_N(1:nrg,0:ntg,0:npg) = fac2*dlog(alph(1:nrg,0:ntg,0:npg))
00018 log_N(0, 0:ntg,0:npg) = -40.0d0
00019
00020 end subroutine compute_psi2wme
