Added a routine to calculate phase at a given time t - right now it contains test time values and can be used with either of the arrival time .dat files
Made a typo, corrected now