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>
&INPUTPP<br>
prefix = A' ,<br>
outdir = '' ,<br>
filplot = 'Apsi' ,<br>
plot_num = 7,<br>
kpoint=1,<br>
kband=200,<br>
/<br>
&PLOT<br>
nfile = 1 ,<br>
filepp(1) = 'A',<br>
weight(1) = 1.0,<br>
fileout = 'Apsi.dat' <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>