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?