24小时热门版块排行榜    

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

linda19918

银虫 (小有名气)

[求助] 第一次用VASP,学了几天了,怎么用脚本来找合适的ENCUT,按照说明说总是出错 已有2人参与

#!/bin/sh
rm WAVECAR
for i in 150 200 250 300 350 400
do
cat > INCAR <<!
SYSTEM = Si-Diamond
ENCUT = $i
ISTART = 0 ; ICHARG = 2
ISMEAR = -5
PREC = Accurate
!
echo "ENCUT = $i eV" ; time vasp
E='grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }''
echo $i   $E >>comment
done

然后运行该脚本,该脚本的名字为new,输入的sh new过后出现一下提示:
: No such file or directory
rm: cannot remove `WAVECAR\r': No such file or directory
'ew: line 4: syntax error near unexpected token `do
'ew: line 4: `do
我该怎么去修改或者我的操作是错误的,
PS:INCAR,POSCAR,POTCAR,KPONTS都有
回复此楼

» 猜你喜欢

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

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

linda19918

银虫 (小有名气)

引用回帖:
3楼: Originally posted by 甲斐之虎 at 2014-03-28 09:14:43
注意,那里是`,而不是'
'grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }''
楼上给的那个是对的

我都改过来了,但是必须要用楼上那种脚本才能运行么?我那种不对么,还有我用的SSH连接到大型计算机上,不知道怎么运行脚本,今天又试着用了一下脚本还是不断出错
5楼2014-04-06 19:21:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 6 个回答

490105179

木虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流。。。 2014-04-07 21:59:49
我的脚本,希望对你有用。
#!/bin/sh
APP_NAME=intelg_mid
NP=12
NP_PER_NODE=12
RUN="RAW"
CURDIR=$PWD

VASP=~/bin/vasp.5.2.11-mp
source /home-gg/compiler/intel/composer_xe_2011_sp1/bin/compilervars.sh intel64
source /home-gg/compiler/mpi/openmpi-1.4.4-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' ',')

rm WAVECAR
for i in 400 450 500 550 600

do
cat > INCAR <<!                             
SYSTEM=SiC
ENCUT=$i
ISTART=0;ICHARG=2
ISMEAR=-5
PREC=high
!
echo "ENCUT= $i eV"
mpirun -np $NP -machinefile $PBS_NODEFILE vasp
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'`
KP=`grep "irreducible" OUTCAR | tail -1 | awk '{printf"%12.4f \n", $5}'`
echo $i $E >> comment
#
#
echo -n "start time  " > time
date >> time
mpirun -np $NP -machinefile $CURDIR/nodelist $VASP  &> vasp.out
echo -n "end   time  " >> time ; date >> time
cat time >> vasp.out
done
#
#
exit
2楼2014-03-27 16:51:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

甲斐之虎

银虫 (正式写手)

【答案】应助回帖

★ ★
感谢参与,应助指数 +1
franch: 金币+2, 谢谢回帖交流。。。 2014-04-07 21:59:56
注意,那里是`,而不是'
'grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5 }''
楼上给的那个是对的
别人笑我太疯颠,我笑他人看不穿。不见五陵豪杰墓,无花无酒锄作田。
3楼2014-03-28 09:14:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

linda19918

银虫 (小有名气)

引用回帖:
2楼: Originally posted by 490105179 at 2014-03-27 16:51:47
我的脚本,希望对你有用。
#!/bin/sh
APP_NAME=intelg_mid
NP=12
NP_PER_NODE=12
RUN="RAW"
CURDIR=$PWD

VASP=~/bin/vasp.5.2.11-mp
source /home-gg/compiler/intel/composer_xe_2011_sp1/bin ...

这个脚本要怎么运行呢??
4楼2014-04-06 19:19:28
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 哪位老哥知道今年的国自然具体哪一天放榜? +5 Ldrop2023 2026-08-13 5/250 2026-08-14 18:48 by ssxclkj
[基金申请] 国基金的申报应该改成非等额制,评价高的钱多评价低的钱少,但是增加资助率 +8 a089 2026-08-07 8/400 2026-08-14 18:46 by beefly
[基金申请] 欢迎发来filecode的Mz6后的代码验证其规律 +22 医学老男孩 2026-08-13 48/2400 2026-08-14 18:33 by wuguocheng
[基金申请] filecode +6 cratir 2026-08-14 10/500 2026-08-14 18:29 by 笑叹辞穷
[基金申请] 咱们一起用铁证分析2026国家社科基金中标与否 +6 启萌科技 2026-08-12 18/900 2026-08-14 18:18 by xueboshi
[基金申请] 大家玩filecode,现在的申请与受理打不开了 +3 dcqxinyang 2026-08-14 6/300 2026-08-14 18:03 by 690092147
[基金申请] 奇怪,两个人的filecode固定段从头到尾一模一样 +8 布布和一二 2026-08-10 11/550 2026-08-14 14:58 by Equinoxhua
[基金申请] 有时候,自然基金真的不能太认真 (我的申报经验) +6 majunge000 2026-08-11 8/400 2026-08-14 14:44 by 医学老男孩
[基金申请] 应该是93bebmhtak前后十一个字符比较关键 +23 Lanmanbaby 2026-08-09 37/1850 2026-08-14 13:40 by Equinoxhua
[基金申请] 关于Filecode分析方法 +9 majunge000 2026-08-10 12/600 2026-08-13 23:20 by iwuli
[文学芳草园] 阿姨 +4 汪汪锅 2026-08-09 4/200 2026-08-13 19:43 by arzu_hma
[硕博家园] 一作与独作在应聘高校教师时区别大吗 +3 mbygzh 2026-08-08 4/200 2026-08-13 19:31 by 龙-樱
[基金申请] 静等基金结果 +8 gjjjzhong 2026-08-10 21/1050 2026-08-13 17:56 by 且听虎啸
[基金申请] 不应该看fileCode +7 且听虎啸 2026-08-12 9/450 2026-08-13 14:27 by flydreamws
[基金申请] Filecode 又变了,巨变 +3 WH3796 2026-08-12 4/200 2026-08-13 14:13 by 小木虫6752397
[基金申请] 为什么网上很多人说本周 12号出结果 +6 瞬息宇宙 2026-08-10 7/350 2026-08-11 19:25 by Tide man
[基金申请] 什么时候出结果,有咨询渠道??? +3 Tide man 2026-08-11 3/150 2026-08-11 17:54 by kudofaye
[基金申请] 据悉今年马上要出结果了 +7 瞬息宇宙 2026-08-10 8/400 2026-08-10 12:42 by Vivilian
[基金申请] 面上项目filecode邪修 +5 西山十月 2026-08-09 7/350 2026-08-10 07:32 by 仁砚薪传
[基金申请] fileCode有新解读? +10 Tide man 2026-08-08 18/900 2026-08-09 12:55 by 仁砚薪传
信息提示
请填处理意见