| 查看: 931 | 回复: 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 |
回复此楼» 收录本帖的淘贴专辑推荐
 科研 |
» 猜你喜欢
with editor 两个月了，什么原因？
已经有12人回复
-
博士白读了
已经有42人回复
-
Chemical Engineering Journal投稿3周了，一直显示With editor状态。这是送审了吗？
已经有11人回复
-
"颜宁:基础研究应顶天立地"能做到基础研究同时顶天立地的才是牛人
已经有4人回复
-
普通院校药学硕士，做合成的，感觉找不到工作
已经有15人回复
-
申博求助
已经有10人回复
-
Nature一直在编辑手里，考虑好几天了，是悬了吗
已经有12人回复
-
刚博士毕业，入职普通一本高校
已经有19人回复
-
催稿后秒拒
已经有12人回复
-
取博导收留
已经有8人回复
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