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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 分子模拟 » LAMMPS/DL_POLY » 【求助】关于DL-POLY报错的问题
121/2返回列表12下一页
查看: 941  |  回复: 11
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前主题已经存档。

baoer3808

银虫 (小有名气)

[交流] 【求助】关于DL-POLY报错的问题

运行DL-POLY时，报错
At line 651 of file vdw_module.f
Fortran runtime error: Bad value during floating point read

TABLE是自己做的，就是再开始读TABLE 时读不下去了，请问这是什么问题啊？怎样解决，请高人指点，非常感谢！
回复此楼

» 猜你喜欢
1楼 2010-01-14 21:18:52
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★
lei0736(金币+1,VIP+0):谢谢 1-15 20:59
baoer3808(金币+2,VIP+0):谢谢 1-17 19:44
文件格式写错了。明显的。
赞一下(11人) 回复此楼
2楼2010-01-15 16:34:58
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baoer3808

银虫 (小有名气)
可是始终找不到错误在哪里
赞一下 回复此楼
3楼2010-01-17 19:45:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baoer3808

银虫 (小有名气)
可是始终找不到错误在哪里
赞一下 回复此楼
4楼2010-01-17 19:45:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★ ★
lei0736(金币+2,VIP+0):谢谢 1-17 21:30
baoer3808(金币+2,VIP+0): 1-18 10:23
你不把它贴出来，谁也不知道怎么回事。你可以把你的表格文件贴出来，并说明你的势函数类型。大量重复的行你把它删掉即可。
赞一下(10人) 回复此楼
5楼2010-01-17 21:19:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baoer3808

银虫 (小有名气)
DL_POLY UN Potentials  in internal units
  4.0064103E-03  10.0000000E+00      2500
