24小时热门版块排行榜    

查看: 1575  |  回复: 4
本帖产生 2 个 1ST强帖 ,点击这里进行查看
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

dingfengbo

新虫 (小有名气)

[求助] encut测试是否需要修改要提交的作业脚本

我想做encut测试这是我的encut脚本:
#!/bin/sh
rm WAVECAR
for a in 400 450 500 550 600 650 700 750
do
cat >INCAR < SYSTEM=CGa
LPLANE=.TRUE.
NPAR=8
LREAL = .TRUE.
PREC    = Accurate
EDIFF   =1E-6
EDIFFG  = -0.03      
IALGO   = 48         
NELMIN  = 4         
ISYM    = 0         
ISPIN = 1
NBANDS = 240
LCHARG  = .TRUE.
LWAVE   = .TRUE.     
LVTOT   = .FALSE.
NBLOCK  = 1         
NSW     = 10
IBRION  = 2
POTIM=0.3
ISMEAR = 0   
SIGMA = 0.05
ENCUT=$a
ISTART=0;ICHARG=2
!
echo "ENCUT=$a eV"
vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo $a $E >>comment
done
测试的时候我想将之提交,下面是我的作业脚本:
#!/bin/sh
#PBS -l nodes=2:ppn=8
#PBS -u XX
#PBS -N YY

JOBINFOR=$PBS_JOBID
MASTERNODE=`hostname`
SCRATCHDIR=$PBS_JOBID
NCPU=`wc -l < $PBS_NODEFILE`

SERVER=$PBS_O_HOST
WORKDIR=$PBS_O_WORKDIR
MKDIR=/bin/mkdir
RSH=/usr/bin/rsh
CP=/bin/cp
cat $PBS_NODEFILE > $WORKDIR/mpihosts
LAUNCH="mpirun --hostfile  $WORKDIR/mpihosts -np $NCPU"
PROGRAMEXEC="/opt/vasp/vasp"

