| 查看: 1444 | 回复: 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 |
回复此楼
木冰天化
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3259.9
- 散金: 180
- 帖子: 809
- 在线: 87.2小时
- 虫号: 4302290
- 注册: 2015-12-21
- 专业: 金属材料表面科学与工程
3楼2019-04-18 14:17:01
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
2楼2019-04-18 08:20:51
