code: def great_circle_Npoints(lonlat1r, lonlat2r, N):
      
  
      ratio = np.linspace(0.0,1.0, N).reshape(-1,1)
  
  
      xyz1 = lonlat2xyz(lonlat1r[0], lonlat1r[1])
      xyz2 = lonlat2xyz(lonlat2r[0], lonlat2r[1])
  
      mids = ratio * xyz2 + (1.0-ratio) * xyz1
      norm = np.sqrt((mids**2).sum(axis=1))
      xyzN = mids / norm.reshape(-1,1)
  
      lonlatN = xyz2lonlat( xyzN[:,0], xyzN[:,1], xyzN[:,2])
  
      return lonlatN
docs: N points along the line joining lonlat1 and lonlat2