Dear everyone,<br>         I am now facing a problem. My charge density contour has little interactions. Although I changed the # of level, I found little difference. Can anybody have a look at my script and point out my mistake?  I really want to get a contour with interactions, in which they should be there. <br>
<br> 1)<a href="http://rho_pp.in">rho_pp.in</a><br><br>        &inputpp<br>    prefix  = 'cntcu'<br>    filplot = 'cntcucharge'<br>    plot_num= 0<br> /<br> &plot<br>    nfile = 1<br>    filepp(1) = 'cntcucharge'<br>
    weight(1) = 1.0<br>    iflag = 2<br>    output_format = 2<br>    fileout = 'cntcu.rho.dat'<br>    e1(1) =1.0, e1(2)=1.0, e1(3) = 0.0,<br>    e2(1) =0.0, e2(2)=0.0, e2(3) = 1.0,<br>    nx=56, ny=40<br> /<br><br>
2) The interactive process:<br><br>Input file > cntcu.rho.dat<br>r0   :   0.0000  0.0000  0.0000<br>tau1 :   1.0000  0.0000  0.0000<br>tau2 :   0.0000  1.0000  0.0000<br>read  25 atomic positions<br>output file >  <a href="http://test.rho.ps">test.rho.ps</a><br>
Read  56 * 40  grid<br>Logarithmic scale (y/n)? > n<br>Bounds:     0.000004    1.403000<br>min, max, # of levels > 0 1.41 6<br><br>3)the <a href="http://rho.ps">rho.ps</a> as attached.<br><br><br>