generates a function factory for usage with sedrate_to_multiadm. At h[i],
the sedimentation rate is given by a gamma distribution with parameters
shapes[i] and rates[i]. Between those values, sedimentation rate is linearly interpolated
Outside of the range of h, behavior is determined by the argument rule
which is passed to approxfun. See there for details.
See also
sedrate_to_multiadm()for estimating age-depth models using the outputssed_rate_from_matrix()for other means of defining sedimentation ratessed_rate_gen_from_bounds()for sed. rate generator based on bounds on sedimentation rates.