24小时热门版块排行榜    

北京石油化工学院2026年研究生招生接收调剂公告
查看: 1088  |  回复: 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的回帖

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的回帖

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的回帖

magicmonk

至尊木虫 (著名写手)

马克

[ 发自手机版 http://muchong.com/3g ]
我们最深的恐惧不是我们能力不够,我们最深的恐惧是能力超越了极限。
6楼2014-04-09 15:47:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 linda19918 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿北化085600材料专硕275|有文章专利|求调剂 +5 Micky11223 2026-03-25 5/250 2026-03-28 15:02 by 果果妈咪
[考研] 一志愿厦门大学化学学硕307求调剂 +10 y7czhao 2026-03-26 10/500 2026-03-28 14:23 by 唐沐儿
[考研] 346求调剂 一志愿070303有机化学 +3 萝卜炖青菜 2026-03-28 3/150 2026-03-28 14:11 by 唐沐儿
[考研] 0703化学/290求调剂/本科经历丰富/工科也可 +9 丹青奶盖 2026-03-26 10/500 2026-03-28 07:45 by barnett0632
[考研] 0703化学求调剂,各位老师看看我!!! +5 祁祺祺 2026-03-25 5/250 2026-03-27 21:44 by 东方猪猪
[考研] 一志愿 西北大学 总分282 英语一62 求调剂 +7 18419759900 2026-03-25 8/400 2026-03-27 16:38 by 18419759900
[考研] 材料求调剂 +5 .m.. 2026-03-25 5/250 2026-03-27 11:08 by 不吃魚的貓
[考研] 一志愿吉大071010,316分求调剂 +3 xgbiknn 2026-03-27 3/150 2026-03-27 10:36 by guoweigw
[考研] 276求调剂。有半年电池和半年高分子实习经历 +10 材料学257求调剂 2026-03-23 11/550 2026-03-27 10:13 by YCIT- LHL
[考研] 284求调剂 +11 junqihahaha 2026-03-26 12/600 2026-03-27 04:37 by wxiongid
[考研] 329求调剂 +7 钮恩雪 2026-03-25 7/350 2026-03-27 04:28 by wxiongid
[考研] 机械学硕310分,数一英一,一志愿211本科双非找调剂信息 +3 @357 2026-03-25 3/150 2026-03-26 16:34 by by.MENG
[考研] 297求调剂 +6 田洪有 2026-03-26 6/300 2026-03-26 15:55 by 不吃魚的貓
[考研] 263求调剂 +6 yqdszhdap- 2026-03-22 10/500 2026-03-26 13:11 by 公瑾逍遥
[考研] 打过很多竞赛,085406控制工程300分,求调剂 +3 askeladz 2026-03-26 3/150 2026-03-26 09:08 by 给你你注意休息
[考研] 一志愿武理085500机械专业总分300求调剂 +3 an10101 2026-03-24 7/350 2026-03-25 00:00 by 山鬼0-
[考研] 300分,材料,求调剂,英一数二 +5 超赞的 2026-03-24 5/250 2026-03-24 21:07 by 星空星月
[考研] 085404电子信息284分求调剂 +4 13659058978 2026-03-24 4/200 2026-03-24 12:15 by syl20081243
[考研] 336化工调剂 +4 王大坦1 2026-03-23 5/250 2026-03-23 18:32 by allen-yin
[考研] 275求调剂 +6 shansx 2026-03-22 8/400 2026-03-22 15:27 by barlinike
信息提示
请填处理意见