U    N     
  1.0816685E+08  1.0749292E+08  1.0682319E+08  1.0615764E+08
  1.0549623E+08  1.0483894E+08  1.0418575E+08  1.0353663E+08
  1.0289155E+08  1.0225049E+08  1.0161343E+08  1.0098033E+08
  1.0035118E+08  9.9725952E+07  9.9104617E+07  9.8487152E+07
  9.7873535E+07  9.7263741E+07  9.6657746E+07  9.6055527E+07
  9.5457060E+07  9.4862322E+07  9.4271289E+07  9.3683939E+07
  9.3100248E+07  9.2520194E+07  9.1943753E+07  9.1370905E+07
  9.0801625E+07  9.0235892E+07  8.9673685E+07  8.9114979E+07
  8.8559755E+07  8.8007991E+07  8.7459663E+07  8.6914753E+07
    ....             ....         ....            ....
  3.9957845E+01  3.9664879E+01  3.9373760E+01  3.9084476E+01
  3.8797015E+01  3.8511366E+01  3.8227518E+01  3.7945459E+01
  3.7665179E+01  3.7386666E+01  3.7109910E+01  3.6834899E+01
  3.6561624E+01  3.6290072E+01  3.6020234E+01  3.5752098E+01
  3.5485655E+01  3.5220893E+01  3.4957803E+01  3.4696374E+01
  3.4436595E+01  3.4178457E+01  3.3921948E+01  3.3667061E+01
  3.3413783E+01  3.3162105E+01  3.2912018E+01  3.2663510E+01
  3.2416574E+01  3.2171198E+01  3.1927374E+01  3.1685091E+01
  3.1444340E+01  3.1205111E+01  3.0967395E+01  3.0731184E+01
  3.0496466E+01  3.0263234E+01  3.0031477E+01  2.9801188E+01
  2.9572356E+01  2.9344972E+01  2.9119029E+01  2.8894516E+01
  2.8671424E+01  2.8449746E+01  2.8229473E+01  2.8010594E+01
  2.7793103E+01  2.7576990E+01  2.7362247E+01  2.7148866E+01
  2.6936837E+01  2.6726152E+01  2.6516804E+01  2.6308784E+01
  2.6102084E+01  2.5896694E+01  2.5692609E+01  2.5489818E+01
  2.5288315E+01  2.5088091E+01  2.4889139E+01  2.4691450E+01
  2.4495017E+01  2.4299831E+01  2.4105886E+01  2.3913173E+01
  2.3721685E+01  2.3531414E+01  2.3342353E+01  2.3154494E+01
  2.2967830E+01  2.2782353E+01  2.2598056E+01  2.2414932E+01
  2.2232973E+01  2.2052173E+01  2.1872523E+01  2.1694018E+01
  2.1516649E+01  2.1340410E+01  2.1165293E+01  2.0991293E+01
  2.0818401E+01  2.0646611E+01  2.0475917E+01  2.0306311E+01
  2.0137787E+01  1.9970337E+01  1.9803956E+01  1.9638636E+01
  1.9474372E+01  1.9311156E+01  1.9148983E+01  1.8987844E+01
  1.8827736E+01  1.8668649E+01  1.8510580E+01  1.8353520E+01
  1.8197465E+01  1.8042407E+01  1.7888341E+01  1.7735260E+01
  1.7583158E+01  1.7432030E+01  1.7281869E+01  1.7132669E+01
  1.6984424E+01  1.6837128E+01  1.6690776E+01  1.6545362E+01
  1.6400879E+01  1.6257322E+01  1.6114686E+01  1.5972964E+01
  1.5832151E+01  1.5692241E+01  1.5553229E+01  1.5415109E+01
  1.5277876E+01  1.5141524E+01  1.5006047E+01  1.4871441E+01
  1.4737699E+01  1.4604817E+01  1.4472789E+01  1.4341610E+01
  1.4211274E+01  1.4081777E+01  1.3953113E+01  1.3825277E+01
  1.3698264E+01  1.3572069E+01  1.3446686E+01  1.3322111E+01
  1.3198339E+01  1.3075364E+01  1.2953182E+01  1.2831788E+01
  1.2711177E+01  1.2591344E+01  1.2472284E+01  1.2353993E+01
  1.2236466E+01  1.2119697E+01  1.2003683E+01  1.1888419E+01
  1.1773899E+01  1.1660120E+01  1.1547076E+01  1.1434764E+01
  1.1323179E+01  1.1212315E+01  1.1102170E+01  1.0992738E+01
  1.0884014E+01  1.0775995E+01  1.0668676E+01  1.0562053E+01
  1.0456121E+01  1.0350876E+01  1.0246315E+01  1.0142432E+01
  1.0039223E+01  9.9366849E+00  9.8348127E+00  9.7336025E+00
  9.6330503E+00  9.5331518E+00  9.4339031E+00  9.3353002E+00
  9.2373390E+00  9.1400155E+00  9.0433259E+00  8.9472662E+00
  8.8518324E+00  8.7570208E+00  8.6628274E+00  8.5692484E+00
  8.4762800E+00  8.3839184E+00  8.2921599E+00  8.2010008E+00
  8.1104372E+00  8.0204656E+00  7.9310822E+00  7.8422834E+00
  7.7540656E+00  7.6664252E+00  7.5793586E+00  7.4928622E+00
  7.4069325E+00  7.3215661E+00  7.2367593E+00  7.1525088E+00
  7.0688111E+00  6.9856628E+00  6.9030604E+00  6.8210006E+00
  6.7394800E+00  6.6584952E+00  6.5780430E+00  6.4981201E+00
  6.4187231E+00  6.3398489E+00  6.2614941E+00  6.1836555E+00
  6.1063300E+00  6.0295144E+00  5.9532054E+00  5.8774001E+00
  5.8020951E+00  5.7272875E+00  5.6529742E+00  5.5791520E+00
  5.5058180E+00  5.4329691E+00  5.3606024E+00  5.2887147E+00
  5.2173032E+00  5.1463649E+00  5.0758968E+00  5.0058961E+00
  4.9363599E+00  4.8672852E+00  4.7986692E+00  4.7305091E+00
  4.6628020E+00  4.5955452E+00  4.5287357E+00  4.4623710E+00
  4.3964481E+00  4.3309644E+00  4.2659171E+00  4.2013036E+00
  4.1371211E+00  4.0733670E+00  4.0100385E+00  3.9471332E+00
  3.8846483E+00  3.8225812E+00  3.7609294E+00  3.6996902E+00
  3.6388611E+00  3.5784396E+00  3.5184232E+00  3.4588092E+00
  3.3995953E+00  3.3407790E+00  3.2823577E+00  3.2243290E+00
  3.1666905E+00  3.1094398E+00  3.0525744E+00  2.9960920E+00
  2.9399902E+00  2.8842667E+00  2.8289190E+00  2.7739449E+00
  2.7193420E+00  2.6651080E+00  2.6112407E+00  2.5577377E+00
  2.5045969E+00  2.4518159E+00  2.3993925E+00  2.3473245E+00
  2.2956097E+00  2.2442459E+00  2.1932309E+00  2.1425626E+00
  2.0922388E+00  2.0422574E+00  1.9926162E+00  1.9433131E+00
  1.8943460E+00  1.8457129E+00  1.7974117E+00  1.7494403E+00
  1.7017966E+00  1.6544786E+00  1.6074844E+00  1.5608118E+00
  1.5144589E+00  1.4684238E+00  1.4227043E+00  1.3772987E+00
  1.3322048E+00  1.2874209E+00  1.2429448E+00  1.1987749E+00
  1.1549091E+00  1.1113455E+00  1.0680823E+00  1.0251177E+00
  9.8244964E-01  9.4007643E-01  8.9799621E-01  8.5620713E-01
  8.1470741E-01  7.7349524E-01  7.3256883E-01  6.9192640E-01
  6.5156620E-01  6.1148646E-01  5.7168545E-01  5.3216142E-01
  4.9291266E-01  4.5393746E-01  4.1523412E-01  3.7680093E-01
  3.3863624E-01  3.0073835E-01  2.6310563E-01  2.2573641E-01
  1.8862905E-01  1.5178194E-01  1.1519344E-01  7.8861953E-02
     ...             ...            ...          ...

     ...              ...            ...          ...
   -1.0570950E-03 -1.0557750E-03 -1.0544560E-03 -1.0531390E-03
  -1.0518230E-03 -1.0505090E-03 -1.0491970E-03 -1.0478860E-03
  -1.0465770E-03 -1.0452700E-03 -1.0439640E-03 -1.0426600E-03
  -1.0413580E-03 -1.0400570E-03 -1.0387580E-03 -1.0374600E-03
  -1.0361640E-03 -1.0348700E-03 -1.0335770E-03 -1.0322860E-03
  -1.0309960E-03 -1.0297090E-03 -1.0284220E-03 -1.0271380E-03
  -1.0258550E-03 -1.0245730E-03 -1.0232930E-03 -1.0220150E-03
  -1.0207380E-03 -1.0194630E-03 -1.0181900E-03 -1.0169180E-03
  -1.0156470E-03 -1.0143790E-03 -1.0131120E-03 -1.0118460E-03
  -1.0105820E-03 -1.0093200E-03 -1.0080590E-03 -1.0068000E-03
  -1.0055420E-03 -1.0042860E-03 -1.0030310E-03 -1.0017780E-03
   
