24小时热门版块排行榜    

查看: 975  |  回复: 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文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 3n8v2C8RimXI 2026-09-13 3/150 2026-09-15 05:49 by 5BDX0d0WFp7t
[教师之家] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 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一区文章,我:8.O.551.O.5.4,科目全,可伽急 +5 QUjhNVAcOSff 2026-09-13 6/300 2026-09-15 05:13 by 5BDX0d0WFp7t
[找工作] 售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一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +4 s3fFTmArrBt6 2026-09-14 5/250 2026-09-14 23:57 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文章T0P,我:8O.551.O54,科目全,可十急 +3 LwdutQ8HoqWP 2026-09-13 3/150 2026-09-14 22:24 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
[考博] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +5 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:36 by vZfe6xYu34yj
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O.5.4,科目全,可+急 +3 QUjhNVAcOSff 2026-09-13 5/250 2026-09-14 21:25 by vZfe6xYu34yj
[考博] 上海工程技术大学激光智能制造课题组招收博士研究生 +4 两三岁ss 2026-09-14 4/200 2026-09-14 20:29 by 班主任的同桌
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 s3fFTmArrBt6 2026-09-14 4/200 2026-09-14 17:48 by Aj1rhIDL5ixY
[教师之家] 售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 萨仁其其格
[考博] 西北工业大学材料学院电化学传感与催化课题组招收2027年推免研究生、博士生、博士后 +3 马方园 2026-09-10 4/200 2026-09-14 09:29 by Chem张zz
[教师之家] 浙江师范大学是怎么坑我的 +10 recruit123 2026-09-09 15/750 2026-09-13 10:45 by seaskyy
信息提示
请填处理意见