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

24小时热门版块排行榜    

Znn3bq.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】vasp测试K点收敛性的脚本问题
汕头大学海洋科学接受调剂
211/1返回列表
查看: 8328  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +20 MAX怅惘 2026-04-09 22/1100 2026-04-14 14:57 by 独醉梦孤城
[考研] 085600材料与化工329分求调剂 +24 叶zilin 2026-04-13 25/1250 2026-04-14 09:20 by 试管破裂
[考研] 300分求调剂 （085501机械专硕，本科扬大） +9 xu@841019 2026-04-11 10/500 2026-04-14 08:48 by 木木mumu～
[考研] B区0809 ，数一英一，290 求调剂 +3 泠潍1111 2026-04-12 4/200 2026-04-13 20:35 by 学员JpLReM
[考研] 求调剂 +16 张番茄不炒蛋 2026-04-10 17/850 2026-04-12 13:58 by 熬夜成！
[考研] 一志愿郑州大学 22408 305分求调剂 +5 安小满zzz 2026-04-08 5/250 2026-04-12 00:41 by 蓝云思雨
[考研] 280求调剂 +7 兮兮夜夜 2026-04-09 10/500 2026-04-12 00:33 by 蓝云思雨
[考研] 一志愿厦大0856，306求调剂 +15 Bblinging 2026-04-11 15/750 2026-04-11 22:53 by 314126402
[考研] 药学专硕调剂 +8 ? 一路生?花? 2026-04-10 10/500 2026-04-11 21:21 by zhouxiaoyu
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 283求调剂 086004考英二数二 +17 那个噜子 2026-04-10 18/900 2026-04-11 16:27 by 明月此时有
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 087100初试311求调剂 +4 任雅琴 2026-04-09 4/200 2026-04-11 10:33 by zhq0425
[考研] 工科273调剂 +6 X1999 2026-04-09 7/350 2026-04-11 10:23 by zhq0425
[考研] 311求调剂 +13 xyp想读书 2026-04-10 14/700 2026-04-11 09:41 by 猪会飞
[考研] 337求调剂 +4 研s. 2026-04-10 4/200 2026-04-11 08:57 by zhq0425
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 一志愿中科大070300化学，314分求调剂 +12 wakeluofu 2026-04-09 12/600 2026-04-10 09:57 by liuhuiying09
[考研] 本科211 工科085400 280分求调剂 可跨专业 +3 LZH（等待调剂中 2026-04-09 3/150 2026-04-09 21:29 by wutongshun
[考研] 复试调剂，一志愿郑州大学材料与化工289分 +31 硕星赴 2026-04-08 31/1550 2026-04-09 16:54 by Delta2012
信息提示
关闭
请填处理意见
关闭