Dear Gabriele, <br>
Thank you very much for your help and answer.<br>
Please let me ask you the last question in this regard. This is the input that I use for getting psi^2 on parallel planes:<br>
 &amp;INPUTPP<br>
                      prefix = A&#39; ,<br>
                      outdir = &#39;&#39; ,<br>
                     filplot = &#39;Apsi&#39; ,<br>
                    plot_num = 7,<br>
                      kpoint=1,<br>
                      kband=200,<br>
 /<br>
 &amp;PLOT<br>
                       nfile = 1 ,<br>
                   filepp(1) = &#39;A&#39;,<br>
                   weight(1) = 1.0,<br>
                     fileout = &#39;Apsi.dat&#39; <br>
                       iflag = 2 ,<br>
               output_format = 2 ,<br>
                       e1(1) = 1,<br>
                       e1(2) = 0,<br>
                       e1(3) = 0,<br>
                       e2(1) = 0,<br>
                       e2(2) = 1,<br>
                       e2(3) = 0,<br>
                       x0(1) = 1,<br>
                       x0(2) = 1,<br>
                       x0(3) = 1,<br>
                          nx = 100 ,<br>
                          ny = 100 ,<br>
 /<br>
<br>
I thought by changing the x0(3) from zero to highest value of structure in z direction, I could get psi^2 on parallel planes.<br>
I use the following input for average.x:<br>
1<br>
Apsi<br>
1<br>
1000<br>
3<br>
3<br>
<br>
however when I change x0(3), the output of average remains
the same. Does it make sense? If positive, how could I get psi^2 on
planes?<br>
Thank you very much for your help; I really really appreciate it. <br>
<br>
M Alaaii<br><br>