应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 第一性原理 » Vasp&MedeA » lev00 运行错误
71/1返回列表
查看: 1933  |  回复: 6
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

无所谓109

银虫 (正式写手)

[求助] lev00 运行错误 已有1人参与

lev00安装成功，试用也成功，对TiO2，CO2，C3N4的电荷处理都没问题，但是一处理吸附问题，当POSCAR中原子比较多时候好像就有问题了，不知道你们遇到过没有？请大家帮忙，谢谢啦~！
删除元素等方法都试过了，不起作用
错误提示为：
lev00 运行错误
1jhkkjgjgyu.jpg
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
向金字塔顶攀登
1楼 2014-01-15 00:54:25
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wekow

禁虫 (小有名气)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
fzx2008: 金币+1, 谢谢回帖！ 2014-01-15 23:43:25
无所谓109: 金币+15, ★★★很有帮助, 谢谢您啦！ 2014-01-16 21:17:15
本帖内容被屏蔽

» 本帖已获得的红花（最新10朵）

无所谓109
2楼2014-01-15 15:32:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

无所谓109

银虫 (正式写手)
送红花一朵
引用回帖:
2楼: Originally posted by wekow at 2014-01-15 15:32:36
你打开多大原子数目的文件？做什么处理的时候出错的？详细说一下。

您这么一提醒，我感觉到是不是原子数目太多了。我做表面吸附性质，打开的一个是五十多个原子，八十多个原子，都不行.....原子种树有四种。主要想用其处理电荷密度分布，之前用VESTA做过，但是2D的没有标尺，并且不是很好看。看见侯老师教程里面用origin和gnuplot画的图很好看，所以就想用lev00.
1、编译时按照论坛里面的教程编译的，lev00是最新版的
2、论坛中搜索，有人说是POSCAR/CONTCAR没删除原子名称,但是我删除了也不行

错误信息：

   47        1.95090    3.90159   28.90359    F F F
   48        3.90159    5.85187   21.04785    F F F
   49        3.90159    5.85228   28.90359    F F F
   50        5.85187    3.90159   21.04785    F F F
   51        5.85228    3.90159   28.90359    F F F
______> Atoms in species   4 ( X) <______
   52        3.90159    3.90159   18.98412    F F F
[1]     7.80318    0.00000    0.00000
[2]     0.00000    7.80318    0.00000
[3]     0.00000    0.00000   33.69439              （原子的坐标信息）
Done geometry!
Trying to read k-points from KPOINTS ...
Trying to read k-points from OUTCAR ...
   0.00000    0.00000    0.00000         0.04000
   0.20000    0.00000    0.00000         0.16000
   0.40000    0.00000    0.00000         0.16000
   0.20000    0.20000    0.00000         0.16000
   0.40000    0.20000    0.00000         0.32000
   0.40000    0.40000    0.00000         0.16000
... K-points have been read in successfully from OUTCAR
Trying to read RWIGS from OUTCAR ...
forrtl: severe (24): end-of-file during read, unit 15, file /home/....../......./OUTCAR   (OUTCAR绝对途径)
Image              PC                Routine            Line        Source            
lev00              00000000006D148A  Unknown               Unknown  Unknown
lev00              00000000006CFF86  Unknown               Unknown  Unknown
lev00              0000000000687D60  Unknown               Unknown  Unknown
lev00              00000000006487EE  Unknown               Unknown  Unknown
lev00              0000000000647D2F  Unknown               Unknown  Unknown
lev00              00000000006660FA  Unknown               Unknown  Unknown
lev00              000000000046B43C  Unknown               Unknown  Unknown
lev00              000000000046EEA7  Unknown               Unknown  Unknown
lev00              00000000004031FC  Unknown               Unknown  Unknown
libc.so.6          00000035BEE1D974  Unknown               Unknown  Unknown
lev00              0000000000403109  Unknown               Unknown  Unknown
赞一下 回复此楼
向金字塔顶攀登
3楼2014-01-16 13:08:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wekow

禁虫 (小有名气)
本帖内容被屏蔽
4楼2014-01-16 18:06:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

无所谓109

银虫 (正式写手)
引用回帖:
4楼: Originally posted by wekow at 2014-01-16 18:06:44
不是原子数目的问题，打开几百个原子的文件都可以。也不是vasp5.2新格式的问题，因为lev00现在支持vasp5.2的文件。forrtl: severe (24): end-of-file during read, unit 15 错误一般是读取文件错误，你在编译的时候 ...

谢谢您啦
我按照您的建议弄了一下啊，我自己没有设置RWIGS，这个是按照POTCAR默认的：如果您方便帮忙看看，我也再琢磨一下：
这个是OUTCAR中的
Atomic Wigner-Seitz radii
   RWIGS  =  -1.00 -1.00 -1.00 -1.00
  virtual crystal weights
   VCA    =   1.00  1.00  1.00  1.00
   NELECT =     314.0000    total number of electrons
   NUPDOWN=      -1.0000    fix difference up-down


