| 查看: 1362 | 回复: 5 | |||
信任开花铁杆木虫 (知名作家)
|
[求助]
在vasp官方论坛看到一段如下代码，不知道如何使用，请各位指教
|
|
就是读取二进制wavecar文件到，具体如下： For reading the header of the WAVECAR file, for example, you could do something like this: #!/usr/bin/python import sys def isprintable(char): return 0x256 <= char <= 0x16f def string(filename): data = open("/Working/Eclipse/plotting/WAVECAR", "rb"  .read()count = 0 line = "" for ch in data: if isprintable(ch): count += 1 line = line + ch else: if count > 1 : print line count = 0 line= "" print line 装啦python，就是不知道这几行程序怎么对齐才能得到正确到wavecar到输出，期待帮助。 |
回复此楼» 猜你喜欢
焦虑
已经有8人回复
-
308求调剂
已经有4人回复
-
NSFC申报书里申请人简历中代表性论著还需要在申报书最后的附件里面再上传一遍吗
已经有14人回复
-
材料与化工一志愿南昌大学327求调剂推荐
已经有6人回复
-
化学调剂0703
已经有7人回复
-
327求调剂
已经有11人回复
-
调剂
已经有8人回复
-
梁成伟老师课题组欢迎你的加入
已经有7人回复
-
伙伴们，祝我生日快乐吧
已经有24人回复
-
中科院材料273求调剂
已经有3人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 请教各位高手，VASP算电荷有序需要设哪些参数啊？
已经有10人回复
- vasp并行编译后出现错误，请大家指教
已经有23人回复
- 如何在vasp中设定离子的价态
已经有18人回复
- VASP的源代码是公开的吗？还有其他可以读懂的程序吗？
已经有12人回复
- 如何用VASP计算表面缺陷
已经有4人回复
- VASP中如何在体系中引入电子和空穴
已经有15人回复
- 当体系中包含两种材料时，如何用VASP计算其中一种材料的功函？
已经有11人回复
- VASP计算出来的PDOS怎么画图啊？求指教
已经有5人回复
- 【求助完结】vasp优化结构运行一段时间后没有出现错误，但是停止了！
已经有20人回复
- 【求助】如何使用VASP计算charged defect的形成焓
已经有12人回复
- 【求助】vasp优化结构未结束，log文件出现下列问题停止计算，请高手指教
已经有6人回复
- 【求助】GaN的费米能是多少？要用VASP计算GaN，不知道算的对不对，请知道的给个参考值
已经有4人回复
- 【求助】vasp 源代码怎么学习
已经有12人回复
- 【求助】请教vasp如何使本征值按每条能带输出
已经有18人回复
- 【求助】如何使用VASP画费米面。
已经有13人回复
- 【讨论】讨论使用vasp超软赝势参数生成pwscf超软赝势文件的可行性
已经有6人回复
后天一
木虫 (小有名气)
我们的自由！
- 1ST强帖: 3
- 应助: 9 (幼儿园)
- 金币: 2544.6
- 红花: 4
- 帖子: 204
- 在线: 311小时
- 虫号: 519608
- 注册: 2008-03-06
- 专业: 凝聚态物性 II ：电子结构
【答案】应助回帖
★
uuv2010(金币+1): 欢迎参加交流 2011-09-07 20:26:54
信任开花(金币+10): 非常感谢，这段程序我也反复对齐测试，结果就是读不出里面到数据，如果不if逻辑判断，直接打印，就是乱码，如果if逻辑判断，则因为isprintable结果都是false，即不可打印，故什么也不输出。不管怎么样，还是感谢热心帮助 2011-09-08 12:43:19
uuv2010(金币+1): 欢迎参加交流 2011-09-07 20:26:54
信任开花(金币+10): 非常感谢，这段程序我也反复对齐测试，结果就是读不出里面到数据，如果不if逻辑判断，直接打印，就是乱码，如果if逻辑判断，则因为isprintable结果都是false，即不可打印，故什么也不输出。不管怎么样，还是感谢热心帮助 2011-09-08 12:43:19
|
#!/usr/bin/python import sys def isprintable(char): return 0x256 <= char <= 0x16f def string(filename): data = open("/Working/Eclipse/plotting/WAVECAR", "rb".read() count = 0 line = "" for ch in data: if isprintable(ch): count += 1 line = line + ch elseif count > 1 : print line count = 0 line= "" print line 试试吧  不确定对，汗~ |
2楼2011-09-07 16:48:57
bigseahorse
禁虫 (正式写手)
|
本帖内容被屏蔽 |
3楼2011-09-07 17:02:19
4楼2011-09-07 20:54:56
dxcharlary
专家顾问 (职业作家)
-
专家经验: +83 - 1ST强帖: 2
- 应助: 251 (大学生)
- 金币: 10417.6
- 散金: 155
- 红花: 65
- 帖子: 3137
- 在线: 757.4小时
- 虫号: 936718
- 注册: 2010-01-04
- 专业: 凝聚态物性I:结构、力学和
- 管辖: 第一性原理
5楼2011-09-08 11:17:46
信任开花
铁杆木虫 (知名作家)
- 应助: 0 (幼儿园)
- 贵宾: 0.02
- 金币: 7745.2
- 散金: 707
- 帖子: 7943
- 在线: 212.6小时
- 虫号: 934954
- 注册: 2009-12-29
- 性别: MM
6楼2011-09-08 12:43:25
