24小时热门版块排行榜    

查看: 3487  |  回复: 13

tina张

新虫 (小有名气)

[求助] lammps读data文件求助! 已有1人参与

求助:
        在lammps上用read_data命令,读入多个文件时,出现以下错误:
      Error :Expected floating point parameter in input or data file .  . /read_data.cpp 164
回复此楼

» 猜你喜欢

不忘初心
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖支持 ( 显示支持度最高的前 50 名 )

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

感谢参与,应助指数 +1
data文件的格式不对,在应该出现浮点数的地方不知道出现了啥,发上来看看

» 本帖已获得的红花(最新10朵)

2楼2016-09-25 17:24:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tina张

新虫 (小有名气)

引用回帖:
4楼: Originally posted by lx_PICO at 2016-09-26 09:47:44
貌似是data文件里少了盒子边界的定义,我加上

-1.7046925403980140e-01 8.0636966986120024e+01 xlo xhi
-1.8855463730960764e-01 4.0308151492031038e+01 ylo yhi
-2.5465666718138280e-01 6.25546566666397 ...

十分感谢!!!我还有一个问题,就是多次读入data文件,atom_ID设置在哪里比较合适,我如果用add  offset的命令后,老出现非法定义的错误,我如果要读入data文件是另一个元素,应该在第一个read_data命令,加上keyword extra/atom/types  9(总共九个元素)
这个data文件设置的对么

LAMMPS data file from restart file: timestep=10000,procs=1

16679 atoms

1 atom types

-10e-01  32.517e+01  xlo xhi
-7e-01  14.85e+01 ylo yhi
-10-10e-01  20.244e+01 zlo zhi

Masses

1  15.999400 # o*
   2   1.007970 # h*
   3   63.546  # Cu
   4   63.546  # Cu
   5   63.546  # Cu
   6   63.546  # Cu
   7   63.546  # Cu
   8   63.546  # Cu
   9   12  # C

Atoms

1 9 3.56 14.24 24.92 0 0 0
2 9 3.56 14.24 28.48 0 0 0
3 9 3.56 14.24 32.04 0 0 0
4 9 3.56 14.24 35.6 0 0 0
不忘初心
5楼2016-09-26 15:42:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

tina张

新虫 (小有名气)

送红花一朵
引用回帖:
2楼: Originally posted by lx_PICO at 2016-09-25 17:24:22
data文件的格式不对,在应该出现浮点数的地方不知道出现了啥,发上来看看

你看,这是一部分data文件,但是要读入一个data文件的话,data文件并没有任何问题啊
LAMMPS data file from restart file: timestep = 100000, procs = 8 (30 25 30)

16679 atoms

9 atom types
1 bond types
   1 angle types

   
   Masses

   1  15.999400 # o*
   2   1.007970 # h*
   3   63.546  # Cu
   4   63.546  # Cu
   5   63.546  # Cu
   6   63.546  # Cu
   7   63.546  # Cu
   8   63.546  # Cu
   9   12  # C



ATOMS

1 7 3.56 14.24 24.92 0 0 0
2 7 3.56 14.24 28.48 0 0 0
3 7 3.56 14.24 32.04 0 0 0
4 7 3.56 14.24 35.6 0 0 0
5 7 3.56 17.8 21.36 0 0 0
6 7 3.56 17.8 24.92 0 0 0
7 7 3.56 17.8 28.48 0 0 0
8 7 3.56 17.8 32.04 0 0 0
9 7 3.56 17.8 35.6 0 0 0
10 7 3.56 17.8 39.16 0 0 0
11 7 3.56 21.36 17.8 0 0 0
12 7 3.56 21.36 21.36 0 0 0
13 7 3.56 21.36 24.92 0 0 0
14 7 3.56 21.36 28.48 0 0 0
15 7 3.56 21.36 32.04 0 0 0
16 7 3.56 21.36 35.6 0 0 0
17 7 3.56 21.36 39.16 0 0 0
18 7 3.56 21.36 42.72 0 0 0
19 7 3.56 24.92 17.8 0 0 0
20 7 3.56 24.92 21.36 0 0 0
21 7 3.56 24.92 24.92 0 0 0
22 7 3.56 24.92 28.48 0 0 0
23 7 3.56 24.92 32.04 0 0 0
24 7 3.56 24.92 35.6 0 0 0
25 7 3.56 24.92 39.16 0 0 0
26 7 3.56 24.92 42.72 0 0 0
27 7 3.56 28.48 17.8 0 0 0
28 7 3.56 28.48 21.36 0 0 0
29 7 3.56 28.48 24.92 0 0 0
30 7 3.56 28.48 28.48 0 0 0
31 7 3.56 28.48 32.04 0 0 0
32 7 3.56 28.48 35.6 0 0 0
不忘初心
3楼2016-09-26 09:03:15
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
tina张: 金币+10, ★★★★★最佳答案, 十分感谢!!! 2016-09-26 15:43:37
引用回帖:
3楼: Originally posted by tina张 at 2016-09-26 09:03:15
你看,这是一部分data文件,但是要读入一个data文件的话,data文件并没有任何问题啊
LAMMPS data file from restart file: timestep = 100000, procs = 8 (30 25 30)

