24小时热门版块排行榜    

查看: 977  |  回复: 7
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

calos818

木虫 (著名写手)

新人菜鸟

[交流] 【求助】脚本出错,没法输出焓和体积,大家帮改改 已有3人参与

我的脚本如下(看侯老师教程上的)
#! /bin/bash
BIN=/bin/vasp
#rm WAVECAR
for i in  550 600; do
cat >INCAR < …………
PSTRESS=$i
…………
!
echo "P= $i" ; $BIN
V=`grep "volume" OUTCAR |tail -1| awk '{printf "12.6f \n",$5)'`
echo $i $V>>SUMMARY.fcc
done
cat SUMMARY.fcc
运行出来就显示P的值550和600,为体积的值没法输出,现在就想揪出OUTCAR里面最后的焓和体积的值,不知脚本哪儿错了,谢谢大家
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xh512

金虫 (正式写手)

★ ★
小木虫(金币+0.5):给个红包,谢谢回帖
youzhizhe(金币+1): 谢谢交流。 2011-04-18 13:04:14
V=`grep "volume" OUTCAR | tail -1| awk '{printf "%12.6f \n", $5}'`
7楼2011-04-18 11:34:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 8 个回答

gavinliu7390

木虫 (著名写手)

叶落鹰飞


zxzj05(金币+1,VIP+0):讨论奖励! 12-1 21:44
首先,你的bin下面需要有vasp。
而且你最好再done上一行加上
cp OUTCAR OUTCAR_$i
cp CONTCAR CONTCAR_$i
真理是一点点接近的!
3楼2009-12-01 19:50:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

calos818

木虫 (著名写手)

新人菜鸟

引用回帖:
Originally posted by gavinliu7390 at 2009-12-1 19:50:
首先,你的bin下面需要有vasp。
而且你最好再done上一行加上
cp OUTCAR OUTCAR_$i
cp CONTCAR CONTCAR_$i

bin下面有vasp模块,这个可以通过,就是没法正确找出体积和焓,似乎那个grep命令不起作用,你帮我看看
4楼2009-12-01 20:01:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

donggaomin

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖
引用回帖:
Originally posted by calos818 at 2009-12-01 12:58:24:
我的脚本如下(看侯老师教程上的)
#! /bin/bash
BIN=/bin/vasp
#rm WAVECAR
for i in  550 600; do
cat >INCAR <<!
…………
PSTRESS=$i
…………
!
echo "P= $i" ; $BIN
V=`grep ...

calos ,您好!能否传我一份“侯老师计算焓的教程”?hfdong@semi.ac.cn   
您的问题是如何解决的?谢谢您!
好好学习
5楼2011-04-18 11:24:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 06:01 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:50 by 5BDX0d0WFp7t
[公派出国] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 4/200 2026-09-15 05:50 by 5BDX0d0WFp7t
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 6/300 2026-09-15 05:38 by 5BDX0d0WFp7t
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:26 by 5BDX0d0WFp7t
[博后之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:01 by mibUvS8DDCwf
[公派出国] 售SCI一区T0P文章,我:8.O55.1.O.54,科目全,可十急 +6 QUjhNVAcOSff 2026-09-13 8/400 2026-09-15 05:01 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 23jxep3nCNZb 2026-09-14 3/150 2026-09-15 01:21 by mibUvS8DDCwf
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:25 by mibUvS8DDCwf
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:21 by mibUvS8DDCwf
[教师之家] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-15 00:12 by vZfe6xYu34yj
[教师之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 3/150 2026-09-14 23:21 by vZfe6xYu34yj
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-13 5/250 2026-09-14 22:45 by vZfe6xYu34yj
[公派出国] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 LwdutQ8HoqWP 2026-09-13 6/300 2026-09-14 22:37 by vZfe6xYu34yj
[找工作] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +5 3n8v2C8RimXI 2026-09-13 5/250 2026-09-14 22:13 by vZfe6xYu34yj
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 QUjhNVAcOSff 2026-09-13 4/200 2026-09-14 21:37 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
[找工作] 浙江师范大学是怎么坑我的 +6 recruit123 2026-09-09 6/300 2026-09-14 12:00 by 萨仁其其格
[考博] 27秋季材料博士申请求助 10+3 Cahal1 2026-09-11 4/200 2026-09-14 08:50 by 北京莱茵编辑
信息提示
请填处理意见