24小时热门版块排行榜    

Znn3bq.jpeg
北京石油化工学院2026年研究生招生接收调剂公告
查看: 794  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料307分求大佬组收留 +11 Hll胡 2026-04-07 11/550 2026-04-08 11:02 by l_paradox
[考研] 0703化学调剂325分 +13 15771691647 2026-04-04 14/700 2026-04-08 10:28 by screening
[考研] 368化学求调剂 +4 wwwwabcde 2026-04-07 5/250 2026-04-07 22:51 by wwwwabcde
[考研] 286求调剂 +20 Faune 2026-04-06 20/1000 2026-04-07 11:33 by 诗与自由
[考研] 材料调剂 +13 一样YWY 2026-04-05 14/700 2026-04-07 09:51 by piklet
[考研] 求调剂 +4 电气小神童 2026-04-04 6/300 2026-04-07 00:14 by guanxin1001
[考研] 285求调剂 +5 mapmath 2026-04-06 6/300 2026-04-06 17:18 by 蓝云思雨
[考研] 材料334求调剂 +19 Eecho# 2026-04-03 19/950 2026-04-06 08:37 by 小小树2024
[考研] 一志愿哈尔滨工业大学085600英一数二337分求调剂 +12 lyz0427 2026-04-03 12/600 2026-04-06 06:37 by houyaoxu
[考研] 262求调剂 +7 天下第一文 2026-04-04 8/400 2026-04-05 21:31 by 激流勇渡
[考研] 22408 总分320,一篇论文二作,两个国三,求调剂 +3 Leomulufu 2026-04-04 5/250 2026-04-05 19:04 by chongya
[考研] 材料专硕(0856) 339分求调剂 +10 哈哈哈鹅哈哈哈 2026-04-04 10/500 2026-04-05 18:51 by 蓝云思雨
[考研] 085600调剂 +9 东照照照 2026-04-04 9/450 2026-04-05 13:44 by ujn_zhuj
[考研] 材料调剂 +9 革微桂 2026-04-04 9/450 2026-04-05 08:27 by 544594351
[考研] 294求调剂 +6 Grey_Ey 2026-04-02 9/450 2026-04-04 22:07 by hemengdong
[考研] 085701求调剂 +7 龚禹铭 2026-04-04 8/400 2026-04-04 13:49 by 小小树2024
[考研] 081200-11408-276学硕求调剂 +6 崔wj 2026-04-02 6/300 2026-04-03 10:19 by 蓝云思雨
[考研] 085801 总分275 本科新能源 求调剂 +19 bradoner 2026-04-01 23/1150 2026-04-03 10:07 by linyelide
[考研] 材料专硕322分 +11 哈哈哈吼吼吼哈 2026-04-01 11/550 2026-04-02 10:52 by lnilvy
[考研] 0710生物学,325求调剂 +3 mkkkkkl 2026-04-01 3/150 2026-04-02 09:48 by Jaylen.
信息提示
请填处理意见