<div dir="ltr"><div>Hi<br><br></div>I forgot to notice that you also need to modify the line<br><div class="gmail_extra"><div class="gmail_quote">   USE cell_base, ONLY : ibrav, celldm<br></div><div class="gmail_quote">as follows:<br><div class="gmail_quote">   USE cell_base, ONLY : ibrav, celldm, at<br><br></div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
I calculated ir spectra of a supercell containing 192 atoms and modified<br>
the dyn file by adding basis vector as follow [...]<br>
The first two numbers do not have a space. Is this a bug?<br></blockquote><div><br></div><div>it is. More exactly: when that code was originally written (20 years ago or so), it was unthinkable to perform such a calculation on more than 99 atoms. Just modify line 39 in PHonon/Gamma/writedyn.f90:<br>        WRITE(iudyn, '(2i3)') na, nb<br></div><div>as<br>        WRITE(iudyn, '(2i5)') na, nb<br></div><div>In order to read your current matrix, if you do not want to manually insert spaces, you may modify, in the routine that reads the dynamical matrix (PHonon/PH/dynmat.f90, routine readmat) the line that reads atomic numbers:<br>           read (1,*) naa, nbb<br></div><div>as<br>           read (1,'(2A3)') naa, nbb<br><br></div><div>Thank you for reporting this<br></div><div>
<br></div><div>Paolo<br></div><div>
_______________________________________________<br>
Pw_forum mailing list<br>
<a href="mailto:Pw_forum@pwscf.org">Pw_forum@pwscf.org</a><br>
<a href="http://pwscf.org/mailman/listinfo/pw_forum" rel="noreferrer" target="_blank">http://pwscf.org/mailman/listinfo/pw_forum</a><br>
</div></div><br><br clear="all"><br>-- <br><div class="gmail_signature"><div dir="ltr"><div><div dir="ltr"><div>Paolo Giannozzi, Dip. Scienze Matematiche Informatiche e Fisiche,<br>Univ. Udine, via delle Scienze 208, 33100 Udine, Italy<br>Phone +39-0432-558216, fax +39-0432-558222<br><br></div></div></div></div></div>
</div></div>