Dear users.<br>     I want to mix the pseudopotential used the  virtual.x . but  it can't finish.I read the  code of virtual.f90 in the uptools,then I  think the codes which is related  with the problem  is following<br>
<br>"<br>       if (mesh(1).ne.mesh(2) ) then<br>          interpolate = .true.<br>

      end if<br>     capel = 0.d0<br>   do i=1,upf_mesh<br>

     capel = capel + abs(r(i,1)-r(i,2)) + abs(rab(i,1)-rab(i,2))<br>

  end do<br>

  if (capel.gt.1.d-6) then<br>

     write (*,*) " pseudopotentials have different mesh "<br>

     interpolate = .true.<br>

  end if<br>

  write (*,*) "INTERPOLATE =", interpolate<br>

  if (interpolate) call errore ("virtual", &<br>

                  "grid interpolation is not working yet",1) <br><br> "<br><br>       If mesh(1) isn't equal mesh(2) ,the value of interpolate is set  to True. if   interpolate =true, the error must  happen. This seem showed  that  two pseudopotentials can't  mix  if they have different mesh grid. But the message of error which "grid interpolation is not working yet"  seems indicates that a subroutine which is used to interpolate should be called before the error. So I am  confused. If someone know how to deal with it , Please tell 
me.Thanks a lot!              <br><br><br>-- <br>Yanchao Wang<br>National lab of Superhard Materials,JiLin Univ.Changchun 130012,P.R.China<br><a href="mailto:Email:wyanchao@gmail.com">Email:wyanchao@gmail.com</a>