Python implementation of the Heinselman and Nicolls Bayesian reconstruction algorithm used to derive resolved ion drift velocity vectors.