| 查看: 951 | 回复: 2 | ||||
| 【悬赏金币】回答本帖问题，作者木冰天化将赠送您 180 个金币 | ||||
木冰天化新虫 (正式写手)
|
[求助]
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 |
回复此楼» 收录本帖的淘贴专辑推荐
 科研 |
» 猜你喜欢
博士毕业找工作选择！
已经有13人回复
-
2024国家自然科学基金资助额度
已经有5人回复
-
Optics Express一直卡在Final Decision
已经有6人回复
-
审稿人一收一拒，编辑给了修
已经有18人回复
-
特聘教授，算高级职称吗？青年基金还没结题的情况下还能报面上吗？
已经有8人回复
-
校级优秀本科毕业论文，比例多少？看指导教师资历吗？
已经有6人回复
-
求推荐博导
已经有4人回复
-
想请问下IEEE Trans最后一轮，你们等了多久？
已经有7人回复
-
等盲审
已经有16人回复
-
25年春秋博士申请
已经有5人回复
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25617.7
- 散金: 602
- 红花: 166
- 帖子: 3765
- 在线: 873.7小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
2楼2019-04-18 08:20:51
木冰天化
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3259.9
- 散金: 180
- 帖子: 809
- 在线: 87.2小时
- 虫号: 4302290
- 注册: 2015-12-21
- 专业: 金属材料表面科学与工程
3楼2019-04-18 14:17:01
5