write/read VTK

testcases/periodicTest/3D/grains.py

@@ -58,14 +58,14 @@ outf = 5                                        # number of iterations between o
 #
 p = scontact.ParticleProblem(3)
 # Initialise particles boundary
-# p.load_msh_boundaries("mesh.msh", ["Bottom", "Top", "X", "Z"])
-# p.add_particle((r,0.5,0.5), r, 4/3*r**3*np.pi*rhop)
-# p.add_particle((L-2*r,0.5,0.5), r, 4/3*r**3*np.pi*rhop)
-# p.velocity()[1,0] = 0.5
+p.load_msh_boundaries("mesh.msh", ["Bottom", "Top", "X", "Z"])
+p.add_particle((r,0.5,0.5), r, 4/3*r**3*np.pi*rhop)
+p.add_particle((L-2*r,0.5,0.5), r, 4/3*r**3*np.pi*rhop)
+p.velocity()[1,0] = 0.5
 
-# p.write_vtk(outputdir,0,0)
-p.read_vtk(outputdir,0)
 p.write_vtk(outputdir,0,0)
+p.read_vtk(outputdir,0)
+# p.write_vtk(outputdir,0,0)
 
 #Initial time and iteration
 t = 0