Well that's a bit of a problem as you can't linearly interpolate the latlong coordinates to get a straight path. AFAIK you have to use some kind of scheme like spherical interpolation using quaternions.

(PS. And remember that my latlongtovec3 function needs coordinates in radians )