| 查看: 1367 | 回复: 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到输出，期待帮助。 |
回复此楼» 猜你喜欢
281求调剂（0805）
已经有9人回复
-
293求调剂
已经有14人回复
-
0703化学336分求调剂
已经有6人回复
-
304求调剂
已经有9人回复
-
环境领域全国重点实验室招收博士1-2名
已经有3人回复
-
材料专硕306英一数二
已经有10人回复
-
301求调剂
已经有6人回复
-
一志愿天津大学化学工艺专业（081702）315分求调剂
已经有7人回复
-
302求调剂
已经有6人回复
-
26博士申请
已经有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 (大学生)
- 金币: 10418.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
