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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】vasp测试K点收敛性的脚本问题
211/1返回列表
查看: 9006  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +4 Aj1rhIDL5ixY 2026-09-14 4/200 2026-09-16 07:48 by i0IGWqQuchto
[硕博家园] 售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急 +4 40ms4Wlo4umh 2026-09-14 4/200 2026-09-16 06:27 by LH5NkK5Ud4bh
[考研] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-16 05:51 by LH5NkK5Ud4bh
[博后之家] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 05:35 by LH5NkK5Ud4bh
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 5/250 2026-09-16 05:03 by LH5NkK5Ud4bh
[博后之家] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-16 04:47 by LH5NkK5Ud4bh
[硕博家园] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 6/300 2026-09-16 03:50 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +4 s3fFTmArrBt6 2026-09-13 4/200 2026-09-16 03:38 by DgNGHc3h5tPl
[找工作] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +7 QUjhNVAcOSff 2026-09-13 8/400 2026-09-16 02:23 by DgNGHc3h5tPl
[考博] 申博发邮件 +4 Lmengk 2026-09-15 4/200 2026-09-16 00:05 by 木心想继续深造
[考研] 售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急 +3 Aj1rhIDL5ixY 2026-09-14 3/150 2026-09-15 20:14 by zNcEhTcH7ihq
[考博] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 17:39 by CgyNCDVNhGVg
[考研] 售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[博后之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:50 by CgyNCDVNhGVg
[考博] 售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-14 7/350 2026-09-15 16:23 by CgyNCDVNhGVg
[博后之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +6 QUjhNVAcOSff 2026-09-13 7/350 2026-09-15 14:13 by VMDqgnlKimtZ
[教师之家] 售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 07:50 by C79jjtjAKjEn
[考研] 售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-15 06:49 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[考研] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +5 LwdutQ8HoqWP 2026-09-13 5/250 2026-09-14 22:25 by vZfe6xYu34yj
信息提示
关闭
请填处理意见
关闭