int eventspercycle = 20; // # events per sphere per cycle int N = 1000; // number of spheres double initialpf = 0.01; // initial packing fraction double maxpf = 0.05; // max packing fraction double temp = 0.2; // initial temp., use 0 for zero velocities double growthrate = 0.001; // growth rate double maxpressure = 100.; // max pressure char* readfile = new // if new, creates new packing char* writefile = write.dat // after mp2 char* datafile = stats.dat // data file up to mp2