“第10届量子化学波函数分析与Multiwfn程序培训班将于5月4-8日于北京举办,这是一次性完整、系统学习波函数分析的各种理论知识和全面掌握强大的Multiwfn波函数分析程序使用的最不可错过的机会!请点击此链接查看详情和报名方式,欢迎参加!

“第18届北京科音分子动力学与GROMACS培训班” 将于5月23-26日于北京举办。这是一次性全面、系统学习分子动力学模拟知识和最流行的分子动力学程序GROMACS的关键机会!报名正在进行中,请点击此链接查看详情,欢迎参加!

计算化学公社

 找回密码 Forget password
 注册 Register
Views: 6769|回复 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-4-24 03:02 , Processed in 0.174825 second(s), 20 queries , Gzip On.

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