00001 subroutine compute_alps2alph(pot,psi) 00002 use phys_constant, only : long 00003 use grid_parameter, only : nrg, ntg, npg 00004 implicit none 00005 real(long), pointer :: pot(:,:,:), psi(:,:,:) 00006 integer :: irg, itg, ipg 00007 do irg = 0, nrg 00008 do itg = 0, ntg 00009 do ipg = 0, npg 00010 pot(irg, itg, ipg) = pot(irg, itg, ipg)/psi(irg, itg, ipg) 00011 end do 00012 end do 00013 end do 00014 end subroutine compute_alps2alph