Use rmatch_augfun also in the US case.

git-svn-id: http://qeforge.qe-forge.org/svn/q-e/trunk/espresso@5540 c92efa57-630b-4861-b058-cf58834340f0
This commit is contained in:
dalcorso 2009-05-01 12:44:37 +00:00
parent 235f4099fc
commit e36788bbbe
1 changed files with 4 additions and 2 deletions

View File

@ -9,7 +9,8 @@ SUBROUTINE pseudo_q (qfunc, qfuncl)
USE kinds, ONLY : DP
USE io_global, ONLY : stdout
USE ld1_parameters, ONLY : nwfsx
USE ld1inc, ONLY : rcut, lls, grid, ndmx, lmx2, nbeta, ikk, ecutrho
USE ld1inc, ONLY : rcut, lls, grid, ndmx, lmx2, nbeta, ikk, ecutrho, &
rmatch_augfun
IMPLICIT NONE
!
REAL(DP), INTENT(IN) :: qfunc(ndmx,nwfsx,nwfsx)
@ -34,7 +35,8 @@ do ns=1,nbeta
! Find the matching point
!
ik=0
rmatch=min(rcut(ns),rcut(ns1))
! rmatch=min(rcut(ns),rcut(ns1))
rmatch=rmatch_augfun
do n=1,mesh
if (grid%r(n)>rmatch) then
ik=n