In the example 05,
cat > si.chdens.in << EOF
&input
nfile = 1
filepp(1) = 'sicharge'
weight(1) = 1.0
iflag = 2
plot_out = 1
output_format = 2
fileout = 'si.rho.dat'
e1(1) =1.0, e1(2)=1.0, e1(3) = 0.0,
e2(1) =0.0, e2(2)=0.0, e2(3) = 1.0,
nx=56, ny=40
is there anyone who can tell which direction it calculated?