| 查看: 3168 | 回复: 17 | |||||
| 当前主题已经存档。 | |||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||||
[交流]
【求助】VASP脚本程序怎么运行？
|
|||||
|
刚学VASP，对其还不是很熟悉；现在正在试着运行脚本文件（按侯老师的例子），准备好INCAR、KPOINTS、POTCAR文件后，应该准备POSCAR文件，那脚本程序run_a0怎样运行？脚本run_a0为： #！/bin/sh rm WAVECAR for i in 3.80 3.85 3.90 3.95 4.00 4.05 4.10 4.15 4.20 4.25 4.30 do cat >POSCAR < Al-fcc $i 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 Direct 0.0 0.0 0.0 ! echo "a =$i angstrom";time vasp E=‘grep "TOTEN" OUTCAR |tail -1|awk ‘{printf "%12.6f \n",$5}'‘ V=‘grep "volume" OUTCAR |tail -1|awk ‘{printf "%12.4f \n",$5}'‘ echo $V $E >>EtVo.dat done VASP运行需四个文件，这里难道用run_a0替代POSCAR文件？但那样的话VASP运行不了，但又不能单独运行脚本程序（出不了结果），到底该怎麽办？哪位能告诉我具体的运行步骤？谢谢！ [ Last edited by aylayl08 on 2010-4-21 at 09:01 ] |
回复此楼» 收录本帖的淘帖专辑推荐
 Photochemistry | 计算-vasp |
» 猜你喜欢
一志愿华中农业071010，320求调剂
已经有19人回复
-
304求调剂
已经有7人回复
-
求博导｜生物质基多孔碳/超级电容方向，已有相关成果，寻能源材料/碳材料方向老师
已经有3人回复
-
二苯甲酮酸类衍生物
已经有6人回复
-
接受任何调剂
已经有4人回复
-
一志愿中科大材料与化工，353分还有调剂学校吗
已经有12人回复
-
320求调剂
已经有4人回复
-
295分求调剂
已经有5人回复
-
294求调剂
已经有8人回复
-
山东双非院校考核超级无底线，领导幸灾乐祸，教师遭殃恐
已经有8人回复
sars518
木虫 (著名写手)
小木虫灌水大队队长
- 1ST强帖: 6
- 应助: 29 (小学生)
- 贵宾: 0.119
- 金币: 3360.2
- 散金: 4434
- 红花: 32
- 帖子: 2641
- 在线: 2213.2小时
- 虫号: 539209
- 注册: 2008-04-03
- 性别: GG
- 专业: 科学社会主义
3楼2009-05-09 16:06:30
2楼2009-05-09 10:54:49
tjpm
金虫 (正式写手)
- 1ST强帖: 1
- 应助: 0 (幼儿园)
- 金币: 1025.7
- 散金: 10
- 红花: 1
- 帖子: 367
- 在线: 19.3小时
- 虫号: 717334
- 注册: 2009-03-07
- 专业: 凝聚态物性 II ：电子结构
★ ★
fegg7502(金币+2,VIP+0):thank you very much! 5-10 18:41
fegg7502(金币+2,VIP+0):thank you very much! 5-10 18:41
|
其实shell脚本很简单的，有C语言基础随便看看就会了，当然还需要熟悉常用的linux命令。 简单写个注释吧 ＝＝＝＝ #！/bin/sh ＃每种shell脚本开头需要指定起类型，比如unix下原始的就是sh,linux都标配bash了， Solaris也是用bash了。功能上现在似乎zsh比较好，目前我是把zsh做默认shell. rm WAVECAR 删除WAVECAR，其实这里有点不好，因为很多系统默认把rm 替换为rm -i 也就是说你要删除这个文件时会跳出一个提示问是否要删除。建议改为 rm WAVECAR -rf for i in 3.80 3.85 3.90 3.95 4.00 4.05 4.10 4.15 4.20 4.25 4.30 do 这个sh/bash下典型的for循环，就是让i逐次取后面的值 cat >POSCAR < Al-fcc $i 0.0 0.5 0.5 0.5 0.0 0.5 0.5 0.5 0.0 Direct 0.0 0.0 0.0 ! 这个就是利用管道创建文件POSCAR，两个“！”分别标识开头和结尾。 echo "a =$i angstrom";time vasp 屏幕输出 a的值，并运行vasp, time命令加在前面就是计算其运行时间。 E=‘grep "TOTEN" OUTCAR |tail -1|awk ‘{printf "%12.6f \n",$5}'‘ 搜索包含"TOTEN"的所有行，利用tail -1取最后一行，然后格式化输出第5个字符段。 V=‘grep "volume" OUTCAR |tail -1|awk ‘{printf "%12.4f \n",$5}'‘ 这个同上 echo $V $E >>EtVo.dat 把V，E的值添加到EtVo.dat done 表示结束for循环。 其实这些工具脚本都很简单的，关键还是要明白输出文件，然后就可以自己写一些适合自己用的。 |
4楼2009-05-09 20:20:29
6楼2009-05-09 21:13:11
