This function uses Matlab's stream2 function with MEaSUREs Antarctic ice surface velocities to calculate stream lines or flow paths.
Updated for compatibility with the new version of the MEaSUREs toolbox, which accesses the raw .nc data rather than re-formatted .mat data.
Minor update: Now if flowlines are plotted and only one output is requested, the lone output is a handle of plotted flowlines.
Upload error in previous attempt.
Introduced plotting in map or polar stereographic coordinates, added support for ps71 coordinates as seed locations, reduced the possibility of aliasing errors by now loading full resolution MEaSUREs dataset for flow calculation.
Now returns along-path velocity, distance, and time vectors.
Inspired: How to estimate subglacial water routes