24小时热门版块排行榜    

查看: 790  |  回复: 1

HSE06

新虫 (小有名气)

[求助] phon处理脚本报错求助

vasp计算完gamma点振动频率之后,用phon提供的脚本提取力常数,报错如下:
./phon2vasp.sh: line 17: syntax error near unexpected token `"   1  0.00477204 -0.00477204  0.00477204 "'
./phon2vasp: line 17: ` "   1  0.00477204 -0.00477204  0.00477204 " \'

这是怎么回事啊?多谢!

phon提供的vasp接口文件内容如下
then
# here the call to VASP
~/bin/vasp
mv OUTCAR OUTCAR.$j
fi
#
# use awk to extract forces
#
echo $i >DYNMAT.$j
awk '
/total drift:/ { start=0}
!/----------/  { if ( start ==1) { printf "%14.8f  %14.8f  %14.8f\n",$4,$5,$6  }}
/TOTAL-FORCE/  { start = 1}
' >DYNMAT.$j

done

# restore POSCAR file
mv POSCAR.phon POSCAR

echo $j > FORCES
i=1
while test $i -le $j
do
cat DYNMAT.$i >> FORCES
i=`expr $i + 1`
done

…………
…………

@gavinliu7390
回复此楼

» 猜你喜欢

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

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

2009537029

新虫 (初入文坛)

这个脚本需要修改一下,具体要把前面部分修改,下面是我的修改,
if test -f POSCAR.phon
then
  echo "POSCAR.phon exist remove before starting batchfile"
  exit 1
fi
cp POSCAR POSCAR.phon
#
#  loop over all exited ions
#
j=0
# here put the content of the DISP file
declare -a arr='( \
"  1  0.00099602  0.00000000  0.00000000 "\
"  1  0.00000000  0.00000000  0.00099602 "\
# " 17  0.00099602  0.00000000  0.00000000 "\
# " 17  0.00000000  0.00099602  0.00000000 "\
# " 17  0.00000000  0.00000000  0.00099602 "\
# " 49  0.00099602  0.00000000  0.00000000 "\
# " 49  0.00000000  0.00000000  0.00099602 "\
# " 65  0.00099602  0.00000000  0.00000000 "\
# " 65  0.00000000  0.00000000  0.00099602 "'
n=0
x=0
while [ $x -le $n ];
do
i=${arr[$x]}
x=`expr $x + 1`
j=`expr $j + 1`
echo "run number $j"

#
# use awk to displace atoms (somewhat tricky)
#
awk '
/SUBSIT/  { npos = $2 ;
            x=$3; y=$4 ; z=$5; }
!/SUBSIT/ { line=line+1
  if (line-7 == npos) {
    printf "%12.9f %12.9f %12.9f\n", $1+x, $2+y, $3+z
  }
  else if (line > 7) {
    printf "%12.9f %12.9f %12.9f\n", $1, $2, $3
  }
  else
    print
} ' >POSCAR <<!
SUBSIT $i
2楼2015-05-15 10:56:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 HSE06 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿中国海洋大学,生物学,301分,求调剂 +4 1孙悟空 2026-03-17 4/200 2026-03-18 17:59 by fivewind
[考研] 085601材料工程专硕求调剂 +6 慕寒mio 2026-03-16 6/300 2026-03-18 14:26 by 007_lilei
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +6 Liwangman 2026-03-15 6/300 2026-03-18 13:21 by 尽舜尧1
[考研] 288求调剂,一志愿华南理工大学071005 +4 ioodiiij 2026-03-17 4/200 2026-03-18 12:36 by Linda Hu
[考研] 303求调剂 +4 睿08 2026-03-17 6/300 2026-03-18 11:01 by Iveryant
[考研] 278求调剂 +5 烟火先于春 2026-03-17 5/250 2026-03-18 08:43 by 星空星月
[考研] 293求调剂 +11 zjl的号 2026-03-16 16/800 2026-03-18 08:10 by zhukairuo
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 290求调剂 +6 孔志浩 2026-03-12 11/550 2026-03-17 14:41 by 周舟舟77
[考研] 211本,11408一志愿中科院277分,曾在中科院自动化所实习 +6 Losir 2026-03-12 7/350 2026-03-17 12:09 by danranxie
[考研] 材料工程专硕274一志愿211求调剂 +6 薛云鹏 2026-03-15 6/300 2026-03-17 11:05 by 学员h26Tkc
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 药学383 求调剂 +3 药学chy 2026-03-15 4/200 2026-03-16 20:51 by 元子^0^
[考研] 070300化学学硕求调剂 +6 太想进步了0608 2026-03-16 6/300 2026-03-16 16:13 by kykm678
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 求调剂 +7 18880831720 2026-03-11 7/350 2026-03-13 16:10 by JourneyLucky
[考研] 材料专硕350 求调剂 +4 王金科 2026-03-12 4/200 2026-03-13 16:02 by ruiyingmiao
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考研] 070303一志愿西北大学学硕310找调剂 +3 d如愿上岸 2026-03-13 3/150 2026-03-13 10:43 by houyaoxu
[考研] 333求调剂 +3 152697 2026-03-12 4/200 2026-03-13 07:08 by Iveryant
信息提示
请填处理意见