| 查看: 1760 | 回复: 2 | |||
木冰天化新虫 (正式写手)
|
[求助]
VASP计算E-V曲线脚本求助 已有1人参与
|
|
我用深圳超算，其脚本如下，我想在计算E-V曲线，在此脚本中加入for循环，请问该如何添加，添加部分如下所示，尝试过添加，但是最后OUTCAR中只有a=2.8时的结果 ，请问各位大神，此脚本该这么书写？感谢感谢 1，脚本 VASP=/home-yg/Soft/Vasp5.4/vasp_std source /home-yg/env/intel-12.1.sh source /home-yg/env/openmpi-1.6.5-intel.sh cd $CURDIR #start creating .nodelist rm -rf $CURDIR/nodelist >& /dev/null for i in `echo $LSB_HOSTS` do echo "$i" >> $CURDIR/nodelist done nodelist=$(cat $CURDIR/nodelist | uniq | awk '{print $1}' | tr '\n' ',') # # echo -n "start time " > time date >> time mpirun -np $NP -machinefile $CURDIR/nodelist $VASP &> output echo -n "end time " >> time ; date >> time cat time >> output # # cd .. done exit 2.循环：for a in 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 do rm WAVECAR cat > POSCAR <<! hcp-Co $a 1 0 0 -0.5 0.86602540378444 0 0 0 1.628 2 Direct 0.0 0.0 0.0 0.33333333 0.66666667 0.5 ! echo a=$a E=$(grep "energy without entropy" OUTCAR |tail -1| awk '{printf "%14.8f \n", $5}') V=$(grep "volume" OUTCAR |tail -1| awk '{printf "%14.8f \n", $5}') echo $E $V >> E-V.txt |
回复此楼» 收录本帖的淘帖专辑推荐
 科研 |
» 猜你喜欢
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
寻找合作团队
已经有3人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
同一个分子，一条来自文献，一条来自AI——不告诉你答案，你会选哪条？
已经有9人回复
-
售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急
已经有4人回复
-
售一区SCI文章T0P，我:8O.551.O54,科目全,可十急
已经有4人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
2楼2019-04-18 08:20:51
木冰天化
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3276.5
- 散金: 180
- 帖子: 815
- 在线: 87.6小时
- 虫号: 4302290
- 注册: 2015-12-21
- 专业: 金属材料表面科学与工程
3楼2019-04-18 14:17:01
