# GIVE TIME OF EVOLUTION IN PS AND NUMBER OF POINTS
time=15
snapshots=int(round(time/dt))
steps_to_be_done = int(round(time/dt))
steps_per_snapshot = steps_to_be_done/snapshots
md.init_simulation(nx=nx,
ny=ny,
nz=nz,
a=0.37317,
R=R,
T0=T0,
m=m,
f=f,
L=L,
epsilon=epsilon,
dt=dt)
md.set_particles()
Ts = np.logspace(0.,4.,num=10,endpoint=False)
for T in Ts:
print(T)
md.change_T(T)
md.termalize()
md.simulate(1.,100,10)
md.end()
