我编写了下面的代码来创建一个函数,该函数以.xyz格式创建一个N粒子位置的文件,可以用来显示粒子轨迹。
但是我对某些事情不太了解。
首先,我想创建并编写一个文件,使其显示如此......
Point = 1
p1 x11 y11 z11
p2 x21 y21 z21
其中所有行都是字符串,在这种情况下取i = 1和j = 1,2。
我觉得我绝对错误地编写了我必须附加.xyz文件的代码部分。
我还应该注意到......
co_ord_data = particles[j].str()
返回4个字符串。
(主要代码)
def VMD_init(timestep, filename)
particle_trajectory = open("filename", "a+")
for i in range(timestep):
if i % 5:
particle_trajectory.write("Point = i")
for j in len(particles):
co_ord_data = particles[j].str()
particle_trajectory.write("co_ord_data")
else:
我希望我能够清楚地解释清楚这一点!请随时回复,请我再试一次。
非常感谢任何帮助!
谢谢: - )