${RSH} ${MASTERNODE}
cd ${WORKDIR}
${MKDIR} ${SCRATCHDIR}
cd ${SCRATCHDIR}
${CP} ${WORKDIR}/* .
${LAUNCH}  ${PROGRAMEXEC} > out 2>&1
请问我怎样修改作业脚本才能进行encut测试的计算呢?谢谢
回复此楼

» 收录本帖的淘帖专辑推荐

Photochemistry

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

haowenping

木虫 (著名写手)

引用回帖:
Originally posted by dingfengbo at 2011-04-22 16:37:10:
我想做encut测试这是我的encut脚本:
#!/bin/sh
rm WAVECAR
for a in 400 450 500 550 600 650 700 750
do
cat >INCAR <<!
SYSTEM=CGa
LPLANE=.TRUE.
NPAR=8
LREAL = .TRUE.
PREC    = Accura ...

你好,请问你有没有实现在pbs下优化你的那个截断能脚本?我现在也被这个问题困惑,如果你会了的话麻烦指教下我,谢谢了!
人有多大胆,地有多大产!
3楼2011-05-07 19:07:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

bingmou

金虫 (著名写手)

【答案】应助回帖


youzhizhe(金币+1, 1ST强帖+1): 谢谢指导。 2011-04-23 08:13:11
如果知识测试小体系的话,建议不要提交,直接单机运行。如果是大体系的话,估计一个ecut就可以了,如果全部是金属,400eV足够;如果有O,N可能要700eV。其他的都可以采用赝势中默认值。
如果一定要提交,可以考虑把这句
PROGRAMEXEC="/opt/vasp/vasp"
改为
PROGRAMEXEC='脚本位置”
2楼2011-04-22 22:58:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

dingfengbo

新虫 (小有名气)

引用回帖:
Originally posted by haowenping at 2011-05-07 19:07:03:
你好,请问你有没有实现在pbs下优化你的那个截断能脚本?我现在也被这个问题困惑,如果你会了的话麻烦指教下我,谢谢了!

你好!按照2楼说的就可以实现了,不知道你想问什么?
4楼2011-05-08 12:52:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiao72379

木虫 (正式写手)

宅男甲

【答案】应助回帖


bingmou(金币+1, 1ST强帖+1): 感谢交流! 2011-05-10 16:52:58
引用回帖:
Originally posted by bingmou at 2011-04-22 22:58:54:
如果知识测试小体系的话,建议不要提交,直接单机运行。如果是大体系的话,估计一个ecut就可以了,如果全部是金属,400eV足够;如果有O,N可能要700eV。其他的都可以采用赝势中默认值。
如果一定要提交,可以考虑 ...

个人意见:PROGRAMEXEC="/opt/vasp/vasp" 改为 PROGRAMEXEC='脚本位置”是不行的,最后一行其实就是一个运行并行的vasp程序的指令。如果改成脚本位置,那运行的程序也就是vasp都不能调用。只是楼主注意下变量PROGRAMEXEC="/opt/vasp/vasp"设置在你的机器上是正确的就行。还有就是我觉得这个脚本写得罗里罗嗦, 我整合个简单点的 看着清楚的 供楼主参考下:
#!/bin/bash
#PBS -N 5-zmos2nrbig
#PBS -j oe
#PBS -q vasp
#PBS -l walltime=9999:00:00
#PBS -l nodes=2:ppn=8

rm WAVECAR
for a in 400 450 500 550 600 650 700 750
do
cat >INCAR < SYSTEM=CGa
LPLANE=.TRUE.
NPAR=8
LREAL = .TRUE.
PREC    = Accurate
EDIFF   =1E-6
EDIFFG  = -0.03      
IALGO   = 48         
NELMIN  = 4         
ISYM    = 0         
ISPIN = 1
NBANDS = 240
LCHARG  = .TRUE.
LWAVE   = .TRUE.     
LVTOT   = .FALSE.
NBLOCK  = 1         
NSW     = 10
IBRION  = 2
POTIM=0.3
ISMEAR = 0   
SIGMA = 0.05
ENCUT=$a
ISTART=0;ICHARG=2
!
echo "ENCUT=$a eV"
mpirun XXX(此处是mpi的设置参数,根据自己的机器设置) ./vasp > ./out(此处将vasp程序 vasp必须的文件和脚本放同一目录,要嫌麻烦就在脚本里加几个变量把路径写进去)
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }'`
echo $a $E >>comment
done
天行健君子以自强不息
5楼2011-05-10 11:04:26
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] filecode,4个jtjc了 +10 ziyangfang 2026-08-19 11/550 2026-08-19 17:19 by wdc28
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +53 医学老男孩 2026-08-13 117/5850 2026-08-19 17:12 by 医学老男孩
[教师之家] 为什么余额宝的年化利率越来越低?主要原因有哪些? +4 瞬息宇宙 2026-08-15 4/200 2026-08-19 16:23 by allmydream
[基金申请] 今天放榜吗? +12 布布和一二 2026-08-19 13/650 2026-08-19 14:33 by 陈晨晨陈啊
[基金申请] 今天基金会出结果吗?20260819 +8 kkkl_v 2026-08-19 8/400 2026-08-19 12:34 by Vivilian
[基金申请] 朋友圈看到的 +6 wangzilk 2026-08-18 8/400 2026-08-19 10:55 by Haru815
[基金申请] 时间戳变了,能看出什么问题? +10 基诺咪客 2026-08-17 14/700 2026-08-19 10:34 by 基诺咪客
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +8 启萌科技 2026-08-12 30/1500 2026-08-19 10:22 by 启萌科技
[基金申请] 快农历七夕节了,轻松一下,男人悄悄话,女施主请不要进来。 +6 Tide man 2026-08-14 7/350 2026-08-19 09:56 by ZJTJZ
[基金申请] 什么时候开奖? +6 CrisMessi 2026-08-18 6/300 2026-08-19 08:06 by Equinoxhua
[论文投稿] 投稿咨询 +4 wwm09 2026-08-17 6/300 2026-08-18 15:36 by wwm09
[基金申请] 时间戳又变了8-15 +14 archvillain 2026-08-15 26/1300 2026-08-18 13:37 by phantomgost
[基金申请] 今天维护系统维护 祝所有人 高中 +8 gjjjzhong 2026-08-18 9/450 2026-08-18 13:01 by 家与远方
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +13 Ldrop2023 2026-08-13 16/800 2026-08-18 12:25 by 淀粉搬运工
[基金申请] 93BebMhtakh前后11位开头都是大写 +4 且听虎啸 2026-08-17 5/250 2026-08-18 00:49 by 蔡棒棒菂
[基金申请] 感觉是下周放榜了 +6 angus9576 2026-08-17 11/550 2026-08-17 23:57 by angus9576
[基金申请] 今天系统多次维护,明天很可能放榜! +8 zju2000 2026-08-16 9/450 2026-08-17 12:20 by lmz0216
[精细化工] 招聘 金属平磨液,抛光液研发工程师 +3 小天0311 2026-08-14 3/150 2026-08-16 07:31 by H9PLUS
[基金申请] 各位道友,我要去昆明玩几天,回来见。 +7 Tide man 2026-08-14 8/400 2026-08-15 01:11 by arzu_hma
[基金申请] 重要来源:本周末出结果 +10 瞬息宇宙 2026-08-12 10/500 2026-08-13 15:46 by likettle
信息提示
请填处理意见