Problem in pp.py script in Pb example folder
Posted: Sun Aug 28, 2016 10:37 pm
Dear moderator,
I would like to point out one probable error I encountered while using "python pp.py" before performing scf and nscf calculations in the following example folder for Lead (pb): espresso-5.4.0/EPW/examples/pb/woSOC/phonons
After running ../../../../bin/ph.x < ph.in >& ph.out &, I obtain following files: pb.dyn1, pb.dyn2, pb.dyn3 ....... pb.dyn16.
To copy these files into the save folder before proceeding to epw calculations, I use the script pp.py (located in the same folder as the above *.dyn files). But it resulted in error something like "cp: cannot stat `pb.dyn1.xml': No such file or directory". This error is self-explanatory as the files to be copied into the save folder are named as pb.dyn.iqpt (like pb.dyn1, pb.dyn2 and so on) without ".xml" extension.
Currently the following command is being used in pp.py script.
I think the ".xml" part should be removed from the script.
Please correct me if I am wrong here.
Secondly, I noticed that in Fig. 4 on the documentation page (http://epw.org.uk/Documentation/FCCLead), transport coupling strength \lambda_{tr} is thin dashed blue line (according to the figure caption) but in the figure it is located above the red line. Although it does not matter much, but looks a big ambiguous in the figure.
Thanks and best wishes
Ankit
I would like to point out one probable error I encountered while using "python pp.py" before performing scf and nscf calculations in the following example folder for Lead (pb): espresso-5.4.0/EPW/examples/pb/woSOC/phonons
After running ../../../../bin/ph.x < ph.in >& ph.out &, I obtain following files: pb.dyn1, pb.dyn2, pb.dyn3 ....... pb.dyn16.
To copy these files into the save folder before proceeding to epw calculations, I use the script pp.py (located in the same folder as the above *.dyn files). But it resulted in error something like "cp: cannot stat `pb.dyn1.xml': No such file or directory". This error is self-explanatory as the files to be copied into the save folder are named as pb.dyn.iqpt (like pb.dyn1, pb.dyn2 and so on) without ".xml" extension.
Currently the following command is being used in pp.py script.
Code: Select all
os.system('cp '+prefix+'.dyn'+str(iqpt)+'.xml save/'+prefix+'.dyn_q'+label+'.xml')
I think the ".xml" part should be removed from the script.
Please correct me if I am wrong here.
Secondly, I noticed that in Fig. 4 on the documentation page (http://epw.org.uk/Documentation/FCCLead), transport coupling strength \lambda_{tr} is thin dashed blue line (according to the figure caption) but in the figure it is located above the red line. Although it does not matter much, but looks a big ambiguous in the figure.
Thanks and best wishes
Ankit