<div dir="ltr">Dear Xcrysden community,<div><br></div><div>I want to render the Fermi surface of a system with a dense kmesh (no. of irreducible kpoints &gt; 10000) from WIEN2k. There was an error saying that the &#39;NKPT is too small&#39;.  So I followed a post on WIEN2k mailing list (<a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg15168.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg15168.html</a>) and did the following:</div><div><br></div><div>1. copied the new source code in Prof. Blaha&#39;s attachment in the above thread into $WIENROOT/SRC_kgen and recompile. The generated case.outputkgen file now shows the kpoint properly.  </div><div>2. In Xcrysden&#39;s source code, i change the NKP parameter in param.inc to 89000 (for NKP larger than this value, i got error during compilation).</div><div>3. replace wn_readbakgen.f with the one provided in this thread: </div><div><a href="http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg15219.html">http://www.mail-archive.com/wien@zeus.theochem.tuwien.ac.at/msg15219.html</a>.</div><div>4. rebuild the excutables by &quot;make all&quot;.</div><div><br></div><div>Then I go back to my calculation and run x kgen and render the fermi surface in Xcrysden. However, there was an error when executing wn_readbakgen which says</div><div><br></div><div><div>&quot;At line 136 of file wn_readbakgen.f (unit = 8, file = &#39;Case.outputkgen&#39;)</div><div>Fortran runtime error: Bad integer for item 1 in list input&quot;</div></div><div><br></div><div><br></div><div>Could someone please help?</div><div>Thanks in advance.</div><div><br></div><div>Best,</div><div>Fermin</div><div><br></div><div>------------------------------</div><div><br></div><div><br></div></div>