Hello, I have searched the archives and other sources, but cannot find an R function that would calculate the bearing of a pair of points (lat, lon). It wouldn't be difficult to write such a function, but has this already been done in a package? Thanks in advance, -- Sebastian P. Luque