Include phase in direct lightcurve array
Include phase in direct lightcurve array

The phase and phase error are now a part of the array returned by the
direct integration method.

TODO: Include similar changes to the arrays returned by the other
methods

JALusk/SNoBoL,JALusk/SuperBoL,JALusk/SNoBoL