| 查看: 2592 | 回复: 22 | ||||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | ||||
wo8165827金虫 (小有名气)
|
[交流]
【求助】帮忙编个用于VASP小程序 已有6人参与
|
|||
|
帮忙编个用于VASP小程序，我的NSW设置为10，算完一次后 cp CONTCAR POSCAR 在算第二次，知道到达精度为止。往往都要 输上好几次。因此想编个小程序来替代重复的工作。 要求是开始计算时就只有INCAR POSCAR POTCAR KPOINTS ，输入命令VASP>log1& 计算，当计算完一次后输入命令 cp CONTCAR POSCAR 和VASP>log2& 再算第二次，直到输出文件logi提示： reached required accuracy - stopping structural energy minimisation 计算结束。 同时每次计算输出文件logi将保留。并输入一个Summary 记录出每一步算出的能量。 先谢谢大家了。。。。。 |
回复此楼» 收录本帖的淘帖专辑推荐
 vasp |
» 猜你喜欢
307求调剂
已经有13人回复
-
071010 323 分求调剂
已经有3人回复
-
材料调剂
已经有9人回复
-
材料专硕调剂
已经有10人回复
-
一志愿南航 335分 | 0856材料化工 | GPA 4.07 | 有科研经历
已经有8人回复
-
求调剂
已经有8人回复
-
284求调剂
已经有15人回复
-
085602化工求调剂（331分）
已经有3人回复
-
一志愿南昌大学324求调剂
已经有8人回复
-
材料科学与工程调剂
已经有5人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于VASP里面一些linux下小程序的使用
已经有4人回复
- VASP会用到一些编程程序么？？？
已经有9人回复
- VASP中脚本程序问题
已经有14人回复
- 大家帮忙看下这个labview小程序哪里有问题？
已经有5人回复
- 求一个简单的读写文件的小程序。。。
已经有4人回复
- 测试POSCAR的脚本不运行
已经有10人回复
- 求处理vasp能带的小程序
已经有14人回复
- 求助用matlab编一个小程序
已经有8人回复
- 编译vasp的lib时就出错，求指点
已经有22人回复
- 求助：vasp优化出错！！！
已经有14人回复
- 求助：vasp结构优化不能收敛
已经有18人回复
- VASP的源代码是公开的吗？还有其他可以读懂的程序吗？
已经有12人回复
- 编译VASP报错，求解释！
已经有16人回复
- 写了一个fortran90的小程序，编译通不过，请大侠帮忙
已经有59人回复
- 【求助】VASP计算的KPOINTS 设置，求助
已经有10人回复
- 【求助】VASP5.2 编译出错make: *** [vasp] 错误 1,请大家帮忙分析分析,谢谢
已经有7人回复
- 【求助成功】求ubuntu中vasp程序安装最后一步出错！！郁闷！！！
已经有12人回复
- 【求助】用mathematic编个二分法的程序
已经有3人回复
- 【求助】帮忙VASP写个小程序
已经有21人回复
tjyl
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3218.1
- 红花: 2
- 帖子: 576
- 在线: 156.6小时
- 虫号: 765184
- 注册: 2009-05-07
- 专业: 生物无机化学
- 管辖: 程序语言
★
gavinliu7390(金币+1):长见识了！ 2010-04-24 23:47
wo8165827(金币+7):谢谢 2010-04-25 08:59
wo8165827(金币+5):谢谢， 2010-04-25 09:02
gavinliu7390(金币+1):长见识了！ 2010-04-24 23:47
wo8165827(金币+7):谢谢 2010-04-25 08:59
wo8165827(金币+5):谢谢， 2010-04-25 09:02
|
不用VASP，不过写这个还是比较容易的。 #!/bin/bash OK=0 try=1 declare -a energy char='reached required accuracy - stopping structural energy minimisation' while [ $OK -eq 0 ] do vasp > log${try} & wait energy[$try]=$(......) #自己写提取出能量值吧 OK=$(grep $char log${try} |wc -l ) cp CONTCAR POSCAR try=$(( try+1 )) done i=1 while [ $i -lt $try ] do echo -e "$i \t ${energy[$i]} " i=$(( i+1 )) done =========================== 自己写提取总能的那点吧 grep，cut,sed,tail,head,awk 用这几个总是足够了 [ Last edited by tjyl on 2010-4-24 at 23:46 ] |
3楼2010-04-24 23:44:07
y1ding
铁杆木虫 (著名写手)
- 1ST强帖: 1
- 应助: 61 (初中生)
- 贵宾: 0.33
- 金币: 5959.3
- 散金: 1
- 红花: 21
- 帖子: 1884
- 在线: 491.1小时
- 虫号: 142265
- 注册: 2005-12-21
- 专业: 凝聚态物性 II ：电子结构
2楼2010-04-24 21:10:49
gavinliu7390
木虫 (著名写手)
叶落鹰飞
- 1ST强帖: 3
- 应助: 51 (初中生)
- 贵宾: 0.481
- 金币: 1564.7
- 散金: 2975
- 红花: 34
- 帖子: 1278
- 在线: 241.4小时
- 虫号: 675997
- 注册: 2008-12-17
- 性别: GG
- 专业: 凝聚态物性I:结构、力学和
4楼2010-04-24 23:44:19
tjyl
金虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3218.1
- 红花: 2
- 帖子: 576
- 在线: 156.6小时
- 虫号: 765184
- 注册: 2009-05-07
- 专业: 生物无机化学
- 管辖: 程序语言
5楼2010-04-24 23:47:44
