first go at abundance_dist.py