24小时热门版块排行榜    

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

gleerat

木虫 (正式写手)

[交流] 【求助成功】使用脚本做sigma优化,输出为空文件

我把周柱峰的sigma脚本复制下来进行sigma优化,结果comment中没有TS的值:
引用回帖:
0.02
0.04
0.06
0.08

但是直接使用命令,在终端中可以显示TS值:
引用回帖:
hpang@hpang-desktop:~/vasp/sigma$ grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'
   -0.034324

我的sigma脚本的内容为:
引用回帖:
#!/bin/sh
rm WAVECAR
for i in 0.02 0.04 0.06 0.08 0.10 0.12 0.14 0.16 0.18 0.20 0.22 0.24 0.26 0.28
do
cat > INCAR < System = CaFe2As2 afm1-opt
PREC = Accurate
ISPIN = 2
MAGMOM = 4*0 3 3 -3 -3 3 -3 -3 3 8*0
ISMEAR = 0
SIGMA = $i
ISTART = 0
EDIFF = 10E-5
ICHARG = 2
ENCUT = 310

LORBIT = 11
LREAL = Auto
!
echo "SIGMA = $i eV" ; time vasp
TS='grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }''
echo $i $TS >>comment
done

在执行脚本之后,到要输出TS值的那一步会出现%12.6f \n,: not found的错误提示。
我用一个脚本test测试:
引用回帖:
#!/bin/sh
TS='grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }''
echo $TS >>out

终端显示信息与sigma脚本一样的:
CODE:
hpang@hpang-desktop:~/vasp/sigma$ ./test
./test: 2: %12.6f \n,: not found

请教了。谢谢

[ Last edited by gleerat on 2010-8-13 at 18:35 ]
回复此楼

» 收录本帖的淘帖专辑推荐

科研

» 猜你喜欢

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

明察、慎思、笃行
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

gleerat

木虫 (正式写手)

引用回帖:
Originally posted by fanchen021 at 2010-08-13 17:01:17:
侯柱峰                        

罪过罪过,我给打错字了
明察、慎思、笃行
3楼2010-08-13 18:17:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 4 个回答

fanchen021

铁杆木虫 (正式写手)

gleerat(金币+2):谢谢提醒,犯下大错了 2010-08-13 18:23:31
gleerat(金币+8):我经常用prettyswan,所以系统不让评价,金币都给你吧 2010-08-13 18:36:09
侯柱峰                        
2楼2010-08-13 17:01:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Prettyswan

木虫 (小有名气)


zzy870720z(金币+1):谢谢回复,楼主和你认识吗?呵呵 2010-08-14 08:09:30
引用回帖:
#!/bin/sh
TS='grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }''
echo $TS >>out

改为
引用回帖:
#!/bin/sh
TS=`grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo $TS >>out

最外边的是反引号`,
4楼2010-08-13 18:21:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 297求调剂 +3 喜欢还是不甘心 2026-03-20 3/150 2026-03-21 18:33 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 279分求调剂 一志愿211 +14 chaojifeixia 2026-03-19 15/750 2026-03-21 13:24 by zhukairuo
[考研] 330求调剂0854 +3 assdll 2026-03-21 3/150 2026-03-21 13:01 by 搏击518
[考研] 306求调剂 +4 chuanzhu川烛 2026-03-18 4/200 2026-03-21 08:25 by laoshidan
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 299求调剂 +6 △小透明* 2026-03-17 6/300 2026-03-21 02:42 by JourneyLucky
[考研] 324分 085600材料化工求调剂 +4 llllkkkhh 2026-03-18 4/200 2026-03-21 01:24 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 311求调剂 +5 冬十三 2026-03-18 5/250 2026-03-21 00:16 by JourneyLucky
[考研] 一志愿南京理工大学085701资源与环境302分求调剂 +4 葵梓卫队 2026-03-18 6/300 2026-03-20 23:02 by JourneyLucky
[考研] 一志愿中南化学(0703)总分337求调剂 +8 niko- 2026-03-19 9/450 2026-03-20 21:57 by luoyongfeng
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 319求调剂 +3 小力气珂珂 2026-03-20 3/150 2026-03-20 19:47 by JourneyLucky
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 085601专硕,总分342求调剂,地区不限 +5 share_joy 2026-03-16 5/250 2026-03-18 14:48 by haxia
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
信息提示
请填处理意见