<div>i use dynmat.x to add lo-to split ,and want to get the new frequencies and modes, and&nbsp; only the frequencies is present .</div>
<div>my infile now:&nbsp;<br><br>&amp;input<br>fildyn=&#39;/raid/xyh/pwscf/fourringwater/pho_scf_relax_fromlqf_via_vaspz/*.dynG&#39;,<br>asr=&#39;one-dim&#39;,<br>q(1)=0,q(2)=0,q(3)=12,<br>filout= &#39;dynmat.out&#39;/<br>&nbsp;</div>

<div>&nbsp;</div>
<div>&nbsp;</div>
<div>but there is not dynmat.out .only the file like&nbsp;pwscf.sh.o8645<br clear="all"></div>
<div>and in pwscf.sh.o8645</div>
<div>&nbsp;</div>
<div><font color="#ff6666">there is :</font></div>
<p>ue Nov 20 17:29:58 CST 2007</p>
<p>&nbsp;&nbsp;&nbsp;&nbsp; Reading Dynamical Matrix from file /raid/xyh/pwscf/fourringwater/pho_scf_relax_fromlqf_via_vaspz/fourringwater.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>
<div>#&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>&nbsp;&nbsp; 12&nbsp;&nbsp;&nbsp; 257.93&nbsp;&nbsp;&nbsp; 7.7324&nbsp;&nbsp;&nbsp; 0.6958<br>&nbsp;&nbsp; 13&nbsp;&nbsp;&nbsp; 564.27&nbsp;&nbsp; 16.9163&nbsp;&nbsp;&nbsp; 0.4540<br>&nbsp;&nbsp; 14&nbsp;&nbsp;&nbsp; 631.40&nbsp;&nbsp; 18.9287&nbsp;&nbsp;&nbsp; 0.1605<br>&nbsp;&nbsp; 15&nbsp;&nbsp;&nbsp; 632.64&nbsp;&nbsp; 18.9659&nbsp;&nbsp;&nbsp; 1.7223<br>&nbsp;&nbsp; 16&nbsp;&nbsp;&nbsp; 633.03&nbsp;&nbsp; 18.9777&nbsp;&nbsp;&nbsp; 1.8109<br>&nbsp;&nbsp; 17&nbsp;&nbsp;&nbsp; 
705.34&nbsp;&nbsp; 21.1455&nbsp;&nbsp;&nbsp; 1.3172<br>&nbsp;&nbsp; 18&nbsp;&nbsp;&nbsp; 750.19&nbsp;&nbsp; 22.4900&nbsp;&nbsp;&nbsp; 3.3584<br>&nbsp;&nbsp; 19&nbsp;&nbsp;&nbsp; 750.80&nbsp;&nbsp; 22.5084&nbsp;&nbsp;&nbsp; 3.3466<br>&nbsp;&nbsp; 20&nbsp;&nbsp;&nbsp; 793.99&nbsp;&nbsp; 23.8030&nbsp;&nbsp;&nbsp; 0.0001<br>&nbsp;&nbsp; 21&nbsp;&nbsp;&nbsp; 854.32&nbsp;&nbsp; 25.6117&nbsp;&nbsp;&nbsp; 3.6181</div>
<div>&nbsp;</div>
<div>...........</div>
<div><br>&nbsp;</div>
<p>&nbsp;</p>
<div><font color="#ff6666">no mode</font>&nbsp;! .i am puzzled .and where is my error ?</div>
<div>&nbsp;</div>
<div><br>-- <br>Xu Yuehua<br>physics Department of Nanjing university<br>China </div>