changed spt axis labels, added legend
run multiple fits on a single object and compare the outputs

At the moment, runs JHK and full SpeX fits when given a BD name and
model filename (so, handles everything I was scripting previously).
Saves the full-page plot for each run then plots a compare_results
plot at the end.

Currently untested!