Modified:   utilities.py. pick_random_sn() and get_sn() have zero-point
phase at r-band maximum.