<div dir="ltr"><div>Thank you. I also had noticed that upf2plotcore wasn't working any longer, but I lacked the courage to try to find out why. There is by the way another code that does the same (upftools/extract_core.x) but it seems to me that  it does strange things (that is: what it extracts is not what one expects)</div><div><br></div><div>Paolo</div><div><br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Nov 25, 2019 at 1:33 PM Aldo Ugolotti <<a href="mailto:a.ugolotti@campus.unimib.it">a.ugolotti@campus.unimib.it</a>> wrote:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">Dear QE developers,<br>
<br>
I recently had to generate a core hole pseudopotential, from which I <br>
needed to extract the core state wavefunction. As usual, after running <br>
the atomic code, I tried to run upf2plotcore.sh - both as included in <br>
the 6.4.1 version of QE - but the output of the latter was empty with <br>
the exception of a single line "# number of core states 1 =  1 0;". <br>
After some checks, I noticed that the awk commands within the script <br>
itself were looking for a slightly different pattern than that included <br>
in the UPF file, and the string manipulation to extract some parameters <br>
was not working as intended.<br>
<br>
Hoping it could be of help, I am attaching a copy of the script I edited <br>
to fix it; the changes are summarized here as well:<br>
<br>
1) in reading the size of the radial grid (mesh_size, nris) the regex <br>
patterns in gsub functions have been substituted as:<br>
<br>
     gsub ( /.*mesh_size="/, "", nris )<br>
<br>
     gsub ( /[^0-9]*".*/, "", nris )<br>
<br>
2) in reading the radial grid, the awk patterns to be searched have been <br>
changed as "/PP_R /" -> "<PP_R>" and "<PP_R"  -> "<PP_R>"<br>
<br>
As for 2), I have not tested it with previous versions, but I noticed <br>
that v6.0 was reporting the "old" format of the tags; about point 1) I <br>
cannot guess why it wouldn't work for me.<br>
<br>
Bests,<br>
<br>
-- <br>
Aldo Ugolotti<br>
<br>
PhD student<br>
Materials Science Dept. U5,<br>
Università degli Studi di Milano-Bicocca<br>
via Cozzi 55,<br>
20125 Milano (MI)<br>
Italy<br>
e-mail: <a href="mailto:a.ugolotti@campus.unimib.it" target="_blank">a.ugolotti@campus.unimib.it</a><br>
<br>
<br>
_______________________________________________<br>
developers mailing list<br>
<a href="mailto:developers@lists.quantum-espresso.org" target="_blank">developers@lists.quantum-espresso.org</a><br>
<a href="https://lists.quantum-espresso.org/mailman/listinfo/developers" rel="noreferrer" target="_blank">https://lists.quantum-espresso.org/mailman/listinfo/developers</a><br>
</blockquote></div><br clear="all"><br>-- <br><div dir="ltr" class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>