lyj714 发表于 2021-7-14 20:31
那就完全没必要自己写，直接g_traj（低版本gmx）或者gmx traj（高版本gmx）提取不就完事了？-mol选项就是 ...

好的，谢谢，我试试

khu910 发表于 2021-7-14 19:53
在水体系中，提取各个水分子和吲哚环在每个时刻的位置，并计算各个水分子每个时刻的质心，以及在各个时刻 ...

那就完全没必要自己写，直接g_traj（低版本gmx）或者gmx traj（高版本gmx）提取不就完事了？-mol选项就是提取的分子属性(提供分子编号index文件)，-ox， -of分别可提取坐标和力，-com提取的是分子质心。

何况你自己写的也不是很对，看你写的没有考虑周期性，盒子边界分子完整。
自己写一般可不是你这样的，直接以trjconv -exec的外挂命令执行处理每一帧数据才是高效的。

trr文件可以直接读取吗?我打开是乱码。

lyj714 发表于 2021-7-14 19:39
我觉得你还是说清楚你什么体系，目的是做什么比较好。看你写的脚本属实有些臃肿，多次调用外部命令本来就会 ...

在水体系中，提取各个水分子和吲哚环在每个时刻的位置，并计算各个水分子每个时刻的质心，以及在各个时刻吲哚对各个水分子的力等。

我觉得你还是说清楚你什么体系，目的是做什么比较好。看你写的脚本属实有些臃肿，多次调用外部命令本来就会导致耗时成倍增加，在我看来完全一个awk命令就能完成的所有工作反倒是让你写的无比臃肿和费事难懂。