24小时热门版块排行榜    

CyRhmU.jpeg
查看: 1243  |  回复: 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到输出,期待帮助。
回复此楼

» 猜你喜欢

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

从不温柔,从不体贴,从不讲理
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dxcharlary

专家顾问 (职业作家)


youzhizhe(金币+1): 谢谢交流。 2011-09-08 16:47:00
二楼应该没问题,python是以对行来表示子程序。就像C++里面的{}
5楼2011-09-08 11:17:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

后天一

木虫 (小有名气)

我们的自由!

【答案】应助回帖


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
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

信任开花

铁杆木虫 (知名作家)

非常感谢,这段程序我也反复对齐测试,结果就是读不出里面到数据,如果不if逻辑判断,直接打印,就是乱码,如果if逻辑判断,则因为isprintable结果都是false,即不可打印,故什么也不输出。不管怎么样,还是感谢热心帮助
从不温柔,从不体贴,从不讲理
6楼2011-09-08 12:43:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见