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

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】vasp测试K点收敛性的脚本问题
51/1返回列表
查看: 8027  |  回复: 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的回帖

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的回帖
查看全部 21 个回答

江山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的回帖
查看全部 21 个回答
普通表情 龙 兔 虎 猫 高级回复(可上传附件)
信息提示
关闭
请填处理意见
关闭