16679 atoms

9 atom types
1 bond ty ...

貌似是data文件里少了盒子边界的定义,我加上

-1.7046925403980140e-01 8.0636966986120024e+01 xlo xhi
-1.8855463730960764e-01 4.0308151492031038e+01 ylo yhi
-2.5465666718138280e-01 6.2554656666639701e+01 zlo zhi

以后是可以读的,另外你的read_data格式写对了吗?

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : in.read
  • 2016-09-26 09:46:33, 601 bytes
  • 附件 2 : try1.dat
  • 2016-09-26 09:46:33, 596 bytes
  • 附件 3 : try2.dat
  • 2016-09-26 09:46:33, 590 bytes
4楼2016-09-26 09:47:44
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

引用回帖:
5楼: Originally posted by tina张 at 2016-09-26 15:42:31
十分感谢!!!我还有一个问题,就是多次读入data文件,atom_ID设置在哪里比较合适,我如果用add  offset的命令后,老出现非法定义的错误,我如果要读入data文件是另一个元素,应该在第一个read_data命令,加上key ...

这是定义的问题,仔细看手册再慢慢试吧。

» 本帖已获得的红花(最新10朵)

6楼2016-09-26 19:37:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tina张

新虫 (小有名气)

送红花一朵
引用回帖:
6楼: Originally posted by lx_PICO at 2016-09-26 19:37:30
这是定义的问题,仔细看手册再慢慢试吧。...

那为什么我用你的命令怎么还说我非法定义
variable                     sphere string data.dat
variable                     waterbe string  water1.data
variable                     waterup string  water2.data
read_data                   ${waterbe}
read_data                   water1.data
Reading data file ...
  orthogonal box = (-1 -1 -1) to (325.17 148.5 202.44)
  4 by 4 by 1 MPI processor grid
  reading atoms ...
  37200 atoms
  scanning bonds ...
  2 = max bonds/atom
  scanning angles ...
  1 = max angles/atom
  reading bonds ...
  24800 bonds
  reading angles ...
  12400 angles
Finding 1-2 1-3 1-4 neighbors ...
  2 = max # of 1-2 neighbors
  1 = max # of 1-3 neighbors
  1 = max # of 1-4 neighbors
  2 = max # of special neighbors
read_data                    ${waterup}   add      0
read_data                    water2.data   add      0
ERROR: Illegal read_data command (../read_data.cpp:164)
这是第一个data文件
LAMMPS data file. msi2lmp v3.9.6 / 11 Sep 2014 / CGCMM for Sketch

  37200 atoms
  24800 bonds
  12400 angles
      0 dihedrals
      0 impropers

   9 atom types
   1 bond types
   1 angle types

  
-10e-01  32.517e+01  xlo xhi
-10e-01  14.85e+01 ylo yhi
-10e-01  20.244e+01 zlo zhi
Masses

   1  15.999400 # o*
   2   1.007970 # h*
   3   63.546  # Cu
   4   63.546  # Cu
   5   63.546  # Cu
   6   63.546  # Cu
   7   63.546  # Cu
   8   63.546  # Cu
   9   12  # C
   
Bond Coeffs # harmonic

   1   540.6336     0.9600 # o*-h*

Angle Coeffs # harmonic

   1    50.0000   104.5000 # h*-o*-h*

Atoms # full

      1      1   1 -0.820000     8.018998451     2.540890217     0.955299609   0   0   0 # o*
      2      1   2  0.410000     8.609850421     3.119996727     0.452578854   0   0   0 # h*
      3      1   2  0.410000     7.374236306     2.239286862     0.299946918   0   0   0 # h*
      4      1   1 -0.820000    39.118998451     2.540890217     0.955299609   0   0   0 # o*
      5      1   2  0.410000    39.709850421     3.119996727     0.452578854   0   0   0 # h*
这是第二个data文件

LAMMPS        data        file.        msi2lmp        v3.9.6        /        11        Sep        2014        /        CGCMM        for        Sketch

        37200        atoms                                                                               
        24800        bonds                                                                               
        12400        angles                                                                               
        0        dihedrals                                                                               
        0        impropers                                                                               

        9        atom        types                                                                       
        1        bond        types                                                                       
        1        angle        types       
       
       
-10e-01        3.25E+02        xlo        xhi                                                                       
        -1.00E+00        1.49E+02        ylo        yhi                                                               
        -1.00E+00        2.02E+02        zlo        zhi                                                               
Masses                                                                                               

        1        15.9994        #        o*                                                               
        2        1.00797        #        h*                                                               
        3        63.546        #        Cu                                                               
        4        63.546        #        Cu                                                               
        5        63.546        #        Cu                                                               
        6        63.546        #        Cu                                                               
        7        63.546        #        Cu                                                               
        8        63.546        #        Cu                                                               
        9        12        #        C                                                               
       
Bond        Coeffs        #        harmonic                                                                       

        1        540.6336        0.96        #        o*-h*                                                       

