| 查看: 1848 | 回复: 5 | ||
牧雪长空铜虫 (正式写手)
|
[求助]
用侯博士的脚本SIGMA测试出的TS都为0！求大侠指点
|
|
为什么我用侯博士的脚本测试的SIGMA值都是零呢，华丽丽的一串零。。如果进入OUTCAR看，会发现，第一个TS其实不是零的。也就是说，应该是我的脚本有错误，可是我又不知道错在哪里。我是在超算上运行的。所以，跟侯博士的原脚本有些差异。下面是我的脚本。 #!/bin/bash rm WAVECAR for i in 0.1 0.12 0.14 do cat > INCAR < SYSTEM = CaCO3 ENCUT = 500 ISTART = 0;ICHARG = 2 NPAR = 4; LREAL = Auto ISMEAR = 0; SIGMA = $i PREC = Accurate LREAL = Auto ; NPAR = 4 ! echo " SIGMA = $i eV "; time /public/software/vasp/vasp-openmpi-ifort/5.3.2/vasp5.3.k TS=`grep "EENTRO" OUTCAR | tail -1 | awk '{printf "%12.6f \n", $5}'` echo $i $TS >> comment done 我在运行的时候，是用qsub提交的作业任务，在这个任务脚本里，路径则是这个程序run_sigma的当下路径了。 |
回复此楼» 猜你喜欢
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有16人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有14人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有15人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有17人回复
-
售SCI-T0P文章，我:8O.5.5.1.O.54,科目齐全,可+急
已经有21人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.5.4,科目全,可+急
已经有21人回复
-
售SCI一区文章，我:8O5.5.1.O5.4,科目全,可伽急
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有8人回复
-
售SCI一区T0P文章，我:8.O55.1.O.54,科目全,可十急
已经有8人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
zbai
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 3466.9
- 红花: 1
- 帖子: 93
- 在线: 104.7小时
- 虫号: 2263643
- 注册: 2013-01-27
- 性别: GG
- 专业: 凝聚态物性 II ：电子结构
2楼2013-10-31 11:34:28
牧雪长空
铜虫 (正式写手)
- 应助: 1 (幼儿园)
- 金币: 1003.8
- 散金: 59
- 帖子: 373
- 在线: 212.7小时
- 虫号: 1465005
- 注册: 2011-10-28
- 专业: 实验地球化学和计算地球化
3楼2013-11-07 19:59:58
490105179
木虫 (正式写手)
- 应助: 30 (小学生)
- 金币: 3400.5
- 散金: 3104
- 红花: 31
- 帖子: 812
- 在线: 332.8小时
- 虫号: 1773039
- 注册: 2012-04-23
- 性别: GG
- 专业: 材料物理化学
4楼2013-11-08 14:42:37
司伟静
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 1661.3
- 散金: 20
- 红花: 1
- 帖子: 143
- 在线: 149.4小时
- 虫号: 2908519
- 注册: 2014-01-02
- 性别: MM
- 专业: 金属非晶态、准晶和纳米晶
5楼2015-06-15 16:06:02
阝东lalala
新虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 94.7
- 散金: 37
- 帖子: 87
- 在线: 53.2小时
- 虫号: 4042245
- 注册: 2015-08-29
- 性别: GG
- 专业: 物理化学
6楼2017-05-09 20:06:28
50