| 查看: 1660 | 回复: 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 |
回复此楼» 收录本帖的淘帖专辑推荐
 科研 |
» 猜你喜欢
什么时候开奖？
已经有6人回复
-
朋友圈看到的
已经有6人回复
-
不应该看fileCode
已经有10人回复
-
欢迎发来filecode的Mz6后的代码验证其规律
已经有113人回复
-
结合人工智能，周易传统文化，filecode打分制来了，3分以上希望很大。
已经有5人回复
-
重要消息，中午系统在维护
已经有11人回复
-
咱们一起用铁证分析2026国家社科基金中标与否
已经有27人回复
-
为什么余额宝的年化利率越来越低？主要原因有哪些？
已经有3人回复
-
2027广东省杰青
已经有9人回复
-
明天放榜？
已经有5人回复
valenhou001
至尊木虫 (职业作家)
- 1ST强帖: 13
- 应助: 241 (大学生)
- 金币: 25701.7
- 散金: 602
- 红花: 166
- 帖子: 3782
- 在线: 873.8小时
- 虫号: 1007127
- 注册: 2010-04-27
- 专业: 凝聚态物性 II ：电子结构
2楼2019-04-18 08:20:51
木冰天化
新虫 (正式写手)
- 应助: 0 (幼儿园)
- 金币: 3271.5
- 散金: 180
- 帖子: 814
- 在线: 87.5小时
- 虫号: 4302290
- 注册: 2015-12-21
- 专业: 金属材料表面科学与工程
3楼2019-04-18 14:17:01
