add a fix to distinguish between centred and edge grids when computing number of points over which we should loop and phasevols

also added a warning in the const_dt spacing function to catch the edge case where tmax or max_evolution_time are shorter than the first timestep, which would lead to an empty interpolation table