#test 2: check time average of L = 32 is 53.9
p1 = 0.5
p2 = 0.5
for l in [16, 32]:
i = 0
tc = []
h1 = []
t = []
sandbox = OS.box(l, N, p1, p2)
while i <= N:
#drive
sandbox.drive()
sandbox.relax()
if sandbox.steady() == True:
h = 0
tc.append(i)
for jg in range(M): # the number of realisations
#open a file named l
fileout = open(f"{l}_{jg}", "wb")
i = 0
h1 = []
t = []
tc = []
sandbox = OS.box(l, N, p1, p2)
while i <= N:
#drive
sandbox.drive()
sandbox.relax()
tt = i
if sandbox.steady() == True:
tc.append(i)
