m = Basemap(projection='cyl')
p0_ll = -73.98, 40.78
p1_ll = 0.08, 51.53
m.drawgreatcircle(p0_ll[0], p0_ll[1], p1_ll[0], p1_ll[1], linewidth=2,
    color='b')
m.drawcoastlines()
m.fillcontinents()
plt.show()
Plotting curved line in Python Basemap