  We present ver. 3 of GLISSANDO, a versatile Monte-Carlo generator for
Glauber-like models of the initial stages of ultra-relativistic heavy-ion
collisions. The present version incorporates the wounded quark model, which is
phenomenologically successful in reproducing multiplicities of particle
production at RHIC and the LHC. Within this model, one can study the nucleon
substructure fluctuation effects, recently explored in p-A collisions. In
addition, the code includes the possibility of investigating collisions of
light nuclei, such as $^{3}$He and $^{3}$H, or the $\alpha$-clustered
$^{7,9}$Be, $^{12}$C, and $^{16}$O, where the deformation of the intrinsic wave
function influences the transverse shape of the initial state. The current
version, being down-compatible, retains the functionality of the previous
releases, such as incorporation of various variants of Glauber-like models, a
smooth $NN$ inelasticity profile in the impact parameter obtained from a
parametrization of experimental data, fluctuating strength of the entropy
deposition, or realistic nuclear distributions of heavy nuclei with
deformation. The code can provide output in the format containing the
event-by-event source location, which may be further used in modeling the
intermediate evolution phase, e.g., with hydrodynamics or transport models. The
interface is simplified, such that in the control input file the user may
supply only the very basic information, such as the collision energy, the mass
numbers of the colliding nuclei, and the model type. GLISSANDO 3 is integrated
with the CERN ROOT platform. The package includes numerous illustrative and
useful ROOT scripts to compute and plot various results.
