<div style="line-height:1.7;color:#000000;font-size:14px;font-family:arial">I checked the code of bands_FS.f90, and found like this:<br><span style="color: rgb(0, 0, 255);"><span style="color: rgb(0, 0, 255);">        open(12,file=<span style="color: rgb(255, 0, 0);">'input_FS'</span>)</span><br><br><span style="color: rgb(0, 0, 255);">        read(12,*) n_start, n_last</span><br><span style="color: rgb(0, 0, 255);">        read(12,*) E_fermi</span><br><span style="color: rgb(0, 0, 255);">        read(12,*) sysname</span><br><span style="color: rgb(0, 0, 255);">        read(12,*) na,nb, nc</span><br><span style="color: rgb(0, 0, 255);">        read(12,*) x(1),x(2),x(3)</span><br><span style="color: rgb(0, 0, 255);">        read(12,*) y(1),y(2),y(3)</span><br><span style="color: rgb(0, 0, 255);">        read(12,*) z(1),z(2),z(3)</span></span><br>So I am just wondering whether it can work well, if you change your input file name as "input_FS", and run only by the command bands_FS.x < CWrium.fs.in > CWrium.fs.out, where CWrium.fs.in is a nscf calculation input file with the K_Points obtained by kvec_FS.x.<br>The more ditailed process is in example 08.<br><div>--<br>GAO Zhe<br>CMC Lab, Materials Science & Engineering Department,<br>Seoul National University, South Korea<br>        
 </div><div id="divNeteaseMailCard"></div><br>At 2012-03-02 11:08:12,"Chan-Woo Lee" <cwandtj@gmail.com> wrote:<br> <blockquote id="isReplyContent" style="PADDING-LEFT: 1ex; MARGIN: 0px 0px 0px 0.8ex; BORDER-LEFT: #ccc 1px solid"><div class="WordSection1"><p class="MsoNormal">Dear QE developers and users,<o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">I am trying to get Fermi surface using kvecs_FS.x and bands_FS.x. I have followed example08 and have succeeded to get K-points using kvecs_FS.x. However, I got following error message when I execute bands_FS.x:<o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">At line 230 of file bands_FS.f90<o:p></o:p></p><p class="MsoNormal">Fortran runtime error: Bad real number in item 13 of list input<o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">Line 230 of bands_FS.f90 is about reading spin-up energiesĄ­. But I have no idea why the code is not working. FYI, I added my input files for kvecs_FS.x and bands_FS.x:<o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">kvecs_FS.in:<o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">0.317093  0.00  0.00<o:p></o:p></p><p class="MsoNormal">0.00  0.317093  0.00<o:p></o:p></p><p class="MsoNormal">0.00  0.00  0.317093<o:p></o:p></p><p class="MsoNormal">12  12  12<o:p></o:p></p><p class="MsoNormal">CWrium<o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">input_FS.in:<o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">10 14<o:p></o:p></p><p class="MsoNormal">22.8158<o:p></o:p></p><p class="MsoNormal">CWrium<o:p></o:p></p><p class="MsoNormal">12 12 12<o:p></o:p></p><p class="MsoNormal">0.317093  0.00  0.00<o:p></o:p></p><p class="MsoNormal">0.00  0.317093  0.00<o:p></o:p></p><p class="MsoNormal">0.00  0.00  0.317093<o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">Here, 10 14 are min. and max. band no. (no. of bands for my SCF calc. is 20), and 22.8158 is Fermi energy. <o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">Again, your help will truly be appreciated. I ccĄŻed Dr. Eyvaz Isaev who is the author of the two codes. <o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">Sincerely,<o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal">Chan-Woo<o:p></o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal"><o:p> </o:p></p><p class="MsoNormal"><i><span style="font-size:10.0pt;font-family:"Arial","sans-serif"">-------<br>Chan-Woo Lee, Ph.D.<o:p></o:p></span></i></p><p class="MsoNormal"><i><span style="font-size:10.0pt;font-family:"Arial","sans-serif"">Postdoctoral Research Associate<o:p></o:p></span></i></p><p class="MsoNormal"><i><span style="font-size:10.0pt;font-family:"Arial","sans-serif""><o:p> </o:p></span></i></p><p class="MsoNormal"><i><span style="font-size:10.0pt;font-family:"Arial","sans-serif"">Department of Chemistry<br>University of Pennsylvania<br>231 South 34th Street<br>Philadelphia, PA 19104-6323 <br>Phone: 1-215-898-3564 (Office)<br><br></span></i><span style="font-size:10.0pt;font-family:"Arial","sans-serif""><o:p></o:p></span></p><p class="MsoNormal"><o:p> </o:p></p></div></blockquote></div><br><br><span title="neteasefooter"><span id="netease_mail_footer"></span></span>