# Input file for part3d

# maximum number of particles to place (fewer will be used if they don't fit)
particles	10000
# all particles have this radius and mass
rad		0.005
mass		0.005

# to create a volume of particles: number, xstart, xend, ystart, yend, zstart, zend
#block		10000  0.8 1.0  0.0 1.0  0.0 1.0
#block		100000  -0.5 0.5  -0.5 0.5  -0.5 0.5
# same block, but force particles to be non-intersecting
#niblock	10000  0.8 1.0  0.0 1.0  0.0 1.0

# or, create a block of randomly-placed particles in a sphere (num, center, radius)
nisphere	10000  0.0 0.0  1.0
# or, same but force non-intersecting
#nisphere	10000  0.5 0.5 0.5  0.2

# to read in a list of particle positions from a file: filename,
#    translate x, y, z, and velocity add u, v, w
#read_part	output100.part	0.0 0.0 0.0  -1.0 0.0 0.0

# to create a strand of spaghetti/cable/rope: length per particle, number of particles, E, xstart, xend, ystart, yend, zstart, zend
#strand 0.0025 1000 3.0e+10 0.45330993288 0.5406039187 0.0025  0.42679000056 0.55287886642 0.95
# we'd also reccommend a small dt for this contact sim:
# dt		0.000001

# set simulation dt lower for higher rk, E
# changed coeff. of friction, from 
dt		0.00005
end_time	1000.0
output_dt	0.01		write output files this often
#output_dt	0.025		write output files this often

# octree parameters (20, 10 are decent settings for most sims)
max_ppc		20		maximum number of particles per level
max_levels	20		maximum number of levels in octree, root=1

# apply r^-2 attraction force between all particles
use_self_grav	yes
G		0.01		gravitation constant
delta		0.005		gravitation anti-discontinuity addend
theta		1.5		B&H parameter
direct		no		compare to direct method?

# apply earth-like, constant-vector gravity
use_uniform_grav yes
g		0.0	-1.0

# use contact forces? set stiffness (rk) and damping (rc)
use_contact	yes
rk		1.0e+3
rc		1.0e+2
contact_per_grav	200
# for strands, part3d will compute rk and b(ending) from E and r

# the outer bounds of the computational domain, particles outside of
#   this bound will be removed from the simulation
comp_domain	-1.2 1.2  -1.2 1.2
boundary	wall

density		512		resolution of the density image
image_size	256		resulution of the point image

outfile_root	output
write_pgm	no
write_gif	no
write_png	yes
write_rad	no
write_part	yes