PYESSENCE evolves linearly perturbed coupled quintessence models with multiple (cold dark matter) CDM fluid species and multiple DE (dark energy) scalar fields, and can be used to generate quantities such as the growth factor of large scale structure for any coupled quintessence model with an arbitrary number of fields and fluids and arbitrary couplings.
https://ui.adsabs.harvard.edu/abs/2016arXiv160800910L and https://ui.adsabs.harvard.edu/abs/2017PhRvD..95l3519L