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

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

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

&nbsp; end do<br>

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

&nbsp;&nbsp;&nbsp;&nbsp; write (*,*) &quot; pseudopotentials have different mesh &quot;<br>

&nbsp;&nbsp;&nbsp;&nbsp; interpolate = .true.<br>

&nbsp; end if<br>

&nbsp; write (*,*) &quot;INTERPOLATE =&quot;, interpolate<br>

&nbsp; if (interpolate) call errore (&quot;virtual&quot;, &amp;<br>

&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; &quot;grid interpolation is not working yet&quot;,1) <br><br>&nbsp;&quot;<br><br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; If mesh(1) isn&#39;t equal mesh(2) ,the value of interpolate is set&nbsp; to True. if&nbsp;  interpolate =true, the error must&nbsp; happen. This seem showed&nbsp; that&nbsp; two pseudopotentials can&#39;t&nbsp; mix&nbsp; if they have different mesh grid. But the message of error which &quot;grid interpolation is not working yet&quot;&nbsp; seems indicates that a subroutine which is used to interpolate should be called before the error. So I am&nbsp; confused. If someone know how to deal with it , Please tell 
me.Thanks a lot! &nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;  <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>