版块导航: 正在加载中...

登录注册

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

24小时热门版块排行榜

>论坛更新日志 (6964)
>考研 (2002)
>导师招生 (1633)
>文献求助 (519)
>考博 (316)
>虫友互识 (307)
>硕博家园 (168)
>休闲灌水 (154)
>博后之家 (129)
>基金申请 (118)
>招聘信息布告栏 (114)
>论文投稿 (99)
>第一性原理 (79)
>教师之家 (74)
>公派出国 (50)
>外文书籍求助 (32)

返回列表

liuchenhan

新虫 (小有名气)

应助: 0 (幼儿园)
金币: 1536.5
散金: 4
帖子: 170
在线: 108小时
虫号: 1508714
注册: 2011-11-25
专业: 传热传质学

[求助] linux集群使用问题已有1人参与

本人需要指定输出能量到comment文件，所用计算资源为北京超级计算中心 intel处理器。
#!/bin/sh
rm WAVECAR
for i in 150 200 250 300 350 400
do
cat > INCAR <<!
SYSTEM = Si-Diamond
ENCUT = $i
ISTART = 0
ICHARG = 2
ISMEAR = -5
PREC = Accurate
!
echo "ENCUT = $i eV"
./vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }’` ### 第15行
echo $i $E >>comment ### 第16行
done

错误如下：
./do: command substitution: line 15: unexpected EOF while looking for matching `''
./do: command substitution: line 16: syntax error: unexpected end of file

还有输出的总能量是系统的总能量绝对值吗？就是有没有减去POTCAR中EATOM，E如果减去了，那么若系统中有2个原子，那就减去了2*EATOM？

回复此楼

» 猜你喜欢

一志愿武理085500机械专业总分300求调剂已经有6人回复
材料学求调剂已经有5人回复
考研调剂已经有4人回复
281求调剂已经有4人回复
0805 316求调剂已经有6人回复
085601求调剂总分293英一数二已经有3人回复
08工学调剂已经有17人回复
340求调剂已经有4人回复
311求调剂已经有3人回复
食品专硕一志愿双一流 328 已经有4人回复

» 本主题相关价值贴推荐，对您同样有帮助:

1楼 2014-12-13 10:50:57

已阅回复此楼关注TA 给TA发消息送TA红花 TA的回帖

未岸亡灵

铁杆木虫 (著名写手)

暂时离开小木虫

1ST强帖: 2
应助: 622 (博士)
金币: 6519.3
散金: 1221
红花: 111
帖子: 1345
在线: 517.7小时
虫号: 1730773
注册: 2012-04-01
性别: GG
专业: 理论和计算化学

【答案】应助回帖

★ ★ ★ ★ ★
liuchenhan: 金币+5, ★★★很有帮助 2014-12-17 08:45:02

第15行引号有问题，可能是反引号和引号，中文引号和英文的关系。你这里不是编辑器，不能显示颜色，我只能告诉你是$5后面两个引号的问题。修改如下：
#!/bin/sh
rm WAVECAR
for i in 150 200 250 300 350 400
do
cat > INCAR <<!
SYSTEM = Si-Diamond
ENCUT = $i
ISTART = 0
ICHARG = 2
ISMEAR = -5
PREC = Accurate
!
echo "ENCUT = $i eV"
./vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo $i $E >>comment
done

你粘贴过去试试。
另外，我建议你改脚本为grep “energy without entropy”而不是grep “TOTEN”

赞一下