L point: (0.5, 0.5, 0.5)
M point: (-0.5, 0.5, 0)
G point: (0., 0.0, 0.0)
vector 1: G-->M
vector 2: G-->L
The two vectors can form a rectangle. You sample a dense k-grid on this rectangle and then do the band structure calculations.
You have to define the other direction that is perpendicular to [111]. And then change the setup in the following lines.
------------
kx = numpy.linspace(0,-0.5,Nk)
ky = numpy.linspace(0.0,0.5,Nk)
kz=numpy.linspace(0.0, 0.5,Nk)
kpoints = []
for i in range(Nk):
for j in range(Nk):
kpoints.append([kx,ky, kz[j]])
----------