How to extract the distance and transport matrices from Scipy's wasserstein_distance?

Question

The scipy.stats.wasserstein_distance function only returns the minimum distance (the solution) between two input distributions, p and q