24小时热门版块排行榜    

查看: 977  |  回复: 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
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

valenhou001

至尊木虫 (职业作家)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
木冰天化(漫天飘雪代发): 金币+20, 谢谢交流 2019-04-18 08:21:27
首先,第二个循环最后少了"done"。

其次,你写了两个循环,前一个循环,好像你的目的是想在每个node上都run一个vasp计算。
你应该把第一去掉。只用第二个,把vasp运行作业的命令放到第二个循环里适当的位置上。
2楼2019-04-18 08:20:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 木冰天化 的主题更新
不应助 确定回帖应助 (注意:应助才可能被奖励,但不允许灌水,必须填写15个字符以上)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[硕博家园] 要不要读博 +7 王乔木 2024-05-24 7/350 2024-05-26 09:59 by 梦燕园
[基金申请] 科研之友阅读量近一周增加了200多。 +11 hdzw9071 2024-05-24 12/600 2024-05-26 09:38 by wanghuawei
[硕博家园] 人生 +15 暮色恋伊人 2024-05-22 15/750 2024-05-26 08:23 by elainzai
[基金申请] 又有两个中科院三、四区计算机领域SCI/SSCI期刊爆出问题了! +6 瞬息宇宙 2024-05-19 8/400 2024-05-26 07:47 by zju2000
[硕博家园] 文科博在木虫上存在感好低呀 +4 hahamyid 2024-05-25 6/300 2024-05-26 07:26 by hahamyid
[教师之家] 被惯着的学生终究要吃大亏 +20 535743368 2024-05-24 21/1050 2024-05-26 00:05 by 张开翅膀的鱼
[论文投稿] 真是奇怪的编辑部? +5 jjdg 2024-05-23 5/250 2024-05-25 21:57 by cqu_zzh
[教师之家] 经常觉得挺累的 +13 zylfront 2024-05-22 19/950 2024-05-25 21:08 by hjc404
[硕博家园] 研0 +4 控制调剂yl 2024-05-25 9/450 2024-05-25 18:21 by 鱼翔浅底1
[论文投稿] 通讯作者一定要放在署名的最后一个吗? 5+3 轨迹永远 2024-05-24 8/400 2024-05-25 17:52 by 盲牛
[硕博家园] 博士复试,申请成绩复核,有机会翻盘吗? +20 长海二声笑 2024-05-21 27/1350 2024-05-25 16:21 by james-james
[考博] 24年还有机会申博吗,求博导收留 +3 鄙人康康 2024-05-23 5/250 2024-05-24 23:11 by 鄙人康康
[电化学] 2200mA锂离子电池设计 +3 小乐987 2024-05-23 3/150 2024-05-24 16:59 by 克拉妮nin
[基金申请] 基金评审 +5 阿呆不呆 2024-05-20 5/250 2024-05-23 20:57 by 3001160025
[论文投稿] Neurocomputing 外审结束 +5 mollyzhang_2003 2024-05-23 5/250 2024-05-23 12:03 by nono2009
[电化学] 如何证明电极上镀上金了? 10+4 刻印时光 2024-05-22 4/200 2024-05-23 11:31 by 牵着小猫跑步
[论文投稿] 因为邮箱故障,没有收到论文接受后进入生产的邮件 5+3 VCOISTNT 2024-05-22 5/250 2024-05-23 10:03 by bobvan
[硕博家园] 超过35岁的海外博士还可以进高校吗? +6 905452934 2024-05-20 7/350 2024-05-23 07:54 by luwangba
[论文投稿] word转成pdf之后公式里面的字体变了,正文字体没变。 +9 1255037206 2024-05-20 11/550 2024-05-23 05:54 by tjushede
[论文投稿] SENSORS AND ACTUATORS B-CHEMICAL 15+4 荣小撇 2024-05-20 4/200 2024-05-21 10:07 by ca0yan9
信息提示
请填处理意见