discretize_depth function

Interpolate over depth to layered profile