Changed `phase` variable to `epoch` for consistency with astronomy and utitlities.py.
Removed code halt when SN r curve is bad. Starting to switch to list for Pymatrix.
Coded some more test for k-correction. Not eyt a solution for error on K.