  We present an update of the CLUMPY code for the calculation of the
astrophysical J-factors (from dark matter annihilation/decay) for any Galactic
or extragalactic dark matter halo including substructures: halo-to-halo
concentration scatter may now be enabled, boost factors can include several
levels of substructures, and triaxiality is a new option for dark matter
haloes. This new version takes advantage of the CFITSIO and HEALPix libraries
to propose FITS output maps using the HEALPix pixelisation scheme. Skymaps for
$\gamma$-ray and $\nu$ signals from generic annihilation/decay spectra are now
direct outputs of CLUMPY. Making use of HEALPix routines, smoothing by a
user-defined instrumental Gaussian beam and computing the angular power
spectrum of the maps is now possible. In addition to these improvements, the
main novelty is the implementation of a Jeans analysis module, to obtain dark
matter density profiles from kinematic data in relaxed spherical systems (e.g.,
dwarf spheroidal galaxies). The code is also interfaced with the GreAT toolkit
designed for Markov Chain Monte Carlo analyses, from which probability density
functions and credible intervals can be obtained for velocity dispersions, dark
matter profiles, and J-factors.
