24小时热门版块排行榜

返回列表

当前主题已经存档。

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zongyan

铁杆木虫 (正式写手)

应助: 0 (幼儿园)
金币: 7064.2
散金: 430
红花: 1
帖子: 600
在线: 190.9小时
虫号: 287555
注册: 2006-10-21
专业: 凝聚态物性 II ：电子结构

[交流] [求助]利用Shell脚本进行VASP计算

我根据侯柱峰写的教程中的脚本进行截断能的测试，但是一运行以后就停在done处，各位是否经历过，请帮忙解决一下。以下是脚本：

#!/bin/bash
BIN= ~/vasp
rm WAVECAR
for i in 150 200 250 300 350 400
do
cat > INCAR < System = Si-Diamond
ISTART = 0 ; ICHARG = 2
ENCUT = $i
ISMEAR = -5
PERC = Accurate
!
echo "ENCUT = $i eV"; time vasp
E = 'grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'
echo $i $E >> comment
done

[ Last edited by wuchenwf on 2009-6-22 at 22:13 ]

回复此楼

» 猜你喜欢

博士自荐已经有6人回复
博士推荐已经有4人回复
求环氧树脂研发1名已经有10人回复
280求调剂已经有5人回复
什么是人一生最重要的？已经有10人回复
面上可以超过30页吧？已经有13人回复
为什么中国大学工科教授们水了那么多所谓的顶会顶刊，但还是做不出宇树机器人？已经有13人回复
版面费该交吗已经有17人回复
【博士招生】太原理工大学2026化工博士已经有8人回复

细推物理须行乐,何用浮名绊此生!

1楼 2008-12-31 19:28:08

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

qphll

金虫 (正式写手)

应助: 18 (小学生)
金币: 2708.7
散金: 2294
红花: 8
帖子: 554
在线: 227.1小时
虫号: 89654
注册: 2005-08-29
性别: GG
专业: 碳素材料与超硬材料

★ ★
fegg7502(金币+2):thank you very much! 1-23 16:00

我想所谓的问题应该是处在这句:
E = 'grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'
楼主应该在comment文件里面看不到输出的E值
正确的格式应该是:

E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'`
请特别注意最后两个是不同的,一个',另外一个是`