Liionpack: How to simulate a pack of cells with different capacities using the ECM model?

Dear PyBaMM team,

To simulate anomalies on pack level, we have been using liionpack with the PyBaMM ECM cell model due to its simplicity and interpretability.

While our simulations work well for cases such as varying internal resistances, we encounter crashes when simulating cells with different capacities. I have created an issue on the liionpack Github:

Since the liionpack package appears to be in maintenance mode, I’m unsure how frequently issues on GitHub are addressed. I got in touch with Tom Tranter who recommended to get in contact with you over this channel.

Is our issue known to you? Do you have any recommendations on how to approach this?
I would greatly appreciate any insights you can share.

Thank you so much and best regards,
Matthias