PICASO (Planetary Intensity Code for Atmospheric Scattering Observations), written in Python, computes the reflected light of exoplanets at any phase geometry using direct and diffuse scattering phase functions and Raman scattering spectral features.
Please see citation information here: https://natashabatalha.github.io/picaso/credit.html#what-to-cite