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

24小时热门版块排行榜    

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

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的回帖
查看全部 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 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 投票:  有多少人是今天查系统知道结果的？ +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 11:51 by xiangy672
[基金申请] 29号明天会评吗 +4 笨笨唐 2026-08-28 4/200 2026-08-31 09:30 by huixian257
[基金申请] 基金不中，共勉 +12 eulota 2026-08-26 12/600 2026-08-31 08:42 by ZJTJZ
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 麻烦专家们看看评委们的意见（F口面上） +6 gdd2018 2026-08-28 11/550 2026-08-30 08:58 by 超级无敌华子
[基金申请] 国自然评审意见 +13 wangmingqi 2026-08-28 19/950 2026-08-29 10:22 by Poppy1104
[基金申请] 系统查不到 +11 董八千 2026-08-26 11/550 2026-08-28 18:06 by Leogzhya
[基金申请] 基金系统什么内容也没有 30+4 winsaint 2026-08-27 9/450 2026-08-28 11:06 by maolC
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +5 yuleib84 2026-08-26 6/300 2026-08-28 00:02 by yudaoqian88
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 怎么看青基中了没有啊 +5 叶九微 2026-08-26 5/250 2026-08-27 10:35 by l_zh2008
[基金申请] 为什么 国际(地区)合作与交流项目 没有放榜？ 10+3 majunge000 2026-08-26 11/550 2026-08-27 08:42 by 北京莱茵编辑
[基金申请] 我不理解！ +15 Edward_pc 2026-08-26 23/1150 2026-08-26 20:34 by zzuzxg
[基金申请] 出来了 +9 trojank 2026-08-26 9/450 2026-08-26 14:25 by 宝贝虫子
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合可查了 +3 paperzjh 2026-08-26 3/150 2026-08-26 10:41 by LemmonTr
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 牛来！米来！面来！ +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
关闭
请填处理意见
关闭