start on estimating the stream track; add estimate of disruption time
add leading= keyword to decide whether to model a leading or a trailing stream
work on determining the track of the stream
determine stream track in R,vR, etc. rather than lb