Commit dcd83f84 authored by Michel Henry's avatar Michel Henry
Browse files

Merge VTK_periodic into nathan_master_release

parent 9e03773f
Pipeline #8787 passed with stages
in 5 minutes and 19 seconds
......@@ -32,7 +32,8 @@ def genInitialPosition(filename, r, H, lx, rhop) :
y = y.flat
# Add a grain at each centre position
for i in range(len(x)) :
R = r*(1-np.random.random()/2)
# R = r*(1-np.random.random()/2)
R = r # Monodisperse
Xi = x[i] + np.random.random()/10000
p.add_particle((Xi, y[i]), R, R**2 * np.pi * rhop)
p.write_vtk(filename,0,0)
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment