<div dir="ltr"><div><div><div><div><div>Dear all,<br><br></div>I think there is a minor bug in the dynmat.f90 in the above versions of the phonon code. This bug is related to the reading of the crystal axis from "matdyn.xml" file (when one does a spin-orbit/non-collinear phonon calculation) and converting them in units of alat.<br><br></div>In line 170 of the code at(3,3) (read from the .xml file) is divided by a0=celldm(1). However, in the .xml file, "at" is already in units of alat. <br><br></div>This results in printing of wrong lattice vectors, when one tries to see the phonon modes in xsf. Commenting off line no 170 gives correct results. <br><br></div>With regards,<br><br></div>Prasenjit <br clear="all"><div><div><div><div><div><div><br>-- <br><div class="gmail_signature"><div dir="ltr">PRASENJIT GHOSH,<br>IISER Pune,<br>Dr. Homi Bhabha Road, Pashan<br>Pune, Maharashtra 411008, India<br><br>Phone: +91 (20) 2590 8203<br>Fax: +91 (20) 2589 9790</div></div>
</div></div></div></div></div></div></div>