Angle        Coeffs        #        harmonic                                                                       

        1        50        104.5        #        h*-o*-h*                                                       

Atoms        #        full        #        #        #        #        #        #        #                       

        37201        1        1        -0.82        8.018998451        2.540890217        0.955299609        0        0        0        #        o*
        37202        1        2        0.41        8.609850421        3.119996727        0.452578854        0        0        0        #        h*
        37203        1        2        0.41        7.374236306        2.239286862        0.299946918        0        0        0        #        h*
        37204        1        1        -0.82        39.11899845        2.540890217        0.955299609        0        0        0        #        o*
        37205        1        2        0.41        39.70985042        3.119996727        0.452578854        0        0        0        #        h*
不忘初心
7楼2016-09-26 21:15:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

lx_PICO

铁虫 (正式写手)

【答案】应助回帖

引用回帖:
7楼: Originally posted by tina张 at 2016-09-26 21:15:12
那为什么我用你的命令怎么还说我非法定义
variable                     sphere string data.dat
variable                     waterbe string  water1.data
variable                     waterup string  w ...

我那个是为了找错的简化版啊,要针对你自己的输入文件改的。你用我给的try1.dat和try2.dat读取会报错吗?会的话说明你的LAMMPS是旧版的。

» 本帖已获得的红花(最新10朵)

8楼2016-09-26 21:40:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tina张

新虫 (小有名气)

引用回帖:
8楼: Originally posted by lx_PICO at 2016-09-26 21:40:03
我那个是为了找错的简化版啊,要针对你自己的输入文件改的。你用我给的try1.dat和try2.dat读取会报错吗?会的话说明你的LAMMPS是旧版的。...

哦,是这样啊!用你给我的太容易try1.data 和try2.data 读的时候报错了,说是data文件有问题 ,这说明我的lammps 是旧版的么?
不忘初心
9楼2016-09-27 10:51:22
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

tina张

新虫 (小有名气)

送红花一朵
引用回帖:
8楼: Originally posted by lx_PICO at 2016-09-26 21:40:03
我那个是为了找错的简化版啊,要针对你自己的输入文件改的。你用我给的try1.dat和try2.dat读取会报错吗?会的话说明你的LAMMPS是旧版的。...

你好!我想问一下你的lammps的版本是多少?新版的lammps安装包中lib中是否含有cuda软件包?
不忘初心
10楼2016-10-20 11:06:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 tina张 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 070300化学求调剂 +3 苑豆豆 2026-03-20 3/150 2026-03-22 14:48 by ColorlessPI
[考研] 一志愿070300浙大化学358分,求调剂! +3 酥酥鱼.. 2026-03-21 3/150 2026-03-22 11:31 by 杨杨杨紫
[考研] 354求调剂 +7 Tyoumou 2026-03-18 10/500 2026-03-22 11:11 by 人来盛
[考研] 328求调剂,英语六级551,有科研经历 +5 生物工程调剂 2026-03-17 9/450 2026-03-21 23:32 by zhujy1982
[考研] 广西大学材料导师推荐 +3 夏夏夏小正 2026-03-17 5/250 2026-03-21 22:20 by 金昊ML
[考研] 求调剂 +4 要好好无聊 2026-03-21 4/200 2026-03-21 18:57 by 学员8dgXkO
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 材料学硕333求调剂 +3 北道巷 2026-03-18 3/150 2026-03-21 18:17 by 学员8dgXkO
[考研] 299求调剂 +5 shxchem 2026-03-20 7/350 2026-03-21 17:09 by ColorlessPI
[考研] 一志愿天津大学化学工艺专业(081702)315分求调剂 +12 yangfz 2026-03-17 12/600 2026-03-21 03:30 by JourneyLucky
[考研] 初始318分求调剂(有工作经验) +3 1911236844 2026-03-17 3/150 2026-03-21 02:33 by JourneyLucky
[考研] 材料学硕297已过四六级求调剂推荐 +11 adaie 2026-03-19 11/550 2026-03-20 21:30 by laoshidan
[考研] 一志愿 南京航空航天大学大学 ,080500材料科学与工程学硕 +5 @taotao 2026-03-20 5/250 2026-03-20 20:16 by JourneyLucky
[考研] 086500 325 求调剂 +3 领带小熊 2026-03-19 3/150 2026-03-20 18:38 by 尽舜尧1
[考研] 求调剂 +3 暗涌afhb 2026-03-16 3/150 2026-03-20 00:28 by 河南大学校友
[考研] 288求调剂,一志愿华南理工大学071005 +5 ioodiiij 2026-03-17 5/250 2026-03-19 18:22 by zcl123
[考研] 085600材料与化工调剂 324分 +10 llllkkkhh 2026-03-18 12/600 2026-03-19 14:33 by llllkkkhh
[考研] [导师推荐]西南科技大学国防/材料导师推荐 +3 尖角小荷 2026-03-16 6/300 2026-03-16 23:21 by 尖角小荷
[考研] 070303 总分349求调剂 +3 LJY9966 2026-03-15 5/250 2026-03-16 14:24 by xwxstudy
信息提示
请填处理意见