| 查看: 735 | 回复: 2 | ||
chenxiaoqian新虫 (小有名气)
|
[求助]
计算热膨胀系数提交不上去
|
|
我在用qha的方法计算热膨胀系数，用的是学校的集群，用./script-qha.sh -r提交可是作业提交不上去，提交后显示取消。 script-qha.sh文件内容为： [chenq@su01 qha-test]$ cat script-qha.sh #!/bin/bash dim="1 1 1" pa="0 1/2 1/2 1/2 0 1/2 1/2 1/2 0" mp="3 3 3" nunit=32 while getopts ':mrc' opt; do case $opt in m) for i in `/bin/ls QHAPOSCAR*`;do qhaname=`echo $i|sed s/POSCAR//` mkdir $qhaname mv $i $qhaname/POSCAR cd $qhaname pwd phonopy -d --dim="$dim" for j in `/bin/ls POSCAR-*`;do dispname=`echo $j|sed s/POSCAR/disp/` mkdir $dispname mv $j $dispname/POSCAR cp ../{INCAR,KPOINTS,POTCAR} $dispname done mkdir perfect cp ../{INCAR,KPOINTS,POTCAR} perfect cp SPOSCAR perfect/POSCAR cd .. done ;; r) for i in `/bin/ls -d QHA-*`;do cd $i pwd vol=`echo $i|sed s/QHA-//` # For displacements for j in `/bin/ls -d disp-*`;do cd $j pwd num=`echo $j|sed s/disp-//` sed s/num/$num/ ../../calc.sh|sed s/vol/$vol/|qsub cd .. done # For perfect cd perfect sed s/num/per/ ../../calc.sh|sed s/vol/$vol/|qsub cd .. cd .. done ;; c) echo "# cell volume energy of cell other than phonon" > e-v.dat for i in `/bin/ls -d QHA-*`;do cd $i phonopy -f disp-*/vasprun.xml phonopy --mp="$mp" -t --dim="$dim" --pa="$pa" --tmax=2004 --tstep=2 volume=`grep volume perfect/vasprun.xml|tail -n 1|awk -F'<|>' '{printf("%20.13f", $3)}'` energy=`grep e_wo_entrp perfect/vasprun.xml|tail -n 1|awk -F'<|>' '{printf("%20.13f", $3)}'` cd .. echo `echo "$volume/$nunit"|bc -ls` `echo "$energy/$nunit"|bc -ls` >> e-v.dat done ;; esac done  1.png  2.png |
回复此楼» 猜你喜欢
279分求调剂 一志愿211
已经有12人回复
-
297求调剂
已经有9人回复
-
085700资源与环境308求调剂
已经有13人回复
-
材料 336 求调剂
已经有4人回复
-
华东师范大学-071000生物学-293分-求调剂
已经有3人回复
-
280求调剂
已经有8人回复
-
324分 085600材料化工求调剂
已经有4人回复
-
一志愿华中科技大学，080502，354分求调剂
已经有5人回复
-
一志愿武理材料305分求调剂
已经有7人回复
-
271材料工程求调剂
已经有8人回复
chenxiaoqian
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 659.4
- 帖子: 152
- 在线: 76.5小时
- 虫号: 2952102
- 注册: 2014-01-31
- 专业: 原子和分子物理
2楼2017-09-21 21:45:58
3楼2017-10-10 23:37:56
40