Generates N-body models through Eddington inversion, and paints stellar probabilities on top using a distribution-function based approach. See EP20.
nbody.py generates the models. See this link for a brief introduction.
npaint.py paints stars on top. See this link.
nbody2hdf5.py generates Gadget2 compatible HDF5 files from 'nbody.py' models. See this link.
To run on python2, add from __future__ import print_function .