<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>
<div><span class="gmail_quote">2 007/11/19, Paolo Giannozzi &lt;<a href="mailto:giannozz@nest.sns.it">giannozz@nest.sns.it</a>&gt;:</span>
<blockquote class="gmail_quote" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><br>On Nov 19, 2007, at 13:45 , xu yuehua wrote:<br><br>&gt; the mode has been recalculate,but there is only new frequency in
<br>&gt; outfile,<br>&gt; no new vibrational vector present which is my wanted. how can i get<br>&gt; it ?<br><br>see the header of file dynmat.f90<br>--<br>Paolo Giannozzi, Dept of Physics, University of Udine<br>via delle Scienze 208, 33100 Udine, Italy
<br>Phone +39-0432-558216, fax +39-0432-558222<br><br><br><br>_______________________________________________<br>Pw_forum mailing list<br><a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br><a href="http://www.democritos.it/mailman/listinfo/pw_forum">
http://www.democritos.it/mailman/listinfo/pw_forum</a><br></blockquote></div><br><br clear="all"><br>-- <br>lalala lalala