24小时热门版块排行榜    

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

butterflyyh

木虫 (著名写手)

[求助] ENCUT测试脚本

使用VASP的ENCUT测试脚本optencut.sh:
rm WAVECAR
for i in 100 150 200 250 300 350 400 450 500 550 600
do
cat > INCAR < SYSTEM = bcc Fe
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
执行后出错:syntax error: unexpected end of file
optencut.sh脚本是我按照VASP计算前的验证写的,请问有没有人知道出错原因?多谢
回复此楼
nothingisimpossible
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

youzhizhe

荣誉版主 (职业作家)

有志者

优秀版主优秀版主

引用回帖:
1楼: Originally posted by butterflyyh at 2011-08-13 22:05:50:
使用VASP的ENCUT测试脚本optencut.sh:
rm WAVECAR
for i in 100 150 200 250 300 350 400 450 500 550 600
do
cat > INCAR <<!
SYSTEM = bcc Fe
ENCUT = $i
...

问题出在最后两句上 修改结果如下
E=`grep "TOTEN" OUTCAR | tail -1 | awk '{printf "%12.6f\n",$5}'`
echo $i $E >> comment

应该就没有问题了。
微博http://weibo.com/10986069
11楼2011-08-19 16:36:47
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 18 个回答

leehb870605

银虫 (小有名气)

【答案】应助回帖

★ ★ ★
butterflyyh(金币+15): 2011-08-17 18:01:40
youzhizhe(金币+3): 谢谢交流。 2011-08-17 19:56:40
1,脚本第一行的
#!/bin/sh
这句话不知道你写了没有,

2,如果你是在WINDOWS环境下的写的脚本,然后再在LINUX环境下运行,可能就会出如下问题,下面是我转的

Shell 尴尬一: syntax error:unexpected end of file
2010-05-23 21:20
在windows下写好了Shell 但是在linux下用
>> sh -n [filesName]  检查语法总是出一个错误 syntax error:unexpected end of file
原因如下:
dos文件传输到unix系统时,会在每行的结尾多一个^M,在vi的时候,当你用如下命令:
>>vi dos.txt
>>:set fileformat=unix
>>:w
就会看到这些存在于每行结尾的^M符号,这个就是产生syntax error:unexpected end of file的原因
解决方案:
在vi下把这些^M都删除后, 搞定!

3,E=$(grep "TOTEN" OUTCAR | tail -1| awk '{printf "%12.6f \n", $5}') 这一句在我那里曾经出过问题,我把它更改过,再每次循环时ECHO一个总能,只是没有这句话后面那些控制输出格式的东西罢了,也可以尝试把这句删掉再试试(个人意见)
2楼2011-08-16 14:57:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

butterflyyh

木虫 (著名写手)

★ ★
youzhizhe: 建议你采用“引用回复”。 2011-08-17 19:57:07
mazuju028(金币+2): 谢谢交流 2011-08-19 08:26:07
1. 第一行的#!/bin/sh加了;

2.是在LINUX下直接编写的,没有每行结尾的^M符号现象;

3.不知道你有没遇到VASP没运算完就读取i 循环的现象?
nothingisimpossible
3楼2011-08-17 18:04:49
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w76990480

荣誉版主 (文坛精英)

super-emuch

优秀版主文献杰出贡献文献杰出贡献

★
youzhizhe(金币+1): 谢谢提示。 2011-08-18 18:00:46
引用回帖:
3楼: Originally posted by butterflyyh at 2011-08-17 18:04:49:
1. 第一行的#!/bin/sh加了;

2.是在LINUX下直接编写的,没有每行结尾的^M符号现象;

3.不知道你有没遇到VASP没运算完就读取i 循环的现象?

如果是运行这个脚本的话,应该没有出现3的现象,我是没有遇见过
专业文献求助
4楼2011-08-17 22:39:29
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[公派出国] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 cqwQDCxMcL3I 2026-09-29 4/200 2026-09-30 13:35 by C3vABopYbZJS
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +10 zju2000 2026-09-26 13/650 2026-09-30 10:26 by yexuqing
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[硕博家园] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:03 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[考博] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +6 GDBe8tDZqE8z 2026-09-28 6/300 2026-09-29 14:41 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +6 IDs3scOF0tjC 2026-09-28 7/350 2026-09-29 14:07 by etmYJ6d2rquH
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 09:59 by JzYBbHIXWSrW
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[考研] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:42 by ZvyPK8n6Nfki
[找工作] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 06:28 by ZvyPK8n6Nfki
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[考博] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 6/300 2026-09-28 22:15 by ez6fGg9abYaj
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-28 19:08 by 9lS3ad5oOymn
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见