<div>hi Paolo:</div>
<div>&nbsp;</div>
<div>&nbsp;</div>
<div>oh .there is a file named&nbsp;&nbsp;dynmatrix.f90 .and i got confused .</div>
<div>i list my in file and my out file here :</div>
<div>&amp;input fildyn=&#39;/raid/xyh/pwscf/*/*.dynG&#39;,asr=&#39;one-dim&#39;,q(1)=0,q(2)=0,q(3)=12/<br>AND THE OUTFILE HERE:</div>
<div>
<p>Thu Oct 25 11:06:00 CST 2007</p>
<p>&nbsp;&nbsp;&nbsp;&nbsp; Reading Dynamical Matrix from file /raid/xyh/pwscf/*/*.dynG<br>asr rotation axis in 1D system=&nbsp;&nbsp;&nbsp; 3<br>&nbsp;&nbsp;&nbsp;&nbsp; Acoustic Sum Rule: || Z*(ASR) - Z*(orig)|| =&nbsp;&nbsp;&nbsp; 0.815659E-01<br>&nbsp;&nbsp;&nbsp;&nbsp; Acoustic Sum Rule: ||dyn(ASR) - dyn(orig)||=&nbsp;&nbsp;&nbsp; 
0.697493E-01</p>
<p>&nbsp;&nbsp;&nbsp;&nbsp; Polarizability (A^3 units)<br>&nbsp;&nbsp;&nbsp;&nbsp; multiply by 0.971817 for Clausius-Mossotti correction<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 5.625414&nbsp;&nbsp;&nbsp; 0.000658&nbsp;&nbsp;&nbsp; 0.000133<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000658&nbsp;&nbsp;&nbsp; 5.624128&nbsp;&nbsp;&nbsp; 0.000114<br>&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.000133&nbsp;&nbsp;&nbsp; 0.000114&nbsp;&nbsp;&nbsp; 
8.860593</p>
<p>&nbsp;&nbsp;&nbsp;&nbsp; IR cross sections are in (D/A)^2/amu units</p>
<p>#&nbsp; mode&nbsp;&nbsp; [cm-1]&nbsp;&nbsp;&nbsp;&nbsp; [THz]&nbsp;&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; IR<br>&nbsp;&nbsp;&nbsp; 1&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp; 0.0000<br>&nbsp;&nbsp;&nbsp; 2&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp; 0.0000<br>&nbsp;&nbsp;&nbsp; 3&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp; 0.0000<br>&nbsp;&nbsp;&nbsp; 4&nbsp;&nbsp;&nbsp;&nbsp;&nbsp; 0.00&nbsp;&nbsp;&nbsp; 0.0000&nbsp;&nbsp;&nbsp; 0.0000<br>&nbsp;&nbsp;&nbsp; 5&nbsp;&nbsp;&nbsp;&nbsp; 49.16&nbsp;&nbsp;&nbsp; 1.4738
 &nbsp;&nbsp;&nbsp; 0.0064<br>&nbsp;&nbsp;&nbsp; 6&nbsp;&nbsp;&nbsp;&nbsp; 63.88&nbsp;&nbsp;&nbsp; 1.9151&nbsp;&nbsp;&nbsp; 0.0043<br>&nbsp;&nbsp;&nbsp; 7&nbsp;&nbsp;&nbsp;&nbsp; 68.56&nbsp;&nbsp;&nbsp; 2.0554&nbsp;&nbsp;&nbsp; 0.0016<br>&nbsp;&nbsp;&nbsp; 8&nbsp;&nbsp;&nbsp; 106.29&nbsp;&nbsp;&nbsp; 3.1865&nbsp;&nbsp;&nbsp; 0.0002<br>&nbsp;&nbsp;&nbsp; 9&nbsp;&nbsp;&nbsp; 156.91&nbsp;&nbsp;&nbsp; 4.7040&nbsp;&nbsp;&nbsp; 0.0002<br>&nbsp;&nbsp; 10&nbsp;&nbsp;&nbsp; 252.93&nbsp;&nbsp;&nbsp; 7.5826&nbsp;&nbsp;&nbsp; 0.2436<br>&nbsp;&nbsp; 11&nbsp;&nbsp;&nbsp; 254.66
 &nbsp;&nbsp;&nbsp; 7.6344&nbsp;&nbsp;&nbsp; 0.4708<br></p><br>&nbsp;and&nbsp; there is not the vector at all!!</div>
<div>and the job script:</div>
<div>
<p>#!/bin/bash<br>#<br>#$ -cwd<br>#$ -j y<br>#$ -S /bin/bash<br>#</p>
<p>MPI_DIR=/home/software/mpich/intel/bin</p>
<p>date</p>
<p>$MPI_DIR/mpirun -np $NSLOTS -machinefile $TMPDIR/machines -nolocal /home/xyh/intel_espresso-3.2/bin/dynmat.x&lt;IR.in</p>
<p>date<br></p></div>
<div>and what is wrong ? i do not know</div>
<div>please give me a detailed advise </div>
<div>thank you </div><br clear="all"><br>-- <br>Xu Yuehua<br>physics Department of Nanjing university<br>China