So I've been playing with a little N-body code in 3D for gravitational problems and have made an OpenGL visualization to go along with it. I have been generating initial conditions (positions and velocities) using explicit formulas. I was wondering of anyone knew of any resources for getting more "interesting" initial conditions. Right now all I generate is a uniform distribution or a spherical distribution. Is there any way to get positions and velocities for something that looks like a real galaxy or something? Does NASA have some huge CSV file for the milky way or something crazy? Google has failed me thus far. Thanks for any direction!