24小时热门版块排行榜    

查看: 2334  |  回复: 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的回帖

obaica

禁言 (著名写手)

本帖内容被屏蔽

2楼2016-08-27 20:23:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)

引用回帖:
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的回帖

magicmonk

至尊木虫 (著名写手)

引用回帖:
5楼: Originally posted by 从心而觅 at 2016-08-29 15:17:33
恩恩,确实是这样打出来的。
想问您一下,这个脚本建立后,是直接提交呢,还是跟incar等文件放到一起提交?另外提交的话是qsub xx.sh吗?...

这个脚本建立后,是直接提交呢,还是跟incar等文件放到一起提交?
跟incar等文件放到一起提交。

另外提交的话是qsub xx.sh吗?
如果是pbs的作业提交系统,当然是这个指令
我们最深的恐惧不是我们能力不够,我们最深的恐惧是能力超越了极限。
6楼2016-08-29 16:57:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

从心而觅

新虫 (小有名气)

引用回帖:
6楼: Originally posted by magicmonk at 2016-08-29 16:57:55
这个脚本建立后,是直接提交呢,还是跟incar等文件放到一起提交?
跟incar等文件放到一起提交。

另外提交的话是qsub xx.sh吗?
如果是pbs的作业提交系统,当然是这个指令...

您好,安装您说的方法,提交出错了。错误提示:
mpiexec_c0106: cannot connect to local mpd (/tmp/mpd2.console
1. no mpd is running on this host
2. an mpd is running but was started without a "console" (-n option)

有一点我很疑惑,感觉mpirun -np 24对于我这边调用不了节点,于是我用我提交作业脚本里面的命令
mpirun -v -genv I_MPI_DEVICE rdma -np $NPROCS /opt/software/vasp.5.3/vasp 还是有点问题。
最后我干脆把我原来提交作业的脚本跟这个脚本嵌套在一起了,作业交上去没问题,正在算,希望能正常吧。
不知道您对我这个问题有没有好的建议呢。
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
7楼2016-08-29 19:07:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ZC鑫

金虫 (初入文坛)

【答案】应助回帖

★ ★
漫天飘雪: 金币+2, 谢谢交流 2016-08-31 21:38:20
#PBS -N Si
#PBS -j oe
#PBS -q batch
#PBS -l nodes=1:ppn=20
#PBS -l walltime=1:00:00

# use openmpi 1.6.5

core_num=$(cat $PBS_NODEFILE|wc -l)

cd $PBS_O_WORKDIR

cat $PBS_NODEFILE

rm WAVECAR
for i in 5.1 5.2 5.3 5.4 5.5 5.6 5.7 5.8
do
cat >POSCAR<<!
Si
$i
0.0 0.5 0.5
0.5 0.0 0.5
0.5 0.5 0.0
Si
2
Direct
0 0 0
0.25 0.25 0.25
!
echo "a=$i"
mpirun -machinefile $PBS_NODEFILE -np ${core_num} --mca btl self,sm,openib --bind-to-core  ~/bin/vasp533_neb180
E=`grep "TOTEN" OUTCAR |tail -1|awk '{printf "%12.6f \n",$5}'`;echo $i $E >>comment
done
cat comment
这是我的感觉你少了vasp位置那一行,去你原来的job.pbs里面找找
8楼2016-08-30 20:15:21
已阅   回复此楼   关注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的回帖

从心而觅

新虫 (小有名气)

引用回帖:
9楼: Originally posted by ZC鑫 at 2016-08-30 20:16:54
#!/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/ ...

我感觉确实是由于提交任务那一行有问题,因为可以生成POSCAR文件。
你这个/public/software/vasp/5.3.3/openmpi-intel/vasp533没有提交任务的命令也可以?
一切有为法,如梦幻泡影!如露亦如电,应作如是观!
10楼2016-08-31 08:05:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 从心而觅 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 材料工程专硕调剂 +5 204818@lcx 2026-03-17 6/300 2026-03-18 22:55 by 204818@lcx
[考研] 085700资源与环境308求调剂 +3 墨墨漠 2026-03-18 3/150 2026-03-18 22:35 by bingxueer79
[考研] 328求调剂,英语六级551,有科研经历 +3 生物工程调剂 2026-03-16 10/500 2026-03-18 20:41 by Wangjingyue
[考研] 一志愿985,本科211,0817化学工程与技术319求调剂 +7 Liwangman 2026-03-15 7/350 2026-03-18 20:08 by walc
[考研] 311求调剂 +11 冬十三 2026-03-15 12/600 2026-03-18 14:36 by 星空星月
[考研] 331求调剂(0703有机化学 +7 ZY-05 2026-03-13 8/400 2026-03-18 14:13 by 007_lilei
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
[考研] 0703化学求调剂 总分331 +3 ZY-05 2026-03-13 3/150 2026-03-18 10:58 by macy2011
[考研] 环境工程调剂 +8 大可digkids 2026-03-16 8/400 2026-03-18 09:36 by zhukairuo
[考研] 梁成伟老师课题组欢迎你的加入 +8 一鸭鸭哟 2026-03-14 10/500 2026-03-17 15:07 by 一鸭鸭哟
[考研] 333求调剂 +3 文思客 2026-03-16 7/350 2026-03-16 18:21 by 文思客
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 326求调剂 +3 mlpqaz03 2026-03-15 3/150 2026-03-16 07:33 by Iveryant
[考研] 080500,材料学硕302分求调剂学校 +4 初识可乐 2026-03-14 5/250 2026-03-14 21:08 by peike
[考研] 复试调剂 +4 z1z2z3879 2026-03-14 5/250 2026-03-14 16:30 by JourneyLucky
[考研] 330求调剂 +3 ?酱给调剂跪了 2026-03-13 3/150 2026-03-14 10:13 by JourneyLucky
[考研] 材料工程调剂 +9 咪咪空空 2026-03-12 9/450 2026-03-13 22:05 by 星空星月
[考研] 304求调剂 +7 7712b 2026-03-13 7/350 2026-03-13 21:42 by peike
[考研] 工科278分求调剂 +5 周慢热啊 2026-03-12 7/350 2026-03-13 15:49 by JourneyLucky
[考博] 福州大学杨黄浩课题组招收2026年专业学位博士研究生,2026.03.20截止 +3 Xiangyu_ou 2026-03-12 3/150 2026-03-13 09:36 by duanwu655
信息提示
请填处理意见