应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】vasp测试K点收敛性的脚本问题
211/1返回列表
查看: 8228  |  回复: 20
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

mazuju028

木虫 (正式写手)

[交流] 【求助成功】vasp测试K点收敛性的脚本问题

在用VASP进行K点收敛测试时，我基本上就是照抄侯柱锋VASP指南上的那个脚本，不过添加了个变量，后面就出现了语法错误，瞎搞了一个多小时都没弄好，请指教。

------run_k脚本如下----------------------------------------------

#!/bin/sh
rm WAVECAR
for i in 3 5 7 9 11
j=`expr "$i" "+" 2`
do
cat > KPOINTS < Automatic generation
0
Monhkorst-Pack
$i $j $i
0.0 0.0 0.0
!
echo "k mesh = $i*$j*$i" ; mpirun -np 8 vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'`
KP=`grep "irreducible" OUTCAR | tail -1 | awk '{printf "%5i \n", $2}'`
echo $i $KP $E >> result
done
----------------------------------------------------------------------------------------
--------log---------------------
./run_k: line 4: syntax error near unexpected token `j=`expr "$i" "+" 2`'

-----------------------------------------------------------------------------------------

谢谢。

[ Last edited by mazuju028 on 2011-4-16 at 16:59 ]
回复此楼

» 收录本帖的淘帖专辑推荐

VASP VASP 第一性原理

» 猜你喜欢

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

» 抢金币啦！回帖就可以得到:

查看全部散金贴
1楼 2011-03-09 09:47:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

g_wcai

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
mpirun -np 8 vasp
请问这句代表什么意思？！
赞一下(1人) 回复此楼
20楼2015-12-02 10:02:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

江山2010

木虫 (小有名气)
mazuju028(金币+1): 不行啊。 2011-03-09 14:45:37
let "j=expr $i + 2" 试一下
赞一下 回复此楼
2楼2011-03-09 11:57:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)
引用回帖:
Originally posted by 江山2010 at 2011-03-09 11:57:27:
let "j=expr $i + 2" 试一下

呵呵，试了，还是不行。
赞一下 回复此楼
3楼2011-03-09 14:45:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

guitar2033

至尊木虫 (职业作家)
★
gzqdyouxia(金币+1): 谢谢提示 2011-03-11 19:02:04
mazuju028(金币+2): 谢谢提示。 2011-03-13 09:53:24
do所在的第5行与j所在的第4行换位置～～
赞一下(1人) 回复此楼
4楼2011-03-11 15:02:19
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

babaleo

木虫 (正式写手)
★
sunyang1988(金币+1): 谢谢交流 2011-03-12 11:00:03
mazuju028(金币+2): 谢谢提示。 2011-03-13 09:53:37
LS正解. 或用
.....
do
j=$[$i + 2]
.....
赞一下(1人) 回复此楼
5楼2011-03-12 09:20:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yshl7237

铜虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by mazuju028 at 2011-03-09 09:47:55:
在用VASP进行K点收敛测试时，我基本上就是照抄侯柱锋VASP指南上的那个脚本，不过添加了个变量，后面就出现了语法错误，瞎搞了一个多小时都没弄好，请指教。

------run_k脚本如下----------------------------- ...

请问你进行K点测试时。KPOINT 文件怎么写呢？
赞一下(1人) 回复此楼
6楼2011-04-15 15:26:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)
★
youzhizhe(金币+1): 谢谢提示。 2011-04-18 08:33:38
引用回帖:
Originally posted by yshl7237 at 2011-04-15 15:26:27:
请问你进行K点测试时。KPOINT 文件怎么写呢？

就是写个run_k脚本啊，我帖子中就写有啊
赞一下(1人) 回复此楼
7楼2011-04-15 22:21:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yshl7237

铜虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
引用回帖:
Originally posted by mazuju028 at 2011-04-15 22:21:35:
就是写个run_k脚本啊，我帖子中就写有啊

那输入文件中就写 INCAR，POSCAR，POTCAR，run-k 这四个就可以是吗？
要是测试ecut的话，就写poscar，potcar，kpoints，run-e。我这样理解对吗？
赞一下(1人) 回复此楼
8楼2011-04-16 10:54:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)
★
youzhizhe(金币+1): 谢谢提示。 2011-04-18 08:33:48
引用回帖:
Originally posted by yshl7237 at 2011-04-16 10:54:52:
那输入文件中就写 INCAR，POSCAR，POTCAR，run-k 这四个就可以是吗？
要是测试ecut的话，就写poscar，potcar，kpoints，run-e。我这样理解对吗？

对啊
回复此楼
9楼2011-04-16 16:58:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

nwx_ky

木虫 (小有名气)
★
小木虫(金币+0.5):给个红包，谢谢回帖交流
qsub: script file 'run' cannot be loaded - No such file or directory
我按上面说的做了。为何会出现这一说？验证encut 时，请问不写INCAR文件吗？
赞一下(1人) 回复此楼

» 本帖已获得的红花（最新10朵）

zhouhongpeng
10楼2011-04-17 10:42:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

mazuju028

木虫 (正式写手)
★
franch(金币+1): 谢谢回帖交流 2011-04-22 14:25:29
引用回帖:
Originally posted by nwx_ky at 2011-04-17 10:42:25:
qsub: script file 'run' cannot be loaded - No such file or directory
我按上面说的做了。为何会出现这一说？验证encut 时，请问不写INCAR文件吗？

你可能没有设置权限吧，要吧run_k设置成可执行。
chmod +x run_k
赞一下(1人) 回复此楼
11楼2011-04-17 13:09:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)
★ ★
小木虫(金币+0.5):给个红包，谢谢回帖
mazuju028(金币+1): 谢谢交流 2011-05-15 11:36:05
j=`expr "$i" "+" 2`

这一行的位置放错了，要放在do所在行的下一行。
赞一下(2人) 回复此楼
12楼2011-05-15 11:32:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tider

木虫 (正式写手)
引用回帖:
Originally posted by yshl7237 at 2011-04-15 15:26:27:
请问你进行K点测试时。KPOINT 文件怎么写呢？

已经包含了
回复此楼
13楼2011-05-17 07:57:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)
请问你的脚本如何使用的？是把脚本和另外三个文件放到一个文件夹里吗？还是别的？求指导
赞一下 回复此楼
14楼2013-10-21 16:04:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

sungjen

银虫 (著名写手)
放到同一个文件夹里后再怎么操作呢？遇到这个问题很久了，一直没弄会，求指导啊
赞一下 回复此楼
15楼2013-10-21 16:06:08
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

计算化学新人

铜虫 (小有名气)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
12楼: Originally posted by valenhou001 at 2011-05-15 11:32:50
j=`expr "$i" "+" 2`

这一行的位置放错了，要放在do所在行的下一行。

请问，如果我想写一个测试k点的脚本，三个参数都可能不一样，请问应该怎么写呢？能不能给一个实例范本啊？谢谢了！
赞一下 回复此楼
16楼2015-03-06 09:00:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
送红花一朵
本帖仅楼主可见
赞一下
17楼2015-12-01 15:45:18
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

mazuju028

木虫 (正式写手)
引用回帖:
17楼: Originally posted by zhouhongpeng at 2015-12-01 15:45:18
请问 你解决这个问题了吗 能不能告诉我一下  我也 遇到了相同的问题...

好好看看回复，上面都有解答
赞一下 回复此楼
18楼2015-12-01 22:06:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

匿名

用户注销 (正式写手)
本帖仅楼主可见
19楼2015-12-02 09:35:49
已阅   申请1ST强帖   回复此楼   编辑   查看我的主页

heechao

至尊木虫 (著名写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
应该是 let j=$i+2
可以在命令窗口里面输入 i=1
                                  let j=$i+2
                                  echo $j
结果会显示 3的
赞一下(1人) 回复此楼
21楼2015-12-20 20:06:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 mazuju028 的主题更新
211/1返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 085600 英一数二272求调剂 5+6 vida_a 2026-03-01 16/800 2026-03-02 19:13 by zhukairuo
[考研] 268求调剂 +7 简单点0 2026-03-02 8/400 2026-03-02 19:07 by zhukairuo
[考研] 材料284求调剂，一志愿郑州大学英一数二专硕 +14 想上岸的土拨鼠 2026-02-28 14/700 2026-03-02 19:06 by zhukairuo
[考研] 299求调剂 +4 kkcoco25 2026-03-02 4/200 2026-03-02 18:04 by barlinike
[考研] 0703 总分319求调剂 +3 Xinyuu 2026-03-02 3/150 2026-03-02 17:20 by houyaoxu
[考研] 085600求调剂 +3 LRZZZZZZ 2026-03-02 3/150 2026-03-02 16:58 by zeng1010
[考研] 材料化工调剂 +12 今夏不夏 2026-03-01 14/700 2026-03-02 16:09 by 今夏不夏
[考研] 材料学硕318求调剂 +9 February_Feb 2026-03-01 9/450 2026-03-02 13:31 by njzyff
[考研] 291 求调剂 +3 化工2026届毕业� 2026-03-02 3/150 2026-03-02 12:55 by houyaoxu
[基金申请] 成果系统访问量大，请15分钟后再尝试。由此给您造成的不便，敬请谅解。 +5 xhuama 2026-03-02 5/250 2026-03-02 12:34 by stidwellNK
[考研] 265分求调剂不调专业和学校有行学上就 +6 礼堂丁真258 2026-02-28 9/450 2026-03-02 12:04 by 52hz~~
[考研] 哈工大计算机刘劼团队招生 +4 hit_aiot 2026-03-01 6/300 2026-03-02 11:53 by 一声问好
[考研] 264求调剂 +4 巴拉巴拉根556 2026-02-28 4/200 2026-03-02 10:48 by yuchj
[考研] 调剂 +3 13853210211 2026-03-02 4/200 2026-03-02 10:16 by 13853210211
[考研] 材料学调剂 +10 提神豆沙包 2026-02-28 12/600 2026-03-02 09:26 by 李老师！
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考研] 0856材料求调剂 +11 hyf hyf hyf 2026-02-28 12/600 2026-03-01 18:57 by 18137688336
[基金申请] 刚录用，没有期刊号，但是在线可看的论文可以放为代表作吗 10+3 arang1 2026-03-01 3/150 2026-03-01 16:43 by babero
[考研] 302材料工程求调剂 +4 Doleres 2026-03-01 5/250 2026-03-01 11:52 by liqiongjy
[论文投稿] 求助coordination chemistry reviews 的写作模板 10+3 ljplijiapeng 2026-02-27 4/200 2026-03-01 09:07 by babero
信息提示
关闭
请填处理意见
关闭