24小时热门版块排行榜    

查看: 3450  |  回复: 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的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

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

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

Tea_Candy

金虫 (小有名气)

引用回帖:
4楼: Originally posted by w76990480 at 2011-08-17 22:39:29:
如果是运行这个脚本的话,应该没有出现3的现象,我是没有遇见过

您好 请教一下 我在运行这个脚本的时候就出现了3中所述的情况
time vasp这句我觉得值得商榷 因为我是在超算上算的
我把time vasp 这句改成了超算上提交任务的命令
结果运行发现 执行上述命令后 任务进入了排队状态 但脚本还是继续执行 结果就出现了3的情况
不知道您有没有办法解决这个问题 谢谢
Life is more than Chemistry.
5楼2011-08-18 20:31:21
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w76990480

荣誉版主 (文坛精英)

super-emuch

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


youzhizhe(金币+1): 谢谢交流。 2011-08-19 16:33:41
引用回帖:
5楼: Originally posted by Tea_Candy at 2011-08-18 20:31:21:
您好 请教一下 我在运行这个脚本的时候就出现了3中所述的情况
time vasp这句我觉得值得商榷 因为我是在超算上算的
我把time vasp 这句改成了超算上提交任务的命令
结果运行发现 执行上述命令后 任务进入了排 ...

我没有加time vasp这句
专业文献求助
6楼2011-08-18 21:30:07
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tea_Candy

金虫 (小有名气)

引用回帖:
6楼: Originally posted by w76990480 at 2011-08-18 21:30:07:
我没有加time vasp这句

不加这句??那怎么执行计算的任务呢?
Life is more than Chemistry.
7楼2011-08-18 23:00:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

w76990480

荣誉版主 (文坛精英)

super-emuch

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

引用回帖:
7楼: Originally posted by Tea_Candy at 2011-08-18 23:00:42:
不加这句??那怎么执行计算的任务呢?

你不加试一下,我也不知道这个time vasp是什么意思
专业文献求助
8楼2011-08-18 23:04:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Tea_Candy

金虫 (小有名气)

引用回帖:
8楼: Originally posted by w76990480 at 2011-08-18 23:04:33:
你不加试一下,我也不知道这个time vasp是什么意思

我试过了 这样不能执行计算任务。。。
能把您的计算脚本贴出来看看吗
Life is more than Chemistry.
9楼2011-08-18 23:14:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

butterflyyh

木虫 (著名写手)

引用回帖:
5楼: Originally posted by Tea_Candy at 2011-08-18 20:31:21:
您好 请教一下 我在运行这个脚本的时候就出现了3中所述的情况
time vasp这句我觉得值得商榷 因为我是在超算上算的
我把time vasp 这句改成了超算上提交任务的命令
结果运行发现 执行上述命令后 任务进入了排 ...

我也是吧time vasp换成提交任务的命令,出现3的情况。现在还没攻克问题。。。谢谢
nothingisimpossible
10楼2011-08-19 09:30:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 butterflyyh 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿吉林大学材料学硕321求调剂 +10 Ymlll 2026-03-18 13/650 2026-03-20 10:19 by 制度的
[考研] 295复试调剂 +4 简木ChuFront 2026-03-19 4/200 2026-03-20 09:33 by yuncha
[考研] 材料080500调剂求收留 +5 一颗meteor 2026-03-13 5/250 2026-03-20 09:14 by xingguangj
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
[考研] 梁成伟老师课题组欢迎你的加入 +9 一鸭鸭哟 2026-03-14 11/550 2026-03-19 17:22 by !本暗一次!
[考研] 321求调剂 +8 何润采123 2026-03-18 10/500 2026-03-19 16:46 by 何润采123
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +11 yangfz 2026-03-17 11/550 2026-03-19 15:06 by houyaoxu
[考研] 287求调剂 +3 晨昏线与星海 2026-03-19 4/200 2026-03-19 12:32 by peike
[考研] 本科郑州大学物理学院,一志愿华科070200学硕,346求调剂 +4 我不是一根葱 2026-03-18 4/200 2026-03-19 09:11 by 浮云166
[考研] 一志愿武理材料305分求调剂 +5 想上岸的鲤鱼 2026-03-18 6/300 2026-03-18 17:53 by 无际的草原
[考研] 297求调剂 +8 戏精丹丹丹 2026-03-17 8/400 2026-03-18 14:30 by laoshidan
[考研] 收复试调剂生 +4 雨后秋荷 2026-03-18 4/200 2026-03-18 14:16 by elevennnne
[考研] 299求调剂 +5 △小透明* 2026-03-17 5/250 2026-03-18 11:49 by 尽舜尧1
[考研] 0703化学336分求调剂 +6 zbzihdhd 2026-03-15 7/350 2026-03-18 09:53 by zhukairuo
[考研] 301求调剂 +4 A_JiXing 2026-03-16 4/200 2026-03-17 17:32 by ruiyingmiao
[考研] 考研化学学硕调剂,一志愿985 +4 张vvvv 2026-03-15 6/300 2026-03-17 17:15 by ruiyingmiao
[考研] 278求调剂 +3 Yy7400 2026-03-13 3/150 2026-03-17 08:24 by laoshidan
[考研] 0854控制工程 359求调剂 可跨专业 +3 626776879 2026-03-14 9/450 2026-03-16 17:42 by 626776879
[考研] 一志愿211 0703方向310分求调剂 +3 努力奋斗112 2026-03-15 3/150 2026-03-16 16:44 by houyaoxu
[考研] 26考研一志愿中国石油大学(华东)305分求调剂 +3 嘉年新程 2026-03-15 3/150 2026-03-15 13:58 by 哈哈哈哈嘿嘿嘿
信息提示
请填处理意见