计算化学公社

 找回密码 Forget password
 注册 Register
Views: 6688|回复 Reply: 0
打印 Print 上一主题 Last thread 下一主题 Next thread

[GROMACS] 增加窗口用rerun重新得到dhdl.xvg文件,算能量时报错

[复制链接 Copy URL]

12

帖子

0

威望

97

eV
积分
109

Level 2 能力者

各位大神,我用rerun重新得到的dhld.xvg后,用alchemical_analysis.py计算结合能,老是提示我“File "/home/test2/K80-test/AFEG/Bina/solvent/dHdl_files/alchemical-analysis-master/alchemical_analysis/alchemical_analysis.py", line 1296, in <module>    main()
  File "/home/test2/K80-test/AFEG/Bina/solvent/dHdl_files/alchemical-analysis-master/alchemical_analysis/alchemical_analysis.py", line 1225, in main
    nsnapshots, lv, dhdlt, u_klt = parser_gromacs.readDataGromacs(P)
  File "/home/test2/K80-test/AFEG/Bina/solvent/dHdl_files/alchemical-analysis-master/alchemical_analysis/parser_gromacs.py", line 324, in readDataGromacs
    nsnapshots[nf,nf] += unixlike.wcPy(f.filename) - f.skip_lines - 1*bLenConsistency
IndexError: index 17 is out of bounds for axis 1 with size 17”


rerun的命令如下:“gmx_mpi mdrun -stepout 100 -s prod.tpr -rerun prod.xtc -dhdl dhdl
大神们帮我看看该如何解决呀?多谢多谢

本版积分规则 Credits rule

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-23 04:12 , Processed in 0.146867 second(s), 20 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list