ion and plot it as a big blue star: 

[intersectlat,intersectlon] = pathcrossingps71(glat,glon,tracklat,tracklon); 
plotps(intersectlat,intersectlon,'bp','markersize',20,'markerfacecolor','blue')

%% Citing AMT
% If this function or any other part of Antarctic Mapping Tools is useful for you, please
% cite the paper that describes AMT.  
% 
% Greene, C. A., Gwyther, D. E., & Blankenship, D. D. Antarctic Mapping Tools for Matlab. 
% _Computers & Geosciences_. 104 (2017) pp.151-157. <http://dx.doi.org/10.1016