This method takes a velocity model and creates a vector containing
slowness-depth layers that, hopefully, adequately sample both
slowness and depth so that the travel time as a function of
distance can be reconstructed from the theta function.
Checks to make sure no slowness layer spans more than maxRangeInterval
and that the (estimated) error due to linear interpolation is less than
maxInterpError.