24小时热门版块排行榜

返回列表

信任开花

铁杆木虫 (知名作家)

应助: 0 (幼儿园)
贵宾: 0.02
金币: 7745.2
散金: 707
帖子: 7943
在线: 212.6小时
虫号: 934954
注册: 2009-12-29
性别: MM

[求助] 在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人回复
博士读完未来一定会好吗已经有23人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有11人回复
读博已经有4人回复
JMPT 期刊投稿流程已经有4人回复
心脉受损已经有5人回复
Springer期刊投稿求助已经有4人回复
小论文投稿已经有3人回复
申请2026年博士已经有6人回复

» 本主题相关价值贴推荐，对您同样有帮助:

请教各位高手，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人回复

从不温柔，从不体贴，从不讲理

1楼2011-09-07 15:25:35

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

后天一

木虫 (小有名气)

我们的自由！

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

#!/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

试试吧

不确定对，汗~