这个是按照您建议编译的lev00
Trying to read RWIGS from OUTCAR ...
forrtl: severe (24): end-of-file during read, unit 15, file /home/export/base/shangjian/sj/BOC/BOC_dos/BOC001/OUTCAR
Image              PC                Routine            Line        Source            
lev00              000000000058377A  Unknown               Unknown  Unknown
lev00              0000000000582276  Unknown               Unknown  Unknown
lev00              000000000053A050  Unknown               Unknown  Unknown
lev00              00000000004FAADE  Unknown               Unknown  Unknown
lev00              00000000004FA01F  Unknown               Unknown  Unknown
lev00              00000000005183EA  Unknown               Unknown  Unknown
lev00              0000000000431983  get_rwigs_                244  read_vasp_geom.f90
lev00              0000000000433E3C  MAIN__                     90  lev00.f90
lev00              00000000004031FC  Unknown               Unknown  Unknown
libc.so.6          0000003ADA41D974  Unknown               Unknown  Unknown
lev00              0000000000403109  Unknown               Unknown  Unknown

这个是错误中的244  read_vasp_geom.f90的位置：
!
!........ default values for RWIGS from pseudopotentials reading from OUTCAR
!   RWIGS=1.0
!
   write(*,'(a)')'Trying to read RWIGS from '//trim(nameout)//' ...'
   open (15,file=trim(nameout),status='old',FORM='FORMATTED',err=30)
   do isp=1,NSPEC
1     read(15,'(a)') line
      if(index(line,'RWIGS').eq.0) go to 1
      call CutStr(line,NumLin,LinPos,LinEnd,0,0,iErr)
      if(NumLin.ge.6) then
         read(line(LinPos(6):LinEnd(6)),*,err=2,end=2) RWIGS(isp)
      else
         go to 2
      end if
   end do
   write(*,'(a/(10(f10.5,x)))') '... default values for RWIGS: ',(RWIGS(isp),isp=1,NSPEC)
   go to 3
赞一下 回复此楼
向金字塔顶攀登
5楼2014-01-16 21:16:45
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

liwei2888379

铜虫 (小有名气)
你好，我也遇到这个问题了，请问你已经找到解决办法了吗？
赞一下 回复此楼
o.o
6楼2016-05-09 17:17:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

wjnnjw

新虫 (初入文坛)
引用回帖:
5楼: Originally posted by 无所谓109 at 2014-01-16 21:16:45
谢谢您啦
我按照您的建议弄了一下啊，我自己没有设置RWIGS，这个是按照POTCAR默认的：如果您方便帮忙看看，我也再琢磨一下：
这个是OUTCAR中的
Atomic Wigner-Seitz radii
   RWIGS  =  -1.00 -1.00 -1.00 -1. ...

你好~我也遇到了一样的问题 请问lz解决了吗？
赞一下 回复此楼
7楼2017-09-12 17:02:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 Janus2024 的主题更新
71/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 2026国自然函评费到账 +18 羊腰板 2026-08-21 21/1050 2026-08-25 10:24 by ma0526
[基金申请] 2026年的国家社科基金项目通讯评审的新规则与新动向、新挑战 +5 process2012 2026-08-23 7/350 2026-08-25 09:42 by huixian257
[基金申请] 我面上完蛋了 +13 且听虎啸 2026-08-20 14/700 2026-08-25 09:10 by mrkang
[基金申请] 明天应该可查了！？ +5 chengyan1220 2026-08-23 5/250 2026-08-24 23:28 by 我4大白菜
[基金申请] 没有任何消息-是不是就凉了 +7 图啦图啦 2026-08-24 8/400 2026-08-24 22:00 by maomao_da
[基金申请] 能否退出参与的面上项目解除限项 +21 koalala 2026-08-24 24/1200 2026-08-24 19:25 by 家与远方
[基金申请] filecode，4个jtjc了 +14 ziyangfang 2026-08-19 17/850 2026-08-24 18:37 by 哈哈蛤？
[基金申请] 今日不放榜？网传国自然预计 8 月 27 日可查结果 +17 医学老男孩 2026-08-20 21/1050 2026-08-24 14:21 by refreshing11
[基金申请] 估计是周四 +4 archvillain 2026-08-18 4/200 2026-08-24 13:53 by zzuzxg
[基金申请] 朋友圈看到的 +7 wangzilk 2026-08-18 9/450 2026-08-24 10:50 by cmrandy
[基金申请] 让我中一个面上吧！ +13 大萍1987 2026-08-20 16/800 2026-08-24 10:23 by 太傻了
[基金申请] 放榜前的不淡定 20+4 snowwithsea 2026-08-19 17/850 2026-08-24 10:20 by echo8914667
[基金申请] 科研孤儿太难了 +18 我4大白菜 2026-08-20 19/950 2026-08-24 09:47 by 凯恩广盛催化分�
[教师之家] 跳槽后在研项目怎么办？ +5 简单化xn 2026-08-22 10/500 2026-08-23 12:38 by 简单化xn
[基金申请] 时间戳今天，20号变了 +5 archvillain 2026-08-20 5/250 2026-08-22 06:12 by hui_daxiao
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 时间戳又变了 +13 wuchongjun 2026-08-20 19/950 2026-08-21 17:21 by 紫杉醇
[基金申请] 应该是下周三26日公布了吧？ +4 哈哈蛤？ 2026-08-21 4/200 2026-08-21 10:58 by Vivilian
[基金申请] 今天放榜没戏了吧 +9 yuleib84 2026-08-19 11/550 2026-08-21 10:06 by gltch
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