24小时热门版块排行榜    

查看: 983  |  回复: 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的回帖

calos818

木虫 (著名写手)

新人菜鸟

大家说说啊
2楼2009-12-01 17:12:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

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的回帖

donggaomin

金虫 (正式写手)


小木虫(金币+0.5):给个红包,谢谢回帖
好像是您的脚本“'{printf "12.6f \n",$5)'`有点小问题,右边的括号变了,所以拿不出volumn
好好学习
6楼2011-04-18 11:28:24
已阅   回复此楼   关注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的回帖

enola

捐助贵宾 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
请问计算摸个压强下的焓,设置PSTRESS后要做结构优化吗?
再做静态计算得到焓,谢谢
timeflies..
8楼2012-05-17 19:59:06
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 calos818 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 6F5UbRU2I5hL 2026-09-14 3/150 2026-09-18 00:28 by cITaGg3p5edr
[硕博家园] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-13 7/350 2026-09-17 23:28 by cuPZTDXS3VOd
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-13 8/400 2026-09-17 23:28 by cuPZTDXS3VOd
[硕博家园] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +7 s3fFTmArrBt6 2026-09-13 9/450 2026-09-17 23:05 by cuPZTDXS3VOd
[考研] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +8 LwdutQ8HoqWP 2026-09-13 8/400 2026-09-17 22:53 by cuPZTDXS3VOd
[硕博家园] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +6 3n8v2C8RimXI 2026-09-13 6/300 2026-09-17 22:29 by cuPZTDXS3VOd
[考博] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +5 40ms4Wlo4umh 2026-09-14 5/250 2026-09-17 14:05 by ObVzqQQh1rrE
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 23jxep3nCNZb 2026-09-14 5/250 2026-09-17 13:53 by ObVzqQQh1rrE
[考研] 售SCI一区T0P文章,我:8O.55.1.O.54,科目全,可伽急 +6 23jxep3nCNZb 2026-09-14 6/300 2026-09-17 13:29 by ObVzqQQh1rrE
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 vZfe6xYu34yj 2026-09-14 3/150 2026-09-17 05:16 by k06BKNblrGRH
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +6 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-17 03:04 by Tql5LQhh5rLK
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +7 3n8v2C8RimXI 2026-09-13 8/400 2026-09-17 01:38 by Tql5LQhh5rLK
[考研] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 23jxep3nCNZb 2026-09-14 4/200 2026-09-16 21:03 by MHf40yA0Uhyb
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +5 6F5UbRU2I5hL 2026-09-14 6/300 2026-09-16 17:27 by i3KBRD19ERkG
[考博] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-16 17:00 by i3KBRD19ERkG
[考博] 上海工程技术大学激光智能制造课题组|2027级博士研究生招生公告 +7 水士口 2026-09-16 8/400 2026-09-16 13:47 by kkklens
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-16 02:14 by DgNGHc3h5tPl
[考研] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +5 s3fFTmArrBt6 2026-09-14 5/250 2026-09-15 16:59 by CgyNCDVNhGVg
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +4 6F5UbRU2I5hL 2026-09-14 4/200 2026-09-15 08:02 by C79jjtjAKjEn
[教师之家] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 LwdutQ8HoqWP 2026-09-13 4/200 2026-09-14 16:13 by Aj1rhIDL5ixY
信息提示
请填处理意见