From 72d24a35ef19b2ddf0075a664f7952df4a071e63 Mon Sep 17 00:00:00 2001 From: imran-marwat Date: Thu, 16 Mar 2017 13:45:49 +0000 Subject: [PATCH] Added MSD_Time list --- Main | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/Main b/Main index 6b1e4ea..229131a 100644 --- a/Main +++ b/Main @@ -155,13 +155,16 @@ If CalculateEnergy == 1.0: Total_KE = [] Total_E = [] -#Uses the velocity verlet method to update position, velocity and force. Outpouts positions to a .xyz file +#Uses the velocity verlet method to update position, velocity and force. Outputs positions to a .xyz file SDTemp = [] SDContain = [] time_last_saved = 0 T = dt*5 l = 0 +MSD_Time_List = [] +MSD_Time = 0 + while t < maxTime: oldForce, PE, KE = totalForce(particles) @@ -191,6 +194,8 @@ while t < maxTime: oldForce = newForce t = t+dt l = l+1 + MSD_Time = MSD_Time + T + MSD_Time_List.append(MSD_Time) #Calls MSD_Calculator method from Particle 3D