24小时热门版块排行榜    

查看: 558  |  回复: 2

yongleli

木虫 (正式写手)

[交流] 【分享】lev00处理vasp5.2输出文件问题 已有2人参与

最近做了点VASP计算,用lev00处理结果。
但是,一旦根据提示输入POSCAR,程序就提示我
FATAL! File POSCAR is bad or absent!

怎么回事呢?
打开read_vasp_geom.f90,赫然发现:
good=.true.
    write(*,*)'Reading geometry from '//trim(nameg)//' ...'
    open (15,file=trim(nameg),status='old',FORM='FORMATTED',err=200)
    read (15,*)
    read (15,*) vscale
    read (15,*,err=200,end=200) DIRC(1,1), DIRC(1,2), DIRC(1,3)
    read (15,*,err=200,end=200) DIRC(2,1), DIRC(2,2), DIRC(2,3)
    read (15,*,err=200,end=200) DIRC(3,1), DIRC(3,2), DIRC(3,3)
    DIRC=DIRC*vscale
!
!_________ reciprocal vectors BCELL without 2*pi
    call bastr(DIRC,BCELL,VOLC,0)
!_________ reciprocal vectors RECC with 2*pi
    call bastr(DIRC,RECC,VOLC,1)
!_________ numbers of atoms NspN in every species  
!
??????????????????????????????????????????
    read(15,*,err=200) (NspN(i),i=1,NSPEC)

因为VASP5.0的CONTCAR会自动增加一行原子名称
(我这步是单点,用的上一步CONTCAR做成的POSCAR),
而lev00没有“与时俱进”,读完了三个晶格矢量紧接着就去找原子个数了。
所以,这里只要在???????????????????????这一行,增加一行
read(15,*,err=200)
重新编译,就可以了。

另外,lev00输出的CHGCAR跟原来的格式有出入,
想要成功使用VMD看CHGCAR的话,
需要手动修改前边6行,即title、scalefactor、晶胞格矢、原子个数、Cart/Direct
部分,另外需要在原子坐标后边增加一行空行。

但愿本信息对用VASP的虫子们有点用处!

[ Last edited by yongleli on 2010-12-6 at 17:05 ]
回复此楼

» 猜你喜欢

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

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

cenwanglai

荣誉版主 (知名作家)

老和山猥琐派九段

优秀版主


zzy870720z(金币+1): 2011-04-15 14:31:01
绝对有用,呵呵。
2楼2010-12-06 17:25:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youzhizhe

荣誉版主 (职业作家)

有志者

优秀版主优秀版主

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖交流
zzy870720z(金币+1): 谢谢提示 2011-04-15 14:31:48
目前已经有最新的lev00 3.22版本支持VASP5.2了  可以选择最新的版本
微博http://weibo.com/10986069
3楼2011-04-15 11:25:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 yongleli 的主题更新
普通表情 高级回复 (可上传附件)
信息提示
请填处理意见