The generation of conformer ensembles and conformational analysis is an essential component of computational chemistry. While CREST performs these functions very well for single molecules and for clusters of rigid molecules, it is not always up to the task for clusters of flexible molecules. Thus, we have developed a new algorithm which extends CREST to such systems. We have termed this new algorithm LEDE-CREST, i.e., a Low-Energy, Diversity-Enhanced variant of CREST. While, in principle, this algorithm could be implemented within CREST, for the present, we have written a script to automate the implementation of our algorithm. This script is designed to run on Unix computing clusters using Slurm schedulers, with CREST 2 (v. 2.12) and xtb 21 (v. 6.5.0) installed, but could be adapted to function on other architectures, as desired. The script, along with instructions, is freely available on Github at . The script is written in such a way that the user is free to adjust many parameters, such as the values of kpush and alpha used, the value of the ratio for conformer selection, and simulation length, among others. As noted, the numerical parameters discussed here were used successfully for a set of seventeen model asphaltene dimers and trimers including microhydrated versions thereof, as well as for the twenty-four cases described here (five complexes occur in both sets), but one may have to do some testing and adjustment for the specific flexible monomers of interest. Our results show that one must make a judicious choice between LEDE-CREST and the various runtypes of CREST, depending on the exact system under study. For detailed information on the thermodynamics of clusters of flexible monomers, LEDE-CREST is to be preferred, but NCI-CREST is better for clusters of rigid monomers, and CREST is better for single molecules.