I am sorry about previous e-mail. The problem was caused by my <span class="enfont">careless~ I should set the parameter </span><span class="enfont">filepp(1) as a file name instead of 1~<br></span><div></div><br>At 2011-04-26 11:47:44,"GAO&nbsp;Zhe"&nbsp;&lt;flux_ray12@163.com&gt; wrote:<br> <blockquote id="isReplyContent" style="border-left: 1px solid rgb(204, 204, 204); margin: 0px 0px 0px 0.8ex; padding-left: 1ex;">Dear QE developers and users:<br>I used pp.x to show the charge density of my structure: TiC (100) surface with one layer of Ni.<br>Of course, this model may not good enough, but, after pw.x running, pp.x gave an error and stopped.<br>The error is:<br><font color="#ff0000">Calling punch_plot, plot_num =&nbsp;&nbsp; 0<br>Writing data to file&nbsp; edensity<br>Reading data from file&nbsp; 1<br>%%%%%%%%%%%%%%%%%%<br>from plot_io : error #&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 2<br>opening file 1<br>%%%%%%%%%%%%%%%%%%</font><br>I also calculated for TiC cell without vacuum slab, every thing went right in this case. The difference is that the charge density file we obtained from filplot = '' in surface case was much larger (3.0MiB) than the one in cell case (smaller than 500KiB).<br>The version of QE I am using is 4.3, this problem occurred both in Red Hat, Mandriva and Ubuntu system.<br>My calculation script:<br><font color="#0000ff">#!/bin/sh<br><br>cat &gt;TiC.scf.in &lt;&lt;EOF<br>&amp;control<br>&nbsp;&nbsp; calculation = 'scf' ,<br>&nbsp;&nbsp; prefix = 'TiC_Ni' ,<br>&nbsp;&nbsp; pseudo_dir = '$PseudoDIR/' ,<br>&nbsp;&nbsp; outdir = '$TempDIR/' ,<br>&nbsp;&nbsp; disk_io = 'low' ,<br>/<br>&amp;system<br>&nbsp;&nbsp; ibrav = 6 ,<br>&nbsp;&nbsp; celldm(1) = 8.17873408 ,<br>&nbsp;&nbsp; celldm(3) = 2.0 ,<br>&nbsp;&nbsp; nat = 14 ,<br>&nbsp;&nbsp; ntyp = 3 ,<br>&nbsp;&nbsp; ecutwfc = 25.52 ,<br>&nbsp;&nbsp; ecutrho = 255.2 ,<br>&nbsp;&nbsp; occupations = 'smearing' ,<br>&nbsp;&nbsp; smearing = 'mp' ,<br>&nbsp;&nbsp; degauss = 0.015 ,<br>/<br>&amp;electrons<br>&nbsp;&nbsp; startingpot = 'atomic' ,<br>&nbsp;&nbsp; startingwfc = 'atomic+random' ,<br>&nbsp;&nbsp; diagonalization = 'david' ,<br>/<br>ATOMIC_SPECIES<br>&nbsp; Ti&nbsp; 47.900&nbsp; Ti.pbe-sp-van_ak.UPF<br>&nbsp;&nbsp; C&nbsp; 12.011&nbsp; C.pbe-van_ak.UPF<br>&nbsp; Ni&nbsp; 58.690&nbsp; Ni.pbe-nd-rrkjus.UPF<br>ATOMIC_POSITIONS crystal<br>&nbsp; Ti&nbsp; 0.00&nbsp; 0.00&nbsp; 0.00<br>&nbsp; Ti&nbsp; 0.50&nbsp; 0.50&nbsp; 0.00<br>&nbsp; Ti&nbsp; 0.50&nbsp; 0.00&nbsp; 0.25<br>&nbsp; Ti&nbsp; 0.00&nbsp; 0.50&nbsp; 0.25<br>&nbsp; Ti&nbsp; 0.00&nbsp; 0.00&nbsp; 0.50<br>&nbsp; Ti&nbsp; 0.50&nbsp; 0.50&nbsp; 0.50<br>&nbsp;&nbsp; C&nbsp; 0.50&nbsp; 0.00&nbsp; 0.00<br>&nbsp;&nbsp; C&nbsp; 0.00&nbsp; 0.50&nbsp; 0.00<br>&nbsp;&nbsp; C&nbsp; 0.00&nbsp; 0.00&nbsp; 0.25<br>&nbsp;&nbsp; C&nbsp; 0.50&nbsp; 0.50&nbsp; 0.25<br>&nbsp;&nbsp; C&nbsp; 0.00&nbsp; 0.50&nbsp; 0.50<br>&nbsp;&nbsp; C&nbsp; 0.50&nbsp; 0.00&nbsp; 0.50<br>&nbsp; Ni&nbsp; 0.00&nbsp; 0.00&nbsp; 0.75<br>&nbsp; Ni&nbsp; 0.50&nbsp; 0.50&nbsp; 0.75<br>K_POINTS automatic<br>5 5 1 0 0 0<br>EOF<br>echo "&nbsp; SCF processing...\c"<br>mpirun -n $1 pw.x -npool $2 &lt;TiC.scf.in &gt;TiC.scf.out<br>echo "done"<br><br>cat &gt;TiC.pp.in &lt;&lt;EOF<br>&amp;inputpp<br>&nbsp;&nbsp; prefix = 'TiC_Ni' ,<br>&nbsp;&nbsp; outdir = '$TempDIR/' ,<br>&nbsp;&nbsp; filplot = 'edensity' ,<br>&nbsp;&nbsp; plot_num = 0 ,<br>/<br>&amp;plot<br>&nbsp;&nbsp; nfile = 1 ,<br>&nbsp;&nbsp; filepp(1) = 1 ,<br>&nbsp;&nbsp; iflag = 2 ,<br>&nbsp;&nbsp; output_format = 3 ,<br>&nbsp;&nbsp; fileout = 'TiC_Ni_cd.xsf' ,<br>&nbsp;&nbsp; e1(1)=0.0,e1(2)=1.0,e1(3)=0.0,<br>&nbsp;&nbsp; e2(1)=0.0,e2(2)=0.0,e2(3)=1.0,<br>&nbsp;&nbsp; x0(1)=0.5,x0(2)=0.0,x0(3)=0.0,<br>&nbsp;&nbsp; nx=52,ny=52,<br>/<br>EOF<br>echo "&nbsp; PP processing...\c"<br>mpirun -n $1 pp.x -npool $2 &lt;TiC.pp.in &gt;TiC.pp.out<br>echo "done"</font><br>Is this problem cased by huge size of charge density file? How could I overcome this problem?<br>Best Regards.<br>        
<br><br><div>--<br>GAO Zhe<br>CMC Lab, MSE, SNU, Seoul, S.Korea<br>        
</div><br><br><span title="neteasefooter"><span id="netease_mail_footer"><hr>
<a href="http://mail.163.com/html/110414_attachment/att1.htm" target="_blank">体验网易邮箱2G超大附件,轻松发优质大电影、大照片,提速3倍!</a>
</span></span></blockquote><br><br><span title="neteasefooter"><span id="netease_mail_footer"><hr/>
<a href="http://mail.163.com/html/110414_attachment/att1.htm" target="_blank">体验网易邮箱2G超大附件,轻松发优质大电影、大照片,提速3倍!</a>
</span></span>