I will try to do it later on when I get back home. For one simulation
with many particles it is easy because you just set the number of
partilcles in gendata.py, then output the initial conditions to
"data", then run with "num_fieldlines=1000" in setup.