Generation of composite galaxies in dynamic equilibrium

We present a program to construct disk galaxies in dynamic equilibrium using the moment-based approach. The method is highly modular, allowing components such as bulges, halos, point-particle black holes, and interstellar clouds to be easily added or removed. The gravitational potential and its derivatives for a given galactic configuration are stored on a high-resolution linear/exponential grid in the R-z plane and can be scaled based on component mass. This allows one grid to be used for different sets of initial conditions, drastically reducing computing time. The method is used to initialize and simulate two coaxial, counterrotating disk galaxies undergoing a vertical merger, among other simulations.