<div dir="ltr">Thanks.<div><br></div><div>But it looks like that charge density is only plotted on a two dimensional plane. Is there a way to plot the 3D wave function (abs and the sign as phase)?</div><div><br></div><div><br></div></div><div class="gmail_extra"><br><div class="gmail_quote">On Fri, Jan 15, 2016 at 12:53 PM, Gavin Abo <span dir="ltr"><<a href="mailto:gsabo@crimson.ua.edu" target="_blank">gsabo@crimson.ua.edu</a>></span> wrote:<br><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">You might have a look at $WIENROOT/SRC_lapw7/<a href="http://read.me" rel="noreferrer" target="_blank">read.me</a> in WIEN2k 14.2.<br>
<br>
The text file mentions three source code files (w1gpl.f, w2gpl.f, and WtoGP.f, which are in $WIENROOT/SRC_lapw7/DOC_psink).<br>
<br>
I have not tried it, but it is my understanding that you can compile them in the DOC_psink directory with something like:<br>
<br>
ifort -o w1gpl w1gpl.f<br>
ifort -o w2gpl w2gpl.f<br>
ifort -o WtoGP WtoGP.f<br>
<br>
The executables are then copied from the DOC_psink directory to the $WIENROOT directory:<br>
<br>
cp w1gpl $WIENROOT/w1gpl<br>
cp w2gpl $WIENROOT/w2gpl<br>
cp WtoGP $WIENROOT/WtoGP<br>
<br>
The programs should then be useable in your case directory to convert case.psink to case.gpl:<br>
<br>
w1gpl < case.psink > case.gpl<br>
<br>
or<br>
<br>
w2gpl < case.psink > case.gpl<br>
<br>
or<br>
<br>
WtoGP < case.psink > case.gpl<br>
<br>
Finally, the data in case.gpl is plotted with gnuplot:<br>
<br>
gnuplot case.gpl<br>
<br>
It looks like an example is given at:<br>
<br>
<a href="http://www.affinity-science.com/wien2k/faq.html#q4" rel="noreferrer" target="_blank">http://www.affinity-science.com/wien2k/faq.html#q4</a><br>
<br>
The website is in Japanese, but Google Translate [ <a href="https://translate.google.com/" rel="noreferrer" target="_blank">https://translate.google.com/</a> ] can be used to get it in English or another available language.<span class=""><br>
<br>
On 1/14/2016 7:20 PM, Yundi Quan wrote:<br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
Hi,<br>
I used lapw7 to generate case.psink which contains the wave function on a grid. Is there a program available to plot the data? I know that the wien2wannier has wplot2xsf routine which converts the case.psink and case.psiarg files into case.xsf compatible with xcrysden.<br>
</blockquote></span>
_______________________________________________<br>
Wien mailing list<br>
<a href="mailto:Wien@zeus.theochem.tuwien.ac.at" target="_blank">Wien@zeus.theochem.tuwien.ac.at</a><br>
<a href="http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien" rel="noreferrer" target="_blank">http://zeus.theochem.tuwien.ac.at/mailman/listinfo/wien</a><br>
SEARCH the MAILING-LIST at:  <a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html" rel="noreferrer" target="_blank">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/index.html</a><br>
</blockquote></div><br></div>