势函数的表达式  U (r )=库伦相互作用项+（bi+bj）e[（ai+aj-rij）/（bi+bj）]  (core
repulsion)+Dij{ .... }(Morse-type)
请老虎大王帮忙看一下吧
赞一下 回复此楼
6楼2010-01-18 10:23:34
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+3,VIP+0):谢谢 1-18 17:36
1. 查格式。 TABLE 文件，要求实型数的格式为“E15.8”， 整型数格式为“i10"。
2. 你有2500个点，也就是前面有1250行是势能值，然后有1250行是对应的力的数值，你是不是这样搞的？
3. 你的TABLE文件是不是跟你的FIELD文件相对应？怎么你这里只有 U-N，没有其它的？ U 和N是你定义好的原子类型吗？
赞一下(11人) 回复此楼
7楼2010-01-18 17:27:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baoer3808

银虫 (小有名气)
1.一直没明白实型数的格式为“E15.8”， 整型数格式为“i10"是什么意思？
2.对的，势能值和对应的力的数值各1250行
3.与FIELD文件是相对应的，现在我把U-U与N-N的作用也加上了
还是原来那样地提示错误
赞一下 回复此楼
8楼2010-01-18 22:44:33
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

老虎大王

木虫 (著名写手)
★ ★ ★ ★ ★
小木虫(金币+0.5):给个红包，谢谢回帖交流
lei0736(金币+4,VIP+0):谢谢 1-19 09:04
E15.8就是，采用科学计数法（也就是你现在用的这种格式），一个数字共宽度15位，其中小数点后8位。也就是这样：**.*******E-03，小数点前两位（如果是负数，负号也占一位），小数点一位，小数点后八位，“E”字占一位，后边一个正/负号占一位，指数部分占两位。

I10就是，这个整数一共占10位，不足部分，前面以空格补足。

你可以打开一个例子，看看人家的TABLE文件是不是这样的格式。也许不一定要在小数点后保留8位，但一个实型数一共要占15位，这是应该的。也就是说，第二个数，要从第16格以后开始写，第三个数，要从第31格以后开始写。整型数也是这样，第二个数从第21格开始写。
赞一下(11人) 回复此楼
9楼2010-01-18 23:35:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

baoer3808

银虫 (小有名气)
万分感谢！对我帮助很大，我试一下。
赞一下 回复此楼
10楼2010-01-19 12:39:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 baoer3808 的主题更新
121/2返回列表12下一页
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
信息提示
关闭
请填处理意见
关闭