24小时热门版块排行榜    

查看: 568  |  回复: 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 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[文学芳草园] 伙伴们,祝我生日快乐吧 +15 myrtle 2026-03-10 24/1200 2026-03-15 21:16 by 苏州_逗号
[考研] 材料工程专硕274一志愿211求调剂 +5 薛云鹏 2026-03-15 5/250 2026-03-15 20:38 by Logic2024
[考研] 311求调剂 +3 26研0 2026-03-15 3/150 2026-03-15 09:12 by JourneyLucky
[考研] 268求调剂 +5 一定有学上- 2026-03-14 6/300 2026-03-14 22:20 by 运气yunqi
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 297一志愿上交085600求调剂 +5 指尖八千里 2026-03-14 5/250 2026-03-14 17:26 by a不易
[考研] 复试调剂 +3 呼呼?~+123456 2026-03-14 3/150 2026-03-14 16:53 by WTUChen
[考研] 085600材料与化工 326 求调剂 +5 热爱生活ing 2026-03-09 5/250 2026-03-14 02:39 by JourneyLucky
[考研] 材料与化工 一志愿山大 321分 求调剂 +7 每天散步 2026-03-09 8/400 2026-03-14 02:18 by JourneyLucky
[考研] 279求调剂 +3 抓着星星的女孩 2026-03-10 3/150 2026-03-13 23:47 by userper
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 求调剂 +3 程雨杭 2026-03-12 3/150 2026-03-13 15:06 by JourneyLucky
[考研] 0856化学工程280分求调剂 +4 shenzxsn 2026-03-11 4/200 2026-03-13 11:55 by ymwdoctor
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 化工学硕306求调剂 +9 42838695 2026-03-12 9/450 2026-03-13 10:16 by houyaoxu
[考研] 293求调剂,一志愿陕师大生物学 +3 ??????.?.??? 2026-03-09 3/150 2026-03-11 10:02 by 学员8dgXkO
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 0856材料与化工353分求调剂 +11 NIFFFfff 2026-03-09 11/550 2026-03-10 18:36 by suyuanhai
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[考研] 家人们 调剂不迷路 看这里 +8 likeihood 2026-03-09 13/650 2026-03-10 08:09 by likeihood
信息提示
请填处理意见