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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » 【求助成功】vasp测试K点收敛性的脚本问题
211/1返回列表
查看: 8932  |  回复: 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返回列表
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +14 爱看书的可乐 2026-08-26 16/800 2026-08-26 18:28 by skyjudy
[基金申请] 我不理解！ +9 Edward_pc 2026-08-26 14/700 2026-08-26 18:03 by hsding2026
[基金申请] 国自然面上复盘~欢迎讨论 +3 晴天加油 2026-08-26 4/200 2026-08-26 17:59 by Tide man
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +6 process2012 2026-08-23 9/450 2026-08-26 15:33 by process2012
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 只有每年这种时候来逛逛小木虫 +26 yaoyewhu2008 2026-08-20 28/1400 2026-08-26 13:49 by zzuzxg
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了，中了面上 +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 怎么看青基中了没有啊 +4 叶九微 2026-08-26 4/200 2026-08-26 10:52 by xiacongxu
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 22/1100 2026-08-25 15:36 by 医学老男孩
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 人气不行了 +11 fansofjerry 2026-08-21 11/550 2026-08-25 11:04 by 孤独的英雄6
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +3 苏东坡二世 2026-08-23 3/150 2026-08-25 10:35 by shisan1313
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