24小时热门版块排行榜    

查看: 2314  |  回复: 12
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

从心而觅

新虫 (小有名气)

[求助] vasp优化构型提交脚本出错 已有2人参与

不论是侯老师还是手册优化构型的时候都提到一个脚本,首先说明我采用的是网上一位前辈的方法。http://emuch.315955.com/html/201105/3234667.html
按照流程,我进行了一系列操作,但是无法运行脚本。具体如下:
1、操作(请帮忙看下是否有误):按照上述前辈的方法,我建了一个脚本rvasp.pbs.sh,通过ssh上传(只上传这一个文件),然后qsub rvasp.pbs.sh,结果出错。
2:脚本内容
# !/bin/sh
rm WAVECAR
for i in  1.00 1.01  1.02 0.99 0.98
do
cat >POSCAR <<!
SYSTEM=Si-C4-BZ
$i
  15.648963576330972       5.532744131682456       9.582993941352564
  -0.000000000000001      16.598232395047376       9.582993941352568
   0.000000000000000       0.000000000000000      19.165987882705135
H C Si
16 40 2
Direct
   0.1692056544865040   0.1692056544865030  -0.0405409025403431            
………………            
   !
echo "a= $i" ;mpirun -np 24  vasp
E=`tail -1 OSZICAR`; echo $i $E >>SUMMARY  
cp CONTCAR CONTCAR$i  
done                                 !第七十六行
3、错误信息(91、92行啥都没有,不知道什么意思)
rm: cannot remove `WAVECAR': No such file or directory
/var/spool/torque/mom_priv/jobs/5584.cluster.hpc.org.SC: line 91: warning: here-document at line 5 delimited by end-of-file (wanted `!')
/var/spool/torque/mom_priv/jobs/5584.cluster.hpc.org.SC: line 92: syntax error: unexpected end of file
4、疑问:脚本里这一行 E=`tail -1 OSZICAR`  分明没有OSZICAR,不知道贴主前辈为啥写这个。
请各位帮忙看下,感激不尽!
回复此楼
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZC鑫

金虫 (初入文坛)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★
从心而觅: 金币+8, ★★★很有帮助 2016-08-31 08:03:52
#!/bin/sh
for alat in $(seq 5.60 0.02 5.80)
do
cat>POSCAR<<eof
Ge-bulk
$alat
0.5 0.5 0.0
0.0 0.5 0.5
0.5 0.0 0.5
Ge
2
Direct
0.0 0.0 0.0
0.25 0.25 0.25
eof
/public/software/vasp/5.3.3/openmpi-intel/vasp533
Etot=$(grep TOTEN OUTCAR | tail -1 | awk '{print $5}')
echo $alat $Etot >> EAA
done
这是另外一个
9楼2016-08-30 20:16:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

从心而觅

新虫 (小有名气)

引用回帖:
2楼: Originally posted by obaica at 2016-08-27 20:23:36
ding!

3楼2016-08-28 13:17:57
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

magicmonk

至尊木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与,应助指数 +1
蛋蛋小童鞋: 金币+2, 谢谢参与,小木虫鼓励积极回答问题者 2016-08-29 16:54:05
从心而觅: 金币+8, ★★★很有帮助, 忘给你金币啦,抱歉。 2016-08-31 07:47:55
键盘左上角,通常是在ESC下面的波浪线那个键,这是正确的反引号
我们最深的恐惧不是我们能力不够,我们最深的恐惧是能力超越了极限。
4楼2016-08-29 14:46:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)

引用回帖:
4楼: Originally posted by magicmonk at 2016-08-29 14:46:41
键盘左上角,通常是在ESC下面的波浪线那个键,这是正确的反引号

恩恩,确实是这样打出来的。
想问您一下,这个脚本建立后,是直接提交呢,还是跟incar等文件放到一起提交?另外提交的话是qsub xx.sh吗?
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
5楼2016-08-29 15:17:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +5 pnpwoqbg8f 2026-02-16 5/250 2026-02-18 18:25 by 7pdbc5bdo2
[博后之家] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 18:12 by 7pdbc5bdo2
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 18:10 by 7pdbc5bdo2
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 17:57 by 7pdbc5bdo2
[基金申请] 情人节自我反思:在爱情中有过遗憾吗? +5 瞬息宇宙 2026-02-15 6/300 2026-02-18 12:51 by 月下雪林
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 11:09 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 i3cz6qj6l2 2026-02-17 3/150 2026-02-18 10:39 by lqtl9djx19
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:53 by lqtl9djx19
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 08:38 by lqtl9djx19
[基金申请] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:40 by lotyj5cz79
[考研] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:38 by lotyj5cz79
[硕博家园] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +4 pnpwoqbg8f 2026-02-16 4/200 2026-02-18 07:23 by lotyj5cz79
[论文投稿] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-18 00:40 by tk2gfblvuz
[找工作] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 4/200 2026-02-18 00:23 by tk2gfblvuz
[公派出国] 售SCI一区文章,我:8 O5 51O 54,科目齐全,可+急 +3 pnpwoqbg8f 2026-02-17 3/150 2026-02-17 23:40 by tk2gfblvuz
[基金申请] 基金正文30页指的是报告正文还是整个申请书 +3 successhe 2026-02-16 4/200 2026-02-17 20:56 by successhe
[基金申请] 今年春晚有几个节目很不错,点赞! +5 瞬息宇宙 2026-02-16 6/300 2026-02-17 12:49 by jymy19840415
[微米和纳米] 球磨粉体时遇到了大的问题,请指教! 10+3 6sbiam 2026-02-12 15/750 2026-02-16 15:03 by tgzxzqj
[基金申请] 过年走亲戚时感受到了所开私家车的鄙视链 +3 瞬息宇宙 2026-02-15 5/250 2026-02-16 14:23 by aspect3000
[硕博家园] 江汉大学解明教授课题组招博士研究生/博士后 +3 cleverlyy 2026-02-12 3/150 2026-02-12 21:02 by qsdf1
信息提示
请填处理意见